Natural selection on coding and noncoding DNA sequences is associated with virulence genes in a plant pathogenic fungus

Optimizing bacteriophage treatment of resistant Pseudomonas

Phage therapy is increasing in relevance as an alternative treatment to combat antibiotic resistant bacteria. Phage cocktails are the state-of-the-art method of administering phages in clinical settings, preferred over monophage treatment because of their ability to eliminate multiple bacterial strains and reduce resistance formation. In our study, we compare monophage applications and phage cocktails to our chosen method of phage sequential treatments. To do so, we isolated four novel bacteriophages capable of infecting Pseudomonas alcaligenes T3, a close relative of P. aeruginosa, and characterized them using sequencing and transmission electron microscopy. While investigating monophage treatments, we observed that different phage concentrations had a strong impact on the timing and amount of resistance formation. When using phage cocktails, we observed that P. alcaligenes were capable of forming resistance in the same timespan it took them to become resistant to single phages. We isolated mutants resistant to each single phage as well as mutants exposed to phage cocktails, resulting in bacteria resistant to all four phages at once. Sequencing these mutants showed that different treatments yielded unique single nucleotide polymorphism mutation patterns. In order to combat resistance formation, we added phages one by one in intervals of 24 h, thus managing to delay resistance development and keeping bacterial growth significantly lower compared to phage cocktails.IMPORTANCEWHO declared antimicrobial resistance a top threat to global health; while antibiotics have stood at the forefront in the fight against bacterial infection, the increasing number of multidrug-resistant bacteria highlights a need to branch out in order to address the threat of antimicrobial resistance. Bacteriophages, viruses solely infecting bacteria, could present a solution due to their abundance, versatility, and adaptability. For this study, we isolated new phages infecting a fast-mutating Pseudomonas alcaligenes strain capable of forming resistance within 30 h. By using a sequential treatment approach of adding one phage after another, we were able to curb bacterial growth significantly more compared to state-of-the-art phage cocktails.

Population Genomics Provide Insights into the Global Genetic Structure of Colletotrichum graminicola, the Causal Agent of Maize Anthracnose

Understanding the genetic diversity and mechanisms underlying genetic variation in pathogen populations is crucial to the development of effective control strategies. We investigated the genetic diversity and reproductive biology of Colletotrichum graminicola isolates which infect maize by sequencing the genomes of 108 isolates collected from 14 countries using restriction site-associated DNA sequencing (RAD-seq) and whole-genome sequencing (WGS). Clustering analyses based on single-nucleotide polymorphisms revealed three genetic groups delimited by continental origin, compatible with short-dispersal of the pathogen and geographic subdivision. Intra- and intercontinental migration was observed between Europe and South America, likely associated with the movement of contaminated germplasm. Low clonality, evidence of genetic recombination, and high phenotypic diversity were detected. We show evidence that, although it is rare (possibly due to losses of sexual reproduction- and meiosis-associated genes) C. graminicola can undergo sexual recombination. Our results support the hypotheses that intra- and intercontinental pathogen migration and genetic recombination have great impacts on the C. graminicola population structure. IMPORTANCE Plant pathogens cause significant reductions in yield and crop quality and cause enormous economic losses worldwide. Reducing these losses provides an obvious strategy to increase food production without further degrading natural ecosystems; however, this requires knowledge of the biology and evolution of the pathogens in agroecosystems. We employed a population genomics approach to investigate the genetic diversity and reproductive biology of the maize anthracnose pathogen (Colletotrichum graminicola) in 14 countries. We found that the populations are correlated with their geographical origin and that migration between countries is ongoing, possibly caused by the movement of infected plant material. This result has direct implications for disease management because migration can cause the movement of more virulent and/or fungicide-resistant genotypes. We conclude that genetic recombination is frequent (in contrast to the traditional view of C. graminicola being mainly asexual), which strongly impacts control measures and breeding programs aimed at controlling this disease.

Maize Anthracnose Stalk Rot in the Genomic Era

