we report here five newly sequenced and assembled penicillium genomes , which we compared with previously available data [ 1619 ] .
the full dataset included the genome sequences of ten penicillium species , six of which are either used as industrial inocula for cheese making ( penicillium roqueforti and penicillium camemberti ) or occur as contaminants in cheeses ( figure  1 ; table s1 ) .
camemberti is only found in cheese and is thought to include a single clonal lineage originating from selection programs at the end of the 19 century from the blue - gray cheese molds used at that time , i.e. , penicillium biforme and penicillium fuscoglaucum [ 20 , 21 ] .
by contrast , p.  roqueforti also occurs in habitats other than cheese , such as silage or wood , and displays substantial genetic diversity [ 22 , 23 ] . for reconstructing a rooted phylogeny of these ten penicillium species , we used four aspergillus species as an outgroup .
we concatenated alignments of 3,089 single - copy genes shared by at least ten species and reconstructed a fully resolved and well - supported maximum - likelihood phylogeny ( figures 1a and 1b ) .
we used this rooted phylogeny to investigate the occurrence of horizontal gene transfers ( hgts ) between penicillium species .
as hgts ( also known as xenology ) result in incongruences between gene genealogies and species trees , we compared all individual gene genealogies with the species tree . for this goal , we used the notung software [ 2527 ] to infer the number of duplication , loss , and hgt events that reconciled the gene genealogies with the species tree .
notung is conservative regarding the inference of hgts because it tests their temporal feasibility , assumes that hgts occur with a low probability , and forces the poorly supported nodes to follow the species tree .
only orthologous groups with at least one homolog in at least eight genomes were analyzed , further rendering our estimates of hgts a lower bound .
notung inferred the highest number of hgts relative to branch length in the clade encompassing p.  camemberti , p.  biforme , and their common ancestor ( figure  1b ) .
8 of the 21 horizontally acquired genes detected in p.  roqueforti were inferred to come from p.  camemberti , p.  biforme , or their common ancestor , indicating recent transfers from species sharing the same ecological niche ( figure  s1 ) .
only five of these eight genes could be assigned putative functions , i.e. , a protein kinase , two transcription factors , a cation transporter , and an integrase - like protein .
cation transport seems particularly relevant for growth in cheeses as several ions are limiting in this medium ( e.g. , iron ions ) , and ph drastically drops during cheese maturation .
another line of evidence for the abundance of hgts in penicillium fungi came from the finding of multiple large genomic islands that were almost 100% identical at the nucleotidic level between distant species , while being absent from closely related species ( figure  1a ) .
the only substitutions detected in these genomic islands corresponded to repeat - induced point mutations , i.e. , c - to - t transitions induced by a specific fungal defense mechanism against transposable elements ( tes ) that can substitute multiple base pairs in a single meiosis . in p.  roqueforti , for example , seven genomic islands larger than 10 kb and displaying above 97% nucleotide identity with multiple other species were found .
such a high level of identity suggests that these genomic islands correspond to recent horizontally transferred regions ( htrs ) , although they could alternatively be recent introgressions .
two lines of evidence , however , support the htr hypothesis rather than introgression : ( 1 ) the presence of several of these regions at non - homologous locations in the different penicillium genomes ( figure  s2 ; ) and ( 2 ) the low mean genome sequence identity between the penicillium species sharing these regions , being less than 90% , an identity level at which no successful interspecific crosses have ever been reported in fungi .
interestingly , these htrs were flanked in p.  roqueforti by copies of tes from a particular family that were rare elsewhere in the genomes ( figures 1a and s3 ) , the i non - ltr retrotransposons .
the other abundant tes ( e.g. , mariner dna transposons and copia retrotransposons ) were in contrast scattered in the genomes ( figure  1a ) .
the genes present in these genomic islands of high sequence similarity partially overlapped with the hgts detected by notung . in p. 
roqueforti for example , 17% of the hgts inferred by notung were located in the seven large htrs ( figure  1a ) .
the fact that not all genes in htrs were detected by notung mainly results from the filter of this analysis where we used only orthologous groups with homologs in at least eight species .
further , 11% of the inferred hgts in p.  roqueforti clustered within 50 kb of the htrs ( figure  1a ) , suggesting that these genomic regions may be prone to integrate foreign dna .
this is consistent with the previous finding that the genomic region where wallaby is inserted in p.  roqueforti carries other species - specific genes in each of p.  camemberti , p.  rubens , and p.  roqueforti .
the identification of multiple very recent htrs , with almost 100% identity in multiple species ( figure  1a ) , together with notung inferences of horizontally transferred genes occurring also elsewhere in the genomes ( at least 104 in total among penicillium fungi ) , indicates that hgts occur frequently among penicillium .
the clustering of i elements in the flanking regions of htrs suggests that they may be involved in the mechanism of horizontal transfers .
it has been shown that tes can pass across species boundaries   and that they promote genomic rearrangements and recombination [ 3234 ] .
the capacity for mycelia fusions may also facilitate the exchange of genetic material in fungi .
five of the seven large htrs detected in p.  roqueforti were shared between penicillium strains isolated from cheese , p.  camemberti carrying four of them ( figure  1a ) .
wallaby   carries a gene encoding an antifungal protein , known for inhibiting the growth of competitors .
the second largest htr was found at terminal edges of scaffolds and was therefore named cheesyter .
this 80-kb region , found as a single block in all the genomes studied , carried 37 putative genes , among which two had relevant putative functions for adaptation to cheese , i.e. , lactose permease and beta - galactosidase ( figure  2a ) .
lactose is present during the first few days of cheese maturation , and it acts as a primary carbon source , being rapidly consumed by lactic acid bacteria .
the two lactose metabolism genes present in cheesyter were among the most strongly expressed in p.  camemberti during the first step of cheese rind maturation , during which lactose is available ( figure  2b ; table s2 ; supplemental experimental procedures ) , indicating a role in the use of the cheese substrate .
we investigated the presence of cheesyter and wallaby by  pcr in 416 strains from 65 fungal species from various environments ( figure  1c ; table s3 ; supplemental experimental procedures ) .
the presence of the transfers was found highly significantly associated with dairy environment both among species (   = 55.7 ; degrees of freedom [ df ]  = 1 ; p value  = 8.571e14 ) and among strains within species ( 2  = 45.2 ; df  = 1 , p value  = 1.774e11 ) .
amplicons were actually obtained only for penicillium species that are frequently isolated in the dairy environment , with the only exception of p.  rubens , the penicillin - producer fungus , in which 16 strains out of 20 carried either one of the two htrs .
the cheesyter and wallaby fragments obtained by pcr showed zero substitution among all strains from all species , including synonymous sites and non - coding regions , as previously found for wallaby .
this result confirms that the presence of cheesyter and wallaby is not ancestral in penicillium species and that these genomic islands have instead been acquired very recently .
this also indicates that the two htrs have spread in several species through recent selective sweeps .
p.  roqueforti was found polymorphic for the presence of wallaby and cheesyter , with all tested strains carrying either both of these regions or neither of them ( table s3 ) . within p. roqueforti , these regions were present only in strains isolated from the cheese environment , suggesting a role in adaptation to the cheese environment .
we therefore investigated whether strains carrying wallaby and cheesyter showed higher fitness in terms of growth on cheese substrate or for competitiveness .
we set up three experiments , focusing on p.  roqueforti , because a large collection of strains was available , isolated from various environments , and including strains carrying both wallaby and cheesyter ( hereafter named w+c+ ) and strains lacking them ( hereafter named wc ) .
we first compared the growth of 50 p.  roqueforti strains on a cheese medium and on a minimal medium ( 26 w+c+ and 24 wc ; table s4 , tab a ) .
neither the presence of wallaby and cheesyter , as a main effect independent of the medium , nor the origin of the strain ( i.e. , cheese versus other environments ) significantly influenced the growth of p.  roqueforti ( table s1 ) .
by contrast , the effect on growth of the medium and its interaction with the presence of the two genomic islands were significant ( figure  3 ; table s1 ) : w+c+ strains had a growth advantage on cheese medium but a slower growth on minimal medium .
second , we investigated whether p.  roqueforti strains carrying wallaby and cheesyter had a higher ability to exclude competitors .
we measured the growth of three fungal strains belonging to species commonly found in cheese but lacking wallaby and cheesyter ( p.  nalgiovense fm193 , p.  biforme lcp05529 , and geotrichum candidum fm074 ) on plates covered with lawns of p.  roqueforti either w+c+ ( n  = 11 ) or wc ( n  = 12 ) strains ( table s4 , tab b ) .
these experiments were carried out on minimal , cheese , and malt agar media . no difference in growth
was detected for the yeast g.  candidum between lawns of w+c+ or wc p.  roqueforti strains ( table s1 ) . by contrast , w+c+ p.  roqueforti strains significantly impaired the growth of the two penicillium challengers on the cheese and malt media ( figure  3b ) .
this was not the case on the minimal medium : the interaction between the presence of the transfers and the medium was significant ( table s1 ) . using the same experimental design , we then investigated the effect of the two genomic islands when present in the challengers . for this goal , we inoculated on p.  roqueforti lawns ( w+c+ , n  = 2 , or wc ,
n  = 2 ) , on cheese medium , different strains of species displaying a polymorphism in the wallaby and/or cheesyter presence ( table s4 , tab c ) .
we used as challengers different strains of p.  camemberti ( w+c , n  = 1 , or w+c+ , n  = 3 ) , p.  biforme ( wc , n  = 2 , or w+c+ , n  = 2 ) , and p.  rubens ( wc , n  = 1 , or w+c ,
, we found that the  p.  roqueforti lawns significantly inhibited the growth of challengers and significantly more so when the p.  roqueforti lawn carried wallaby and cheesyter .
interestingly , the presence of either cheesyter or wallaby in the challengers allowed better  growth on wc p. 
roqueforti lawns while neither had significant effect on the growth on w+c+ p.  roqueforti lawns ( table s1 ) .
third , we investigated competition among p.  roqueforti strains carrying ( w+c+ , n  = 8) or lacking ( wc , n  = 11 ) wallaby and cheesyter .
we grew p.  roqueforti strains on cheese medium as pairwise face - to - face confrontations , and we measured the deviations from symmetrical growth ( table s4 , tab d ; figure  3c ; supplemental experimental procedures ) . for the w+c+ versus w+c+ confrontations ,
the mean growth deviation from a boundary in the exact middle of the petri dish was not significantly different from zero ( t test , t  = 1.5 ; df  = 36 ; p value  = 0.14 ) .
similar results were obtained for the wc versus wc confrontations ( t test , t  = 0.25 ; df  = 70 ; p value  = 0.80 ) .
for the w+c+ versus wc confrontations , deviations were measured by taking the w+c+ strain as the focal strain ; the mean growth deviation was significantly different from zero and positive , the w+c+ strains thus growing farther than the wc strain ( t  test , t  = 12.32 ; df  = 90 ; p value  < 0.0001 ) .
the mean deviations were significantly higher in the wc versus w+c+ confrontations than in the w+c+ versus w+c+ or wc versus wc confrontations ( tukey - kramer test , p value  < 0.0001 ) , while the means between these two latter were not significantly different .
this experiment shows that the competitive advantage of w+c+ strains against wc strains also holds within the species p.  roqueforti .
altogether , these experimental results strongly support the existence of fitness advantages for the penicillium strains carrying the horizontally transferred genomic islands , both in the use of cheese substrate and in competition with fungal competitors .
the two cheese species studied  here underwent parallel adaptation to the cheese medium , and this involved the transfers of identical regions across species boundaries .
hgt events have been reported in fungi [ 35 , 38 , 39 ] , particularly in environments created by humans , in domesticated yeasts and fungal pathogens of crops [ 10 , 40 , 41 ] .
the extent and timing of gene transfers and the number of species having received the same htrs are here particularly striking .
furthermore , we provide experimental evidence of fitness advantages for strains carrying these htrs on a human - made medium . these findings altogether are potentially useful for guiding modern strain improvement programs . indeed , together with the protocol for inducing sex in p.  roqueforti [ 22 , 42 ] , the identification here of several key candidate genes important for cheese metabolism and competition may allow further selecting interesting traits for cheese industry using the great genetic variability present in p.  roqueforti strains without wallaby or cheesyter .
in addition , our results suggest that caution is required concerning the introduction of genes into microorganisms , as these genes could readily be transferred to other species in the food environment .
indeed , the rapid spread of wallaby and cheesyter into many species of the dairy environment , even when occurring only as contaminants , indicates that transgenes may readily cross species boundaries in the food chain . finally , the findings here of rapid adaptation through frequent horizontal gene transfers among distant species under selection in novel , human - made media contribute to our understanding of the evolutionary genomic mechanisms allowing rapid adaptation to environmental changes in eukaryotes .
; investigation , a.b . , m.l .- v . , r.c.r.d.l.v . , and j.r .
 summarydomestication is an excellent model for studies of adaptation because it involves recent and strong selection on a few , identified traits [ 15 ] . 
 few studies have focused on the domestication of fungi , with notable exceptions [ 611 ] , despite their importance to bioindustry [ 12 ] and to a general understanding of adaptation in eukaryotes [ 5 ] . 
 penicillium fungi are ubiquitous molds among which two distantly related species have been independently selected for cheese making  
 p .  
 roqueforti for blue cheeses like roquefort and p.  camemberti for soft cheeses like camembert . 
 the selected traits include morphology , aromatic profile , lipolytic and proteolytic activities , and ability to grow at low temperatures , in a matrix containing bacterial and fungal competitors [ 1315 ] . by comparing the genomes of ten penicillium species , 
 we show that adaptation to cheese was associated with multiple recent horizontal transfers of large genomic regions carrying crucial metabolic genes . 
 we identified seven horizontally transferred regions ( htrs ) spanning more than 10 kb each , flanked by  specific transposable elements , and displaying nearly 100% identity between distant penicillium species . 
 two htrs carried genes with functions involved in the utilization of cheese nutrients or competition and were found nearly identical in multiple strains and species of cheese - associated penicillium fungi , indicating recent selective sweeps ; they were experimentally associated with faster growth and greater competitiveness on cheese and contained genes highly expressed in the early stage of cheese maturation . 
 these findings have industrial and food safety implications and improve our understanding of the processes of adaptation to rapid environmental changes . 