Magnaporthales is an order in Ascomycota within Fungi that contains over 200 species, including pathogens of cereals and grasses as well as saprobes on submerged wood[@b1][@b2][@b3]. The best-studied species in Magnaporthales is the rice blast fungus *Pyricularia* *oryzae* (= *Magnaporthe oryzae*), which is one of the most devastating threats to food security worldwide[@b4][@b5]. This pathogen is better known for its *Pyricularia* asexual state ([fig. 1](#f1){ref-type="fig"}), which has been used for decades as a paradigm for fungal leaf infection studies[@b6][@b7]. The rice blast fungus also has the recently discovered root-infecting capability[@b8]. *Gaeumannomyces graminis* (take-all pathogen) and *Magnaporthiopsis poae* (summer patch pathogen), two other model species in Magnaporthales, are root-infecting pathogens of cereals and grasses. In addition to the plant pathogens, several saprotrophic wood-inhabiting genera (e.g., *Pseudohalonectria* and *Ophioceras*) are also placed in Magnaporthales[@b9].

The studies of Magnaporthales date back to the late 19th century when the rice stem rot fungus *Nakataea oryzae* (= *Magnaporthe salvinii*) and the crabgrass pathogen *Pyricularia grisea* were first recorded[@b10][@b11]. Cannon[@b12] established the family Magnaporthaceae to accommodate these grass-associated fungi centered on *Magnaporthe*. Historically, Magnaporthaceae was placed in various orders in Ascomycota, such as Diaporthales[@b13], Phyllachorales, and Xylariales[@b14], due to a lack of robust phylogenetic evidence and paucity of convincing phenotypic characters. A new order Magnaporthales was recently proposed for these fungi[@b2]. Molecular phylogenetic studies indicated that Magnaporthales belongs to the subclass Sordariomycetidae (Sordariomycetes, Ascomycota). However, consensus regarding additional phylogenetic affinities has not yet been reached. For example, in Zhang and Blackwell\'s[@b15] 18S rDNA tree, Magnaporthales formed a sister clade with Ophiostomatales. Huhndorf et al.[@b16] generated a 28S rDNA tree that suggested a close relationship of Magnaporthales with Sordariales, Chaetosphaeriales and Boliniales. In the 18S rDNA tree of Thongkantha et al.[@b2], Magnaporthales formed a sister clade to Diaporthales and Ophiostomatales. In the same paper, however, the 28S rDNA tree grouped it with Chaetosphaeriales and Sordariaceae. Zhang et al.[@b17] generated a four-gene phylogeny and found the rice blast fungus to be most closely related to Diaporthales; however, no other Magnaporthales taxa were included in the analysis. Poor taxon sampling due to a lack of genome data from Magnaporthales and phylogenetic heterogeneity among gene loci[@b18][@b19] likely explain the conflicting tree topologies in the studies described above. Indeed, the current public availability of genome sequences from Magnaporthales is limited to the three model species: *Pyricularia oryzae*, *Magnaporthiopsis poae*, and *Gaeumannomyces graminis* var. *tritici*[@b7]([www.broadinstitute.org](http://www.broadinstitute.org)).

Therefore, despite their pathogenic diversity and economic importance, several key questions regarding the evolution of Magnaporthales remain unanswered. 1) Does Magnaporthales constitute a monophyletic group? 2) What is the phylogenetic position of Magnaporthales within Sordariomycetes? 3) What are the phylogenetic relationships among key taxa in Magnaporthales? and, 4) How have pathogenicity and phenotypes evolved among Magnaporthales taxa? In order to address these issues, we reconstructed the Magnaporthales phylogeny using a phylogenomic approach. Phylogenomics can often overcome problems associated with single locus or multi-locus phylogenetic analyses. To this end, we generated transcriptome data (\~420 M paired-end reads) from 21 species of Magnaporthales, of which five were chosen for draft genome sequencing. Our work resulted in the largest collection of genome data to date for this poorly studied group and offers the opportunity to resolve the Magnaporthales phylogeny and to understand the evolution of its constituent lineages.

Results and Discussion
======================

Transcriptome sequencing from 21 Magnaporthales species
-------------------------------------------------------