Anthracnose stalk rot (ASR) of maize results in millions of dollars in losses annually in the United States. ASR, together with anthracnose leaf blight and anthracnose top dieback, is caused by the fungus Colletotrichum graminicola. Current ASR management recommendations emphasize host resistance and reduction of plant stressors (e.g., drought, heat, low fertility, or soil acidity). Stress reduction may be more difficult to achieve in the future due to more high-intensity production protocols and climate change. Moreover, cultural and chemical management practices may conflict with other important goals, including environmental sustainability and maximization of yield potential. Thus, future ASR management may rely more heavily on host resistance, for which there are relatively few highly effective sources. The last comprehensive review of C. graminicola and maize anthracnose was written over two decades ago. The genomic age has brought important new insights into mechanisms governing the host-pathogen interaction from the application of molecular and cytological technologies. This review provides a summary of our current model of maize anthracnose etiology, including how increased knowledge of molecular and cellular events could contribute to better ASR management. Improved understanding of C. graminicola taxonomy has confirmed that the fungus is specific to Zea mays, and that it colonizes living maize tissues via a critical biotrophic phase. Successful biotrophic establishment relies on an array of secreted protein effectors and secondary metabolites produced at different stages of infection and dispersed to multiple locations. These molecules could provide therapeutic targets for the next generation of transgenic or gene-edited ASR-resistant hybrids.

Proteinaceous Effector Discovery and Characterization in Plant Pathogenic Colletotrichum Fungi

Anthracnose caused by plant pathogenic Colletotrichum fungi results in large economic losses in field crop production worldwide. To aid the establishment of plant host infection, Colletotrichum pathogens secrete numerous effector proteins either in apoplastic space or inside of host cells for effective colonization. Understanding these effector repertoires is critical for developing new strategies for resistance breeding and disease management. With the advance of genomics and bioinformatics tools, a large repertoire of putative effectors has been identified in Colletotrichum genomes, and the biological functions and molecular mechanisms of some studied effectors have been summarized. Here, we review recent advances in genomic identification, understanding of evolutional characteristics, transcriptional profiling, and functional characterization of Colletotrichum effectors. We also offer a perspective on future research.

Fungal Pathogenesis-Related Cell Wall Biogenesis, with Emphasis on the Maize Anthracnose Fungus Colletotrichum graminicola

The genus Colletotrichum harbors many plant pathogenic species, several of which cause significant yield losses in the field and post harvest. Typically, in order to infect their host plants, spores germinate, differentiate a pressurized infection cell, and display a hemibiotrophic lifestyle after plant invasion. Several factors required for virulence or pathogenicity have been identified in different Colletotrichum species, and adaptation of cell wall biogenesis to distinct stages of pathogenesis has been identified as a major pre-requisite for the establishment of a compatible parasitic fungus-plant interaction. Here, we highlight aspects of fungal cell wall biogenesis during plant infection, with emphasis on the maize leaf anthracnose and stalk rot fungus, Colletotrichum graminicola.

Comparative Genomics and Gene Pool Analysis Reveal the Decrease of Genome Diversity and Gene Number in Rice Blast Fungi by Stable Adaption with Rice

Magnaporthe oryzae caused huge losses in rice and wheat production worldwide. Comparing to long-term co-evolution history with rice, wheat-infecting isolates were new-emerging. To reveal the genetic differences between rice and wheat blast on global genomic scale, 109 whole-genome sequences of M. oryzae from rice, wheat, and other hosts were reanalyzed in this study. We found that the rice lineage had gone through stronger selective sweep and fewer conserved genes than those of Triticum and Lolium lineages, which indicated that rice blast fungi adapted to rice by gene loss and rapid evolution of specific loci. Furthermore, 228 genes associated with host adaptation of M. oryzae were found by presence/absence variation (PAV) analyses. The functional annotation of these genes found that the fine turning of genes gain/loss involved with transport and transcription factor, thiol metabolism, and nucleotide metabolism respectively are major mechanisms for rice adaption. This result implies that genetic base of specific host plant may lead to gene gain/loss variation of pathogens, so as to enhance their adaptability to host. Further characterization of these specific loci and their roles in adaption and evaluation of the fungi may eventually lead to understanding of interaction mechanism and develop new strategies of the disease management.

Divergence of a genomic island leads to the evolution of melanization in a halophyte root fungus

Understanding how organisms adapt to extreme living conditions is central to evolutionary biology. Dark septate endophytes (DSEs) constitute an important component of the root mycobiome and they are often able to alleviate host abiotic stresses. Here, we investigated the molecular mechanisms underlying the beneficial association between the DSE Laburnicola rhizohalophila and its host, the native halophyte Suaeda salsa, using population genomics. Based on genome-wide Fst (pairwise fixation index) and Vst analyses, which compared the variance in allele frequencies of single-nucleotide polymorphisms (SNPs) and copy number variants (CNVs), respectively, we found a high level of genetic differentiation between two populations. CNV patterns revealed population-specific expansions and contractions. Interestingly, we identified a ~20 kbp genomic island of high divergence with a strong sign of positive selection. This region contains a melanin-biosynthetic polyketide synthase gene cluster linked to six additional genes likely involved in biosynthesis, membrane trafficking, regulation, and localization of melanin. Differences in growth yield and melanin biosynthesis between the two populations grown under 2% NaCl stress suggested that this genomic island contributes to the observed differences in melanin accumulation. Our findings provide a better understanding of the genetic and evolutionary mechanisms underlying the adaptation to saline conditions of the L. rhizohalophila-S. salsa symbiosis.

