# Duplications and losses in gene families of rust pathogens highlight putative effectors

Amanda L. Pendleton, Katherine E. Smith, Nicolas Feau, Francis Martin, Igor V. Grigoriev, Richard Hamelin, C. Dana Nelson, J. Gordon Burleigh, John M. Davis

► To cite this version:

Amanda L. Pendleton, Katherine E. Smith, Nicolas Feau, Francis Martin, Igor V. Grigoriev, et al. Duplications and losses in gene families of rust pathogens highlight putative effectors. Frontiers in Plant Science, 2014, 5, 13 p. 10.3389/fpls.2014.00299. hal-01268958

HAL Id: hal-01268958

https://hal.science/hal-01268958

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Duplications and losses in gene families of rust pathogens highlight putative effectors

Amanda L. Pendleton<fnref n="1" />, Katherine E. Smith<fnref n="2" />, Nicolas Feau<fnref n="3" />, Francis M. Martin<fnref n="4" />, Igor V. Grigoriev<fnref n="5" />, Richard Hamelin<fnref n="3" />, C. Dana Nelson<fnref n="2" />, J. Gordon Burleigh<fnref n="1" />,<fnref n="6" />,<fnref n="7" /> and John M. Davis<fnref n="1" />,<fnref n="7" />,<fnref n="8" />*

**Edited by:**

Sébastien Duplessis, INRA, France

**Reviewed by:**

Dag Ahren, Lund University,
Sweden
Yao-Cheng Lin, VIB/Ghent
University, Belgium

**\*Correspondence:**

John M. Davis, School of Forest
Resources and Conservation,
University of Florida, 365
Newins-Ziegler Hall, Gainesville,
FL 32611, USA
e-mail: jmdavis@ufl.edu

Rust fungi are a group of fungal pathogens that cause some of the world's most destructive diseases of trees and crops. A shared characteristic among rust fungi is obligate biotrophy, the inability to complete a lifecycle without a host. This dependence on a host species likely affects patterns of gene expansion, contraction, and innovation within rust pathogen genomes. The establishment of disease by biotrophic pathogens is reliant upon effector proteins that are encoded in the fungal genome and secreted from the pathogen into the host's cell apoplast or within the cells. This study uses a comparative genomic approach to elucidate putative effectors and determine their evolutionary histories. We used OrthoMCL to identify nearly 20,000 gene families in proteomes of 16 diverse fungal species, which include 15 basidiomycetes and one ascomycete. We inferred patterns of duplication and loss for each gene family and identified families with distinctive patterns of expansion/contraction associated with the evolution of rust fungal genomes. To recognize potential contributors for the unique features of rust pathogens, we identified families harboring secreted proteins that: (i) arose or expanded in rust pathogens relative to other fungi, or (ii) contracted or were lost in rust fungal genomes. While the origin of rust fungi appears to be associated with considerable gene loss, there are many gene duplications associated with each sampled rust fungal genome. We also highlight two putative effector gene families that have expanded in Cqf that we hypothesize have roles in pathogenicity.

**Keywords:** effectors, rust pathogens, secretome, genome evolution, comparative genomics

## INTRODUCTION

Rust fungi are plant infecting filamentous fungi in the order Pucciniales (Basidiomycota) that are unified by obligate biotrophy (Voegele and Mendgen, 2011). This form of pathogenicity requires a live host to establish a parasitic relationship. This is accomplished through the establishment of a molecularly intimate interaction at the host-pathogen interface characterized by the secretion of an arsenal of proteins from the pathogen that suppress host defense mechanisms and promote the acquisition of essential nutrients by the pathogen (Dodds et al., 2009; Stergiopoulos and de Wit, 2009). Such proteins, termed effectors, are thought to establish and maintain a compatible interaction between the pathogen and host. The processes that drive evolution of effector diversity are of great interest because pathogen's effector genes and host resistance genes are the interacting "gene-for-gene" pairs that drive coevolution in these pathosystems (Jones and Dangl, 2006; Stergiopoulos and de Wit, 2009).

Secreted proteins can be identified from whole genome sequences through the utilization of bioinformatic tools to isolate proteins with N-terminal secretion signals. Bioinformatic pipelines can then be used to narrow predicted secreted protein sets to putative effectors. These proteins contain features of known effectors such as elevated cysteine content (greater than 2%), that would enable the formation of stabilizing disulfide bridges (Stergiopoulos and de Wit, 2009), and protein domains associated with pathogenicity. Length is a criteria used to identify small secreted proteins (SSPs) from within putative effector protein sets, as SSPs are effector-like proteins with lengths less than 300 amino acids. Sequence comparisons alone do not provide a reliable means to identify putative effectors since some known effectors are lineage-specific while others are conserved across taxa (Rep, 2005; Saunders et al., 2012; Giraldo and Valent, 2013). Candidate effectors, a further distinction, are putative effectors that have additional support for roles in pathogenicity

<footnote n="1">¹ Plant Molecular and Cellular Biology Program, University of Florida, Gainesville, FL, USA</footnote>

<footnote n="2">² Southern Research Station, Southern Institute of Forest Genetics, USDA Forest Service, Saucier, MS, USA</footnote>

<footnote n="3">³ Department of Forest Sciences, University of British Columbia, Vancouver, BC, Canada</footnote>

<footnote n="4">⁴ Interactions Arbres/Microorganismes, Laboratoire d'Excellence ARBRE, INRA-Nancy, UMR Institut National de la Recherche Agronomique - Université de Lorraine, Champenoux, France</footnote>

<footnote n="5">⁵ US Department of Energy, Joint Genome Institute, Walnut Creek, CA, USA</footnote>

<footnote n="6">⁶ Biology Department, University of Florida, Gainesville, FL, USA</footnote>

<footnote n="7">⁷ Genetics Institute, University of Florida, Gainesville, FL, USA</footnote>

<footnote n="8">⁸ School of Forest Resources and Conservation, University of Florida, Gainesville, FL, USA</footnote>
(i.e., induced transcription or elevated expression in planta). Genetic evidence for functional redundancy of effectors, presumably due to multigene families of effector proteins, whose members share similar functions, has been reported in several pathogens (Kamper et al., 2006; Rafiqi et al., 2012; Saitoh et al., 2012; Giraldo and Valent, 2013). This suggests it would be useful to characterize families of proteins with effector-like characteristics so as to identify families that have expanded during evolution in association with the acquisition of pathogenic life history characteristics. Examining the evolutionary history of protein families across a set of diverse fungal taxa should help identify lineage-specific, putative effector protein families, families that may have evolved similar functions in more distantly related taxa, and families that may exhibit functional redundancy.

*Cronartium quercuum* f. sp. *fusiforme* (*Cqf*) is a rust pathogen that has a complex life cycle with five spore types and exhibits alternation between two hosts, oak (*Quercus* spp.) and southern pines (*Pinus* spp.). The fungus incites fusiform rust disease on southern pines, leading to significant economic losses to the forest products industry. The impact of the disease on pine production has motivated extensive research on the genetic interaction between *Cqf* and pine. The objective of this study is to identify putative effector gene families in the *Cqf* genome through comparative genomic analyses between *Cqf* and 15 other fungal taxa, including two other rust pathogens. We have identified families that have expanded in *Cqf* that we hypothesize are involved in conditioning stem gall phenotypes observed on the pine host. Our analyses provide a more thorough perspective on *Cqf* and rust pathogen evolution and also highlight the evolutionary patterns of putative effector families that *Cqf* employs to establish disease on two taxonomically diverse host species.

# MATERIALS AND METHODS

## GENE FAMILY CONSTRUCTION

Complete proteomes were downloaded from the public databases of the National Center for Biotechnology Information (www.ncbi.nlm.nih.gov/genome), U.S. Department of Energy's Joint Genome Institute (jgi.doe.gov/fungi), and the Broad Institute (www.broadinstitute.org). Sixteen proteomes were obtained: (Basidiomycota) *Cronartium quercuum* f.sp. *fusiforme* G11 version 1.0 (*Cqf*; unpublished, jgi.doe.gov/Cronartium), *Melamspora larici-populina* version 1.0 (*Mlp*; Duplessis et al., 2011a,b), *Puccinia graminis* f.sp. *tritici* CRL 75-36-700-3 race SCCL (*Pgt*; Duplessis et al., 2011a,b), *Mixia osmundae* IAM 14324 version 1.0 (*Mos*; Toome et al., 2014), *Sporobolomyces roseus* version 1.0 (*Sro*; with permission; jgi.doe.gov/fungi), *Rhodotorula graminis* strain WP1 version 1.1 (*Rgr*; with permission; jgi.doe.gov/fungi), *Ustilago maydis* strain 521 (*Uma*; Kamper et al., 2006), *Malasezzia globosa* CBS 7966 (*Mgl*; Xu et al., 2007), *Pisolithus tinctorius* Marx 270 version 1.0 (*Pti*; with permission; jgi.doe.gov/fungi), *Phanerochaete chrysosporium* version 2.0 (*Pch*; Martinez et al., 2004), *Heterobasidion irregulare* version 2.0 (*Hir*; Olson et al., 2012), *Serpula lacrymans* S7.3 version 2.0 (*Sla*; Eastwood et al., 2011), *Agaricus bisporus* var. *bisporus* H97 version 2.0 (*Abi*; Morin et al., 2012), *Laccaria bicolor* version 2.0 (*Lbi*; Martin et al., 2008), *Amanita muscaria* Koide version 1.0 (*Amu*; with permission; jgi.doe.gov/fungi), and (Ascomycota) *Saccharomyces cerevisiae* S288C (*Sce*; Goffeau