We generated \~420 M paired-end reads (totaling \~65 Gbp) from the 21 targeted Magnaporthales taxa ([tables S1 and S2](#s1){ref-type="supplementary-material"}). These data were assembled into a total of 315 K contigs with 13--31 K contigs from each taxon (N50 = 1.3--2.1 Kbp). To aid detection of potential gene duplications among Magnaporthales (see Methods), five of these taxa were subjected to complete genome sequencing, resulting in 35--43 Mbp genome assemblies that encoded \~12 K genes from each lineage ([table S3](#s1){ref-type="supplementary-material"}). The mRNA-seq data have been submitted to NCBI (SRP050514, Bioproject number: PRJNA269089), whereas the results of the genome sequences and assembled contigs are presented at the project ([dblab.rutgers.edu/FunGI/](http://dblab.rutgers.edu/FunGI/)). In this paper, we focused on reconstruction of the Magnaporthales tree of life based on a set of highly conserved "core" genes. To this end, we selected 309 individual CEGMA (Core Eukaryotic Genes Mapping Approach; Parra et al.[@b38]) gene markers that are largely free of gene duplication in Magnaporthales and other Ascomycota (see Methods).

Phylogenomics resolves the Magnaporthales tree of life
------------------------------------------------------

To address the phylogenetic position of Magnaporthales among Sordariomycetes, a maximum likelihood (ML) tree was generated using 83,616 aligned amino acids from 226 single-copy genes within six Magnaporthales species and 15 non-Magnaporthales species representing the major lineages of Pezizomycotina with two Saccharomycetes as outgroup ([fig. 1](#f1){ref-type="fig"}, [table S4](#s1){ref-type="supplementary-material"}). The ML tree provided strong bootstrap support (90--100%) for all interior nodes in the tree except for the union of *Cryphonectria parasitica*, *Grosmannia clavigera*, and Magnaporthales (41%). Bayesian inference (BI) using the same dataset showed a congruent and highly supported topology (see [fig. S1](#s1){ref-type="supplementary-material"}). The topology of a well-supported coalescent model-based tree using this dataset of 23 species ([fig. S2](#s1){ref-type="supplementary-material"}) was also consistent with the results shown in [fig. 1](#f1){ref-type="fig"}, except for the swapping of branches containing *Cryphonectria parasitica* and the *Neurospora crassa*-*Chaetomium globosum* clade. Therefore, the phylogeny of the major Sordariomycetidae lineages appears to be largely resolved, except for the relationships among Sordariales, *Cryphonectria parasitica*, and *Grosmannia clavigera.*

To study in detail the relationships within Magnaporthales, a ML tree was built using a 226-protein alignment that included 82,715 amino acids from 24 Magnaporthales and five outgroup taxa. This ML tree provided robust bootstrap support (98--100%) for the monophyly of the Magnaporthales and most interior nodes therein, except for the branch uniting *Gaeumannomyces radicicola* and *G. graminis* species (55%; [fig. 2](#f2){ref-type="fig"}). The BI tree and the coalescent model-based tree using this dataset ([fig. S3 and S4](#s1){ref-type="supplementary-material"}) resulted in very similar topologies. The minor differences between the multi-protein ML tree and the coalescent tree include: 1) The poorly supported paraphyletic relationship between *Gaeumannomyces graminis* var. *graminis* and *G. radicicola* in the ML tree ([fig. 2](#f2){ref-type="fig"}), that was replaced with a monophyletic relationship in the coalescent tree ([fig. S4](#s1){ref-type="supplementary-material"}); and 2) The positions of *Cryphonectria parasitica* and *Neurospora crassa* were swapped in the two trees.

In summary, we inferred well-supported phylogenies of Sordariomycetes and Magnaporthales using \>200 conserved single-copy gene markers. Notably, the consistent topologies generated using multi-protein concatenation and the concatenation-free coalescent model suggests that our results are not biased by phylogenetic heterogeneity across genetic loci that may mislead multi-protein phylogeny construction (e.g., Song et al.[@b20]).

The phylogenetic position of Magnaporthales within Sordariomycetes
------------------------------------------------------------------

Ascomycota, the largest group within Fungi includes three subphyla, Taphrinomycotina, Saccharomycotina, and Pezizomycotina, which are differentiated by the ascus dehiscence mechanism and fruiting body structure[@b21][@b22]. The order Magnaporthales, together with other filamentous ascomycetes, comprise the class Sordariomycetes that belongs to subphylum Pezizomycotina. Our phylogenomic analysis resolved relationships among the classes in Pezizomycotina included in the study ([fig. 1](#f1){ref-type="fig"}). The Orbiliomycetes-Pezizomycetes clade (*Monacrosporium haptotylum* and *Tuber melanosporum*) represents an early diverging lineage in Pezizomycotina. The affinity between these two classes received stronger support from our analyses when compared to a previous study[@b21]. This result is consistent with shared Pezizomycotina characters such as the basic ascus type, the simple spore dehiscence mechanism, and the exposed hymenium[@b21][@b23]. The classes Eurotiomycetes (e.g., *Aspergillus niger*), Lecanoromycetes (*Cladonia grayi*), and Dothideomycetes (*Phaeosphaeria nodorum*) form a strongly supported clade. Their monophyletic relationship was recovered in earlier phylogenetic studies, however with lower bootstrap support than reported here[@b21][@b22][@b24][@b25]. Sordariomycetes and Leotiomycetes (e.g., *Blumeria graminis*) are the later diverging lineages and their monophyletic relationship is consistent with previous studies[@b18][@b19][@b21][@b24]. The shared characters between these two lineages include unitunicate and inoperculate asci. The latter is regarded as the advanced spore dispersal mechanism in Pezizomycotina[@b22][@b24]. Given these results, we surmise that single-copy CEGMA genes are useful phylogenetic markers that appear to recapitulate and consolidate previously suggested relationships amongst Pezizomycotina.

Although Magnaporthales formed a strongly supported monophyletic clade in Sordariomycetes ([fig. 1](#f1){ref-type="fig"}), the phylogenetic position of this lineage relative to Ophiostomatales and Diaporthales is less well resolved using different phylogenetic methods (i.e., compare [figs. 1](#f1){ref-type="fig"}, [S1, S2](#s1){ref-type="supplementary-material"}). This was also the case in previous analyses of morphological and molecular data[@b2][@b16][@b21]. Nonetheless, our analyses support a sister group relationship between Magnaporthales and Ophiostomatales; the latter is represented here by the mountain pine beetle-associated blue stain fungus *Grosmannia clavigera*. This result is consistent with the shared characteristics between the two orders of non-stromatic perithecia and a hyphomycetous asexual state. In contrast, Diaporthales typically have stromatic perithecia and a coelomycetous asexual state ([fig. 1](#f1){ref-type="fig"}). The grouping of Magnaporthales and Ophiostomatales (*Gromannia clavigera* or *Ophiostoma piliferum*) was also reported in two previous phylogenomic analyses that included fewer taxa[@b19][@b26].

Phylogenetic relationships among Magnaporthales taxa
----------------------------------------------------

In our phylogenomic analyses, three well-supported major clades are recognized within Magnaporthales ([fig. 2](#f2){ref-type="fig"}). These three clades correspond to Magnaporthaceae, Pyriculariaceae and Ophioceraceae, the recently classified families in Magnaporthales that were identified using a 2-locus phylogeny[@b27]. Most species in clade C are associated with grass roots, use a hyphopodium as the penetration structure, and have a phialophora-like asexual state. The exceptions are *Buergenerula spartinae*, *Bussabanomyces longisporus*, *Nakataea oryzae*, and *Omnidemptus affinis* which are associated with the above-ground parts of host plants. The clade B species that cause blast diseases in Poaceae and other monocot plants use an appressorium as the host penetration structure and have a *Pyricularia* asexual state. The saprotrophic *Ophioceras* and *Pseudohalonectria* species form the basal clade A, which are associated with submerged wood tissue. Currently the asexual state of the majority of these species has not yet been observed[@b3].

Recent phylogenetic studies support the idea that terrestrial habits and a saprotrophic lifestyle are plesiomorphies of the Sordariomycetes fungi[@b17][@b28][@b29]. Our results indicate that the clade A in Magnaporthales shifted from a terrestrial to an aquatic environment but maintained the ancestral saprotrophic lifestyle. In addition to the Wood clade in Magnaporthales, aquatic lineages occur in several other orders in Sordariomycetes, including Halosphaeriales, Sordariales, Diaporthales, and Xylariales[@b17]. Therefore, such a habit shift has apparently occurred several times in Sordariomycetes[@b17][@b28].

Aerial and root infections are two strategies employed by plant pathogens to obtain host nutrients[@b30]. Species in clade B in Magnaporthales have adapted to aerial infection. The reproductive hyphae of the *Pyricularia* asexual state are long and erect, which facilitates conidium dispersal. After conidia attach to the host, a specialized structure known as an appressorium develops from conidia and penetrates the leaf cuticle using high turgor pressure, which is essential for fungal invasion of the host aerial parts ([fig. 2](#f2){ref-type="fig"})[@b30]. However, most clade C members appear to lack aerial infection ability and have adapted to root infection. Their asexual states are phialophora-like, which is characterized by simple and short reproductive hyphae[@b12][@b31]. These taxa form lobed and swollen hyphal structures named hyphopodia to penetrate the host root tissue and obtain nutrients ([fig. 2](#f2){ref-type="fig"})[@b30].

Traditionally, ascospore morphology is the basis for genus and species delimitation in Magnaporthales and in many other ascomycetes. For example, in Magnaporthales, species that produce spindle-shaped ascospores are placed in *Magnaporthe*, whereas those that produce filiform ascospores are in *Gaeumannomyces*. The results of our study indicate that morphology-based generic concepts are unreliable as systematic markers; e.g., *Magnaporthe* and *Gaeumannomyces* are polyphyletic, whereas ecological and pathogenicity features correspond better to the phylogenetic history of Magnaporthales. The clear separation of the Blast clade and the *Magnaporthe* type species, *M. salvinii*, indicates that the rice blast fungus does not belong to the genus *Magnaporthe*, which was mistakenly classified in the 1970s based on ascospore morphology[@b14]. *Pyricularia*, the older asexual generic name for the rice blast fungus is informative and legitimate[@b32], and therefore is suggested as the correct scientific name for this model taxon[@b27][@b33][@b34]. Taxonomic revision of *Magnaporthe*, *Gaeumannomyces*, and other Magnaporthales is underway or already published in other papers[@b35][@b36].

In conclusion, our results demonstrate that Magnaporthales is a robustly supported monophyletic order within Sordariomycetes with its sister group provisionally being Ophiostomatales. The early diverging lineage of Magnaporthales (clade A) is mostly aquatic and saprotrophic on wood substrates. The pathogenicity of the remaining lineages is likely a derived feature with the B clade evolved to infect host aerial parts using appressoria and clade C adapted to plant root infection with hyphopodia. Our results significantly enrich sequence data from the order Magnaporthales and, more importantly, establish a framework for future comparative genomic and functional studies to address the evolution and ecology of this important lineage.

Methods
=======

Fungal strains
--------------

Based on previously generated six-locus phylogenies[@b3][@b35][@b36], we targeted 21 Magnaporthales species ([table S1](#s1){ref-type="supplementary-material"}) for transcriptome sequencing, of which five (*Magnaporthe salvinii*, *Magnaporthiopsis incrustans*, *Magnaporthiopsis rhizophila*, *Ophioceras dolichostomum*, and *Pseudohalonectria lignicola*) were also used for whole-genome sequencing. The sequencing, annotation, and analyses of the genome data will be presented elsewhere (see [table S3](#s1){ref-type="supplementary-material"} for general information). The species sampled were obtained from the American Type Culture Collection, Manassas, Virginia (ATCC), the Centraalbureau voor Schimmelcultures Fungal Biodiversity Centre, the Netherlands (CBS) or the Rutgers Mycological Herbarium, New Brunswick, New Jersey (RUTPP).

Transcriptome sequencing
------------------------

The transcriptomes from 21 Magnaporthales taxa ([table S1](#s1){ref-type="supplementary-material"}) were generated using Illumina sequencing. Total RNA was extracted using the RNeasy Plant Mini Kit (Qiagen, CA) from 1--2 week old mycelium grown in four different conditions: PDA under light, PDA in dark, minimal medium agar (MM; 10 g sucrose, 1 g Ca(NO~3~)~2~·4H~2~O, 0.2 g KH~2~PO~4~, 0.25 g MgSO~4~·7H~2~O, 0.15 g NaCl, 15 g agar in 1 liter distilled water) under light, and MM in dark. The TruSeq RNA Sample Preparation Kit (Illumina, CA) was used to construct mRNA libraries based on standard PCR methods. The libraries were sequenced using 300 cycle paired-end reagents on the MiSeq (Illumina) following the manufacturer\'s instructions.

Transcriptome assembly and protein prediction
---------------------------------------------

Transcriptome data from the five Magnaporthales taxa for which we generated whole-genome data ([table S3](#s1){ref-type="supplementary-material"}) were trimmed to remove adaptors and assembled using Trinity with guidance of the corresponding in-house genome assemblies[@b37]. The mRNA-seq data from each of the 16 remaining Magnaporthales taxa ([table S1 and S2](#s1){ref-type="supplementary-material"}) were cleaned and then assembled using CLC workbench 7 (<http://www.clcbio.com/>) under the default settings. Briefly, reads were cleaned to remove low-quality regions followed by adaptor sequence trimming. The resulting reads were used for *de novo* assembly. The reads were mapped to resulting contigs and the consensus contig sequences were optimized during the mapping procedure. The contig sequences were split at regions with \<2X coverage and the consensus sequences (length ≥ 200 bp) were exported for further analyses.

Protein sequences from the five Magnaporthales species with whole-genome data were predicted using a combination of three gene prediction programs (i.e., SNAP, Augustus, and the self-training GeneMark-ES). The predicted protein sequences were used to build a local database that included complete proteomes from 20 other Ascomycota that are available in the public domain ([table S4](#s1){ref-type="supplementary-material"}). The assembled contig sequences for the remaining 16 taxa were used as a query to search the local database using BLASTx (*e*-value ≤ 1e-10). The translated protein sequence from each query contig was derived from the BLASTx output.

Selection of gene markers
-------------------------

We chose the CEGMA (Core Eukaryotic Genes Mapping Approach) gene set[@b38] as markers to reconstruct the phylogeny. CEGMA proteins represent 458 KOGs (eukaryotic orthologous groups) that are conserved among eukaryotes[@b39]. To avoid misleading effects caused by gene duplication (e.g., Qiu et al.[@b40]), we limited the analyses to single-copy CEGMA genes. To this end, we downloaded the seven whole-proteome datasets that are annotated in the KOG database (<http://www.ncbi.nlm.nih.gov/COG/>), including those from nematode (*Caenorhabditis elegans*), common fruit fly (*Drosophila melanogaster*), human (*Homo sapiens*), mouse-ear cress (*Arabidopsis thaliana*), baker\'s yeast (*Saccharomyces cerevisiae*), fission yeast (*Schizosaccharomyces pombe*) and microsporidian parasite (*Encephalitozoon cuniculi*). We mapped a total of 25 complete proteomes, including those from the five in-house Magnaporthales taxa ([table S3](#s1){ref-type="supplementary-material"}) and 20 other Ascomycota taxa ([table S4](#s1){ref-type="supplementary-material"}), to these seven KOG-annotated proteomes based on the top hit from BLASTp (*e*-value ≤ 1e-5) searches. KOGs with more than one gene in two or more mapped Ascomycota taxa were removed, leaving 309 single-copy CEGMA KOG genes that were used for downstream analyses.

Construction of single-protein alignments
-----------------------------------------

We identified the genes corresponding to the 309 single-copy CEGMA genes in each of the tested taxa using a reciprocal-best-hit strategy. Ascomycota proteomes or proteins translated from EST contigs were used as query to search a local database comprising all CEGMA proteins using BLASTp (*e*-value ≤ 1e-5) and *vice versa*. For each CEGMA gene, sequences from different taxa were collected and aligned using MAFFT version 6.923b (\--auto)[@b41]. The alignments were trimmed using Gblocks (-b4 = 5, -b5 = h)[@b42]. The resulting alignments were then processed with T-COFFEE version 9.03[@b43] to remove poorly aligned residues (conservation score ≤ 5) within conserved sequence blocks. Columns with missing data (\>80%) were removed and short alignments (\<100 amino acids) were discarded. The remaining alignments were used for phylogenetic tree reconstruction. We generated a Pezizomycotina dataset (226 genes) including six Magnaporthales species and 15 other species representing major lineages in Pezizomycotina with two Saccharomycetes as outgroup ([fig. 1](#f1){ref-type="fig"}), and a Magnaporthales-focused dataset (another 226 genes) including 24 Magnaporthales taxa and five closely related Sordariomycetes species ([fig. 2](#f2){ref-type="fig"}).

Phylogenetic reconstruction using concatenated multi-protein alignments
-----------------------------------------------------------------------

For each multi-protein alignment, we first inferred the best protein evolutionary model using Prottest version 3.2[@b44]. The LG + Γ + F model was identified as the top model for both multi-protein alignments (i.e., those used to infer the trees shown in [figs. 1](#f1){ref-type="fig"} and [2](#f2){ref-type="fig"}). Maximum likelihood trees were built using RAxML version 7.2.8[@b45] under the selected model with branch support inferred using 1000 bootstrap replicates. Bayesian trees were built using MrBayes version 3.2.2[@b46]. Because MrBayes does not include LG model[@b47], we modified the source code 'model.c' to incorporate it and then compiled the MPI-version executable. For each alignment, the analyses were performed with two independent runs each with one million generations. Trees were sampled every 100 generations. Consensus trees were built and the posterior probability values for all branches were calculated after removing 25% of the trees as burn-in ([figs. S1 and S3](#s1){ref-type="supplementary-material"}).

Tree reconstruction using the coalescent model
----------------------------------------------

For the Pezizomycotina dataset, we reconstructed a maximum pseudo-likelihood tree ([fig. S2](#s1){ref-type="supplementary-material"}) using a coalescent model that does not require the concatenation of multiple single-gene alignments[@b48]. Statistical support was estimated using 100 multilocus replicates generated using Seo\'s method[@b49]. For each replicate, we randomly sampled 226 genes with replacement. For each sampled gene, a pseudo-alignment was generated by random sampling of amino acid site from the corresponding original alignment with replacement. A ML tree was built for each pseudo-alignment using IQtree[@b50] under the best amino acid evolutionary model selected on the fly. The ML trees were rooted using *Yarrowia lipolytica* or *Monacrosporium haptotylum* (in case of missing data in *Y. lipolytica*) as outgroup. The rooted trees were then used for maximum pseudo-likelihood tree construction using MP-EST[@b48] under the default settings. The resulting 100 maximum pseudo-likelihood replicate trees were summarized following a majority rule using 'Consense' function in the Phylip package (<http://evolution.genetics.washington.edu/phylip.html>). The coalescent model-based tree for the Magnaporthales-focused dataset was generated following a similar procedure ([fig. S4](#s1){ref-type="supplementary-material"}). The single gene trees were rooted using *Verticillium dahliae* or *Botrytis cinerea* (in case of missing data in *V. dahliae*).

Author Contributions
====================

J.L. and G.C. conducted most of the experiments in collaboration with N.E.W. who was in charge of Illumina DNA sequencing. H.Q. and G.C. completed data analysis. N.Z. and D.B. conceived and directed the project. J.L., H.Q., G.C., D.B. and N.Z. wrote and reviewed the manuscript.

Supplementary Material {#s1}
======================

###### Supplementary Information

Supplementary Information

This work was supported by a grant from the National Science Foundation (DEB 1145174) awarded to N. Zhang and D. Bhattacharya. We thank Udi Zelzion (Rutgers University) for assistance in assembling the transcriptome data.

![Maximum likelihood tree of 21 Pezizomycotina species and two Saccharomycetes used as outgroup species based on 83,616 amino acid positions derived from 226 genes.\
The morphology of the sexual and asexual states for Magnaporthales and closely related species is shown (right panel). The supporting values for each node were estimated using 1000 bootstrap replicates. The strain number and genome size for each species are provided. The schematic drawings were prepared by co-author J. Luo.](srep09448-f1){#f1}

![Maximum likelihood tree of 24 Magnaporthales species and five Sordariomycetes used as outgroup species based on 82,715 amino acid positions derived from 226 genes.\
The infection sites and the morphology of sexual and asexual states for the three major Magnaporthales clades are illustrated (right panel). The supporting values for each node were estimated using 1000 bootstrap replicates. The strain number and genome size (if available) for each species are provided. \**Buergenerula spartinae*, *Bussabanomyces longisporus*, *Nakataea oryzae*, and *Omnidemptus affinis* are associated with above-ground parts of host plants, which is an exception in clade C. *Pyricularia oryzae* is associated with both the leaf and root of the host plant, which is an exception in clade B. Drawings were done by author J. Luo. The photos were prepared by co-authors J. Luo and N. Zhang.](srep09448-f2){#f2}

[^1]: These authors contributed equally to this work.