Peroxiredoxin alleviates the fitness costs of imidacloprid resistance in an insect pest of rice

Chemical insecticides have been heavily employed as the most effective measure for control of agricultural and medical pests, but evolution of resistance by pests threatens the sustainability of this approach. Resistance-conferring mutations sometimes impose fitness costs, which may drive subsequent evolution of compensatory modifier mutations alleviating the costs of resistance. However, how modifier mutations evolve and function to overcome the fitness cost of resistance still remains unknown. Here we show that overexpression of P450s not only confers imidacloprid resistance in the brown planthopper, Nilaparvata lugens, the most voracious pest of rice, but also leads to elevated production of reactive oxygen species (ROS) through metabolism of imidacloprid and host plant compounds. The inevitable production of ROS incurs a fitness cost to the pest, which drives the increase or fixation of the compensatory modifier allele T65549 within the promoter region of N. lugens peroxiredoxin (NlPrx) in the pest populations. T65549 allele in turn upregulates the expression of NlPrx and thus increases resistant individuals' ability to clear the cost-incurring ROS of any source. The frequent involvement of P450s in insecticide resistance and their capacity to produce ROS while metabolizing their substrates suggest that peroxiredoxin or other ROS-scavenging genes may be among the common modifier genes for alleviating the fitness cost of insecticide resistance.

Mosaicism in Human Health and Disease

Mosaicism refers to the occurrence of two or more genomes in an individual derived from a single zygote. Germline mosaicism is a mutation that is limited to the gonads and can be transmitted to offspring. Somatic mosaicism is a postzygotic mutation that occurs in the soma, and it may occur at any developmental stage or in adult tissues. Mosaic variation may be classified in six ways: (a) germline or somatic origin, (b) class of DNA mutation (ranging in scale from single base pairs to multiple chromosomes), (c) developmental context, (d) body location(s), (e) functional consequence (including deleterious, neutral, or advantageous), and (f) additional sources of mosaicism, including mitochondrial heteroplasmy, exogenous DNA sources such as vectors, and epigenetic changes such as imprinting and X-chromosome inactivation. Technological advances, including single-cell and other next-generation sequencing, have facilitated improved sensitivity and specificity to detect mosaicism in a variety of biological contexts.

No Evidence for Recent Selection at FOXP2 among Diverse Human Populations

FOXP2, initially identified for its role in human speech, contains two nonsynonymous substitutions derived in the human lineage. Evidence for a recent selective sweep in Homo sapiens, however, is at odds with the presence of these substitutions in archaic hominins. Here, we comprehensively reanalyze FOXP2 in hundreds of globally distributed genomes to test for recent selection. We do not find evidence of recent positive or balancing selection at FOXP2. Instead, the original signal appears to have been due to sample composition. Our tests do identify an intronic region that is enriched for highly conserved sites that are polymorphic among humans, compatible with a loss of function in humans. This region is lowly expressed in relevant tissue types that were tested via RNA-seq in human prefrontal cortex and RT-PCR in immortalized human brain cells. Our results represent a substantial revision to the adaptive history of FOXP2, a gene regarded as vital to human evolution.

The diversity, evolution and epidemiology of plant viruses: A phylogenetic view

During the past four decades, the scientific community has seen an exponential advance in the number, sophistication, and quality of molecular techniques and bioinformatics tools for the genetic characterization of plant virus populations. Predating these advances, the field of Phylogenetics has significantly contributed to understand important aspects of plant virus evolution. This review aims at summarizing the impact of Phylogenetics in the current knowledge on three major aspects of plant virus evolution that have benefited from the development of phylogenetic inference: (1) The identification and classification of plant virus diversity. (2) The mechanisms and forces shaping the evolution of plant virus populations. (3) The understanding of the interaction between plant virus evolution, epidemiology and ecology. The work discussed here highlights the important role of phylogenetic approaches in the study of the dynamics of plant virus populations.

Mind the gap! The mitochondrial control region and its power as a phylogenetic marker in echinoids