et al., 1996), for a total of 200,313 proteins. Gene families were delineated by OrthoMCL v.5.0 software (Li et al., 2003) using default parameters (minimum *e*-value of 1e-05, minimum similarity of 50%).

## SECRETOME PREDICTION

The collective set of secreted proteins, or the secretome, of *Cqf* was identified bioinformatically. Annotation of a secreted protein is determined by signal peptide (SignalP 3.0 and 4.0; Bendtsen et al., 2004; Petersen et al., 2011), protein localization (TargetP 1.1; Emanuelsson et al., 2000), and transmembrane domain (TMHMM 2.0; Krogh et al., 2001) bioinformatics prediction software (Feau et al. in prep.). Proteins predicted by TargetP 1.1 to be targeted for the mitochondrion (with RC values between 1 and 3) were discarded and residual proteins are submitted to TMHMM 2.0. If no TM-domain is identified in the protein, or a TM-domain is predicted in the N-terminal region of the protein (i.e., in the first 70 amino acids), the protein is re-oriented toward SignalP 4.0; in any other case, the protein is discarded. SignalP 4.0 either implements the SignalP-TM network to discriminate between a true signal peptide and an N-terminal trans-membrane region or the SignalP-noTM network if the program does not identify a TM-like domain in the N-terminal region of the protein. In this last case (i.e., if the SignalP-noTM network is implemented by SignalP 4.0), the protein is re-oriented toward SignalP 3.0 and a signal peptide prediction is positive if either both NN and HMM converged in a positive result or if NN D-score returns a positive result with a D-score ≥ 0.5.

## ESTIMATION OF GENE TREES

The protein sequences from each gene family were aligned using MUSCLE (Edgar, 2004). We assembled a collection of amino acid alignments from gene families with at least four sequences. For each of the gene family alignments, we performed a maximum likelihood (ML) search to find the optimal topology using RAxML v.7.2.8 with the PROTCATJTT model (Stamatakis et al., 2005). Gene tree estimates often contain much error and can be improved with knowledge of the underlying species tree (e.g., Rasmussen and Kellis, 2011). We constructed a species tree from a phylogenetic matrix of 2404 single copy genes with sequences from at least eight fungal taxa. We performed a ML search using RAxML v.7.2.8 with the PROTCATJTT model on the concatenated single gene matrix to estimate the species tree. For each of the gene trees, we used TreeFix version 1.1.8 (Wu et al., 2013) to improve on the ML topology given the species tree. TreeFix searches for a statistically equivalent rooted gene tree topology that minimizes the number of duplications and losses implied by the species tree. For 10 of the gene families, the TreeFix runs did not complete in 1 week. For these gene trees, we rooted the ML tree with a root that minimizes the number of implied duplications and losses using the program OptRoot (www.wehe.us). For all of the gene trees output from TreeFix or OptRoot, the locations of the implied duplications and losses were mapped on the species tree using URec version 1.02 (Gorecki and Tiuryn, 2007).

## FUNCTIONAL ANNOTATION OF PROTEINS

