High-mobility group (HMG) B proteins are eukaryotic DNA-binding proteins characterized by the HMG-box functional motif. These transcription factors play a pivotal role in global genomic functions and in the control of genes involved in specific developmental or metabolic pathways. The filamentous ascomycete Podospora anserina contains 12 HMG-box genes. Of these, four have been previously characterized; three are mating-type genes that control fertilization and development of the fruit-body, whereas the last one encodes a factor involved in mitochondrial DNA stability. Systematic deletion analysis of the eight remaining uncharacterized HMG-box genes indicated that none were essential for viability, but that seven were involved in the sexual cycle. Two HMG-box genes display striking features. PaHMG5, an ortholog of SpSte11 from Schizosaccharomyces pombe, is a pivotal activator of mating-type genes in P. anserina, whereas PaHMG9 is a repressor of several phenomena specific to the stationary phase, most notably hyphal anastomoses. Transcriptional analyses of HMG-box genes in HMG-box deletion strains indicated that PaHMG5 is at the hub of a network of several HMG-box factors that regulate mating-type genes and mating-type target genes. Genetic analyses revealed that this network also controls fertility genes that are not regulated by mating-type transcription factors. This study points to the critical role of HMG-box members in sexual reproduction in fungi, as 11 out of 12 members were involved in the sexual cycle in P. anserina. PaHMG5 and SpSte11 are conserved transcriptional regulators of mating-type genes, although P. anserina and S. pombe diverged 550 million years ago. Two HMG-box genes, SOX9 and its upstream regulator SRY, also play an important role in sex determination in mammals. The P. anserina and S. pombe mating-type genes and their upstream regulatory factor form a module of HMG-box genes analogous to the SRY/SOX9 module, revealing a commonality of sex regulation in animals and fungi.
Podospora anserina, a coprophilous fungus, is used extensively as a model organism to address questions of sexual development and mating-type functions. Its mating-type locus contains three HMGB genes that encode transcription factors involved in fertilization and fruit-body development. We present the functional characterization of the remaining HMGB genes, which revealed that 11 of 12 HMGB genes were involved in sexual development. An analysis of the relationships between these genes uncovered a regulatory network governing the expression of mating-type genes. PaHMG5 is a key transcription factor that operates upstream of mating-type genes in this network. A homolog of PaHMG5 performs a similar function in the fission yeast Schizosaccharomyces pombe, which diverged from P. anserina 550 million years ago. The conservation of a regulatory circuit over such a prolonged timeframe is a striking exception to the general observation that sex developmental pathways are highly variable, even across closely related lineages. A module consisting of two HMGB transcription factors (Sry and Sox9) is a key regulator of sex determination in mammals. We propose that the module containing PaHMG5 and mating-type HMGB genes is the fungal counterpart of the mammalian module, revealing a commonality of sex regulation in animals and fungi.
Citation: Ait Benkhali J, Coppin E, Brun S, Peraza-Reyes L, Martin T, Dixelius C, et al. (2013) A Network of HMG-box Transcription Factors Regulates Sexual Cycle in the Fungus Podospora anserina. PLoS Genet 9(7): e1003642. doi:10.1371/journal.pgen.1003642
Editor: Joseph Heitman, Duke University Medical Center, United States of America
Received: December 28, 2012; Accepted: June 3, 2013; Published: July 18, 2013
Copyright: © 2013 Ait Benkhali et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: JAB was funded by the doctoral school GGC 426. TM and CD were funded by the Sida SWE-2005-453 and the Swedish University of Agricultural Sciences. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
High-mobility-group box (HMGB) proteins  include chromatin architectural proteins as well as specific transcription factors that are involved in highly diverse functions ranging from sex determination  to extracellular immune signaling , . All of these functions rely on the HMG-box, an HMGB conserved motif containing approximately 80 amino acids arranged in a distinctive L-shaped three-α-helical fold , . This HMG-box motif sharply bends DNA  and facilitates the assembly of transcriptional complexes that involve other proteins by distorting chromatin . Based on phylogenetic analyses of the HMG-box, the HMGB superfamily can be divided into two families, HMGB-UBF_HMG and SOX/TCF/MATA_HMG . These families were named after the best known representative in each group, namely hUBF , the SOX (Sry-type HMG-box) genes , TCF-1  and MATa-1 . The HMGB-UBF_HMG family is considered to be sequence non-specific, as generic UBF proteins can bind both ribosomal DNA regulatory sequences and sequences across the entire ribosomal DNA repeat . Members of this family are present in plants, fungi and animals. On the other hand, the SOX/TCF/MATA_HMG family contains proteins that bind specific DNA sequences, with a common T/A rich core (reviewed in ). The SOX/TCF/MATA_HMG family is subdivided into the SOX-TCF_HMG and the MATA_HMG subfamilies. The MATA_HMG subfamily includes exclusively fungal proteins, most of which are involved in sexual processes, while SOX-TCF_HMG genes appear to be restricted to animals .
Members of the HMGB-UBF_HMG family function in many processes such as transcription, genomic stability and the three R's (replication, recombination and DNA repair) (reviewed in ). An exhaustive analysis of HMGB-UBF_HMG_box encoding genes was performed in Saccharomyces cerevisiae in a previous study, and it showed HMGB-UBF_HMG-box genes to have highly diversified functions. In brief, the six HMGB-UBF_HMG_box genes from S. cerevisiae are involved in ribosomal DNA transcription , chromatin remodeling complexes , , and mitochondrial metabolism and energy , , . Analysis of HMGB-UBF_HMG-box genes in S. pombe mainly focused on the splsd1 and splsd2 genes. The SpLsd1 protein functions as a histone demethylase, and both SpLsd1/2 are believed to affect the epigenetic state of the cell .
Members of the SOX/TCF/MATA_HMG family are important regulators of differentiation and the sexual process. In mammals, the Y-linked testis-determining factor, SRY, is the founding member of the SOX subfamily . SRY induces male sex determination by regulating Sox9 , , . Unlike Sry, Sox9 is conserved among non-mammalian vertebrates and induces male-to-female sex-reversal when mutated , illustrating that it is ancestral and critical for vertebrate function. HMG-box genes also play a central role in the fungal mating process. Almost all mating-type loci from Ascomycota contain at least one MATA_HMG gene , . The most notable exception to this is S. cerevisiae (Saccharomycotina), which lost one MATA_HMG mating-type gene during evolution  but retained the prototypical α1 protein MATα1p. MATα1p was identified as a member of the MATA_HMG subfamily . ROX1 is the other gene of this subfamily in S. cerevisiae and one of the few fungal MATA_HMG-box genes that are not involved in sexual reproduction. Instead, it represses the expression of hypoxic genes . In S. pombe, mating-type gene transcription is regulated by the MATA_HMG-box gene, SpSte11 . The activation of the sexual process in S. pombe relies on the nuclear accumulation of SpSte11, which is triggered by starvation and pheromone signaling . An allelic system of HMG-box genes that putatively determine mating types was also identified in Zygomycota , , ,  and Microsporidia , , but they were part of the SOX-TCF_HMG and HMGB-UBF_HMG groups . In Basidiomycota, mating-type loci do not contain HMG-box genes, but MATA_HMG-box genes play an essential role in sexual development, as shown in Ustilago maydis , , , Coprinopsis cinerea  and Cryptococcus neoformans , .
A genome-wide systematic deletion analysis of transcription factors performed in the homothallic fungus Fusarium graminearum indicated that several HMG-box genes are involved in sexual development . However, not all HMG-box genes from F. graminearum were deleted in that study and their genetic relationships were not investigated. Here, we report the systematic deletion analysis of HMG-box genes in the heterothallic fungus Podospora anserina as well as an in-depth analysis of their genetic interactions. The P. anserina genome encodes a total of 12 HMG-box genes. Of these, four have been previously characterized, including three mating-type genes (FPR1, FMR1 and SMR2) belonging to the MATA_HMG subfamily , and the nuclear gene encoding a mitochondrion-targeted HMGB-UBF-box protein (mtHMG1, Pa_1_13340)  (Table 1). FPR1 is encoded by the mat+ idiomorph , while FMR1 and SMR2 are both encoded by the mat- idiomorph , . These mating-type genes control fertilization and the development of the fruit-body (reviewed in ). The remaining eight uncharacterized HMG-box genes encode two MATA_HMG-box proteins, four HMGB-UBF_HMG-box proteins, and two proteins that contain an HMG-box with an atypical residue. In the present study, these genes were individually deleted and the phenotypes of mutant strains were carefully assessed. The results revealed that all genes were dispensable for viability. Moreover, six mutant strains displayed sexual reproduction defects, and one of the two remaining mutants displayed vegetative defects. Transcriptional analyses indicated that an ortholog of SpSte11 is at the hub of the HMG-box gene network and that it regulates the expression of FMR1 and FPR1 and fertility.
Identification of HMG-box genes in P. anserina
The Podospora protein database (Materials and Methods) was searched for HMGB proteins using the previously defined HMG-box consensus sequence  as query in a Blastp search . P. anserina genome annotation in Fungal Transcription Factor  and Superfamily databases  was examined to identify possible missing HMGB genes in our analyses. A total of 12 HMG-box genes were identified. Of these, eight had not been previously characterized and were here named PaHMG2 to PaHMG9 (Table 1). CD-Search  identified five MATA_HMG-box proteins and five HMGB-UBF_HMG-box proteins (Table 1), but failed to categorize two HMG-box proteins. These two proteins, PaHMG7 and PaHMG9, contained an atypical charged residue instead of a conserved aliphatic or aromatic amino acid at position 9 (Figure 1). However, they contained the specific aromatic amino acids that anchor the recognition helix of the HMG-box domain to the hydrophobic core at position 8 and 11, confirming that they belong to the HMGB superfamily. Additional domains are shown in Table 1. Except for mtHMG1, which contained a mitochondrial targeting signal, the remaining 11 HMG-box proteins from P. anserina were predicted to localize to the nucleus.
A search for orthologs of P. anserina HMG-box genes with FUNGIpath  in selected fungal species indicated that three had an ortholog in S. cerevisiae, six had an ortholog in S. pombe, and all had an ortholog in Neurospora crassa (Table 1). PaHMG6, PaHMG7 and PaHMG8 had orthologs in Basidiomycota but the functions of these orthologs have not been investigated yet. Prf1, Rop1 and HMG3 from U. maydis , , Mat2 from C. neoformans  and Pcc1 from C. cinerea  encoded HMG-box regulators that are related to the sexual cycle in Basidiomycota. FUNGIpath did not detect any P. anserina orthologs for these HMG-box genes.
Most HMG-box genes localized to chromosome 1, which contains the mating-type locus . Only three out of 12 HMGB genes mapped outside of chromosome 1, suggesting that the distribution of HMG-box genes may be biased to chromosome 1. However, statistical tests indicated this bias to be inconclusive (Material and Methods).
Evolutionary analysis of fungal HMG-box proteins
The HMG-box of the 12 P. anserina HMGB proteins was extracted and grouped in a phylogram with the HMG-box motif of selected plant, animal and fungal species, including S. cerevisiae, S. pombe, N. crassa, Aspergillus nidulans, Cochliobolus heterostrophus and U. maydis (Figure 2). A total of 154 HMG-box domains were clustered in four groups, which overall corresponded to the previously defined HMGB groups : MATα_HMG (groupe A, Figure 2) , MATA_HMG (group B, Figure 2), SOX-TCF_HMG (group C, Figure 2) and HMGB-UBF_HMG (group D, Figure 2).
Clustering of core amino acid sequences using maximum-likelihood and model LG+G . Color labeling: MATα_HMG (A, green), MATA_HMG (B, yellow), SOX-TCF_HMG (C, orange), HMGB-UBF_HMG (D, H, blue), MAT1-1-3 in MATA_HMG (E, white) and STE11 in MATA_HMG (G, red). Other labels: Microsporidia MAT sex locus in HMGB-UBF_HMG (F, grey), Phycomyces blakesleeanus (Zygomycota) sexM (Phybl8) and sexP (Phybl9) (purple) and P. anserina proteins (red dots). LR-ELW values greater than 70% are shown. Abbreviations: Ailuropoda melanoleuca (Ailme); Ajellomyces capsulatus (Ajeca); Alternaria alternata (Altal); Alternaria brassicicola (Altbr); Anopheles gambiae (Anoga); Antonospora locustae (Antlo); Arabidopsis thaliana (Arath); Aspergillus fumigatus (Aspfu); Aspergillus nidulans (Aspni); Botryotinia fuckeliana (Botfu); Bipolaris sacchari (Bipsa); Caenorhabditis elegans (Caeel); Candida albicans (Canal); Cervus elaphus yarkandensis (Cerel); Ciona savignyi (Ciosa); Cochliobolus heterostrophus (Coche); Cochliobolus homomorphus (Cocho); Coprinopsis cinerea (Copci); Cryphonectria parasitica (Crypa); Culex quinquefasciatus (Culqu); Danio rerio (Danre); Dothistroma pini (Dotpi); Drosophila melanogaster (Drome); Enterocytozoon bieneusi (Entbi); Encephalitozoon cuniculi (Enccu); Fusarium acaciae-mearnsii (Fusac); Gibberella fujikuroi (Gibfu); Gibberella zeae (Gibze); Homo sapiens (Homsa); Lachancea thermotolerans (Lacth); Magnaporthe oryzae (Magor); Mycosphaerella graminicola (Mycgr); Podospora anserina (Podan); Neurospora crassa (Neucr); Penicillium marneffei (Penma); Phycomyces blakesleeanus (Phybl); Pneumocystis carinii (Pneca); Pyrenopeziza brassicae (Pyrbr); Pyrenophora teres (Pyrte); Rhynchosporium secalis (Rhyse); Saccharomyces cerevisiae (Sacce); Schizosaccharomyces japonicus (Schja); Schizosaccharomyces pombe (Schpo); Sordaria macrospora (Sorma); Stemphylium sarciniforme (Stesa); Strongylocentrotus purpuratus (Strpu); Takifugu rubripes (Takru); Ustilago maydis (Ustma); Verticillium dahliae (Verda); Xenopus laevis (Xenla); and Zygosaccharomyces rouxii (Zygro). Numbers after species names indicate α1 proteins (1), MATA_HMG (2), MAT1-1-3 (3), SOX (4), HMGB-UBF_HMG (5) and other HMG domains (6–9). When more than one domain was present for the same species, the suffix a, b or c was used. Units indicate the number of amino acid changes per position. Species codes and accession numbers grouped by evolutionary affinity are listed in Table S2.
Group A (Figure 2) included exclusively the mating-type transcription factors with an α1 domain, which was proposed to correspond to a new class of HMG-box , . Group A formed a clade related to the MATA_HMG subfamily and contained FMR1 (Podan1). The other HMG-box mating-type transcription factors FPR1 (Podan2a) and SMR2 (Podan3) were placed in group B (Figure 2), which contained mostly MATA_HMG-box proteins. PaHMG5 (Podan2c) clustered within clade G, with SpSte11 (Schpo2c) of S. pombe  and NCU09387/FMF-1 (Neucr2c) of N. crassa , . The orthology of these three proteins, suggested by their phylogenic position, was independently confirmed with FUNGIpath. Moreover, the relationship of NCU09387 (Neucr2c) with SpSte11 (Schpo2c) was previously reported . PaHMG8 (Podan2b) was the only MATA_HMG-box protein encoded by a gene outside chromosome 1. Its ortholog in S. cerevisiae is ROX1p (Sacce2), a repressor of hypoxic genes, and one of the few MATA_HMG-box proteins not related to mating processes . Earlier evolutionary analysis of the oxygen-responding system in Kluyveromyces lactis and S. cerevisiae suggested that ROX1 was recruited specifically to control this system in S. cerevisiae . Rfg1p (Canal2), the ortholog of ROX1p (Sacce2) in Candida albicans, controls filamentous growth and virulence . Interestingly, the clade defined by PaHMG8 (Podan2b), ROX1p (Sacce2) and Rfg1p (Canal2) included a MAT1-2-1 mating-type protein (Ajeca2, see ‘&’ in clade B, Figure 2) and was placed close to the MAT1-2-1 (Podan2a, Aspni2a, Neucr2a and Sorma2) and MAT1-1-3 mating-type proteins (clade E, Figure 2). The placement of the PaHMG8 (Podan2b)/ROX1p (Sacce2) clade in the MATA_HMG subfamily and our functional analyses of PaHMG8 (Podan2b) (see below) supports the idea that the ancestral ROX1 gene was primarily involved in sexual development before being recruited for various other functions in S. cerevisiae and C. albicans. Surprisingly, group B contained an HMGB-UBF_HMG clade (clade H, Figure 2), which included PaHMG3 (Podan5e) and several fungal orthologs, notably the histone demethylase SpLsd1 (Schpo7a) of S. pombe . All proteins from clade H have SWIRM and amino oxydase domains, which are characteristic of histone demethylases, thus supporting their orthologous relationship.
Five P. anserina HMG-box proteins were classified into group D (Figure 2), which comprised members of the HMGB-UBF_HMG family, except the clade grouping PaHMG7 (Podan6b), NCU06874 (Neucr6b) from N. crassa and NP_595970 (Schpo7b) from S. pombe. These three proteins are characterized by a conserved atypical residue in the HMG-box domain [see the above section for its description in PaHMG7 (Podan6b)] and a SprT-like domain (Interpro accession number IPR006640). The four HMGB-UBF_HMG-box proteins, PaHMG2 (Podan5c), PaHMG4 (Podan5b), mtHMG1 (Podan5d), and PaHMG6 (Podan5a), branched into different clades of group D. PaHMG2 (Podan5c) belonged to a branch that included proteins with a Sterile Alpha Domain (SAM), which is involved in protein-protein and protein-RNA interactions . In some cases, the orthologous relationships, as defined by FUNGIpath, were not in agreement with the HMG-box based phylogeny. The first instance was the mitochondrial protein mtHMG1 (Podan5d) and its ortholog NCU02695 (Neucr5d) in N. crassa, which failed to group into the same clade within the HMGB-UBF_HMG family although FUNGIpath provided a high confidence score for the orthology of mtHMG1 (Podan5d) and NCU02695 (Neucr5d) HMG-box domain. Moreover, mtHMG1 (Podan5d) and NCU02695 (Neucr5d) have a mitochondrial targeting signal and are characterized by a DUF1898 domain, supporting an orthologous relationship. In another instance, Nhp6Ap (Sacce5a) and Nhp6Bp (Sacce5b), two functionally redundant putative inparalogs from S. cerevisiae collectively referred to as Nhp6p , , were placed in distant clades with PaHMG6 (Podan5a) and PaHMG4 (Podan5b), respectively, and may be outparalogs instead. In contrast, FUNGIpath identified Nhp6p (Sacce5a, b) as co-orthologs of PaHMG6 (Podan5a). Similarly, Hmo1p (Sacce5d) from S. cerevisiae belonged to the same clade as PaHMG4 (Podan5b), while FUNGIpath identified Hmo1p (Sacce5d) as an ortholog of PaHMG6 (Podan5a). However, the FUNGIpath confidence score for orthology of PaHMG6 (Podan5a), Nhp6p (Sacce5a, b), and Hmo1p (Sacce5d) was low. Further analyses will be necessary to resolve these phylogenetic ambiguities.
PaHMG9 (Podan6a), along with its N. crassa ortholog NCU07568 (Neucr6a), were placed between group C and D (Figure 2). Accordingly, CD-Search failed to place PaHMG9 (Podan6a) and NCU07568 (Neucr6a) into either the MATA_HMG or HMGB-UBF groups (Table 1). Interestingly, the PaHMG9 HMG-box (Podan6a) is characterized by an atypical residue (see above section) that is conserved in NCU07568 (Neucr6a). Our phylogenic analysis also placed a MATA_HMG-box protein from N. crassa (NCU02326, Neucr2d) between group C and D (Figure 2). Previous work reported that NCU02326 (Neucr2d) is related to SpeSte11 (Schpo2c) and is the closest homolog to NCU09387 (Neucr2c) . In agreement with this report, FUNGIpath identified NCU02326 (Neucr2d) as an ortholog of SpSte11 (Schpo2c, clade G in Figure 2) and an inparalog of NCU09387 (Neucr2c). However, unlike NCU09387 (Neucr2c), NCU02326 (Neucr2d) was not placed within clade G with all other SpSte11 orthologs. Such inconsistency is not unprecedented, as there are several evidences in N. crassa that effective defense against duplicated sequences prevents the maintenance of closely related paralogs in this species . To independently assess which of the two N. crassa inparalogs is conserved in P. anserina, we compared the environment of these genes in the two species. The shared synteny observed upstream and downstream of PaHMG5 (Podan2c) and NCU09387 (Neucr2c) (Figure S1A) confirms that these two genes are orthologs. The absence of conserved organization between the NCU02326 (Neucr2d) locus and its putative counterpart in P. anserina indicates that this gene is absent in P. anserina (Figure S1B).
Phenotypic analysis of HMG-box deletion strains during vegetative growth
The eight additional HMG-box genes identified in this study were inactivated by targeted gene deletion to assess their role in the life cycle of P. anserina. Deletions were verified by Southern blot analysis (Materials and Methods). The eight mutants were examined for growth and for macroscopic mycelium alterations on minimal medium in mat + and mat− context. The ΔPahmg2, ΔPahmg6 and ΔPahmg9 strains displayed reduced growth (Table 2). The eight mutants were also tested for cold- and thermo-sensitive growth at 18°C and 36°C and for sensitivity to caffeine (phosphodiesterase inhibitor), fludioxonil, which acts on osmoregulation, and sodium dodecylsulfate, which acts on the cell membrane, without revealing any additive phenotype.
Microscopic examination of a 4-day-old ΔPahmg9 thallus revealed that the leading edge (apical and sub-apical hyphae) displayed several developmental processes that typically occur in old hyphae in the wild-type strain. For instance, male gametes (spermatia) were observed in the sub-apical portion of a ΔPahmg9 thallus (Figure 3). P. anserina differentiates special cell structures dedicated to breach and to exploit solid cellulosic substrates like cellophane . These structures resemble appressoria, the specialized infective structures of fungal plant pathogens. The P. anserina appressorium-like structures were identical in wild-type and ΔPahmg9 strains (Figure 3A). However, they emerged sooner at the leading edge in the mutant strain than in the wild-type strain. The first appressorium-like structures were found at 1.58±0.31 mm and 4.08±0.20 mm from the hyphal tip in the ΔPahmg9 and wild-type strains, respectively. Finally, the most striking phenotype of the ΔPahmg9 strain was the dramatic increase in hyphal anastomoses. While few anastomoses are usually observed in the sub-apical portion of a wild-type thallus (Figure 3B), sub-apical hyphae from the ΔPahmg9 strain fused with each other, leading to bundles of tightly attached hyphae (Figure 3B). Moreover, apical hyphae often fused with each other in the ΔPahmg9 strain (Figure 3B). We never observed wild-type apical hypha anastomoses, suggesting that the ΔPahmg9 mutation deregulates cell fusion. We therefore named this gene KEF1 (KEep-on-Fusing 1). To our knowledge, KEF1 is the first identified gene to have a repressor function during cell fusion, whereas many genes activate this process in N. crassa ,  and P. anserina ,  (reviewed in ).
(A) Appressorium-like structures in the Δkef1 strain and in the wild type (WT). Development of appressorium-like structures was not affected in the Δkef1 strain. The pictures were taken every 1 µm as Z stacks, with “0” corresponding to hyphae growing onto the cellophane layer (see panel C). Arrowheads indicate needle-like hyphae that penetrated the cellophane layer; p: palm-like structures within the cellophane layer. (B) Comparison of hyphal fusion (anastomosis) in wild-type (WT) and Δkef1 strains. In the wild type, anastomoses were never observed in apical hyphae (1), while they were occasionally observed in the subapical area (2). In Δkef1 strain, anastomoses were profuse in apical (1) and in subapical (2) areas. Most notably, anastomosis occurred between the apex of hyphae at the leading edge and neighbouring apical hyphae (see arrows in Δkef1 in 1). (C) Schematic of appressorium-like development in P. anserina. a: appressorium-like structure including the palm-like structure and the needle-like hyphae; c: cellophane layer on which mycelium is growing. (D) Spermatia and their spermatogonia were present in the subapical area of the Δkef1 strain but were absent from the same area in the wild-type strain. Scale bar = 10 µm in all panels.
Phenotypic analysis of HMG-box deletion strains during sexual development
P. anserina is a heterothallic (self-sterile) fungus that has two mating types, mat+ and mat−. Haploid strains of each mating type initiate the sexual cycle by differentiating male gametes (spermatia) and female organs (protoperithecia). A phenotypic analysis of sexual development separately evaluates the male and female fertility of each mutant. The female organ is a multicellular structure comprising protective maternal hyphae and the ascogonium, which contains the female gametic nuclei. The spermatia are independent cells that can be collected and used to fertilize any strain of opposite mating type. A pheromone/receptor signaling system allows the ascogonium to recognize and fuse with sexually compatible spermatia. Therefore, fertilization can be controlled to initiate the development of the fruit-body. The male gametic nucleus is delivered into the ascogonium, which differentiates into the hymenium. Karyogamy takes place in the hymenium and is immediately followed by meiosis. Subsequently, the haploid nuclei are packed into ascospores. Ultimately, the mature ascospores are forcibly discharged from fruit-bodies (perithecia).
To determine whether HMG-box gene deletions impair sexual reproduction, each mutant strain was tested for male and female fertility in both mating types (mat+ and mat−). The test for male fertility consisted in fertilizing wild-type female organs with spermatia from each HMG-box deletion mutant (Figure 4A). This analysis indicated that all mutants, except ΔPahmg5, produced functional spermatia. Moreover, fertilized female organs differentiated into mature fruit-bodies producing ascospores, indicating that HMG-box deletions in the male nuclei did not affect any developmental steps in the fertilized female organs. For each mutant, reciprocal crosses with mat− and mat+ wild-type strains behaved similarly, indicating that the observed phenotype was not dependent on mating type. We further quantified spermatia produced by each deletion mutant and measured their activity in fertilization assays (Materials and Methods). Two mutants, ΔPahmg4 and Δkef1 produced five and seven times more spermatia, respectively, than the wild-type strain. Spermatia from these two mutants and from other mutants displayed the same fertilization ability as wild-type spermatia. The overproduction of spermatia observed in Δkef1 was in agreement with microscopic observations, illustrating that even young hyphae produce spermatia. The ΔPahmg5 mat+ and ΔPahmg5 mat− strains were found to produce as many spermatia as wild-type strains. Thus, male sterility in these ΔPahmg5 strains can be attributed either to an inability of spermatia to fertilize wild-type protoperithecia or to an arrest of perithecial development shortly after fertilization.
(A) Analysis of male and female fertility of HMG-box gene deletion mutants in crosses with wild-type tester strains. When cultures were confluent, sterile water was poured and dispersed over the surface of the mycelium. As each strain produced spermatia (male cells) and protoperithecia (female organs) regardless of its mating type, reciprocal fertilization of the mutant and wild-type strains took place and indicated whether the mutant was fertile as a male (donor) or as a female (receptor). For ΔPahmg5 some perithecia differentiated at the contact zone where mutant and wild-type mycelia fuse. In the resulting heterokaryotic mycelium, wild-type nuclei complement the male and/or female sterility defect of the mutant. Those perithecia were fertile and expelled numerous asci allowing genetic analysis. (B) Analysis of perithecium distribution in homozygous crosses of HMG-box gene deletion mutants. Fragmented mycelium from mat+ and mat− strains with the same HMG-box deletion were deposited at the center of a Petri dish and incubated until perithecia formed. Typically, the wild-type strains differentiated perithecia within a ring-like area.
To test the female fertility, we examined the formation of fruit-bodies by mycelia from each mutant fertilized with wild-type spermatia (Figure 4A). No fruit-bodies were produced on mutant mycelium in crosses involving ΔPahmg5, ΔPahmg6, ΔPahmg8, Δkef1 and Δmthmg1 strains, demonstrating that these mutants were sterile as female partners. In contrast, ΔPahmg2, ΔPahmg3, ΔPahmg4 and ΔPahmg7 strains produced fruit-bodies that developed normally and ejected as many ascospores as fruit-bodies from a wild-type cross. For each mutant, reciprocal crosses with mat− and mat+ wild-type strains behaved similarly, indicating that the observed phenotype was not dependent on mating type.
We also examined the distribution of perithecia on the mycelium. These assays were based on homozygous crosses and were performed by inoculating a mixture of fragmented mat+ and mat− mycelia in the center of a plate containing minimum medium under constant light illumination (Figure 4B). In this test, a wild-type strain differentiates fruit-bodies mostly in a ring-like area, which was 1 cm wide and located approximately 1 cm away from the inoculation point . Only ΔPahmg2 and ΔPahmg7 differentiated fruit-bodies with the same pattern as the wild-type strain. The alterations were striking for ΔPahmg4, since fruit-bodies were distributed on the entire surface of the culture and the ring-like area was no longer visible. The pattern was also different for ΔPahmg3, which developed a wider ring-like fruit-body area. Furthermore, we observed that this phenotype was exacerbated by the deletion of PaHMG7. This difference was confirmed by quantitative measurement of the ring-like diameter in the ΔPahmg3 ΔPahmg7 double mutant (Table 2). The increase in ring diameter (54.9 mm instead of 53.7 mm) observed in the double mutant suggested that the PaHMG3 and PaHMG7 genes together control the distribution of fruit-bodies. Fruit-bodies were never observed for ΔPahmg5, ΔPahmg8 and Δkef1. The ΔPahmg6 mutant did not produce fruit-bodies as fast as the wild-type strains, but extended incubation revealed that it was weakly female fertile. Fruit-bodies were 50 times less abundant than in a wild-type cross and they displayed an altered shape, with a smaller body and a larger neck than the wild-type fruit-bodies (data not shown). They began to eject ascospores 8 days after fertilization, whereas completion of the sexual process took 4 days in a wild-type cross. Genetic analysis of their progeny confirmed that they were produced by ΔPahmg6 parents and not by contaminating mycelium.
In conclusion, all the eight HMG-box deletion strains analyzed in this study are altered phenotypically. Complementation assays confirmed that the deletion was responsible for the observed phenotype (Materials and Methods). None of the eight HMG-box genes analyzed in this study were essential for viability, although the deletion of three genes (PaHMG2, PaHMG6 and KEF1) resulted in growth alterations. Most HMG-box gene deletions affected sexual reproduction (Table 2). PaHMG2 is the only P. anserina HMG-box gene whose deletion exclusively affected vegetative growth. Deletions of PaHMG3, PaHMG4 and PaHMG7 affected the distribution of fruit-bodies. More importantly, these deletions failed to affect vegetative growth and fruit-body development. Deletion of PaHMG6, PaHMG8 and KEF1 impaired female fertility, but mutant strains remained male fertile. Strikingly, the ΔPahmg5 mutant exhibits both male and female sterility, an uncommon phenotype in P. anserina that was only reported for mating-type mutants .
Heterokaryotic complementation of female sterility caused by HMG-box gene deletions
Microscopic observation of ΔPahmg5, ΔPahmg6, ΔPahmg8 and Δkef1 cultures revealed the presence of protoperithecia, indicating that female sterility was not due to the absence of female reproductive structures. Female sterility of ΔPahmg5, ΔPahmg6, ΔPahmg8 and Δkef1 can be attributed to a defect either in the development of the fruit-body envelope or in the formation of the hymenium, or both. To determine which tissue requires the HMG-box proteins, we performed trikaryon mosaics. In this experiment, mat+ and mat− HMG-box mutant strains were co-cultured with a strain containing a deletion of the mat+ idiomorph (Δmat) . Except for the loss of the mat+ mating-type gene, the Δmat strain had a wild-type genotype. Because the Δmat strain lacked the mating-type gene needed for fertilization, it could not participate in sexual reproduction and, thus, could not rescue hymenium defects. However, the Δmat strain maintained wild-type vegetative characteristics and could act as a helper strain by complementing maternal defects in a contest of sexually compatible mutant strains. Notably, the Δmat strain can provide maternal haploid tissues to form the perithecial envelope .
Fertile fruit-bodies were observed in the three trikaryotic cultures of ΔPahmg6, ΔPahmg8 and Δkef1 (data not shown), indicating that Δmat cells surrounding the hymenia of these HMG-box mutants promote the development of functional perithecia. These results demonstrated that PaHMG6, PaHMG8 and KEF1 are required for making the perithecial envelope and that they are dispensable for the hymenium. We also examined the behaviour of the Δmthmg1 strain in a trikaryon mosaic test. As previously described , most Δmthmg1 cultures died quickly and did not differentiate protoperithecia. A trikaryotic culture involving sexually compatible Δmthmg1 strains and Δmat strain produced perithecia, but no ascospores were formed (data not shown). These results indicate that the Δmat tissue complemented the vegetative requirement for mtHMG1 during protoperithecium formation, while the absence of mtHMG1 in the hymenium led to an arrest of its development. The protein mtHMG1 is therefore required for both protoperithecial and hymenium development.
Trikaryotic cultures of ΔPahmg5 did not produce any perithecia. It is noteworthy that mating-type gene mutations do not affect PaHMG5 transcription , thus excluding the possibility that the Δmat strain is unable to complement the ΔPahmg5 defects because it is itself affected for PaHMG5 expression. Failure to restore fertility of ΔPahmg5 in the trikaryotic test could reflect the need for this gene in both male and female fertility, which can be difficult to be restored simultaneously. Therefore, we designed an assay to separately evaluate the restoration of female and male fertility of ΔPahmg5. To evaluate whether the female sterility of ΔPahmg5 could be rescued by a Δmat PaHMG5+ strain, co-cultures of mat+ ΔPahmg5 and Δmat PaHMG5+ strains were used as female partners in crosses fertilized by mat− wild-type spermatia. Conversely, to evaluate whether Δmat PaHMG5+ could restore the male competency of ΔPahmg5, cultures of the mat− wild-type strain were fertilized with spermatia issued from a co-culture of mat+ ΔPahmg5 and Δmat PaHMG5+ strains (Materials and Methods). Two fruit-bodies were observed on 10 Petri dishes for the female sterility restoration assay, while three fruit-bodies were observed on 10 mat− plates for the male restoration assay. Several thousand perithecia were formed with wild-type strains in similar experiments, indicating that male and female sterility of ΔPahmg5 strain was inefficiently complemented by the Δmat PaHMG5+ strain. Taken together, these data are consistent with a requirement of PaHMG5 for fertilization or an early stage of hymenium development. The arrest in the development of the hymenium at an early stage precludes any conclusion on the role of PaHMG5 in the development of the maternal perithecial envelope.
Defects caused by overexpression of PaHMG5
In the complementation assay of the ΔPahmg5 mutant, introduction of the PaHMG5 wild-type allele restored male fertility but transformants remained female sterile (Materials and Methods). Moreover, the transformants displaying the highest efficiency as male partners were vegetatively altered, displaying a flat vegetative mycelium without aerial hyphae. These data suggested that deregulation of PaHMG5 may be detrimental to the fungus. To examine the consequences of unregulated PaHMG5 expression, a plasmid was constructed to express a fusion of PaHMG5 with the Pagpd (glyceraldehyde-3 phosphate-dehydrogenase) promoter and initiation codon  (Materials and Methods). When this fusion construct was introduced into protoplasts from the wild-type strain, most of the recovered transformants displayed a flat mycelium. In addition, they exhibited female sterility in crosses with the wild-type strain. Genetic analysis (Materials and Methods) demonstrated that the vegetative defect (flat mycelium) and female sterility consistently co-segregated with the Pagpd::PaHMG5 transgene. These observations strongly suggested that female fertility relies on tight regulation of PaHMG5 expression.
Transcriptional expression of mating-related HMG-box genes and downstream target genes
PaHMG5, PaHMG6, PaHMG8, KEF1 and mtHMG1  are involved in the development of male and female organs, raising the question of their genetic interactions. To assess the relationships between these genes, quantitative real-time RT-PCR (RT-qPCR) was used to examine their expression patterns in ΔPahmg5 (Table S3), ΔPahmg6 (Table S4), ΔPahmg8 (Table S5), Δkef1 (Table S6) and Δmthmg1 (Table S7) strains in a mat+ and mat− context. Genes were defined as up-regulated in the mutant strain if the fold change (FC) was >1, with a p-value of <0.05 (see Materials and Methods for FC computation). On the other hand, genes were defined as downregulated in the mutant strain if 0<FC<1 with a p-value of <0.05. FCs with a 95% confidence interval including the value of 1 were not considered significant . The results are summarized in Figure 5. Most of the deletions had similar effects in mat+ and mat− strains. The few exceptions were statistically non-significant results in mat+ or mat− strains, e.g., mtHMG1 had a significant FC (1.35) in mat+ Δkef1 relative to wild-type and a non-significant FC (1.1) in mat− Δkef1 (Table S6). A map of the genetic interactions among PaHMG5, PaHMG6, PaHMG8, KEF1 and mtHMG1 was constructed based on the assumption that deletion of a regulatory gene can affect downstream target genes, while non–affected genes are either upstream regulatory genes or genes acting in an independent pathway (Figure 6). Interestingly, PaHMG6, PaHMG8, KEF1 and mtHMG1 appeared to interact and converge on a key regulator, PaHMG5. We detected an effect of mtHMG1 on the transcription of KEF1 and downstream HMG-box genes, suggesting that this protein is localized to the nucleus although it has been shown to be targeted to the mitochondria . Re-examination of its sequence using PSORTII  revealed the presence of several monopartite and bipartite nuclear localization signals.
FCs from Table S3 to S7 were converted to a heat map using Matrix2png . Squares for FCs with non-significant statistical values are in black. Squares for inapplicable values are in white. (A) Heat map in the mat+ strains. The FCs of HMG-box genes, FPR1 and selected FPR1 target genes are represented as indicated on the scale for the following strains: mat+ ΔPahmg6 (delPahmg6), mat+ Δmthmg1 (delmthmg1), mat+ Δkef1 (delkef1), mat+ ΔPahmg5 (delPahmg5), and mat+ ΔPahmg8 (delPahmg8). Gene number: MFP, Pa_2_2310; PRE2, Pa_4_1380; AOX, Pa_3_1710; PEPCK, Pa_4_3160. (B) Heat map in the mat− strains. The FCs of HMG-box genes, FMR1 and selected FMR1 target genes were represented as indicated on the scale in (A) for the following strains: mat− ΔPahmg6 (delPahmg6), mat− Δmthmg1 (delmthmg1), mat− Δkef1 (delkef1), mat− ΔPahmg5 (delPahmg5), and mat− ΔPahmg8 (delPahmg8). Gene number: MFM, Pa_1_8290; PRE1, Pa_7_9070.
Arrows with heads and blunt ends indicate activation and repression, respectively. The numbers next to the arrows indicate the average FC in gene expression between the wild-type and mutant strains.The relationship between mtHMG1 and PaHMG5 may be mediated by KEF1. Alternatively, mtHMG1 may by-pass KEF1 to repress PaHMG5 directly or indirectly. The consistency in the FCs suggest that PaHMG6 by-passes this cascade to activate PaHMG5 and PaHMG8 either directly or indirectly. The numbers next to the arrows connecting the mating-type genes (FMR1 and FPR1) and the downstream target genes are FCs that were obtained from .
We further examined the transcription levels of mating-type genes and of a selection of specific target genes , including the pheromone and pheromone receptor genes, in the different mutant strains (Table S3 to S7 and Figure 5). We found that the transcription of the mating-type genes FMR1 and FPR1 was reduced in ΔPahmg5 and ΔPahmg8 strains. As expected from the cascade shown in Figure 6, FMR1 was downregulated in mat− ΔPahmg6; in contrast, we failed to identify a significant reduction in FPR1 transcription in mat+ ΔPahmg6. However, all tested FPR1 target genes were either up- or downregulated in this latter strain, suggesting that PaHMG6 by-passes FPR1 to control its target genes. Overall, transcription of the mating-type target genes increased in Δkef1 strains, in agreement with the repressor effect of KEF1 on mating-type genes. Interestingly, transcription of pheromone receptor genes decreased in female sterile HMG-box mutants, even in the Δkef1 strain. This is in agreement with their essential role in female fertility .
PaHMG5 appears to be the major regulator of FPR1 and FMR1 mating-type genes in the network presented in Figure 6. This raises the possibility that constitutive expression of FPR1 and FMR1 could compensate for the absence of PaHMG5 and, thus, rescue sterility of the ΔPahmg5 mutant. Transgenic versions of FMR1 and FPR1 driven by the Angpd constitutive promoter were previously reported to complement loss-of-function of the corresponding gene . Functional Angpd-FMR1 or Angpd-FPR1 transgenes were introduced by genetic crosses in mat− and mat+ ΔPahmg5 mutant strains, respectively. We observed that ΔPahmg5 strains carrying these transgenes remained male and female sterile, indicating that PaHMG5 function is not limited to the transcriptional activation of mating-type genes. It may regulate mating-type target genes as a cofactor of mating-type transcription factors. Alternatively, it may regulate fertility genes that are different from mating-type target genes.
Promoter region analysis of HMG-box genes and target genes of HMG-box transcription factors reveals HMG-box binding sites
The regulatory network shown in Figure 6 consists of a cascade of HMG-box genes, suggesting that each gene may contain a binding site for the upstream regulating HMG-box factor. A search for a conserved binding site using MEME  identified a consensus motif (A/G)ACAAAGAA in KEF1, mtHMG1, PaHMG5, PaHMG8, and the FMR1 and FPR1 mating-type genes (Figure 7A). This consensus motif is very similar to the common core DNA motif A(A/T)CAA(A/T)G that is recognized by HMG-box transcription factors  (reviewed in ). The remaining P. anserina HMG-box genes either contained a sequence that displayed some differences to the A(A/T)CAA(A/T)G core sequence (PaHMG2, PaHMG3, PaHMG4 and PaHMG6) (Figure 7A), or they did not contain any related sequence (PaHMG7 and SMR2). Further analyses of mating-type target genes using MEME revealed that the mat+ pheromone receptor gene (PRE1), alternative oxydase gene (AOX), phospho-enol pyruvate kinase gene (PEPCK), Pa_1_24410 and Pa_6_7350 also contained the core HMG-box binding site ACAAAGA (Figure 7A). Interestingly, the two pheromone genes (MFM and MFP) displayed the same conserved core sequence, ATCAAAG. The mat− pheromone receptor (PRE2), Pa_4_80, Pa_4_3858 and Pa_5_9770 did not contain the core HMG-box binding site, suggesting that these genes are secondary targets of HMG-box genes. A total of eight genes contained the (A/G)ACAAAGAA consensus site. The comparison with the distribution of this site in the P. anserina genome indicated that the consensus site is significantly enriched in the selected set of genes examined here (p-value = 0.016) (Materials and Methods).
(A) Conserved sequences from a number of HMG-box genes and mating-type target genes identified using MEME program . Position −1 is the first nucleotide upstream of the translation initiation codon. (B) Weblogo of the consensus sequence generated by MEME (assembled from sequences listed in A). (C) The P. anserina consensus sequence aligned with the binding site of HMG-box proteins: MATa-1 (N. crassa) , SpSte11 (S. pombe) , Prf2 (U. maydis) ,  and Mat2 (C. neoformans) .
Electrophoretic mobility shift assays
We further tested the binding of the heterologously expressed His-tagged version of the entire PaHMG5 protein (HMG5His) in electrophoretic mobility shift assays (EMSAs) with primers containing the putative binding sites identified above. We first examined the binding of HMG5His to its putative target genes (FMR1, FPR1 and PaHMG8, Figure 6), including PaHMG5. These genes contained the consensus site (A/G)ACAAAGAA. Under the binding conditions used, HMG5His strongly bound to primers containing these sequences and showed decreased binding upon addition of unlabelled competitor (Figure 8A). Moreover, a primer with a scrambled FMR1 sequence did not display any binding of HMG5His, confirming the specificity of the sequence recognition. These experiments support the idea that PaHMG5 directly regulates these target genes as well as its own transcription, as previously demonstrated for SpSte11 . We further examined the affinity of HMG5His for the putative binding sites of the mating-type target genes (Figure S3), which displayed binding sites different from the consensus. Binding of HMG5His to sites that were different from the consensus (A/G)ACAAAGAA was much weaker than to the consensus binding site (see MFM, MFP AOX and PEPCK in Figure S3). The binding to MFM, MFP and PRE1 sites was more carefully investigated (Figure 8, B to D). Reciprocal competition of MFM, MFP and PRE1 sites with the FMR1 consensus site confirmed that HMG5His had a greater affinity for the consensus sequence than for the sites of these mating-type target genes. These results indicate that PaHMG5 recognized the sites in the 5′UTR of mating-type target genes, but may require the mating-type transcription factors to increase the efficiency of binding. We also evaluated the affinity of HMG5His for the putative binding sites of the HMG-box genes that are not regulated by PaHMG5 in our proposed network (Figure 6). HMG5His bound strongly to KEF1 and mtHMG1 sites, which are identical to the consensus sequence, while binding to PaHMG2, PaHMG3, PaHMG4 and PaHMG6 sites was much less efficient (Figure S4). Strikingly, the less efficient binding occurred with sites that were different from the consensus by modification of the central part of the site (see PaHMG2, PahMG3 and PaHMG6 in Figure S4). Reciprocal competition confirmed that HMG5His has much less affinity for the PaHMG2, PaHMG3, PaHMG4 and PaHMG6 sites than for the consensus site (Figure S5). We assume that the binding of HMG5His to these sites corresponded to the recognition of an HMG-box binding site by an HMG-box protein, but this did not demonstrate in vivo regulation by PaHMG5.
(A) Interaction of His tagged PaHMG5 (HMG5His) with probes corresponding to FMR1, FMR1 scrambled sequence (FMR1Sc), FPR1, PaHMG8 and PaHMG5 oligonucleotides. The probes are indicated below each panel. The interaction of HMG5His with probe was analyzed without competitor and in the presence of increasing amounts (given as fold molar excess below the triangles) of competitor (indicated above the triangles). (B) Interaction of the HMG5His protein with the MFM probe and reciprocal competition between FMR1 and MFM oligonucleotides. The probes are indicated below each panel. The competitors are indicated above the triangles. A control of the interaction of HMG5His with the FMR1 probe was performed as in A and included in the assay. HMG5His has a greater affinity for FMR1 than for MFM sequence, as indicated by the efficient exclusion of MFM probe by FMR1 competitor and the very inefficient exclusion of FMR1 probe by MFM competitor. (C) Interaction of the HMG5His protein with the MFP probe and reciprocal competition between FMR1 and MFP oligonucleotides. Legend as in B. HMG5His has a greater affinity for FMR1 than for MFP sequence, as indicated by the efficient exclusion of MFP probe by FMR1 competitor and the inefficient exclusion of FMR1 probe by MFP competitor. (D) Interaction of the HMG5His protein with the PRE1 probe and reciprocal competition between FMR1 and PRE1 oligonucleotides. Legend as in B. HMG5His has a greater affinity for FMR1 than for PRE1 oligonucleotides, as indicated by the efficient exclusion of PRE1 probe by FMR1 competitor and the inefficient exclusion of FMR1 probe by PRE1 competitor.
The Weblogo (Figure 7B) showed a consensus binding site obtained from the entire set of analyzed sequences. This consensus corresponded to the sequence that was recognized with the best efficiency by HMG5His in EMSA. Comparison of this consensus with other HMG-box binding sites (Figure 7C) reveals a striking similarity with the TR-box, which is bound by SpSte11 of S. pombe , . The P. anserina consensus site matched the TR-box in nine consecutive basepairs. The similarity of the HMG5His consensus binding site with the PRE-boxes from U. maydis and C. neoformans was reduced. In particular, the PRE-box of U. maydis had an A instead of a G in position 7 in the P. anserina consensus. This transition correlates with a weak binding of HMG5His (see PaHMG2 and PaHMG6 in Figure S4). Taken together, these data indicate that PaHMG5 binds most efficiently to sites that are almost identical to those that are recognized by SpSte11, its ortholog in S. pombe.
An analysis of the distribution of the (A/G)ACAAAGAA binding site in the P. anserina genome indicated that this motif is preferentially localized on segments of 1000 bp upstream of the predicted translational start sites (p-value of <0.0001) (Materials and Methods). In contrast, a scrambled binding site did not display significant enrichment in these 1000 bp regions corresponding to promoter regions (p-value = 0.08). A total of 502 genes contain the (A/G)ACAAAGAA motif, suggesting that up to 5% of P. anserina genes may be directly controled by MATA_HMG-box transcription factors.
Significance and conservation of the HMG-box gene functions in fungi
The data presented here provide new insights into the role and the relationships of HMG-box genes in the fungus P. anserina (Pezizomycotina). We revealed a network of HMG-box genes upstream of the FMR1 and FPR1 mating-type genes, which are themselves HMG-box genes. PaHMG5 plays a central role in this network by regulating mating-type gene transcription. We identified SpSte11 as the ortholog of PaHMG5 in S. pombe by different phylogenetic analyses, and we demonstrated that both recognize almost identical binding sites. KEF1 is another important member of the network, acting upstream of PaHMG5 as a repressor. KEF1 also appears as a repressor of hyphal anastomoses, which are dramatically deregulated in a Δkef1 strain. Moreover, the control of mating-type target genes by FMR1 and FPR1 was described in a previous report , which, together with the results presented here, provides the first exhaustive view of the regulatory circuits upstream and downstream of mating-type genes in a filamentous Ascomycete.
Strikingly, ten of the 12 HMG-box genes identified in the genome of P. anserina control fertility and sexual development. Although two genes (PaHMG2 and PaHMG7) were not directly involved in sexual reproduction, the deletion of PaHMG7 displayed a synergistic effect with the deletion of PaHMG3 on the perithecium distribution in co-cultures of mat− and mat+ strains. This observation demonstrates that 11 of the HMG-box genes have a direct or indirect function during sexual reproduction. The genome-wide systematic deletion analysis of F. graminearum transcription factors  allowed us to compare the resulting phenotypes with those of P. anserina (Table S8). A total of six out of the 11 HMG-box genes deleted from F. graminearum were involved in perithecium development. One HMG-box gene (FGSG_06760) was not deleted, and the perithecium distribution on the mycelium was not tested in F. graminearum as described here. Taken together, these data implicate the HMGB superfamily in sexual reproduction in Pezizomycotina, although further analyses in other species should validate this finding.
Comparative analyses of sex regulatory pathways indicate that transcription factors have often evolved to accommodate unique rather than conserved functions, even across closely related lineages ,  (reviewed in , ). For example, the pheromone-response pathway is controlled by an HMG-box protein in C. neoformans  and a transcription factor distantly related to the homeodomain family in S. cerevisiae , providing evidence for a dramatic change in a key regulator in these species, despite the presence of strong effector conservation . Our study revealed that PaHMG5 and SpSte11  are striking exceptions to the unusual plasticity of pathways regulating sex. These two orthologous HMG-box proteins positively regulate mating-type gene transcription in P. anserina and S. pombe. Moreover, the transcription factor analysis published by Son et al  suggests that the SpSte11 ortholog in F. graminearum (FGSG_01366, GzHMG010) also shares this conserved function. The phenotype of F. graminearum deleted for GzHMG010 recapitulates the phenotype of mating-type gene deletions (Table S8, see GzHMG010, MAT1-1-1 and MAT1-2-1) and mating-type locus deletion . This feature is expected for genes operating in the same pathway, but molecular analyses will be necessary to determine whether the control of mating-type gene expression by SpSte11 orthologs extends to other fungi. N. crassa contains two co-orthologs of PaHMG5, indicating that control of mating-type gene expression in N. crassa may be more complex than in P. anserina. Based on phylogenetic analyses and analyses of synteny, we identified NCU09387 as the ortholog of PaHMG5, while NCU02326, the inparalog of NCU09387, has no counterpart in P. anserina. Mutations of NCU09387 resulted in strains that could mate, but fruit-body development arrested before ascospore formation , . This phenotype is clearly different from the mating-type deletion, which resulted in a strain that is unable to mate . This observation raises the interesting possibility that NCU02326 may control the expression of mating-type genes during fertilization, while NCU09387 is involved in the regulation of mating-type genes after fertilization.
The HMG-box gene network regulates vegetative and mating-type fertility genes
Sexual reproduction relies on an interplay between vegetative tissues (the mycelium and the maternal hyphae) and sexual tissue (the hymenium). Vegetative tissues critically contribute to optimal conditions for ascogonium formation and to fertility by providing nutrients to growing fruit-bodies . Inside the fruit-body, the hymenium goes through karyogamy and meiosis, and provides signals to the vegetative tissues to sustain nutrient mobilization . This duality raises the question as to whether the alteration that entailed sterility in the mutant strains affects vegetative tissues, or sexual tissues via mating-type genes. The strains deleted for PaHMG6, PaHMG8 and KEF1 showed female sterility which was restored by complementation with the Δmat strain in trikaryotic mosaic tests. Although PaHMG6, PaHMG8 and KEF1 control transcription of mating-type genes, the trikaryotic test revealed that they also control fertility genes that are independent of mating-type genes but, nevertheless, critical for sexual reproduction.
The regulatory circuit shown in Figure 6 and the complete sterility resulting from the PaHMG5 deletion points to this gene being a major regulator of sexual reproduction in mat+ and mat− strains. We demonstrated that this gene controls mating-type genes and the pheromone/receptor systems. Pheromone and pheromone receptor genes were identified as the most critical mating-type targets for male and female fertility, respectively , . The deletion of PaHMG5 had a moderate effect on the transcription of mating-type genes (FPR1, FC = 0.25; FMR1, FC = 0.1; Table S3), but strongly reduced the transcription of pheromone genes (MFP, FC = 0.0004; MFM, FC = 0.004; Table S3). This effect on pheromone genes can explain the male sterility phenotype of ΔPahmg5 strains. Moreover, some genes necessary for the biogenesis of pheromones are regulated by mating-type genes , and their transcription may also be decreased in ΔPahmg5 strains, thereby enhancing male sterility. By contrast, transcription of the PRE2 receptor gene was reduced 4-fold, a value that is unlikely to result in complete female sterility, which was characteristic of the mat+ ΔPahmg5 strain. This observation suggests that the deletion of PaHMG5 affects the transcription of another target gene that is critical for female fertility. Further experiments will be necessary to identify this target gene and to determine if it is a target of mating-type genes or a vegetative critical fertility gene.
Relationship between the stationary phase and the HMG-box gene network
Several lines of evidence indicate that the competence of P. anserina for sexual reproduction is acquired during the stationary phase. Development of reproductive structures takes place during this phase , and the expression of mating-type genes increases up to 1000-fold upon entry into the stationary phase . Our study provides the first evidence of a link between the stationary phase and the HMG-box genes that control the sexual cycle. Deletion of KEF1 resulted in the formation of anastomoses, appressorium-like structures, and spermatia in young hyphae, all hallmarks of the stationary phase. We propose that this gene is a critical repressor of the switch to the stationary phase, which underlies the morphological transitions of this stage. However, this gene only contributed moderately to the increase in mating-type gene expression, as its deletion resulted in a 3-fold increase of FMR1 and FPR1 transcription. Other pathways are likely to connect entry into the stationary phase and sexual competence. Two MAP kinase pathways, PaMpk2 and PaMpk1, have essential roles in establishing the stationary phase in P. anserina , , . No experimental evidence is yet available in P. anserina to support the link between these MAP kinase pathways and HMG-box genes. Exhaustive analyses of MAP kinase pathways in S. cerevisiae and S. pombe revealed numerous connections with HMG-box proteins. For instance, SpSte11 is a direct target of the Spk1 MAP kinase pathway in S. pombe , suggesting that PaHMG5 may be a phosphorylation target of the Spk1 ortholog in P. anserina (PaMpk2). Another interesting connection was uncovered in S. cerevisiae between Nhp6p and the Mpk1p (Slk2p) MAP kinase pathway , which is orthologous to the PaMpk1 MAP kinase cascade in P. anserina. These connections provide future direction to find the molecular pathways linking the stationary phase and the sexual cycle.
A module of two HMG-box genes involved in sexual development is present in animals and fungi
Sex determination is highly variable, in contrast to other developmental systems that are well conserved through evolution. Even within a single kingdom, studies on common laboratory model organisms reveal that the genetic mechanisms of sex determination bear little, if any, resemblance. Strikingly, PaHMG5 and SpSte11 have conserved their function as mating-type regulators in P. anserina and S. pombe , although these two organisms have diverged 550 million years ago . Moreover, the mating-type genes in both fungi are themselves HMG-box-genes, thereby defining an HMG-box module that is conserved in both fungi. An HMG-box module involved in sexual reproduction is also present in the basidiomycete fungus U. maydis. The Rop1 protein directly regulates the transcription of prf1, another HMG-box gene . The Prf1 protein in turn induces the expression of mating-type genes , . HMG-box genes also play a critical role in sex determination in vertebrates. Sry, the mammalian Y-chromosomal testis-determining HMG-box gene is an activator of Sox9. Sox9 is also conserved among non-mammalian vertebrate species and has an ancestral and pivotal role in sex determination . The conservation of a similar regulatory HMG-box module in vertebrates and Dikarya reveals a commonality of sex regulation in animals and fungi. Sry and Sox9 are not orthologs of the fungal HMGB module genes. The two modules are thus analogs, not homologs. However, Martin et al  noted that the SexM protein of Phycomyces blakesleeanus was classified within the SOX-TCF_HMG subfamily. This placement is confirmed in the phylogram presented here (Figure 2: SexM, Phybl8). The discovery of a SOX-TCF_HMG-box regulator of mating-type genes in P. blakesleeanus would indicate an ancestral origin for an HMG box module involved in sex determination in Opisthokonta. Further investigations on HMG-box transcription factors and sex regultation in fungi should provide relevant information about the conservation and evolution of such modules.
Materials and Methods
Strains and media
The genetic and biological features of P. anserina were first described by Rizet and Engelmann  and current culture techniques can be found at http://podospora.igmors.u-psud.fr/methods.php. The strains used in this study were all derived from the S strain , which was used to determine the P. anserina genome sequence .
Identification of HMG-box genes and structural analyses
The P. anserina genome and protein databases are available at http://podospora.igmors.u-psud.fr/. Identification of HMG-box genes was checked using Fungal Transcription Factor (http://ftfd.snu.ac.kr/intro.php)  and the Superfamily databases (http://supfam.cs.bris.ac.uk/SUPERFAMILY/index.html) . For domain identification, Pfam  and CD-searches  were run on http://pfam.sanger.ac.uk/ and http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml, respectively. N. crassa HMG-box genes were obtained from http://www.broadinstitute.org/annotation/genome/neurospora/MultiHome.html. FUNGIpath  was used to identify systemically orthologous groups of genes (http://embg.igmors.u-psud.fr/fungipath/). Subcellular localization was predicted with PSORTII  at http://psort.hgc.jp/form2.html.
Statistical test for HMG-box gene distribution in the P. anserina genome
A total of nine HMG-box genes mapped to chromosome I and three HMG-box genes mapped to other chromosomes. Among the nine genes that mapped to chromosome I, three mapped to the mating-type locus and corresponded to mat+ and mat− idiomorphs (reviewed in ). The idiomorph considered for statistical analysis was the mat+ idiomorph, which contained one HMG-box gene (FPR1). The total number of HMG-box genes on chromosome I was therefore, seven. The sizes of chromosome I to VII were 8,813,526 bp, 5,165,605 bp, 4,712,833 bp, 3,808,397 bp 4,734,309 bp, 4,264,133 bp and 4,087,213 bp, respectively. The total size of the P. anserina genome was 35,686,016 bp. The expected number of HMG-box genes on chromosome I was between two and three. The p-value was calculated on a contingency table http://www.graphpad.com/quickcalcs/contingency1.cfm using Fisher's exact test (p-value = 0.07 for two genes; p-value = 0.18 for three genes). The expected number of HMG-box genes on chromosome I was not significantly different from seven genes.
Sequence acquisition, identification of consensus amino acids and phylogenetic analysis were performed as described previously . The HMG-box domains used to build the phylogenetic tree (Figure 2) and their alignment are in dataset S1 and dataset S2, respectively.
Gene deletion and complementation
To delete the chromosomal copy of the eight HMG-box genes, eight plasmids containing deletion cassettes conferring resistance to hygromycin B [HygR] were constructed according to the N. crassa strategy for high-throughput generation of gene deletion  with modifications aimed at minimizing errors in the 5′ and 3′ flanking regions  (see Table S9 for primer sequences). The deletion cassette was released from the vector by AscI digestion prior to transformation of ΔPaKu70 protoplasts . This transformation assay consistently yielded a high percentage of transformants with the correct deletion (>90%). One or two transformants obtained from each assay were subjected to Southern blot analysis to confirm the deletion (Figure S2), and one transformant with the expected hybridization pattern was selected for further analysis. The eight primary transformants containing corresponding targeted deletions were genetically purified by crossing with a wild-type strain of opposite mating type. This eliminated untransformed nuclei and segregated out the ΔPaKu70 mutation through its phleomycin resistant [PhleoR] phenotype. Screening for [HygR, PhleoS] strains allowed the identification of mat+ and mat− strains containing the HMG-box deletion but lacking ΔPaKu70. These strains constituted the stock of the deletion mutant for subsequent studies.
To ensure that the phenotype(s) observed for the HMG-box gene deletion mutants was actually due to inactivation of the relevant gene, the wild-type allele was reintroduced by transformation into the corresponding mutant. Wild-type alleles were obtained by amplifying fragments encompassing the corresponding gene (see Table S10 for primer sequences) and these were used directly for co-transformation of the mutant strain with the pPable vector , which conferred resistance to phleomycin. A significant number of co-transformants displaying a restored wild-type phenotype were recovered in each assay, demonstrating that phenotypes were not due to additional mutations (Table S11). It should be noted that introduction of the PaHMG5 wild-type allele into the ΔPaHMG5 mutant only rescued the male defect without restoring female fertility, indicating partial complementation.
Measure of spermatium production and activity
Mutants were grown for 7 and 14 days on Petri dishes containing minimal agar medium. Spermatia were recovered by washing the surface of the dish with 1.5 ml of water and were counted with a haemacytometer chamber. A diluted spermatium suspension was used to fertilize a wild-type strain of opposite mating type and perithecia were counted after incubation for 4 days. Typically, 30% to 50% of spermatia from the wild-type strain were fertilizing, giving rise to perithecia. Duplicates were carried out for all strains and the whole experiment was performed twice to confirm the data.
Microscopic observations were made on 4-day-old mycelia growing on cellophane placed on solid M0 medium (minimal medium lacking dextrin as carbon source) in a Petri dish. Small pieces (1 cm2) of cellophane containing mycelium were cut with a scalpel and mounted upside down in water. Pictures were taken with a Leica DMIRE 2 microscope coupled to a 10 MHz Cool SNAPHQ charge-coupled device camera (Roper Instruments). Since penetration of the cellophane through apressorium-like structures occurs perpendicularly to the surface, pictures were obtained at one micrometer increments to capture this process. Stacks of pictures were analyzed with ImageJ (http://rsb.info.nih.gov/ij) and deconvolution was performed with CombineZP (Alan Bradlay; email@example.com). Calculations of the mean and the standard deviation of the distance between emergence of appressorium-like structures and the leading edge of the thallus were made using 20 and 10 individual measurements in wild-type and ΔPahmg9 strains, respectively.
Construction of the ΔPahmg3 ΔPahmg7 double mutant
One heterokaryotic ΔPahmg3 mat+/ΔPahmg7 mat− culture was self-crossed, and mat+ and mat− homokaryotic double mutants were isolated from the progeny. Single and double mutants displayed a [HygR] phenotype and were, thus, undistinguishable. Therefore, a search was performed for asci showing first division segregation of the [HygR] phenotype (i.e., asci with two [HygR] and two [HygS] ascospores). Sensitive ascospore-derived cultures from these asci carried wild-type alleles of both genes; hence, resistant ones harbored mutations in both.
Deregulation of the PaHMG5 gene
The plasmid pBHGSTE11 contains a fusion of the P. anserina gpd promoter and initiation codon  with the coding phase of PaHMG5. The Pagpd promoter was a 0.38 kbp fragment obtained from the pPable plasmid , digested with NcoI, treated with Klenow and digested with XbaI. The PaHMG5 sequence was amplified with Pfu (Promega) from GA0AB103CF05  using 5PSTE11 and 3HindSTE11 primers (see Table S10 for primer sequences) and digested with HindIII. The Pagpd promoter and the PaHMG5 fragment were then ligated into the HindIII and XbaI sites of plasmid pBCHygro  to yield the pBHGSTE11 plasmid. Sequencing of the entire fusion confirmed that the first Pagpd codon was in frame with PaHMG5 and that no mutation altered PaHMG5.The pBHGSTE11 plasmid was introduced into mat+ wild-type protoplasts and 10 [HygR] transformants were phenotypically analyzed in a cross with a mat− wild-type strain. A total of nine transformants showed a flat and female sterile mycelium, but they were fertile as male partners. To determine more precisely the phenotypic effects resulting from the integration of the Pagpd::PaHMG5 fusion, the progeny from three representative transformants were subjected to genetic analysis. Segregation of the Pagpd::PaHMG5 fusion was scored through the [HygR] phenotype. Most unpigmented ascospores did not germinate. Their genotype could nevertheless be deduced from tetrad analysis. For two transformants, the presence of the Pagpd::PaHMG5 fusion was responsible for an ascospore pigmentation and germination defect (although most unpigmented ascospores did not germinate their genotype could nevertheless be deduced from tetrad analysis). However, pigmented ascospores giving rise to a [HygR] mycelium were recovered in the same progeny. These displayed a similar phenotypic vegetative alteration as observed in the primary transformants (flat mycelium and female sterility). In a transformant corresponding to a different integration site (different second division segregation % of [HygS]/[HygR]), the presence of the Pagpd::PaHMG5 fusion did not affect ascospore pigmentation; instead it conferred the vegetative mycelium defect (flat mycelium and female sterility). These data and conclusions were subsequently confirmed by analyzing second generation progeny which were obtained by crossing the purified [HygR] Pagpd::PaHMG5 bearing transformants with the wild-type strain.
Vegetative cultures for RNA preparation were performed on Petri dishes containing minimal medium and covered with a cellophane sheet (Bio Rad Hercules, USA). These cultures were inoculated with nine implants from mat+, mat− or HMG-box mutants of either mating type. Dishes were placed at 27°C under constant light (0 h) and were removed from the incubation room at 96 h, at which time P. anserina was competent for fertilization . Mycelia were harvested and RNAs were extracted as described previously . Purified RNAs were submitted to an additional DNase digestion in solution and cleaned up once more on RNeasy Plant Mini Kit (Qiagen, Hilden, Germany). Total RNAs were reverse transcribed with SuperScript III (Lifes Technologies) according to manufacturer's instructions. Each time the expression of an intronless gene was quantified, a non-reverse transcribed (NRT) control was performed for each biological replicate. For genes with introns, all primers were designed against two consecutive exons (see Table S12 for primer sequences) and an NRT control was systematically performed on a pool of biological replicates. Each RT-qPCR experiment contained at least five biological replicates and each point was performed in technical duplicate. Normalization genes for ΔPahmg5, ΔPahmg8 and Δmthmg1 were selected from a pool of ten housekeeping genes using geNorm  as described previously . geNorm failed to select normalization genes for ΔPahmg6 and ΔPahmg9 strains, probably because metabolism was altered in these strains. A single stable reference gene was identified in these strains using NormFinder . The normalization genes are listed in Table S13. RT-qPCR normalization was performed according to the relative quantification method with kinetic PCR efficiency correction. Standard error and 95% confidence interval calculations, and other statistical analyses were performed using REST 2009 software (Qiagen, Hilden, Germany) . The FC in the expression of a gene of interest was computed as the normalized relative quantity of cDNA in sample relative to that in the control:
FC = relative quantity of cDNA for the gene of interest x (geometric mean of relative quantity of cDNA for the normalization genes)−1
The relative quantity (RQ) of a cDNA was:
RQ = efficiency of amplification (arithmetic mean for WT strain replicates – arithmetic mean for mutant strain replicates)
The efficiency of amplification was above 1.8 for all analyzed genes. Genes were defined as downregulated in the mutant strain if 0<FC<1 with a p-value of <0.05. On the other hand, genes were defined as up-regulated in the mutant strain if FC>1, with a p-value of <0.05. FCs with a 95% confidence interval including the value of 1 were not considered significant .
Consensus motif search and analysis of HMG-box binding motif distribution
Motif searches were conducted using MEME http://meme.nbcr.net/meme/cgi-bin/meme.cgi  on segments of 1000 bp upstream of the predicted translational start sites. A group of core sequences including KEF1, mtHMG1, PaHMG5, PaHMG8, and the FMR1 and FPR1 mating-type genes was first analyzed using MEME to identify the consensus HMG-box binding site. Subsequently, each candidate sequence was included with the core sequences for MEME analysis. Segments of 1500 bp upstream of the predicted translational start sites were analyzed when genes yielded no hit using MEME. The number of occurrences of the (A/G)ACAAAGAA binding site was counted for the P. anserina genome (n = 1149 in 35,686,016 bp). Assuming a random distribution throughout the genome, the (A/G)ACAAAGAA sequence is expected to occur 0.55 time in the promoter regions of the genes selected for MEME analysis (17,000 bp). The observed and expected numbers of (A/G)ACAAAGAA sequence were compared on a contingency table (http://www.graphpad.com/quickcalcs/contingency1.cfm) and a p-value was computed using Fisher's exact test. The (A/G)ACAAAGAA sequence is significantly enriched in the selected set of gene analyzed with MEME (p-value = 0.016). We further analyzed the distribution of the (A/G)ACAAAGAA sequence in the promoter region of all P. anserina genes. The number of occurrences of the (A/G)ACAAAGAA sequene was counted on segments of 1000 bp upstream from the predicted translational start site of each P. anserina gene (n = 502 in 10,460,890 bp). Assuming a random distribution throughout the genome, the (A/G)ACAAAGAA sequence is expected to occur 339 times in the 10,460,890 bp of the promoter region DNA. Thus, the observed distribution of the [(G/A)ACAAAGAA] sequence in the promoter sequences is significantly different from the random distribution (p-value of <0.0001). The scrambled PaHMG5 binding site [A(G/A)AAGAACA] occurred 1398 times in the genome of P. anserina and 457 times in the promoter regions. A random distribution of the scrambled site would result in 413 occurrences in promoter regions. Thus, the observed distribution of the [(G/A)ACAAAGAA] sequence in promoter regions is not significantly different from the random distribution (p-value = 0.08).
Protein expression, purification and electrophoretic mobility shift assays
The full length PaHMG5 cDNA was amplified by PCR with LA Taq (TaKaRa, Shiga, Japan) from reverse transcribed total RNAs with primers Nde13 and HisBam13 (Table S10) according to the manufacturer's instructions. HisBam13 was designed to introduce an His6 tag downstream of the 3′ coding sequence. The PCR products were cloned into the pET28 vector (Novagen) between NdeI and BamHI restriction sites and inserts were sequenced to identify cDNA without mutations. E. coli BL21 (DE3) transformed with the recombinant PET28 vector was grown in 2×YT medium (MP Biomedicals) supplemented with kanamycin at 50 µg/ml. Approximately 800 ml of culture medium was incubated in a shaker at 200 rpm at 37°C until OD600~0.6–0.8. Protein expression was then induced with 0.5 mM isopropyl β-D-thiogalactopyranoside (Sigma) and the cell culture was further incubated at 15°C overnight. Cells were harvested by centrifugation, resuspended in 40 ml of 20 mM Tris pH9.0, 500 mM NaCl, 5 mM β-mercaptoethanol and protease inhibitor cocktail (Roche), and stored at −20°C. Cell lysis was achieved by sonication, and the cell extract was centrifuged at 20000 g for 30 min at 4°C. The His-tagged protein from the soluble fraction was purified on a nickel-nitrilotriacetic acid column (Qiagen Inc.) and eluted with an isocratic imidazole gradient, followed by a cation exchange step on a HiTrap Heparin column (GE Healthcare) equilibrated in 20 mM Tris (pH9.0), 300 mM NaCl, 5 mM β-mercaptoethanol, and 5% glycerol. The protein was eluted using a linear salt gradient. The production of recombinant protein was confirmed by SDS-PAGE. Lysis, soluble and some purified fractions were tested by western blot analysis. Proteins were resolved by SDS-PAGE and transferred onto nitrocellulose membrane (Protran, Whatman). The membrane was blocked by incubation for 1 hour at room temperature with 5% nonfat milk in TBS-T (Tris buffered saline-Tween 20, pH 7.5), incubated with rocking for 1 hour with anti-6His IgG conjugated with peroxidase (1/2,000) (Roche), and developed using the BM Blue POD substrate reagents from Roche.
Complementary primers (Table S14) were annealed to yield double stranded 33 bp oligonucleotides with single base 5′ overhangs consisting of a guanine to promote efficient labeling by T4 polynucleotide kinase . The double stranded oligonucleotides were 5′-end labeled by T4 polynucleotide kinase (Thermo Scientific) and [γ-32P]ATP (222TBq/mmole), according to the manufacturer's instructions. The probe was purified on MicroSpin G-25 columns (GE Healthcare) and further processed as described in .
EMSAs were performed with 340 µg of purified HMG5His in 20 mM Tris (pH9.0), 500 mM NaCl, 5 mM β-mercaptoethanol, 5 mM MgCl2 and 5% glycerol supplemented with 1.25 µg of poly(dI-dC) in a total volume of 11 µl. Samples were incubated for 20 min on ice, labeled probe was added, and incubation was continued for 1 h on ice. For competition experiments, unlabeled double stranded oligonucleotides were added after incubating 20 min with poly(dI-dC) and incubated further for 15 min before adding the labeled probe.
Protein-DNA complexes were separated in polyacrylamide gel (6%) in 0.25× Tris-borate-EDTA buffer at 200 V per gel for ~80 min. Radioactive probes were visualized using a Typhoon laser scanner (GE Healthcare).
Alignment of HMG-box domains used for Figure 2.
Comparative organization around SpeSte11 orthologs in P. anserina and N. crassa. Orthologs were determined using FUNGIpath . The gene names are enclosed in arrowed boxes indicating gene orientation. Orthologous genes are enclosed in boxes of identical color and connected with double arrows. Genes enclosed in white boxes do not have any ortholog in P. anserina or N. crassa. Gene sizes are not to scale. A: synteny in P. anserina and N. crassa for genes upstream and downstream of PaHMG5 and NCU09387. The conserved synteny indicates that PaHMG5 is the ortholog of NCU09387. B: organization of genes upstream and downstream of NCU02326 in N. crassa and search for a synteny in P. anserina. The absence of a conserved synteny confirms that the ortholog of NCU02326 is absent in P. anserina.
Genomic Southern blots of HMG-box mutant strains probed with the hygromycin sequence. 1: DNA from ΔPahmg6 transformants # 1 and # 2 digested with KpnI. Both transformants displayed the expected pattern for homologous recombination. 2: DNA from ΔPahmg2 transformants # 1 and # 2 digested with KpnI. Only transformant # 2 displayed the expected pattern for homologous recombination. 3: DNA from ΔPahmg3 transformants # 1 and # 2 digested with EcoRI. Both transformants displayed the expected pattern for homologous recombination. 4: DNA from ΔPahmg5 transformants # 1 and # 2 digested with EcoRI. Both transformants displayed the expected pattern for homologous recombination. 5: DNA from ΔPahmg6 transformants # 1 and # 2 digested with EcoRI. Both transformants displayed the expected pattern for homologous recombination. 6: DNA from ΔPahmg4 transformants # 1 and # 2 digested with EcoRI. Only transformant #1 displayed the expected pattern for homologous recombination. 7: DNA from Δkef1 transformant # 1 digested with EcoRI. This transformant displayed the expected pattern for homologous recombination. 8: DNA from ΔPahmg7 transformants # 1 and # 2 digested with KpnI. Both transformants displayed the expected pattern for homologous recombination.
Electrophoretic mobility shift assays with PaHMG5 and mating-type target gene oligonucleotides. (A) Interaction of His tagged PaHMG5 (HMG5His) with probes corresponding to mat− mating-type target genes: Pa_1_8290 (MFM), Pa_1_9070 (PRE1) and Pa_1_7350 (7350). The probe is indicated below each panel with the sequence of its core HMG-box binding site. The interaction of HMG5His with probe was analyzed without competitor and in the presence of competitor at 100 fold molar excess. (B) Interaction of His tagged PaHMG5 (HMG5His) with probes corresponding to mat+ mating-type target genes: Pa_2_2310 (MFP), Pa_3_1710 (AOX), Pa_4_3160 (PEPCK) and Pa_1_24410 (24410). Legend as in (A).
Electrophoretic mobility shift assays with PaHMG5 and HMG-box gene oligonucleotides. Interaction of His tagged PaHMG5 (HMG5His) with probes corresponding to HMG-box genes: KEF1, mtHMG1, PaHMG2, PaHMG3, PaHMG4 and PaHMG6. Legend as in Figure S3 (A).
Electrophoretic mobility shift reciprocal competition assays with PaHMG5 and HMG-box gene oligonucleotides. Legend as in Figure 8 (B). (A) Interaction of the HMG5His protein with the PaHMG2 probe and reciprocal competition between FMR1 and PaHMG2 oligonucleotides. HMG5His has a greater affinity for FMR1 than for PaHMG2 oligonucleotides, as indicated by the efficient exclusion of PaHMG2 probe by FMR1 competitor and the inefficient exclusion of FMR1 probe by PaHMG2 competitor. (B) Interaction of the HMG5His protein with the PaHMG3 probe and reciprocal competition between FMR1 and PaHMG3 oligonucleotides. HMG5His has a greater affinity for FMR1 than for PaHMG3 oligonucleotides, as indicated by the efficient exclusion of PaHMG3 probe by FMR1 competitor and the inefficient exclusion of FMR1 probe by PaHMG3 competitor. (C) Interaction of the HMG5His protein with the PaHMG4 probe and reciprocal competition between FMR1 and PaHMG4 oligonucleotides. HMG5His has a greater affinity for FMR1 than for PaHMG4 oligonucleotides, as indicated by the efficient exclusion of PaHMG4 probe by FMR1 competitor and the inefficient exclusion of FMR1 probe by PaHMG4 competitor. (D) Interaction of the HMG5His protein with the PaHMG6 probe and reciprocal competition between FMR1 and PaHMG6 oligonucleotides. HMG5His has a greater affinity for FMR1 than for PaHMG6 oligonucleotides, as indicated by the efficient exclusion of PaHMG6 probe by FMR1 competitor and the inefficient exclusion of FMR1 probe by PaHMG6 competitor.
Color scheme used for Jalview.
Code and accession numbers for proteins shown in Figure 2.
Relative quantification of HMG-box gene and mating-type target gene transcription in ΔPahmg5 (ΔPa_1_13940) and WT strains.
Relative quantification of HMG-box gene and mating-type target gene transcription in ΔPahmg6 (ΔPa_1_14230) and WT strains.
Relative quantification of HMG-box gene and mating-type target gene transcription in ΔPahmg8 (ΔPa_6_4110) and WT strains.
Relative quantification of HMG-box gene and mating-type target gene transcription in Δkef1 (ΔPahmg9, ΔPa_1_7190) and WT strains.
Relative quantification of HMG-box gene and mating-type target gene transcription in Δmthmg1 (ΔPa_1_13340) and WT strains.
Phenotype of F. graminearum strains deleted for HMGB genes (adapted from ).
Oligonucleotide primers used for HMG-box gene deletion.
Oligonucleotide primers used for HMG-box gene amplification.
Complementation of P. anserina mutant strains.
Oligonucleotide primers used for RT-qPCR.
Normalization genes (REF) used for RT-qPCR.
Oligonucleotides primers used for EMSA.
We thank Philippe Silar for valuable discussion and comments, Véronique Contamine for the gift of the Δmthmg1 strains, and Marie-Claude Serre for helpful advices on EMSA experiments. We acknowledge the assistance of Vivien Debuchy for the informatic analysis of binding site distribution.
Conceived and designed the experiments: RD EC CD HvT. Performed the experiments: JAB EC SB TM NL LPR RD. Analyzed the data: JAB EC SB TM NL RD . Wrote the paper: RD EC.
- 1. Bustin M (2001) Revised nomenclature for high mobility group (HMG) chromosomal proteins. Trends Biochem Sci 26: 152–153. doi: 10.1016/s0968-0004(00)01777-1
- 2. Koopman, P. 2010. HMG Domain Superfamily of DNA-bending Proteins: HMG, UBF, TCF, LEF, SOX, SRY and Related Proteins. eLS. DOI: 10.1002/9780470015902.a0002325.pub2
- 3. Bianchi ME, Manfredi AA (2007) High-mobility group box 1 (HMGB1) protein at the crossroads between innate and adaptive immunity. Immunol Rev 220: 35–46. doi: 10.1111/j.1600-065x.2007.00574.x
- 4. Yanai H, Ban T, Wang Z, Choi MK, Kawamura T, et al. (2009) HMGB proteins function as universal sentinels for nucleic-acid-mediated innate immune responses. Nature 462: 99–103. doi: 10.1038/nature08512
- 5. Read CM, Cary PD, Crane-Robinson C, Driscoll PC, Norman DG (1993) Solution structure of a DNA-binding domain from HMG1. Nucleic Acids Res 21: 3427–3436. doi: 10.1093/nar/21.15.3427
- 6. Weir HM, Kraulis PJ, Hill CS, Raine AR, Laue ED, et al. (1993) Structure of the HMG box motif in the B-domain of HMG1. EMBO J 12: 1311–1319.
- 7. Giese K, Cox J, Grosschedl R (1992) The HMG domain of lymphoid enhancer factor 1 bends DNA and facilitates assembly of functional nucleoprotein structures. Cell 69: 185–195. doi: 10.1016/0092-8674(92)90129-z
- 8. Stros M (2010) HMGB proteins: interactions with DNA and chromatin. Biochim Biophys Acta 1799: 101–113. doi: 10.1016/j.bbagrm.2009.09.008
- 9. Soullier S, Jay P, Poulat F, Vanacker JM, Berta P, et al. (1999) Diversification pattern of the HMG and SOX family members during evolution. J Mol Evol 48: 517–527. doi: 10.1007/pl00006495
- 10. Jantzen HM, Admon A, Bell SP, Tjian R (1990) Nucleolar transcription factor hUBF contains a DNA-binding motif with homology to HMG proteins. Nature 344: 830–836. doi: 10.1038/344830a0
- 11. Gubbay J, Collignon J, Koopman P, Capel B, Economou A, et al. (1990) A gene mapping to the sex-determining region of the mouse Y chromosome is a member of a novel family of embryonically expressed genes. Nature 346: 245–250. doi: 10.1038/346245a0
- 12. van de Wetering M, Oosterwegel M, Dooijes D, Clevers H (1991) Identification and cloning of TCF-1, a T lymphocyte-specific transcription factor containing a sequence-specific HMG box. EMBO J 10: 123–132. doi: 10.1093/intimm/3.11.1189
- 13. Staben C, Yanofsky C (1990) Neurospora crassa a mating-type region. Proc Natl Acad Sci U S A 87: 4917–4921. doi: 10.1073/pnas.87.13.4917
- 14. O'Sullivan AC, Sullivan GJ, McStay B (2002) UBF binding in vivo is not restricted to regulatory sequences within the vertebrate ribosomal DNA repeat. Mol Cell Biol 22: 657–668. doi: 10.1128/mcb.22.2.657-668.2002
- 15. Gadal O, Labarre S, Boschiero C, Thuriaux P (2002) Hmo1, an HMG-box protein, belongs to the yeast ribosomal DNA transcription system. Embo J 21: 5498–5507. doi: 10.1093/emboj/cdf539
- 16. Ray S, Grove A (2009) The yeast high mobility group protein HMO2, a subunit of the chromatin-remodeling complex INO80, binds DNA ends. Nucleic Acids Res 37: 6389–6399. doi: 10.1093/nar/gkp695
- 17. Stillman DJ (2010) Nhp6: a small but powerful effector of chromatin structure in Saccharomyces cerevisiae. Biochim Biophys Acta 1799: 175–180. doi: 10.1016/j.bbagrm.2009.11.010
- 18. Kao LR, Megraw TL, Chae CB (1993) Essential role of the HMG domain in the function of yeast mitochondrial histone HM: functional complementation of HM by the nuclear nonhistone protein NHP6A. Proc Natl Acad Sci U S A 90: 5598–5602. doi: 10.1073/pnas.90.12.5598
- 19. Lambert JR, Bilanchone VW, Cumsky MG (1994) The ORD1 gene encodes a transcription factor involved in oxygen regulation and is identical to IXR1, a gene that confers cisplatin sensitivity to Saccharomyces cerevisiae. Proc Natl Acad Sci U S A 91: 7345–7349. doi: 10.1073/pnas.91.15.7345
- 20. Tsaponina O, Barsoum E, Astrom SU, Chabes A (2011) Ixr1 is required for the expression of the ribonucleotide reductase Rnr1 and maintenance of dNTP pools. PLoS Genet 7: e1002061. doi: 10.1371/journal.pgen.1002061
- 21. Lan F, Zaratiegui M, Villen J, Vaughn MW, Verdel A, et al. (2007) S. pombe LSD1 homologs regulate heterochromatin propagation and euchromatic gene transcription. Mol Cell 26: 89–101. doi: 10.1016/j.molcel.2007.02.023
- 22. Barrionuevo F, Bagheri-Fam S, Klattig J, Kist R, Taketo MM, et al. (2006) Homozygous inactivation of Sox9 causes complete XY sex reversal in mice. Biol Reprod 74: 195–201. doi: 10.1095/biolreprod.105.045930
- 23. Foster JW, Dominguez-Steglich MA, Guioli S, Kwok C, Weller PA, et al. (1994) Campomelic dysplasia and autosomal sex reversal caused by mutations in an SRY-related gene. Nature 372: 525–530. doi: 10.1038/372525a0
- 24. Wagner T, Wirth J, Meyer J, Zabel B, Held M, et al. (1994) Autosomal sex reversal and campomelic dysplasia are caused by mutations in and around the SRY-related gene SOX9. Cell 79: 1111–1120. doi: 10.1016/0092-8674(94)90041-8
- 25. Knower KC, Kelly S, Harley VR (2003) Turning on the male-SRY, SOX9 and sex determination in mammals. Cytogenet Genome Res 101: 185–198. doi: 10.1159/000074336
- 26. Debuchy R, Berteaux-Lecellier V, Silar P (2010) Mating systems and sexual morphogenesis in Ascomycetes. In: Borkovich KA, Ebbole DJ, editors. Cellular and Molecular Biology of Filamentous Fungi. Washington, DC: ASM Press. pp. 501–535.
- 27. Martin T, Lu SW, van Tilbeurgh H, Ripoll DR, Dixelius C, et al. (2010) Tracing the origin of the fungal alpha1 domain places its ancestor in the HMG-box superfamily: implication for fungal mating-type evolution. PLoS One 5: e15199. doi: 10.1371/journal.pone.0015199
- 28. Butler G, Kenny C, Fagan A, Kurischko C, Gaillardin C, et al. (2004) Evolution of the MAT locus and its Ho endonuclease in yeast species. Proc Natl Acad Sci U S A 101: 1632–1637. doi: 10.1073/pnas.0304170101
- 29. Kwast KE, Burke PV, Brown K, Poyton RO (1997) REO1 and ROX1 are alleles of the same gene which encodes a transcriptional repressor of hypoxic genes in Saccharomyces cerevisiae. Curr Genet 32: 377–383. doi: 10.1007/s002940050291
- 30. Sugimoto A, Iino Y, Maeda T, Watanabe Y, Yamamoto M (1991) Schizosaccharomyces pombe ste11+ encodes a transcription factor with an HMG motif that is a critical regulator of sexual development. Genes Dev 5: 1990–1999. doi: 10.1101/gad.5.11.1990
- 31. Qin J, Kang W, Leung B, McLeod M (2003) Ste11p, a high-mobility-group box DNA-binding protein, undergoes pheromone- and nutrient-regulated nuclear-cytoplasmic shuttling. Mol Cell Biol 23: 3253–3264. doi: 10.1128/mcb.23.9.3253-3264.2003
- 32. Gryganskyi AP, Lee SC, Litvintseva AP, Smith ME, Bonito G, et al. (2010) Structure, function, and phylogeny of the mating locus in the Rhizopus oryzae complex. PLoS One 5: e15273. doi: 10.1371/journal.pone.0015273
- 33. Idnurm A, Walton FJ, Floyd A, Heitman J (2008) Identification of the sex genes in an early diverged fungus. Nature 451: 193–196. doi: 10.1038/nature06453
- 34. Lee SC, Corradi N, Byrnes EJ 3rd, Torres-Martinez S, Dietrich FS, et al. (2008) Microsporidia evolved from ancestral sexual fungi. Curr Biol 18: 1675–1679. doi: 10.1016/j.cub.2008.09.030
- 35. Wetzel J, Burmester A, Kolbe M, Wostemeyer J (2012) The mating-related loci sexM and sexP of the zygomycetous fungus Mucor mucedo and their transcriptional regulation by trisporoid pheromones. Microbiology 158: 1016–1023. doi: 10.1099/mic.0.054106-0
- 36. Lee SC, Corradi N, Doan S, Dietrich FS, Keeling PJ, et al. (2010) Evolution of the sex-related locus and genomic features shared in microsporidia and fungi. PLoS One 5: e10539. doi: 10.1371/journal.pone.0010539
- 37. Brefort T, Muller P, Kahmann R (2005) The high-mobility-group domain transcription factor Rop1 is a direct regulator of prf1 in Ustilago maydis. Eukaryot Cell 4: 379–391. doi: 10.1128/ec.4.2.379-391.2005
- 38. Hartmann HA, Kahmann R, Bolker M (1996) The pheromone response factor coordinates filamentous growth and pathogenicity in Ustilago maydis. EMBO J 15: 1632–1641.
- 39. Urban M, Kahmann R, Bolker M (1996) Identification of the pheromone response element in Ustilago maydis. Mol Gen Genet 251: 31–37. doi: 10.1007/bf02174341
- 40. Murata Y, Fujii M, Zolan ME, Kamada T (1998) Molecular analysis of pcc1, a gene that leads to A-regulated sexual morphogenesis in Coprinus cinereus. Genetics 149: 1753–1761.
- 41. Kruzel EK, Giles SS, Hull CM (2012) Analysis of Cryptococcus neoformans sexual development reveals rewiring of the pheromone-response network by a change in transcription factor identity. Genetics 191: 435–449. doi: 10.1534/genetics.112.138958
- 42. Lin X, Jackson JC, Feretzaki M, Xue C, Heitman J (2010) Transcription factors Mat2 and Znf2 operate cellular circuits orchestrating opposite- and same-sex mating in Cryptococcus neoformans. PLoS Genet 6: e1000953. doi: 10.1371/journal.pgen.1000953
- 43. Son H, Seo YS, Min K, Park AR, Lee J, et al. (2011) A phenome-based functional analysis of transcription factors in the cereal head blight fungus, Fusarium graminearum. PLoS Pathog 7: e1002310. doi: 10.1371/journal.ppat.1002310
- 44. Dequard-Chablat M, Allandt C (2002) Two copies of mthmg1, encoding a novel mitochondrial HMG-like protein, delay accumulation of mitochondrial DNA deletions in Podospora anserina. Eukaryot Cell 1: 503–513. doi: 10.1128/ec.1.4.503-513.2002
- 45. Debuchy R, Coppin E (1992) The mating types of Podospora anserina: functional analysis and sequence of the fertilization domains. Mol Gen Genet 233: 113–121. doi: 10.1007/bf00587568
- 46. Debuchy R, Arnaise S, Lecellier G (1993) The mat- allele of Podospora anserina contains three regulatory genes required for the development of fertilized female organs. Mol Gen Genet 241: 667–673. doi: 10.1007/bf00279909
- 47. Turgeon BG, Debuchy R (2007) Cochliobolus and Podospora: mechanism of sex determination and the evolution of reproductive lifestyle. In: Heitman J, Kronstad JW, Taylor JW, Casselton LA, editors. Sex in Fungi, molecular determination and evolutionary implications. Washington, DC: ASM Press. pp. 93–121.
- 48. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402. doi: 10.1093/nar/25.17.3389
- 49. Park J, Park J, Jang S, Kim S, Kong S, et al. (2008) FTFD: an informatics pipeline supporting phylogenomic analysis of fungal transcription factors. Bioinformatics 24: 1024–1025. doi: 10.1093/bioinformatics/btn058
- 50. Gough J, Karplus K, Hughey R, Chothia C (2001) Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure. J Mol Biol 313: 903–919. doi: 10.1006/jmbi.2001.5080
- 51. Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, et al. (2011) CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res 39: D225–229. doi: 10.1093/nar/gkq1189
- 52. Grossetete S, Labedan B, Lespinet O (2010) FUNGIpath: a tool to assess fungal metabolic pathways predicted by orthology. BMC Genomics 11: 81. doi: 10.1186/1471-2164-11-81
- 53. Espagne E, Lespinet O, Malagnac F, Da Silva C, Jaillon O, et al. (2008) The genome sequence of the model ascomycete fungus Podospora anserina. Genome Biol 9: R77. doi: 10.1186/gb-2008-9-5-r77
- 54. Jackson D, Lawson T, Villafane R, Gary L (2013) Modeling the structure of yeast MATα1: an HMG-box motif with a C-terminal helical extension. Open J Biophys 3: 1–12. doi: 10.4236/ojbiphy.2013.31001
- 55. Iyer SV, Ramakrishnan M, Kasbekar DP (2009) Neurospora crassa fmf-1 encodes the homologue of the Schizosaccharomyces pombe Ste11p regulator of sexual development. J Genet 88: 33–39. doi: 10.1007/s12041-009-0005-2
- 56. Johnson TE (1979) A neurospora mutation that arrests perithecial development as either male or female parent. Genetics 92: 1107–1120.
- 57. Borkovich KA, Alex LA, Yarden O, Freitag M, Turner GE, et al. (2004) Lessons from the genome sequence of Neurospora crassa: tracing the path from genomic blueprint to multicellular organism. Microbiol Mol Biol Rev 68: 1–108. doi: 10.1128/mmbr.68.1.1-108.2004
- 58. Fang ZA, Wang GH, Chen AL, Li YF, Liu JP, et al. (2009) Gene responses to oxygen availability in Kluyveromyces lactis: an insight on the evolution of the oxygen-responding system in yeast. PLoS One 4: e7561. doi: 10.1371/journal.pone.0007561
- 59. Kadosh D, Johnson AD (2001) Rfg1, a protein related to the Saccharomyces cerevisiae hypoxic regulator Rox1, controls filamentous growth and virulence in Candida albicans. Mol Cell Biol 21: 2496–2505. doi: 10.1128/mcb.21.7.2496-2505.2001
- 60. Kim CA, Bowie JU (2003) SAM domains: uniform structure, diversity of function. Trends Biochem Sci 28: 625–628. doi: 10.1016/j.tibs.2003.11.001
- 61. Costigan C, Kolodrubetz D, Snyder M (1994) NHP6A and NHP6B, which encode HMG1-like proteins, are candidates for downstream components of the yeast SLT2 mitogen-activated protein kinase pathway. Mol Cell Biol 14: 2391–2403. doi: 10.1128/mcb.14.4.2391
- 62. Dowell NL, Sperling AS, Mason MJ, Johnson RC (2010) Chromatin-dependent binding of the S. cerevisiae HMGB protein Nhp6A affects nucleosome dynamics and transcription. Genes Dev 24: 2031–2042. doi: 10.1101/gad.1948910
- 63. Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, et al. (2003) The genome sequence of the filamentous fungus Neurospora crassa. Nature 422: 859–868. doi: 10.1038/nature01554
- 64. Brun S, Malagnac F, Bidard F, Lalucque H, Silar P (2009) Functions and regulation of the Nox family in the filamentous fungus Podospora anserina: a new role in cellulose degradation. Mol Microbiol 74: 480–496. doi: 10.1111/j.1365-2958.2009.06878.x
- 65. Fu C, Iyer P, Herkal A, Abdullah J, Stout A, et al. (2011) Identification and characterization of genes required for cell-to-cell fusion in Neurospora crassa. Eukaryot Cell 10: 1100–1109. doi: 10.1128/ec.05003-11
- 66. Simonin AR, Rasmussen CG, Yang M, Glass NL (2010) Genes encoding a striatin-like protein (ham-3) and a forkhead associated protein (ham-4) are required for hyphal fusion in Neurospora crassa. Fungal Genet Biol 47: 855–868. doi: 10.1016/j.fgb.2010.06.010
- 67. Lalucque H, Malagnac F, Brun S, Kicka S, Silar P (2012) A non-Mendelian MAPK-generated hereditary unit controlled by a second MAPK pathway in Podospora anserina. Genetics 191: 419–433. doi: 10.1534/genetics.112.139469
- 68. Silar P (2011) Grafting as a method for studying development in the filamentous fungus Podospora anserina. Fungal Biol 115: 793–802. doi: 10.1016/j.funbio.2011.06.005
- 69. Read ND, Fleissner A, Roca GM, Glass L, N. (2010) Hyphal fusion. In: Borkovich KA, Ebbole DJ, editors. Cellular and Molecular Biology of Filamentous Fungi. Washington, DC: ASM Press. pp. 260–273.
- 70. Coppin E, Berteaux-Lecellier V, Bidard F, Brun S, Ruprich-Robert G, et al. (2012) Systematic deletion of homeobox genes in Podospora anserina uncovers their roles in shaping the fruiting body. PLoS One 7: e37488. doi: 10.1371/journal.pone.0037488
- 71. Coppin E, Arnaise S, Contamine V, Picard M (1993) Deletion of the mating-type sequences in Podospora anserina abolishes mating without affecting vegetative functions and sexual differentiation. Mol Gen Genet 241: 409–414. doi: 10.1007/bf00284694
- 72. Jamet-Vierny C, Debuchy R, Prigent M, Silar P (2007) IDC1, a pezizomycotina-specific gene that belongs to the PaMpk1 MAP kinase transduction cascade of the filamentous fungus Podospora anserina. Fungal Genet Biol 44: 1219–1230. doi: 10.1016/j.fgb.2007.04.005
- 73. Bidard F, Ait Benkhali J, Coppin E, Imbeaud S, Grognet P, et al. (2011) Genome-wide gene expression profiling of fertilization competent mycelium in opposite mating types in the heterothallic fungus Podospora anserina. PLoS One 6: e21476. doi: 10.1371/journal.pone.0021476
- 74. Ridder R, Osiewacz HD (1992) Sequence analysis of the gene coding for glyceraldehyde-3-phosphate dehydrogenase (gpd) of Podospora anserina: use of homologous regulatory sequences to improve transformation efficiency. Curr Genet 21: 207–213. doi: 10.1007/bf00336843
- 75. du Prel JB, Hommel G, Rohrig B, Blettner M (2009) Confidence interval or p-value?: part 4 of a series on evaluation of scientific publications. Dtsch Arztebl Int 106: 335–339.
- 76. Nakai K, Horton P (1999) PSORT: a program for detecting sorting signals in proteins and predicting their subcellular localization. Trends Biochem Sci 24: 34–36. doi: 10.1016/s0968-0004(98)01336-x
- 77. Coppin E, Debuchy R (2000) Co-expression of the mating-type genes involved in internuclear recognition is lethal in Podospora anserina. Genetics 155: 657–669.
- 78. Bailey TL, Elkan C (1994) Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol 2: 28–36.
- 79. van Beest M, Dooijes D, van De Wetering M, Kjaerulff S, Bonvin A, et al. (2000) Sequence-specific high mobility group box factors recognize 10–12-base pair minor groove motifs. J Biol Chem 275: 27266–27273. doi: 10.1074/jbc.m004102200
- 80. Kunitomo H, Higuchi T, Iino Y, Yamamoto M (2000) A zinc-finger protein, Rst2p, regulates transcription of the fission yeast ste11(+) gene, which encodes a pivotal transcription factor for sexual development. Mol Biol Cell 11: 3205–3217. doi: 10.1091/mbc.11.9.3205
- 81. Aono T, Yanai H, Miki F, Davey J, Shimoda C (1994) Mating pheromone-induced expression of the mat1-Pm gene of Schizosaccharomyces pombe: identification of signalling components and characterization of upstream controlling elements. Yeast 10: 757–770.
- 82. Tsong AE, Tuch BB, Li H, Johnson AD (2006) Evolution of alternative transcriptional circuits with identical logic. Nature 443: 415–420. doi: 10.1038/nature05099
- 83. Tuch BB, Galgoczy DJ, Hernday AD, Li H, Johnson AD (2008) The evolution of combinatorial gene regulation in fungi. PLoS Biol 6: e38. doi: 10.1371/journal.pbio.0060038
- 84. Dyer PS (2008) Evolutionary biology: genomic clues to original sex in fungi. Curr Biol 18: R207–R209. doi: 10.1016/j.cub.2008.01.014
- 85. Tsong AE, Brian BT, Johnson AD (2007) Rewiring transcriptional circuitry: mating-type regulation in Saccharomyces cerevisiae and Candida albicans as a model for evolution. In: Heitman J, Kronstad JW, Taylor JW, Casselton LA, editors. Sex in Fungi, molecular determination and evolutionary implications. Washington, DC: ASM Press. pp. 75–89.
- 86. Wong Sak Hoi J, Dumas B (2010) Ste12 and Ste12-like proteins, fungal transcription factors regulating development and pathogenicity. Eukaryot Cell 9: 480–485. doi: 10.1128/ec.00333-09
- 87. Desjardins AE, Brown DW, Yun SH, Proctor RH, Lee T, et al. (2004) Deletion and complementation of the mating type (MAT) locus of the wheat head blight pathogen Gibberella zeae. Appl Environ Microbiol 70: 2437–2444. doi: 10.1128/aem.70.4.2437-2444.2004
- 88. Ferreira AV, An Z, Metzenberg RL, Glass NL (1998) Characterization of mat A-2, mat A-3 and ΔmatA mating-type mutants of Neurospora crassa. Genetics 148: 1069–1079.
- 89. Coppin E, de Renty C, Debuchy R (2005) The function of the coding sequences for the putative pheromone precursors in Podospora anserina is restricted to fertilization. Eukaryot Cell 4: 407–420. doi: 10.1128/ec.4.2.407-420.2005
- 90. Kicka S, Silar P (2004) PaASK1, a mitogen-activated protein kinase kinase kinase that controls cell degeneration and cell differentiation in Podospora anserina. Genetics 166: 1241–1252. doi: 10.1534/genetics.166.3.1241
- 91. Kicka S, Bonnet C, Sobering AK, Ganesan LP, Silar P (2006) A mitotically inheritable unit containing a MAP kinase module. Proc Natl Acad Sci U S A 103: 13445–13450. doi: 10.1073/pnas.0603693103
- 92. Kjaerulff S, Lautrup-Larsen I, Truelsen S, Pedersen M, Nielsen O (2005) Constitutive activation of the fission yeast pheromone-responsive pathway induces ectopic meiosis and reveals ste11 as a mitogen-activated protein kinase target. Mol Cell Biol 25: 2045–2059. doi: 10.1128/mcb.25.5.2045-2059.2005
- 93. Berbee ML, Taylor JW (2010) Dating the molecular clock in fungi - how close are we? Fungal Biol Rev 24: 1–16. doi: 10.1016/j.fbr.2010.03.001
- 94. Rizet G, Engelmann C (1949) Contribution à l'étude génétique d'un ascomycète tatrasporé: Podospora anserina. Rhem Rv Cytol Biol Veg 11: 201–304.
- 95. Rizet G (1952) Les phénomènes de barrage chez Podospora anserina. I. Analyse génétique des barrages entre souches S and s. Rev Cytol Biol Veg 13: 51–92.
- 96. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, et al. (2012) The Pfam protein families database. Nucleic Acids Res 40: D290–D301. doi: 10.1093/nar/gkr1065
- 97. Colot HV, Park G, Turner GE, Ringelberg C, Crew CM, et al. (2006) A high-throughput gene knockout procedure for Neurospora reveals functions for multiple transcription factors. Proc Natl Acad Sci U S A 103: 10352–10357. doi: 10.1073/pnas.0601456103
- 98. El-Khoury R, Sellem CH, Coppin E, Boivin A, Maas MF, et al. (2008) Gene deletion and allelic replacement in the filamentous fungus Podospora anserina. Curr Genet 53: 249–258. doi: 10.1007/s00294-008-0180-3
- 99. Silar P (1995) Two new easy to use vectors for transformations . Fungal Genet Newslett 42: 73.
- 100. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3: RESEARCH0034.
- 101. Andersen CL, Jensen JL, Orntoft TF (2004) Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res 64: 5245–5250. doi: 10.1158/0008-5472.can-04-0496
- 102. Pfaffl MW, Horgan GW, Dempfle L (2002) Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res 30: e36. doi: 10.1093/nar/30.9.e36
- 103. van Houten V, Denkers F, van Dijk M, van den Brekel M, Brakenhoff R (1998) Labeling efficiency of oligonucleotides by T4 polynucleotide kinase depends on 5′-nucleotide. Anal Biochem 265: 386–389. doi: 10.1006/abio.1998.2900
- 104. Holden NS, Tacon CE (2011) Principles and problems of the electrophoretic mobility shift assay. J Pharmacol Toxicol Methods 63: 7–14. doi: 10.1016/j.vascn.2010.03.002
- 105. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, et al. (2007) Clustal W and Clustal X version 2.0. Bioinformatics 23: 2947–2948. doi: 10.1093/bioinformatics/btm404
- 106. Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ (2009) Jalview Version 2-a multiple sequence alignment editor and analysis workbench. Bioinformatics 25: 1189–1191. doi: 10.1093/bioinformatics/btp033
- 107. Le SQ, Gascuel O (2008) An improved general amino acid replacement matrix. Mol Biol Evol 25: 1307–1320. doi: 10.1093/molbev/msn067
- 108. Pavlidis P, Noble WS (2003) Matrix2png: a utility for visualizing matrix data. Bioinformatics 19: 295–296. doi: 10.1093/bioinformatics/19.2.295
- 109. Philley ML, Staben C (1994) Functional analyses of the Neurospora crassa MT a-1 mating type polypeptide. Genetics 137: 715–722.