BACKGROUND

In Metazoa, mitochondrial markers are the most commonly used targets for inferring species-level molecular phylogenies due to their extremely low rate of recombination, maternal inheritance, ease of use and fast substitution rate in comparison to nuclear DNA. The mitochondrial control region (CR) is the main non-coding area of the mitochondrial genome and contains the mitochondrial origin of replication and transcription. While sequences of the cytochrome oxidase subunit 1 (COI) and 16S rRNA genes are the prime mitochondrial markers in phylogenetic studies, the highly variable CR is typically ignored and not targeted in such analyses. However, the higher substitution rate of the CR can be harnessed to infer the phylogeny of closely related species, and the use of a non-coding region alleviates biases resulting from both directional and purifying selection. Additionally, complete mitochondrial genome assemblies utilizing next generation sequencing (NGS) data often show exceptionally low coverage at specific regions, including the CR. This can only be resolved by targeted sequencing of this region.

RESULTS

Here we provide novel sequence data for the echinoid mitochondrial control region in over 40 species across the echinoid phylogenetic tree. We demonstrate the advantages of directly targeting the CR and adjacent tRNAs to facilitate complementing low coverage NGS data from complete mitochondrial genome assemblies. Finally, we test the performance of this region as a phylogenetic marker both in the lab and in phylogenetic analyses, and demonstrate its superior performance over the other available mitochondrial markers in echinoids.

CONCLUSIONS

Our target region of the mitochondrial CR (1) facilitates the first thorough investigation of this region across a wide range of echinoid taxa, (2) provides a tool for complementing missing data in NGS experiments, and (3) identifies the CR as a powerful, novel marker for phylogenetic inference in echinoids due to its high variability, lack of selection, and high compatibility across the entire class, outperforming conventional mitochondrial markers.

Gene cluster conservation provides insight into cercosporin biosynthesis and extends production to the genus Colletotrichum

Species in the fungal genus Cercospora cause diseases in many important crops worldwide. Their success as pathogens is largely due to the secretion of cercosporin during infection. We report that the cercosporin toxin biosynthesis (CTB) gene cluster is ancient and was horizontally transferred to diverse fungal plant pathogens. Because our analyses revealed genes adjacent to the established CTB cluster with similar evolutionary trajectories, we evaluated their role in Cercospora beticola to show that four are necessary for cercosporin biosynthesis. Lastly, we confirmed that the apple pathogen Colletotrichum fioriniae produces cercosporin, the first case outside the family Mycosphaerellaceae. Other Colletotrichum plant pathogens also harbor the CTB cluster, which points to a wider role that this toxin may play in virulence.

Species in the genus Cercospora cause economically devastating diseases in sugar beet, maize, rice, soy bean, and other major food crops. Here, we sequenced the genome of the sugar beet pathogen Cercospora beticola and found it encodes 63 putative secondary metabolite gene clusters, including the cercosporin toxin biosynthesis (CTB) cluster. We show that the CTB gene cluster has experienced multiple duplications and horizontal transfers across a spectrum of plant pathogenic fungi, including the wide-host range Colletotrichum genus as well as the rice pathogen Magnaporthe oryzae. Although cercosporin biosynthesis has been thought to rely on an eight-gene CTB cluster, our phylogenomic analysis revealed gene collinearity adjacent to the established cluster in all CTB cluster-harboring species. We demonstrate that the CTB cluster is larger than previously recognized and includes cercosporin facilitator protein, previously shown to be involved with cercosporin autoresistance, and four additional genes required for cercosporin biosynthesis, including the final pathway enzymes that install the unusual cercosporin methylenedioxy bridge. Lastly, we demonstrate production of cercosporin by Colletotrichum fioriniae, the first known cercosporin producer within this agriculturally important genus. Thus, our results provide insight into the intricate evolution and biology of a toxin critical to agriculture and broaden the production of cercosporin to another fungal genus containing many plant pathogens of important crops worldwide.

Intraspecific comparative genomics of isolates of the Norway spruce pathogen (Heterobasidion parviporum) and identification of its potential virulence factors

BACKGROUND

Heterobasidion parviporum is an economically most important fungal forest pathogen in northern Europe, causing root and butt rot disease of Norway spruce (Picea abies (L.) Karst.). The mechanisms underlying the pathogenesis and virulence of this species remain elusive. No reference genome to facilitate functional analysis is available for this species.

RESULTS