Functional annotations were obtained from the Joint Genome Institute's (JGI) Mycocosm (jgi.doe.gov/fungi; Grigoriev et al.,
2014) for the 16 organisms included in the phylogenetic and gene family analyses. Protein domains were identified using the online InterPro interface (http://www.ebi.ac.uk/ interpro/; Hunter et al., 2012). Transmembrane domain regions were identified in amino acid sequences of proteins using TMpred (www.ch.embnet.org/software/TMPRED_form.html; Hofmann, 1993). Glycosylphosphatidylinositol (GPI) anchor sites were predicted using big-PI Predictor (http://mendel.imp. ac.at/gpi/gpi_server.html; Eisenhaber et al., 1999).

# RESULTS

## GENE FAMILY ANALYSIS

The OrthoMCL analysis of the proteomes identified 19,489 gene families that contained 152,964 proteins. This protein count was ~76% (152,964/200,313) of the total proteins input into OrthoMCL analysis. Protein counts per gene family ranged from 2 (minimum size for a gene family) to 343 proteins, and the average family size was 7.8 proteins. Approximately 42% of the gene families had proteins encoded from only a single taxon, and families with proteins encoded in two or three taxa were the next most abundant families (Figure 1). Relatively few families contained proteins detected in 4–14 taxa, but more families contained proteins detected in 15 or 16 taxa (~12% of all families; 2,277/19,489). The families broadly conserved across all 16 sampled taxa are likely to contain core essential fungal proteins. The remaining ~24% of input proteins that did not group into families are considered true singletons, as

they lack homologs within their own proteome or in the other taxa.

To highlight gene families specific to the rust pathogen lineage, we compared gene family conservation between four pathogen genomes belonging to the subphylum Pucciniomycotina, which include three rust pathogens (Cqf, Mlp, and Pgt; Pucciniales) and a non-rust fern pathogen, Mixia osmundae (Mos; Mixiales). We selected the 4673 gene families containing proteins from at least one of these four pathogens (and no proteins from other sampled taxa) from the complete OrthoMCL family dataset. These families contained 22,784 proteins and exhibited varying patterns of conservation across the four taxa (Figure 2A). Most prominently, 14,978 of the 22,784 proteins (65.7%) were encoded in only one of the four pathogen genomes, illustrating high levels of species specificity (Figure 2A). Fewer proteins were shared between two or more rust fungi in this subset of families (7512/22,784 or 33.0%) (Figure 2A). Of the 19,485 families determined by OrthoMCL, 656 families (or 3.4%; Figure 2B) consisted of gene models found only in the three rust pathogen genomes, where each of the three rust pathogens had a representative gene model in the family. A total of 3466 proteins (Figure 2A) were ascribed to these rust pathogen-specific families. The largest family contained 249 proteins, and the smallest had 3 proteins. These 656 families represent the “core” rust pathogen protein set. Of the sampled genomes, the two pathogens with the most uniquely shared families are Cqf and Mlp, which have 878 conserved families.

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Gene families are predominantly species-specific in the sampled taxa. The proportions of gene families with proteins encoded in one through 16 fungal taxa genomes (taxa count) are displayed for the 19,489 OrthoMCL gene families.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**FIGURE 2 | Conserved proteins and families within only four**
**Pucciniomycete pathogen genomes are mostly species-specific.** Gene
family (OrthoMCL) conservation within Pucciniomycete pathogens; *Mixia*
*osmundae* (Mos), *Cronartium quercuum* f.sp. *fusiforme* (Cqf), *Melampsora*</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>*larici-populina* (Mlp), and *Puccinia graminis* f.sp. *tritici* (Pgt). The values indicate the total number of (A) gene models or (B) gene families conserved in only these four species and absent in the remaining 12 fungal taxa included in the OrthoMCL analysis.</figcaption></figure>

IDENTIFYING PUTATIVE EFFECTORS

We identified gene families encoding putative effectors in the *Cqf*
genome. To highlight putative effectors, the predicted secretome
 predicted secretome (predicted secreted proteins; see Methods) was analyzed for cys-
teine content and family-level conservation. The *Cqf* secretome
harbors 666 SSPs, which are secreted proteins with fewer than
300 amino acids (aa). The range in protein lengths within the
secretome was 51–1716 aa with a median length of 249 aa.

Analysis of Cqf gene families elucidated the evolutionary his-
tories of secreted putative effectors. To identify putative effector
families within the Cqf genome, we selected gene families with
at least two secreted proteins, as these families would then con-
tain at least two paralogous putative effectors and the family
would have therefore expanded in the Cqf genome. In total, 132
putative effector families were identified. Sixty-five of these fam-
ilies were conserved effector families, with proteins from two or
more fungal taxa. These families had sequences from 6.94 taxa
on average (Table 1) and represent potential effectors with func-
tions that can occur in a wide range of hosts. Alternatively, 67
novel effector families were considered to be evolutionary inno-
vations since the family members consisted of only Cqf proteins
(Table 2). The average family size for conserved effector families
(18.23 proteins) was significantly larger than Cqf-specific fami-
lies (3.54 proteins; t-test, p-value < 0.001). However, there was
no difference in the number of Cqf proteins per family in con-
served (mean = 5.02 proteins) and Cqf-specific families (mean
= 2.4 proteins). Families where all Cqf protein members are pre-
dicted to be secreted were found in both candidate effector family
types and at proportions that were not significantly different from
one another (conserved families = 40/65, Cqf-specific = 44/67;
Tables 1, 2). Evidence for potential sub- and/or neofunctionaliza-
tion was observed in 23 of the 67 (34.3%) Cqf-specific putative
effector families, as only a subset of proteins within these families
received secretion predictions, suggesting distinct biological roles
among family members.

GENE GAINS AND LOSSES

Gene gain and loss was quantified across all 16 sampled fungal
taxa. We mapped the gene trees from gene families with at least

four proteins onto a species tree to determine the patterns of
duplication and loss across the 16 fungal taxa. In total, we exam-
ined 10,371 gene trees containing 131,863 protein sequences.
These gene trees implied a minimum of 49,539 duplications (i.e.,
gene family gains) and 21,789 losses (i.e., gene family contrac-
tions and/or entire family loss). Over 93.9% of the duplications
and 67.9% of the losses are species-specific, occurring in a sin-
gle lineage at the tips of the species tree (Figure 3). The number
of species-specific duplications was positively correlated with the
size of a taxon’s proteome (R² = 0.93), suggesting that gene
duplication is a mechanism for proteomic expansion and diver-
sification for the selected fungal taxa (Figure 4). There was no
obvious relationship between proteome size and species-specific
duplication with life history forms (i.e., symbiotic, pathogenic, or
free-living) (Figure 4). Species-specific losses were not correlated
with the proteome size, but the rust pathogen lineage exhibited
fewer losses than other sampled taxa (Figure 5).

We identified genes that were gained and lost specifically in the rust pathogen clade. There were many gene losses (1217 events within 1148 families) associated origin of the rust pathogen clade within Pucciniomycotina (*Cqf*, *Mlp*, and *Pgt*) compared to the number of gains (248 events, 142 families) (site R in **Figure 3**). The number of taxa represented in these 1148 families range from 2 to 16 species, with the largest proportion of families (10.6%) having representatives from all 16 taxa in the analysis (**Figure 6**). Families lost genes at the origin of the rust fungi appear to occur in few of the sampled fungal lineages than those that had duplications in the rust fungi. Fifty percent of duplicated families contain proteins from 14 or more sampled taxa (**Figure 6**). Though the disproportionate level of gene losses prior to the common ancestor of rust pathogens is striking, each of the three rust fungal species shows evidence of high species-specific rates of duplication (**Figure 3**). In fact, 32.1% of all the duplications across the tree are specific to only one of the rust species (**Figure 3**).

The proteome of the *Cqf* rust pathogen is enriched for novel
proteins whose expansion has presumably contributed to spe-
cialization in its pathosystem. Numerous species-specific dupli-
cations have occurred following within the *Cqf* lineage (2730
duplication events in 549 families; **Figure 3**). Of the 549 families
<table><caption>Table 1 | Conserved putative effector families have broad and narrow taxonomic distributions.</caption><thead><tr><th>OrthoMCL group ID</th><th>Total proteins in family</th><th>Cqf proteins in family</th><th>Proteins in Cqf secretome</th><th>Number of taxa represented in family</th></tr></thead><tbody><tr><td>5485</td><td>12</td><td>7*</td><td>7</td><td>3</td></tr><tr><td>2168</td><td>17</td><td>6</td><td>5</td><td>3</td></tr><tr><td>2725</td><td>16</td><td>13</td><td>5</td><td>3</td></tr><tr><td>1053</td><td>94</td><td>6</td><td>4</td><td>13</td></tr><tr><td>2731</td><td>16</td><td>5</td><td>4</td><td>5</td></tr><tr><td>5853</td><td>11</td><td>5</td><td>4</td><td>5</td></tr><tr><td>6604</td><td>9</td><td>5</td><td>4</td><td>2</td></tr><tr><td>1101</td><td>65</td><td>62</td><td>3</td><td>2</td></tr><tr><td>1281</td><td>33</td><td>3*</td><td>3</td><td>3</td></tr><tr><td>1293</td><td>32</td><td>5</td><td>3</td><td>10</td></tr><tr><td>1397</td><td>27</td><td>3*</td><td>3</td><td>12</td></tr><tr><td>1831</td><td>19</td><td>3*</td><td>3</td><td>6</td></tr><tr><td>2730</td><td>16</td><td>3*</td><td>3</td><td>10</td></tr><tr><td>5067</td><td>13</td><td>4</td><td>3</td><td>7</td></tr><tr><td>6199</td><td>10</td><td>3*</td><td>3</td><td>3</td></tr><tr><td>6608</td><td>9</td><td>3*</td><td>3</td><td>3</td></tr><tr><td>6599</td><td>9</td><td>3*</td><td>3</td><td>5</td></tr><tr><td>7036</td><td>8</td><td>4</td><td>3</td><td>2</td></tr><tr><td>9412</td><td>5</td><td>3*</td><td>3</td><td>2</td></tr><tr><td>10,437</td><td>4</td><td>3*</td><td>3</td><td>2</td></tr><tr><td>1014</td><td>170</td><td>69</td><td>2</td><td>10</td></tr><tr><td>1219</td><td>39</td><td>3</td><td>2</td><td>14</td></tr><tr><td>1382</td><td>28</td><td>2*</td><td>2</td><td>13</td></tr><tr><td>1410</td><td>27</td><td>3</td><td>2</td><td>15</td></tr><tr><td>1444</td><td>26</td><td>2*</td><td>2</td><td>14</td></tr><tr><td>1428</td><td>26</td><td>2*</td><td>2</td><td>10</td></tr><tr><td>1507</td><td>24</td><td>3</td><td>2</td><td>15</td></tr><tr><td>1546</td><td>23</td><td>2*</td><td>2</td><td>13</td></tr><tr><td>1541</td><td>23</td><td>2*</td><td>2</td><td>16</td></tr><tr><td>1593</td><td>22</td><td>2*</td><td>2</td><td>13</td></tr><tr><td>1703</td><td>21</td><td>2*</td><td>2</td><td>8</td></tr><tr><td>1898</td><td>19</td><td>2*</td><td>2</td><td>14</td></tr><tr><td>1829</td><td>19</td><td>2*</td><td>2</td><td>15</td></tr><tr><td>1957</td><td>18</td><td>2*</td><td>2</td><td>14</td></tr><tr><td>2180</td><td>17</td><td>3</td><td>2</td><td>9</td></tr><tr><td>2191</td><td>17</td><td>2*</td><td>2</td><td>10</td></tr><tr><td>2172</td><td>17</td><td>2*</td><td>2</td><td>16</td></tr><tr><td>4359</td><td>15</td><td>2*</td><td>2</td><td>13</td></tr><tr><td>3833</td><td>15</td><td>2*</td><td>2</td><td>12</td></tr><tr><td>3825</td><td>15</td><td>2*</td><td>2</td><td>11</td></tr><tr><td>4560</td><td>14</td><td>2*</td><td>2</td><td>10</td></tr><tr><td>5075</td><td>13</td><td>2*</td><td>2</td><td>7</td></tr><tr><td>5492</td><td>12</td><td>3</td><td>2</td><td>6</td></tr><tr><td>5799</td><td>12</td><td>2*</td><td>2</td><td>10</td></tr><tr><td>6367</td><td>10</td><td>3</td><td>2</td><td>6</td></tr><tr><td>6213</td><td>10</td><td>3</td><td>2</td><td>3</td></tr><tr><td>6206</td><td>10</td><td>2*</td><td>2</td><td>3</td></tr><tr><td>6601</td><td>9</td><td>7</td><td>2</td><td>2</td></tr><tr><td>6606</td><td>9</td><td>3</td><td>2</td><td>7</td></tr><tr><td>6629</td><td>9</td><td>2*</td><td>2</td><td>6</td></tr></tbody></table>

(Continued)

<table><caption>Table 1 | Continued</caption><thead><tr><th>OrthoMCL group ID</th><th>Total proteins in family</th><th>Cqf proteins in family</th><th>Proteins in Cqf secretome</th><th>Number of taxa represented in family</th></tr></thead><tbody><tr><td>7040</td><td>8</td><td>2*</td><td>2</td><td>2</td></tr><tr><td>7590</td><td>7</td><td>2*</td><td>2</td><td>3</td></tr><tr><td>9416</td><td>5</td><td>4</td><td>2</td><td>2</td></tr><tr><td>9414</td><td>5</td><td>4</td><td>2</td><td>2</td></tr><tr><td>9843</td><td>5</td><td>2*</td><td>2</td><td>3</td></tr><tr><td>9462</td><td>5</td><td>2*</td><td>2</td><td>2</td></tr><tr><td>9446</td><td>5</td><td>2*</td><td>2</td><td>2</td></tr><tr><td>9431</td><td>5</td><td>2*</td><td>2</td><td>4</td></tr><tr><td>9424</td><td>5</td><td>2*</td><td>2</td><td>3</td></tr><tr><td>10,479</td><td>4</td><td>3</td><td>2</td><td>2</td></tr><tr><td>10,474</td><td>4</td><td>3</td><td>2</td><td>2</td></tr><tr><td>10,432</td><td>4</td><td>3</td><td>2</td><td>2</td></tr><tr><td>11,968</td><td>3</td><td>2*</td><td>2</td><td>2</td></tr><tr><td>11,958</td><td>3</td><td>2*</td><td>2</td><td>2</td></tr><tr><td>11,908</td><td>3</td><td>2*</td><td>2</td><td>2</td></tr><tr><td>Total</td><td>1185 proteins</td><td>326 proteins</td><td>162 proteins</td><td>-</td></tr><tr><td>Average</td><td>18.23 proteins per family</td><td>5.02 Cqf proteins per family</td><td>2.5 secreted proteins per family</td><td>6.94 taxa per family</td></tr></tbody></table>

Gene families with greater than two Cqf predicted secreted proteins are listed. Data is ranked by the number of Cqf secreted proteins. The total number of proteins in each family is provided as well as the number of proteins belonging to the Cqf secretome (i.e., predicted secreted proteins). Asterisks adjacent to total protein counts indicate families where all members are Cqf secretome members. If no asterisk is present, only a portion of the family received secretion predictions. Family 5485 (bold) will be detailed later in article.

that have undergone Cqf-specific duplications, 248 (or 45.17%) contain proteins not observed in any other analyzed fungal taxa. These 248 novel families comprise 14.5% of the annotated Cqf proteome (2017/13,903 proteins), highlighting the rapid expansion of novel, likely pathogenicity-related gene families. The vast majority (98.8%) of these novel families do not have BLASTp hits in the NCBI non-redundant database or have hits to unknown proteins (minimum e-value of 1e-10; Table 3) and 94.5% do not contain InterPro domains (unpublished, jgi.doe.gov/Cronartium; Hunter et al., 2012). Since the families that are unique to the Cqf lineage are largely uncharacterized, they likely follow the assumptions for putative pathogenicity factors or effectors. Nearly 12% (234/2,017) of proteins encoded in the 248 novel Cqf families are members of the predicted Cqf secretome. This is significantly greater than the ~8% of entire Cqf proteome that also belongs to the secretome (Chi-square = 25.418, p-value = 0.0001). The protein characteristics of these secreted proteins are effector-like, as the average cysteine content is 2.2% and the median protein length is 272 amino acids.

The families that were duplicated in the Cqf lineage and contain sequences from other taxa exhibit patterns of conservation that differ from the families duplicated or depleted in rust
<table><caption>Table 2 | Potential sub- and neo-functionalization within Cqf-specific putative effector families.</caption><thead><tr><th>Gene family ID</th><th>Cqf proteins in family</th><th>Proteins in Cqf secretome within family</th></tr></thead><tbody><tr><td>8514</td><td>6*</td><td>6</td></tr><tr><td><strong>9417</strong></td><td>5*</td><td>5</td></tr><tr><td>9892</td><td>5*</td><td>5</td></tr><tr><td>2663</td><td>17</td><td>4</td></tr><tr><td>6030</td><td>11</td><td>4</td></tr><tr><td>9436</td><td>5</td><td>4</td></tr><tr><td>9897</td><td>5</td><td>4</td></tr><tr><td>7037</td><td>8</td><td>3</td></tr><tr><td>9890</td><td>5</td><td>3</td></tr><tr><td>9891</td><td>5</td><td>3</td></tr><tr><td>10,467</td><td>4</td><td>3</td></tr><tr><td>11,111</td><td>4</td><td>3</td></tr><tr><td>11,876</td><td>3*</td><td>3</td></tr><tr><td>11,879</td><td>3*</td><td>3</td></tr><tr><td>11,934</td><td>3*</td><td>3</td></tr><tr><td>12,788</td><td>3*</td><td>3</td></tr><tr><td>12,794</td><td>3*</td><td>3</td></tr><tr><td>12,804</td><td>3*</td><td>3</td></tr><tr><td>12,812</td><td>3*</td><td>3</td></tr><tr><td>5202</td><td>13</td><td>2</td></tr><tr><td>7921</td><td>7</td><td>2</td></tr><tr><td>7928</td><td>7</td><td>2</td></tr><tr><td>8502</td><td>6</td><td>2</td></tr><tr><td>9880</td><td>5</td><td>2</td></tr><tr><td>10,433</td><td>4</td><td>2</td></tr><tr><td>10,463</td><td>4</td><td>2</td></tr><tr><td>11,112</td><td>4</td><td>2</td></tr><tr><td>11,928</td><td>3</td><td>2</td></tr><tr><td>11,944</td><td>3</td><td>2</td></tr><tr><td>12,777</td><td>3</td><td>2</td></tr><tr><td>12,803</td><td>3</td><td>2</td></tr><tr><td>12,807</td><td>3</td><td>2</td></tr><tr><td>12,814</td><td>3</td><td>2</td></tr><tr><td>14,526</td><td>2*</td><td>2</td></tr><tr><td>14,527</td><td>2*</td><td>2</td></tr><tr><td>14,537</td><td>2*</td><td>2</td></tr><tr><td>14,552</td><td>2*</td><td>2</td></tr><tr><td>14,554</td><td>2*</td><td>2</td></tr><tr><td>14,563</td><td>2*</td><td>2</td></tr><tr><td>14,570</td><td>2*</td><td>2</td></tr><tr><td>14,577</td><td>2*</td><td>2</td></tr><tr><td>14,589</td><td>2*</td><td>2</td></tr><tr><td>14,623</td><td>2*</td><td>2</td></tr><tr><td>15,977</td><td>2*</td><td>2</td></tr><tr><td>15,979</td><td>2*</td><td>2</td></tr><tr><td>15,980</td><td>2*</td><td>2</td></tr><tr><td>15,992</td><td>2*</td><td>2</td></tr><tr><td>16,012</td><td>2*</td><td>2</td></tr><tr><td>16,034</td><td>2*</td><td>2</td></tr><tr><td>16,036</td><td>2*</td><td>2</td></tr><tr><td>16,051</td><td>2*</td><td>2</td></tr></tbody></table>

(Continued)

<table><caption>Table 2 | Continued</caption><thead><tr><th>Gene family ID</th><th>Cqf proteins in family</th><th>Proteins in Cqf secretome within family</th></tr></thead><tbody><tr><td>16,052</td><td>2*</td><td>2</td></tr><tr><td>16,078</td><td>2*</td><td>2</td></tr><tr><td>16,079</td><td>2*</td><td>2</td></tr><tr><td>16,080</td><td>2*</td><td>2</td></tr><tr><td>16,081</td><td>2*</td><td>2</td></tr><tr><td>16,091</td><td>2*</td><td>2</td></tr><tr><td>16,101</td><td>2*</td><td>2</td></tr><tr><td>16,102</td><td>2*</td><td>2</td></tr><tr><td>16,106</td><td>2*</td><td>2</td></tr><tr><td>16,119</td><td>2*</td><td>2</td></tr><tr><td>16,123</td><td>2*</td><td>2</td></tr><tr><td>16,129</td><td>2*</td><td>2</td></tr><tr><td>16,146</td><td>2*</td><td>2</td></tr><tr><td>16,160</td><td>2*</td><td>2</td></tr><tr><td>16,163</td><td>2*</td><td>2</td></tr><tr><td>16,191</td><td>2*</td><td>2</td></tr><tr><td>Total</td><td>237 proteins</td><td>164 proteins</td></tr><tr><td>Average</td><td>3.54 Cqf proteins per family</td><td>2.4 secreted proteins per family</td></tr></tbody></table>

Cqf-specific gene families with greater than two predicted secreted proteins are listed. The total number of proteins in each Cqf-specific family is provided as well as the number of proteins belonging to the Cqf secretome (i.e., predicted secreted proteins) are indicated. Asterisks adjacent to total protein counts indicate families where all members are Cqf secretome members. If no asterisk is present, only a portion of the family received secretion predictions. Family 9417 (bold) will be detailed later in article.

pathogens (site R, **Figure 3**). Instead, these Cqf-specific duplicated families ($n = 549$ families) are predominantly conserved in not only Cqf, but also 2-3 taxa (**Figure 6**).

### Cqf PUTATIVE EFFECTOR GENE FAMILIES—DISTRIBUTION AND EXPANSION

Family 5485 is the largest family represented in the predicted Cqf secretome. The family contains 12 orthologous proteins (7 Cqf, 4 Mlp, and 1 Pgt proteins). Eleven of the 12 proteins have predicted N-terminal signal peptides (SignalP 4.0; Petersen et al., 2011), and all seven members from Cqf are annotated as belonging to the Cqf secretome. Domain architecture and conservation data for Family 5485 proteins helps to predict their biological functions and putative roles in establishing infection. Additionally, 11 of the 12 proteins in this family contained three multicopper oxidase (MCO) domains and the remaining protein (Pgt_20719) contained two of the three domains (**Figure 7A**). The Interpro domains identified include: Cupredoxin domain (IPR008972), Multicopper Oxidase, Type 1 (IPR001117), and Multicopper Oxidase, Type 2 (IPR011706), and Multicopper Oxidase, Type 3 domain (IPR011707) (**Figure 7A**). A Copper-Binding Site domain (IPR002355) was identified in only three Cqf family members. These three proteins have a distinct phylogenetic history from other family members (**Figure 7B**). Generally, the
<figure><img src="image_6.png" /></figure>

**FIGURE 3** | Gene loss and gain in the Basidiomycota fungal lineage highlights shared loss in the rust pathogen lineage and high levels of species-specific gain. Mapping putative gene duplications and losses across 16 fungal taxa. Values in blue are associated with gains/duplications, whereas orange indicates loss.

Outside of parentheses are the number of gain or loss events that have occurred on the branch preceding a node, and within parentheses are the number of gene families associated with duplications or losses. The node denoted with R indicates the last common ancestor of the rust pathogens.

phylogenetic relationships, as well as the genomic colocalization of the proteins in this family mirrors the domain architecture, providing insight into how these proteins evolved (Figure 7A). Several additional families of MCOs are present in the Cqf

genome (i.e., Families 5853, 1542, and 1053), however, by definition, Family 5485 has a distinct evolutionary history from other families as evidenced by distinct family placement by OrthoMCL.
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

Family 9417 is the third largest family in the *Cqf* secretome, with all five of its proteins predicted as secreted. This family contains putative effectors likely involved in the establishment of disease, as all members have signal peptides, short lengths (average 207 aa), and high cysteine content (average 6.5%). All five family members contain at least one fungal extracellular membrane (CFEM) domain (Interpro IPR008427). Five additional proteins encoded in the *Cqf* proteome contain CFEM domains. Two of the five do not belong to a gene family, and the remaining three proteins each were ascribed to different families containing orthologs from multiple fungal taxa, unlike *Cqf*-specific family 9417. Similar to Family 5485, proteins of Family 9417 also colocalize in the genome, as three members are located on scaffold 43 of the *Cqf* assembly and the remaining two proteins are adjacent to one another on scaffold 5 (Figure 8). Protein members of Family 9417 adhere to consensus domain structure and subcellular targeting of previously identified CFEM proteins. Online prediction algorithms detected transmembrane domain regions (Tmpred; Hofmann, 1993) and glycosylphosphatidylinositol (GPI) anchor sites (big-Pi Predictor; Nielsen et al., 1997) in a subset of the family proteins (Figure 8). All proteins, excluding Cqf91696, were predicted to have N-terminal transmembrane helices spanning amino acids 3-23 for both proteins. Two members, Cqf712797 and Cqf651034 had C-terminal GPI anchor sites at amino acids 223 and 302, respectively (p-values 1.25E-04 and 2.10E-04). Only Cqf91696 had no bioinformatic evidence of association with the fungal membrane.

DISCUSSION

This study provides the first detailed analysis of the secretome
of the fusiform rust pathogen, *Cqf*, since the recent assembly

and annotation of a draft reference genome (unpublished,
rgi.doe.gov/Cronartium). Additional criteria used in isolating
putative effectors from within the *Cqf* genome and its corre-
sponding secretome, included proteins exhibiting rust pathogen-
specific and *Cqf*-specific gene family membership. Following
gene family constructions, we highlighted putative effectors with
paralogs (within *Cqf*) or orthologs/paralogs (between *Cqf* and
other taxa) within the *Cqf* secretome. Over half (51%) of pro-
teins considered to be effector-like (small, cysteine-rich, secreted
proteins) belong to gene families. This is comparable to results
found in the hemibiotrophic pathogens *Phytophthora ramorum*
and *P. sojae* where 77% of their secretomes are found in multigene
families (Tyler et al., 2006). These findings demonstrate the value
of an evolutionary perspective for highlighting families harbor-
ing putative *Cqf* effectors. Altogether, the large-scale comparative
genomics analyses in this study help elucidate the unique patterns
of evolution in a rust proteome and its associated secretome.

PUTATIVE EFFECTOR FAMILIES

With the completion of the *Cqf* draft genome, it is important to identify proteins that may be involved in establishing disease, such as effectors, on oak and pine hosts. Based on the evolutionary forces presumed to act on effectors, in combination with a trio of rust pathogen genomes facilitating comparative analyses, we can now do experiments not previously feasible. We suggest this is a reasonable approach to identifying putative effectors that complements more conventional methods. Previous studies searched for effectors within other systems based on the presence of a signal peptide, cysteine richness, and short protein lengths (<300 aa) (Joly et al., 2010; Cantu et al., 2011; Duplessis et al., 2011b;
<figure><img src="image_9.png" /><figcaption>**FIGURE 6** | Families gained and lost in rust pathogens and *Cqf* have varying levels of taxonomic conservation. The proportion of gene families (y-axis) that contain protein members from 1 to 16 fungal taxa (x-axis) among those that have either expanded or contracted in the rust fungi or *Cqf*.</figcaption></figure>

<table><caption>**Table 3** | *Cqf*-specific duplicated gene families contain predominantly uncharacterized proteins.</caption><thead><tr><td>Gene family annotation</td><td>Number of gene families (proteins)</td></tr></thead><tbody><tr><td>No hits</td><td>183 (1510)</td></tr><tr><td>Unknown protein</td><td>57 (430)</td></tr><tr><td>20S Proteasome subunit alpha 6</td><td>2 (9)</td></tr><tr><td>Zinc finger CCHC-type protein</td><td>1 (33)</td></tr><tr><td>HIV-1 retropepsin, polyprotein</td><td>1 (13)</td></tr><tr><td>Polysaccharide lyase family 4</td><td>1 (7)</td></tr><tr><td>Reverse transcriptase</td><td>1 (5)</td></tr><tr><td>MFS transporter, inorganic phosphate transporter</td><td>1 (5)</td></tr><tr><td>CFEM domain containing protein</td><td>1 (5)</td></tr></tbody></table>

Functional annotation of the 248 gene families duplicated only in *Cqf* by BLASTp against the non-redundant NCBI database (minimum e-values of 1e-10). A family was ascribed a function if more than two proteins in the family received the same top annotated BLASTp hit. The number of proteins within families is indicated in parentheses.

Hacquard et al., 2012; Saunders et al., 2012). This study highlights the usefulness of comparative genomic analyses to examine the evolutionary history of each secretome member, and that this approach can also be complemented with structural characteristics of predicted secreted proteins. The rationale behind these comparisons is that effector families conserved in rust fungi and unique to *Cqf* are candidates for conditioning rust pathogen and *Cqf* infection strategies, respectively.

We observed species-specific proteomic gene family gains/duplications in the *Cqf* lineage, a subset of which represents putative effectors. The paralogous nature (i.e., multi-copy) of their protein family members indicates functional redundancy,

which is consistent with other pathogenic fungi (Kamper et al., 2006; Saitoh et al., 2012). We have identified two lines of evidence that point toward neo- and sub-functionalization in *Cqf* putative effector families. First, differential subcellular localization predictions have been observed within putative effector families. In about 34% of *Cqf*-specific families, only a subset of proteins are secreted from the fungal cell, while remaining family members are not predicted for secretion, thus remaining within the fungal cell. This pattern suggests that secreted proteins with effector function may have evolved from non-secreted proteins without an effector function or vice versa. Second, changes in domain architecture of proteins within putative effector families also points to neo- or subfunctionalization. For example family 5485 contains MCO laccase-like enzymes and a single clade of three *Cqf* proteins that have acquired a MCO copper binding site in the evolution of this family. It is possible that these proteins have novel or distinct functions within *Cqf* than their paralogs within the genome. This family is a strong putative effector family because all *Cqf* members belong to the predicted secretome and it has undergone *Cqf*-specific family duplications. Protein members within this family co-localize in the genome, possibly resulting from tandem duplication from non-equal crossing over. Various functions have been ascribed to previously identified fungal MCOs including lignin degradation (Leonowicz et al., 2001; Lundell et al., 2010), melanin synthesis (Langfelder et al., 2003), fruiting body formation (Kues and Liu, 2000), and pathogenicity on hosts (Zhu and Williamson, 2004). This family has expanded in *Cqf*, the first sequenced rust pathogen that forms stem galls in woody tissues, and we hypothesize that these enzymes play a role in gall formation. The most common function for laccases/MCOs in basidiomycete fungi is lignin metabolism (Thurston, 1994; Kües and Rühl, 2011). However, this gene family exhibits a lack of conservation with known
<figure><img src="image_10.png" /></figure>

**FIGURE 7** | Domain structure, phylogenetic relationships, and colocalization of copper oxidases. (A) Phylogenetic relationship estimated with TreeFix (Wu et al., 2013) for the 12 members of Family 5485. Tree reconstructed with Tree of Life viewer (Letunic and Bork, 2011). See methods for additional details for gene tree estimation. Branch lengths are not informative for this tree. Domain architectures are indicated for each protein: N-terminal signal peptides per SignalP 4.0 (blue; Petersen et al., 2011); (MCu1—IPR01117) Multicopper oxidase Type 1 (green); (MCu2—IPR011706) Multicopper oxidase Type 2 (orange); (MCu3—IPR011707) Multicopper oxidase Type 3 (purple); and (CuBS—IPR002355) Multicopper oxidase, copper-binding site as determined by InterPro (red; Hunter et al., 2012). The colocalization of proteins on scaffolds is indicated by branches sharing the identical colors. Thin black branches do not colocalize. (B) Co-localization of proteins within MCO Family 5485 on scaffolds 7, 91, and 114. Family 5485 members are denoted with the ID above each gene model. Secretome members are blue arrows and non-secreted proteins are orange arrows. Gene orientation on the scaffold is indicated with arrows. Note: gene lengths are not to scale.

<figure><img src="image_11.png" /></figure>

**FIGURE 8** | Family 9417 proteins harbor both signal peptides and CFEM domains and colocalize in the Cqf genome. Domain architecture of Family 9417 protein members where signal peptides (blue), CFEM Interpro domains (green), GPI anchor sites (orange), and predicted transmembrane domains (red) are indicated for each protein. Three proteins colocalize together on scaffold 43, the remaining on scaffold 5 (indicated on right).

MCOs of lignin-degrading wood rots (*P. chrysosporium* and *S. lacrymans*), which points to the possibility that these enzymes may be involved in pathogenicity or may metabolize a plant substrate other than lignin. On both hosts, *Cqf* infects primary tissue that lacks high levels of lignification such as spongy mesophyll cells of oak leaves (Mims et al., 1996) and vascular cambium of pine (Gray et al., 1982). If the Family 5485 enzymes are involved in lignin degradation, the enzymatic activity may occur late in gall development on the pine host, where the tissues are more heavily lignified due to secondary wall formation. Though their biochemical targets are unknown *in planta*, we hypothesize that Family 5485 enzymes are secreted during infection and condition the gall phenotype on the pine host. Further studies are required to elucidate their true role in disease.

A second gene family that has expanded in the *Cqf* lineage is Family 9417, which includes five *Cqf*-specific paralogs that co-localize in the genome. Similar to Family 5485, differential domain architecture within this family implies that neo- or subfunctionalization may have occurred. Family 9417 contains putative effectors that harbor conserved, fungal-specific CFEM-domains. These domains exhibit a characteristic cysteine distribution and have a broad taxonomic conservation in fungi (Kulkarni et al., 2003; Martin et al., 2008; Perez et al., 2011). Predicted functions of proteins harboring CFEM domains include critical roles in appressorial development (Choi and Dean, 1997; DeZwaan et al., 1999), signal transducers, adhesion and cell-surface receptors (Kulkarni et al., 2003). In contrast to Family 5485 proteins, which may interact with the host during infection, the molecular target for Family 9417 proteins could be fungal. We hypothesize these proteins are secreted and may play roles during infection of the host.

## EVOLUTION OF GENE GAIN AND LOSS

Patterns of gene family loss and gain for rust fungi highlight major shifts in their proteomes, possibly associated with the rust pathogen's obligate biotrophic lifestyle. The origin of the rust pathogen clade is associated with nearly five times more losses, or family contractions, than duplications. There are many possible mechanisms for gene loss in rust fungi. For this reason, further investigations are required to both identify specific mechanisms and quantify their levels of effects on gene family evolution in rust fungi. However, we hypothesize that the since obligate biotrophy has evolved multiple times in fungi (Spanu, 2012), the skew toward gene loss in the rust pathogen lineage might be associated with the shift from the life history of its ancestral state to that of the obligate biotrophic pathogens we observe today. These lost and/or contracted families exhibit broad taxonomic
conservation and may have been constituents of the ancestral
"core" fungal gene set, suggesting that they are unnecessary for
obligate biotrophic but may be necessary for free-living and sym-
biotic species. For example, enzymes integral to the sulfur and
nitrogen assimilation pathways are missing in *Cqf* (unpublished,
jgi.doe.gov/Cronartium), *Mlp*, and *Pgt* (Duplessis et al., 2011b).
This also suggests that evolution for obligate biotrophy drives
toward an irreversible life history shift (Spanu, 2012).

Although the rusts have undergone considerable gene fam-
ily losses and contractions, they exhibit some of the largest
proteomes in fungi. Much of their proteome size appears
to be due to species-specific duplications. Nearly one-third
(32.1%) of all observed duplications across all of the sam-
pled basidiomycete fungi are rust taxon-specific duplications.
The high levels of species-specific duplication yield disproport-
tionately greater numbers of newly-evolved genes in the rust
pathogen genomes compared to ancient or conserved genes
(genes shared with older lineages) in each proteome. The pres-
ence of so many species-specific duplications suggests that
the rusts have highly labile genomes. This is consistent with
the large (>10%) genomic size variation detected in progeny
from a single *Cqf* cross relative to parental isolates (Anderson
et al., 2010). Such rapid changes, occurring in the span of
a single generation, could facilitate the gene gains and losses
observed in our analyses. The close association with hosts
may foster a labile and diverse genome, enabling the para-
sites to rapidly adapt to the continually evolving host resistance
pathways.

COMPARATIVE ANALYSIS AND GENETIC MAPPING TO VALIDATE
PUTATIVE EFFECTORS

Further characterization of putative effectors in *Cqf* could be
accomplished with analysis of selection potentially arising from
host resistance mechanisms (Allen et al., 2004; Aguileta et al.,
2009; Barrett et al., 2009; Thrall et al., 2012). In addition, expres-
sion analysis can be informative, since secreted proteins with
specific expression profiles during infection are stronger effec-
tor candidates (Ellis et al., 2009). Time-course experiments have
been successful in other rust pathogen systems to elucidate the
effector-like proteins involved in multiple or highly specific stages
during infection (Joly et al., 2010; Duplessis et al., 2011a; Bruce
et al., 2014). Also, resequencing of closely related rust pathogens
such as *Cronartium ribicola*, *C. flaccidum*, and *Peridermium hark-
nessii* (Vogler and Bruns, 1998) would improve precision of gene
family delineations and identification of true singleton *Cqf* effec-
tors, which are likely to be more newly evolved than effectors in
families, and may therefore be products of highly-specific host-
*Cqf* coevolution. Finally, a subset of the predicted effectors are
avirulence proteins and are, by definition, involved in genotype-
specific "gene-for-gene" interactions with hosts. These putative
avirulence effectors can be validated through genetic mapping
to their corresponding host resistance genes, an approach that
has previously been successful in identifying the first avirulence
protein locus in *Cqf* (Kubisiak et al., 2011). Altogether, these val-
idation approaches will yield true members of the *Cqf* secretome
and provide additional insight into the biological functions for
effectors infecting oak and pine.

ACKNOWLEDGMENTS

We acknowledge funding from the USDA Forest Service, Southern Research Station (Agreement 11-CA-11330126-120) (to John M. Davis). This work was partly supported by the French National Research Agency through the Laboratory of Excellence ARBRE (ANR-11-LABX-0002-01) (to Francis M. Martin). The work conducted by the U.S. Department of Energy Joint Genome Institute is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. We thank the following fungal genome project PIs for pre-publication access to genome sequence information: Drs. Sharon Doty (*Rhodotorula graminis*) and Ken Wolfe (*Sporobolomyces roseus*). The *Cqf* genome assembly and annotations can be interactively accessed through the JGI fungal genome portal MycoCosm at jgi.doe.gov/Cronartium (unpublished). We acknowledge Alicia Clym and Andrea Aerts at JGI for their roles in assembling and annotating the *Cqf* reference genome, respectively. Dr. Annegret Kohler at INRA, as well as Dr. Alan Kuo at JGI, are acknowledged for their leading roles in the acquisition of mycorrhizal genome sequence data.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online
at: http://www.frontiersin.org/journal/10.3389/fpls.2014.00299/
abstract

REFERENCES

Aguileta, G., Hood, M., Refregier, G., Giraud, T., Kader, J., and Delseny, M. (2009).
Genome evolution in plant pathogenic and symbiotic fungi. *Adv. Bot. Res.* 49,
151–193. doi: 10.1016/S0065-2296(08)00603-4

Allen, R., Bittner-Eddy, P., Grenvitte-Briggs, L., Meitz, J., Rehmany, A., Rose, L., et al. (2004). Host-parasite coevolutionary conflict between Arabidopsis and downy mildew. *Science* 306, 457, 1957–1960. doi: 10.1126/science.1104022

Anderson, C., Kubisiak, T., Nelson, C., Smith, J., and Davis, J. (2010). Genome size variation in the pine fusiform rust pathogen *Cronartium quercuum* f.sp. fusiforme as determined by flow cytometry. *Mycologia* 102, 1295–1302. doi: 10.3852/10-040

Barrett, L., Thrall, P., Dodds, P., van der Merwe, M., Linde, C., Lawrence, G., et al. (2009). Diversity and evolution of effector loci in natural populations of the plant pathogen *Melampsora lini*. *Mol. Biol. Evol.* 26, 2499–2513. doi: 10.1093/molbev/msp166

Bendtsen, J., Nielsen, H., Von Heijne, G., and Brunak, S. (2004). Improved prediction of signal peptides: SignalP 3.0. *J. Mol. Biol.* 340, 783–795. doi: 10.1016/j.jmb.2004.05.028

Bruce, M., Neugebauer, K. A., Joly, D. L., Migeon, P., Cuomo, C. A., Wang, S., et al. (2014). Using transcription of six *Puccinia triticina* races to identify the effective secretome during infection of wheat. *Front. Plant Sci.* 4:520. doi: 10.3389/fpls.2013.00520

Cantu, D., Govindarajulu, M., Kozik, A., Wang, M., Chen, X., Kojima, K., et al. (2011). Next Generation sequencing provides rapid access to the genome of *Puccinia striiformis* f. sp *tritici*, the causal agent of wheat stripe rust. *PLoS ONE* 6:e24230. doi: 10.1371/journal.pone.0024230

Choi, W., and Dean, R. (1997). The adenylate cyclase gene MAC1 of *Magnaporthe grisea* controls appressorium formation and other aspects of growth and development. *Plant Cell* 9, 1973–1983. doi: 10.1105/tpc.9.11.1973

DeZwaan, T., Carroll, A., Valent, B., and Sweigard, J. (1999). *Magnaporthe grisea* *Pth11p* is a novel plasma membrane protein that mediates appressorium differentiation in response to inductive substrate cues. *Plant Cell* 11, 2013–2030. doi: 10.1105/tpc.11.10.2013

Dodds, P., Rafiqi, M., Gan, P., Hardham, A., Jones, D., and Ellis, J. (2009). Effectors of biotrophic fungi and oomycetes: pathogenicity factors and triggers of host resistance. *New Phytol.* 183, 993–999. doi: 10.1111/j.1469-8137.2009.02922.x
Duplessis, S., Cuomo, C. A., Lin, Y. C., Aerts, A., Tisserant, E., Veneault-Fourrey, C., et al. (2011b). Obligate biotrophy features unraveled by the genomic analysis of rust fungi. *Proc. Natl. Acad. Sci. U.S.A.* 108, 9166–9171. doi: 10.1073/pnas.1019315108

Duplessis, S., Hacquard, S., Delaruelle, C., Tisserant, E., Frey, P., Martin, F., et al. (2011a). *Melampsora larici-populina* transcript profiling during germination and timecourse infection of poplar leaves reveals dynamic expression patterns associated with virulence and biotrophy. *Mol. Plant Microbe Interact.* 24, 808–818. doi: 10.1094/MPMI-01-11-0006

Eastwood, D. C., Floudas, D., Binder, M., Majcherczyk, A., Schneider, P., Aerts, A., et al. (2011). The plant cell wall–decomposing machinery underlies the functional diversity of forest fungi. *Science* 333, 762–765. doi: 10.1126/science.1205411

Edgar, R. (2004). MUSCLE: a multiple sequence alignment method with reduced time and space complexity. *BMC Bioinformatics* 5:113. doi: 10.1186/1471-2105-5-113

Eisenhaber, B., Bork, P., and Eisenhaber, F. (1999). Prediction of potential GPI-modification sites in proprotein sequences. *J. Mol. Biol.* 292, 741–758.

Ellis, J., Rafiqi, M., Gan, P., Chakrabarti, A., and Dodds, P. (2009). Recent progress in discovery and functional analysis of effector proteins of fungal and oomycete plant pathogens. *Curr. Opin. Plant Biol.* 12, 399–405. doi: 10.1016/j.pbi.2009.05.004

Emanuelsson, O., Nielsen, H., Brunak, S., and Von Heijne, G. (2000). Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. *J. Mol. Biol.* 300, 1005–1016. doi: 10.1006/jmbi.2000.3903

Giraldo, M., and Valent, B. (2013). Filamentous plant pathogen effectors in action. *Nat. Rev. Microbiol.* 11, 800–814. doi: 10.1038/nrmicrob3119

Goffeau, A., Barrell, B., Bussey, H., Davis, R., Dujon, B., Feldmann, H., et al. (1996). Life with 6000 genes. *Science* 274, 546–567.

Gorecki, P., and Tiuryn, J. (2007). URec: a system for unrooted reconciliation. *Bioinformatics* 23, 511–512. doi: 10.1093/bioinformatics/btl634

Gray, D., Amerson, H., and Van Dyke, C. (1982). An ultrastructural comparison of monokaryotic and dikaryotic haustoria formed by the fusiform rust fungus *Cronartium-quercuum* f-sp *fusiforme*. *Can. J. Bot.* 60, 2914–2922. doi: 10.1139/b82-352

Grigoriev, I., Nikitin, R., Haridas, S., Kuo, A., Ohm, R., Otillar, R., et al. (2014). MycoCosm portal: gearing up for 1000 fungal genomes. *Nucleic Acids Res.* 42, D699–D704. doi: 10.1093/nar/gkt1183

Hacquard, S., Joly, D., Lin, Y., Tisserant, E., Feau, N., Delaruelle, C., et al. (2012). A comprehensive analysis of genes encoding small secreted proteins identifies candidate effectors in *Melampsora larici-populina* (Poplar Leaf Rust). *Mol. Plant Microbe Interact.* 25, 279–293. doi: 10.1094/MPMI-09-11-0238

Hofmann, K. (1993). TMbase-A database of membrane spanning protein segments. *Biol. Chem.* 374, 166.

Hunter, S., Jones, P., Mitchell, A., Apweiler, R., Attwood, T., Bateman, A., et al. (2012). InterPro in 2011: new developments in the family and domain prediction database. *Nucleic Acids Res.* 40, D306–D312. doi: 10.1093/nar/gkr948

Joly, D., Feau, N., Tanguay, P., and Hamelin, R. (2010). Comparative analysis of secreted protein evolution using expressed sequence tags from four poplar leaf rusts (*Melampsora* spp.). *BMC Genomics* 11: 422. doi: 10.1186/1471-2164-11-422

Jones, J., and Dangl, J. (2006). The plant immune system. *Nature* 444, 323–329. doi: 10.1038/nature05286

Kamper, J., Kahmann, R., Bolker, M., Ma, L., Brefort, T., Saville, B., et al. (2006). Insights from the genome of the biotrophic fungal plant pathogen *Ustilago maydis*. *Nature* 444, 97–101. doi: 10.1038/nature05248

Krogh, A., Larsson, B., Von Heijne, G., and Sonnhammer, E. (2001). Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. *J. Mol. Biol.* 305, 567–580. doi: 10.1006/jmbi.2000.4315

Kubisiak, T., Anderson, C., Amerson, H., Smith, J., Davis, J., and Nelson, C. (2011). A genomic map enriched for markers linked to Avr1 in *Cronartium quer-cuum* f.sp. *fusiforme*. *Fungal Genet. Biol.* 48, 266–274. doi: 10.1016/j.fgb.2010.09.008

Kues, U., and Liu, Y. (2000). Fruiting body production in basidiomycetes. *Appl. Microbiol. Biotechnol.* 54, 141–152. doi: 10.1007/s002530000396

Kües, U., and Rühl, M. (2011). Multiple multi-copper oxidase gene families in basidiomycetes-what for? *Curr. Genomics* 12, 72–94. doi: 10.2174/138920211795564377

Kulkarni, R., Kelkar, H., and Dean, R. (2003). An eight-cysteine-containing CFEM domain unique to a group of fungal membrane proteins. *Trends Biochem. Sci.* 28, 118–121. doi: 10.1038/nature03449

Langfelder, K., Streibel, M., Jahn, B., Haase, G., and Brakhage, A. (2003). Biosynthesis of fungal melanins and their importance for human pathogenic fungi. *Fungal Genet. Biol.* 38, 143–158. doi: 10.1016/S1087-1845(02)00526-1

Leonowicz, A., Cho, N., Luterek, J., Wilkolazka, A., Wojtas-Wasilewska, M., Matuszewska, A., et al. (2001). Fungal laccase: properties and activity on lignin. *J. Basic Microb.* 41, 185–227. doi: 10.1002/1521-4028(200107)41:3/4<185::AID-JOBM185>3.0.CO;2-T

Letunic, I., and Bork, P. (2011). Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy. *Nucleic Acids Res.* 39, W475–W478. doi: 10.1093/nar/gkr201

Li, L., Stoeckert, C., and Roos, D. (2003). OrthoMCL: identification of ortholog groups for eukaryotic genomes. *Genome Res.* 13, 2178–2189. doi: 10.1101/gr.122450

Lundell, T., Makela, M., and Hilden, K. (2010). Lignin-modifying enzymes in filamentous basidiomycetes - ecological, functional and phylogenetic review. *J. Basic Microb.* 50, 5–20. doi: 10.1002/jobm.200900338

Martin, F., Aerts, A., Ahren, D., Brun, A., Danchin, E., Duchaussoy, F., et al. (2008). The genome of *Laccaria bicolor* provides insights into mycorrhizal symbiosis. *Nature* 452, 88–92. doi: 10.1038/nature06556

Martinez, D., Larrondo, L. F., Putnam, N., Gelpke, M. D. S., Huang, K., Chapman, J., et al. (2004). Genome sequence of the lignocellulose degrading fungus *Phanerochaete chrysosporium* strain RP78. *Nat. Biotechnol.* 22, 695–700. doi: 10.1038/nbt967

Mims, C. W., Liljebelke, K. A., and Covert, S. F. (1996). Ultrastructure of telia and teliospores of the rust fungus *Cronartium quercuum* f. sp. *fusiforme*. *Mycologia* 88, 47–56. doi: 10.2307/376078

Morin, E., Kohler, A., Baker, A. R., Foulongne-Oriol, M., Lombard, V., Nagye, L. G., et al. (2012). Genome sequence of the button mushroom Agaricus bisporus reveals mechanisms governing adaptation to a humic-rich ecological niche. *Proc. Natl. Acad. Sci. U.S.A.* 109, 17501–17506. doi: 10.1073/pnas.1206847109

Nielsen, H., Engelbrecht, J., Brunak, S., and von Heijne, G. (1997). Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. *Protein Eng.* 10, 1–6. doi: 10.1093/protein/10.1.1

Olson, Å., Aerts, A., Asiegbu, F., Belbahri, L., Bouzid, O., Broberg, A., et al. (2012). Insight into trade off between wood decay and parasitism from the genome of a fungal forest pathogen. *New Phytol.* 194, 1001–1013. doi: 10.1111/j.1469-8137.2012.04128.x

Perez, A., Ramage, G., Blanes, R., Murgui, A., Casanova, M., and Martinez, J. P. (2011). Some biological features of *Candida albicans* mutants for genes coding fungal proteins containing the CFEM domain. *FEMS Yeast Res.* 11, 273–284. doi: 10.1111/j.1567-1364.2010.00714.x

Petersen, T., Brunak, S., Von Heijne, G., and Nielsen, H. (2011). SignalP 4.0: discriminating signal peptides from transmembrane regions. *Nat. Methods* 8, 785–786. doi: 10.1038/nmeth.1701

Rafiqi, M., Ellis, J., Ludowici, V., Hardham, A., and Dodds, P. (2012). Challenges and progress towards understanding the role of effectors in plant-fungal interactions. *Curr. Opin. Plant Biol.* 15, 477–482. doi: 10.1016/j.pbi.2012.05.003

Rasmussen, M., and Kellis, M. (2011). A Bayesian approach for fast and accurate gene tree reconstruction. *Mol. Biol. Evol.* 28, 273–290. doi: 10.1093/mol-bev/msq189

Rep, M. (2005). Small proteins of plant-pathogenic fungi secreted during host colonization. *FEMS Microbiol. Lett.* 253, 19–27. doi: 10.1016/j.femsle.2005.09.014

Saitoh, H., Fujisawa, S., Mitsuoka, C., Ito, A., Hirabuchi, A., Ikeda, K., et al. (2012). Large-scale gene disruption in *Magnaporthe oryzae* identifies MC69, a secreted protein required for infection by monocot and dicot fungal pathogens. *PLoS Pathog.* 8:e1002711. doi: 10.1371/journal.ppat.1002711

Saunders, D. G., Win, J., Cano, L. M., Szabo, L. J., Kamoun, S., and Raffaele, S. (2012). Using hierarchical clustering of secreted protein families to classify and rank candidate effectors of rust fungi. *PLoS ONE* 7:e29847. doi: 10.1371/journal.pone.0029847

Spanu, P. D. (2012). The genomics of obligate (and nonobligate) biotrophs. *Annu. Rev. Phytopathol.* 50, 91–109. doi: 10.1146/annurev-phyto-081211-173024

Stamatakis, A., Ludwig, T., and Meier, H. (2005). RAxML-III: a program for sequential, parallel and distributed inference of large phylogenetic trees. *ConcurR. Comput. Pract. Exp.* 17, 1705–1723. doi: 10.1093/bioinformatics/bti191
Stergiopoulos, I., and de Wit, P. (2009). Fungal effector proteins. *Annu. Rev. Phytopathol.* 47, 233–263. doi: 10.1146/annurev.phyto.112408.132637

Thrall, P., Laine, A., Ravensdal, M., Nemri, A., Dodds, P., Barretts, L., et al. (2012). Rapid genetic change underpins antagonistic coevolution in a natural host-pathogen metapopulation. *Ecol. Lett.* 15, 425–435. doi: 10.1111/j.1461-0248.2012.01749.x

Thurston, C. (1994). The structure and function of fungal laccases. *Microbiol* 140, 19–26. doi: 10.1099/13500872-140-1-19

Toome, M., Ohm, R. A., Riley, R. W., James, T. Y., Lazarus, K. L., Henrissat, B., et al. (2014). Genome sequencing provides insight into the reproductive biology, nutritional mode and ploidy of the fern pathogen *Mixia osmundae*. *New Phytol.* 202, 554–564. doi: 10.1111/nph.12653

Tyler, B., Tripathy, S., Zhang, X., Dehal, P., Jiang, R., Aerts, A., et al. (2006). *Phytophthora* genome sequences uncover evolutionary origins and mechanisms of pathogenesis. *Science* 313, 1261–1266. doi: 10.1126/science.1128796

Voegele, R., and Mendgen, K. (2011). Nutrient uptake in rust fungi: how sweet is parasitic life? *Euphytica* 179, 41–55. doi: 10.1007/s10681-011-0358-5

Vogler, D., and Bruns, T. (1998). Phylogenetic relationships among the pine stem rust fungi (*Cronartium* and *Peridermium* spp.). *Mycologia* 90, 244–257. doi: 10.2307/3761300

Wu, Y., Rasmussen, M., Bansal, M., and Kellis, M. (2013). TreeFix: statistically informed gene tree error correction using species trees. *Syst. Biol.* 62, 110–120. doi: 10.1093/sysbio/sys076

Xu, J., Saunders, C. W., Hu, P., Grant, R. A., Boekhout, T., Kuramae, E. E., et al. (2007). Dandruff-associated Malassezia genomes reveal convergent and

divergent virulence traits shared with plant and human fungal pathogens. *Proc. Natl. Acad. Sci. U.S.A.* 104, 18730–18735. doi: 10.1073/pnas.07067-56104

Zhu, X., and Williamson, P. (2004). Role of laccase in the biology and virulence of *Cryptococcus neoformans*. *FEMS Yeast Res.* 5, 1–10. doi: 10.1016/j.femsyr.2004.04.004

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Received: 01 April 2014; accepted: 06 June 2014; published online: 26 June 2014.

Citation: Pendleton AL, Smith KE, Feau N, Martin FM, Grigoriev IV, Hamelin R, Nelson CD, Burleigh JG and Davis JM (2014) Duplications and losses in gene families of rust pathogens highlight putative effectors. *Front. Plant Sci.* 5:299. doi: 10.3389/fpls.2014.00299

This article was submitted to Plant-Microbe Interaction, a section of the journal *Frontiers in Plant Science*.

Copyright © 2014 Pendleton, Smith, Feau, Martin, Grigoriev, Hamelin, Nelson, Burleigh and Davis. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.