To better understand the virulence factor at both phenotypic and genomic level, we characterized 15 H. parviporum isolates originating from different locations across Finland for virulence, vegetative growth, sporulation and saprotrophic wood decay. Wood decay capability and latitude of fungal origins exerted interactive effects on their virulence and appeared important for H. parviporum virulence. We sequenced the most virulent isolate, the first full genome sequences of H. parviporum as a reference genome, and re-sequenced the remaining 14 H. parviporum isolates. Genome-wide alignments and intrinsic polymorphism analysis showed that these isolates exhibited overall high genomic similarity with an average of at least 96% nucleotide identity when compared to the reference, yet had remarkable intra-specific level of polymorphism with a bias for CpG to TpG mutations. Reads mapping coverage analysis enabled the classification of all predicted genes into five groups and uncovered two genomic regions exclusively present in the reference with putative contribution to its higher virulence. Genes enriched for copy number variations (deletions and duplications) and nucleotide polymorphism were involved in oxidation-reduction processes and encoding domains relevant to transcription factors. Some secreted protein coding genes based on the genome-wide selection pressure, or the presence of variants were proposed as potential virulence candidates.

CONCLUSION

Our study reported on the first reference genome sequence for this Norway spruce pathogen (H. parviporum). Comparative genomics analysis gave insight into the overall genomic variation among this fungal species and also facilitated the identification of several secreted protein coding genes as putative virulence factors for the further functional analysis. We also analyzed and identified phenotypic traits potentially linked to its virulence.

Ceratocystis cacaofunesta genome analysis reveals a large expansion of extracellular phosphatidylinositol-specific phospholipase-C genes (PI-PLC)

BACKGROUND

The Ceratocystis genus harbors a large number of phytopathogenic fungi that cause xylem parenchyma degradation and vascular destruction on a broad range of economically important plants. Ceratocystis cacaofunesta is a necrotrophic fungus responsible for lethal wilt disease in cacao. The aim of this work is to analyze the genome of C. cacaofunesta through a comparative approach with genomes of other Sordariomycetes in order to better understand the molecular basis of pathogenicity in the Ceratocystis genus.

RESULTS

We present an analysis of the C. cacaofunesta genome focusing on secreted proteins that might constitute pathogenicity factors. Comparative genome analyses among five Ceratocystidaceae species and 23 other Sordariomycetes fungi showed a strong reduction in gene content of the Ceratocystis genus. However, some gene families displayed a remarkable expansion, in particular, the Phosphatidylinositol specific phospholipases-C (PI-PLC) family. Also, evolutionary rate calculations suggest that the evolution process of this family was guided by positive selection. Interestingly, among the 82 PI-PLCs genes identified in the C. cacaofunesta genome, 70 genes encoding extracellular PI-PLCs are grouped in eight small scaffolds surrounded by transposon fragments and scars that could be involved in the rapid evolution of the PI-PLC family. Experimental secretome using LC-MS/MS validated 24% (86 proteins) of the total predicted secretome (342 proteins), including four PI-PLCs and other important pathogenicity factors.

CONCLUSION

Analysis of the Ceratocystis cacaofunesta genome provides evidence that PI-PLCs may play a role in pathogenicity. Subsequent functional studies will be aimed at evaluating this hypothesis. The observed genetic arsenals, together with the analysis of the PI-PLC family shown in this work, reveal significant differences in the Ceratocystis genome compared to the classical vascular fungi, Verticillium and Fusarium. Altogether, our analyses provide new insights into the evolution and the molecular basis of plant pathogenicity.

MM, Prakash Mishra G, Singh D, Ye W, Tyler BM, Tripathy S. EumicrobeDBLite: a lightweight genomic resource and analytic platform for draft oomycete genomes

We have developed EumicrobeDBLite-a lightweight comprehensive genome resource and sequence analysis platform for oomycete organisms. EumicrobeDBLite is a successor of the VBI Microbial Database (VMD) that was built using the Genome Unified Schema (GUS). In this version, GUS has been greatly simplified with the removal of many obsolete modules and the redesign of others to incorporate contemporary data. Several dependences, such as perl object layers used for data loading in VMD, have been replaced with independent lightweight scripts. EumicrobeDBLite now runs on a powerful annotation engine developed at our laboratory, called 'Genome Annotator Lite'. Currently, this database has 26 publicly available genomes and 10 expressed sequence tag (EST) datasets of oomycete organisms. The browser page has dynamic tracks presenting comparative genomics analyses, coding and non-coding data, tRNA genes, repeats and EST alignments. In addition, we have defined 44 777 core conserved proteins from 12 oomycete organisms which form 2974 clusters. Synteny viewing is enabled by the incorporation of the Genome Synteny Viewer (GSV) tool. The user interface has undergone major changes for ease of browsing. Queryable comparative genomics information, conserved orthologous genes and pathways are among the new key features updated in this database. The browser has been upgraded to enable user upload of GFF files for quick view of genome annotation comparisons. The toolkit page integrates the EMBOSS package and has a gene prediction tool. Annotations for the organisms are updated once every 6 months to ensure quality. The database resource is available at www.eumicrobedb.org.

A comparative genomic analysis of putative pathogenicity genes in the host-specific sibling species Colletotrichum graminicola and Colletotrichum sublineola

Background

Colletotrichum graminicola and C. sublineola cause anthracnose leaf and stalk diseases of maize and sorghum, respectively. In spite of their close evolutionary relationship, the two species are completely host-specific. Host specificity is often attributed to pathogen virulence factors, including specialized secondary metabolites (SSM), and small-secreted protein (SSP) effectors. Genes relevant to these categories were manually annotated in two co-occurring, contemporaneous strains of C. graminicola and C. sublineola. A comparative genomic and phylogenetic analysis was performed to address the evolutionary relationships among these and other divergent gene families in the two strains.

Results

Inoculation of maize with C. sublineola, or of sorghum with C. graminicola, resulted in rapid plant cell death at, or just after, the point of penetration. The two fungal genomes were very similar. More than 50% of the assemblies could be directly aligned, and more than 80% of the gene models were syntenous. More than 90% of the predicted proteins had orthologs in both species. Genes lacking orthologs in the other species (non-conserved genes) included many predicted to encode SSM-associated proteins and SSPs. Other common groups of non-conserved proteins included transporters, transcription factors, and CAZymes. Only 32 SSP genes appeared to be specific to C. graminicola, and 21 to C. sublineola. None of the SSM-associated genes were lineage-specific. Two different strains of C. graminicola, and three strains of C. sublineola, differed in no more than 1% percent of gene sequences from one another.

Conclusions

Efficient non-host recognition of C. sublineola by maize, and of C. graminicola by sorghum, was observed in epidermal cells as a rapid deployment of visible resistance responses and plant cell death. Numerous non-conserved SSP and SSM-associated predicted proteins that could play a role in this non-host recognition were identified. Additional categories of genes that were also highly divergent suggested an important role for co-evolutionary adaptation to specific host environmental factors, in addition to aspects of initial recognition, in host specificity. This work provides a foundation for future functional studies aimed at clarifying the roles of these proteins, and the possibility of manipulating them to improve management of these two economically important diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3457-9) contains supplementary material, which is available to authorized users.

Subtle genetic population structure in Pacific halibut Hippoglossus stenolepis

Pacific halibut Hippoglossus stenolepis from 10 sampling locations throughout their range were investigated for signs of population structure. Two genetic data sets were created: (1) all individuals (n = 828) at few anonymous microsatellite markers (number of loci = 16); (2) fewer individuals (n = 435) genotyped at anonymous as well as expressed sequence-tag linked microsatellites (number of loci = 61). A combination of multidimensional scaling plots, discriminant analysis of principal components and pairwise differentiation estimates suggested that samples from the Aleutian Islands, particularly the western Aleutian Islands, were genetically distinct from samples collected in other regions. In addition, outlier analyses found that two markers linked to expressed sequence tags may be under directional selection and could explain the differentiation among samples. These results confirm findings from previous research and suggest that population structure may exist within a current management unit (i.e. International Pacific Halibut Commission Regulatory Area 4B).

Comparative Genomic Analysis among Four Representative Isolates of Phytophthora sojae Reveals Genes under Evolutionary Selection

Comparative genomic analysis is useful for identifying genes affected by evolutionary selection and for studying adaptive variation in gene functions. In Phytophthora sojae, a model oomycete plant pathogen, the related study is lacking. We compared sequence data among four isolates of P. sojae, which represent its four major genotypes. These isolates exhibited >99.688%, >99.864%, and >98.981% sequence identities at genome, gene, and non-gene regions, respectively. One hundred and fifty-three positive selection and 139 negative selection candidate genes were identified. Between the two categories of genes, the positive selection genes were flanked by larger intergenic regions, poorly annotated in function, and less conserved; they had relatively lower transcription levels but many genes had increased transcripts during infection. Genes coding for predicted secreted proteins, particularly effectors, were overrepresented in positive selection. Several RxLR effector genes were identified as positive selection genes, exhibiting much stronger positive selection levels. In addition, candidate genes with presence/absence polymorphism were analyzed. This study provides a landscape of genomic variation among four representative P. sojae isolates and characterized several evolutionary selection-affected gene candidates. The results suggest a relatively covert two-speed genome evolution pattern in P. sojae and will provide clues for identification of new virulence factors in the oomycete plant pathogens.

A highly conserved metalloprotease effector enhances virulence in the maize anthracnose fungus Colletotrichum graminicola

Colletotrichum graminicola causes maize anthracnose, an agronomically important disease with a worldwide distribution. We have identified a fungalysin metalloprotease (Cgfl) with a role in virulence. Transcriptional profiling experiments and live cell imaging show that Cgfl is specifically expressed during the biotrophic stage of infection. To determine whether Cgfl has a role in virulence, we obtained null mutants lacking Cgfl and performed pathogenicity and live microscopy assays. The appressorium morphology of the null mutants is normal, but they exhibit delayed development during the infection process on maize leaves and roots, showing that Cgfl has a role in virulence. In vitro chitinase activity assays of leaves infected with wild-type and null mutant strains show that, in the absence of Cgfl, maize leaves exhibit increased chitinase activity. Phylogenetic analyses show that Cgfl is highly conserved in fungi. Similarity searches, phylogenetic analysis and transcriptional profiling show that C. graminicola encodes two LysM domain-containing homologues of Ecp6, suggesting that this fungus employs both Cgfl-mediated and LysM protein-mediated strategies to control chitin signalling.

Sarcococca Blight: Use of Whole-Genome Sequencing for Fungal Plant Disease Diagnosis

Early and accurate diagnosis of new plant pathogens is vital for the rapid implementation of effective mitigation strategies and appropriate regulatory responses. Most commonly, pathogen identification relies on morphology and DNA marker analysis. However, for new diseases, these approaches may not be sufficient for precise diagnosis. In this study, we used whole-genome sequencing (WGS) to identify the causal agent of a new disease affecting Sarcococca hookeriana (sarcococca). Blight symptoms were observed on sarcococca and adjacent Buxus sempervirens (boxwood) plants in Maryland during 2014. Symptoms on sarcococca were novel, and included twig dieback and dark lesions on leaves and stems. A Calonectria sp. was isolated from both hosts and used to fulfill Koch's postulates but morphology and marker sequence data precluded species-level identification. A 51.4-Mb WGS was generated for the two isolates and identified both as Calonectria pseudonaviculata. A single-nucleotide polymorphism at a noncoding site differentiated between the two host isolates. These results indicate that the same C. pseudonaviculata genotype has the ability to induce disease on both plant species. This study marks the first application of WGS for fungal plant pathogen diagnosis and demonstrates the power of this approach to rapidly identify causal agents of new diseases.

Genus-Wide Comparative Genome Analyses of Colletotrichum Species Reveal Specific Gene Family Losses and Gains during Adaptation to Specific Infection Lifestyles

Members from Colletotrichum genus adopt a diverse range of lifestyles during infection of plants and represent a group of agriculturally devastating pathogens. In this study, we present the draft genome of Colletotrichum incanum from the spaethianum clade of Colletotrichum and the comparative analyses with five other Colletotrichum species from distinct lineages. We show that the C. incanum strain, originally isolated from Japanese daikon radish, is able to infect both eudicot plants, such as certain ecotypes of the eudicot Arabidopsis, and monocot plants, such as lily. Being closely related to Colletotrichum species both in the graminicola clade, whose members are restricted strictly to monocot hosts, and to the destructivum clade, whose members are mostly associated with dicot infections, C. incanum provides an interesting model system for comparative genomics to study how fungal pathogens adapt to monocot and dicot hosts. Genus-wide comparative genome analyses reveal that Colletotrichum species have tailored profiles of their carbohydrate-degrading enzymes according to their infection lifestyles. In addition, we show evidence that positive selection acting on secreted and nuclear localized proteins that are highly conserved may be important in adaptation to specific hosts or ecological niches.

Common protein sequence signatures associate with Sclerotinia borealis lifestyle and secretion in fungal pathogens of the Sclerotiniaceae

Fungal plant pathogens produce secreted proteins adapted to function outside fungal cells to facilitate colonization of their hosts. In many cases such as for fungi from the Sclerotiniaceae family the repertoire and function of secreted proteins remains elusive. In the Sclerotiniaceae, whereas Sclerotinia sclerotiorum and Botrytis cinerea are cosmopolitan broad host-range plant pathogens, Sclerotinia borealis has a psychrophilic lifestyle with a low optimal growth temperature, a narrow host range and geographic distribution. To spread successfully, S. borealis must synthesize proteins adapted to function in its specific environment. The search for signatures of adaptation to S. borealis lifestyle may therefore help revealing proteins critical for colonization of the environment by Sclerotiniaceae fungi. Here, we analyzed amino acids usage and intrinsic protein disorder in alignments of groups of orthologous proteins from the three Sclerotiniaceae species. We found that enrichment in Thr, depletion in Glu and Lys, and low disorder frequency in hot loops are significantly associated with S. borealis proteins. We designed an index to report bias in these properties and found that high index proteins were enriched among secreted proteins in the three Sclerotiniaceae fungi. High index proteins were also enriched in function associated with plant colonization in S. borealis, and in in planta-induced genes in S. sclerotiorum. We highlight a novel putative antifreeze protein and a novel putative lytic polysaccharide monooxygenase identified through our pipeline as candidate proteins involved in colonization of the environment. Our findings suggest that similar protein signatures associate with S. borealis lifestyle and with secretion in the Sclerotiniaceae. These signatures may be useful for identifying proteins of interest as targets for the management of plant diseases.

Nucleotide variation and balancing selection at the Ckma gene in Atlantic cod: analysis with multiple merger coalescent models

High-fecundity organisms, such as Atlantic cod, can withstand substantial natural selection and the entailing genetic load of replacing alleles at a number of loci due to their excess reproductive capacity. High-fecundity organisms may reproduce by sweepstakes leading to highly skewed heavy-tailed offspring distribution. Under such reproduction the Kingman coalescent of binary mergers breaks down and models of multiple merger coalescent are more appropriate. Here we study nucleotide variation at the Ckma (Creatine Kinase Muscle type A) gene in Atlantic cod. The gene shows extreme differentiation between the North (Canada, Greenland, Iceland, Norway, Barents Sea) and the South (Faroe Islands, North-, Baltic-, Celtic-, and Irish Seas) with F_ST > 0.8 between regions whereas neutral loci show no differentiation. This is evidence of natural selection. The protein sequence is conserved by purifying selection whereas silent and non-coding sites show extreme differentiation. The unfolded site-frequency spectrum has three modes, a mode at singleton sites and two high frequency modes at opposite frequencies representing divergent branches of the gene genealogy that is evidence for balancing selection. Analysis with multiple-merger coalescent models can account for the high frequency of singleton sites and indicate reproductive sweepstakes. Coalescent time scales vary with population size and with the inverse of variance in offspring number. Parameter estimates using multiple-merger coalescent models show that times scales are faster than under the Kingman coalescent.

Common protein sequence signatures associate with Sclerotinia borealis lifestyle and secretion in fungal pathogens of the Sclerotiniaceae

Fungal plant pathogens produce secreted proteins adapted to function outside fungal cells to facilitate colonization of their hosts. In many cases such as for fungi from the Sclerotiniaceae family the repertoire and function of secreted proteins remains elusive. In the Sclerotiniaceae, whereas Sclerotinia sclerotiorum and Botrytis cinerea are cosmopolitan broad host-range plant pathogens, Sclerotinia borealis has a psychrophilic lifestyle with a low optimal growth temperature, a narrow host range and geographic distribution. To spread successfully, S. borealis must synthesize proteins adapted to function in its specific environment. The search for signatures of adaptation to S. borealis lifestyle may therefore help revealing proteins critical for colonization of the environment by Sclerotiniaceae fungi. Here, we analyzed amino acids usage and intrinsic protein disorder in alignments of groups of orthologous proteins from the three Sclerotiniaceae species. We found that enrichment in Thr, depletion in Glu and Lys, and low disorder frequency in hot loops are significantly associated with S. borealis proteins. We designed an index to report bias in these properties and found that high index proteins were enriched among secreted proteins in the three Sclerotiniaceae fungi. High index proteins were also enriched in function associated with plant colonization in S. borealis, and in in planta-induced genes in S. sclerotiorum. We highlight a novel putative antifreeze protein and a novel putative lytic polysaccharide monooxygenase identified through our pipeline as candidate proteins involved in colonization of the environment. Our findings suggest that similar protein signatures associate with S. borealis lifestyle and with secretion in the Sclerotiniaceae. These signatures may be useful for identifying proteins of interest as targets for the management of plant diseases.