Population genomics of Fusarium graminearum isolates from the Americas.

Fusarium head blight (FHB) is a major disease of wheat and barley worldwide and is caused by different species in the genus Fusarium, Fusarium graminearum being the most important. We conducted population genomics analyses using SNPs obtained through genotyping by sequencing of over 500 isolates of F. graminearum from the US Upper Midwest, New York, Louisiana, and Uruguay. PCA and STRUCTURE analyses group our isolates into four previously described populations: NA1, NA2, Southern Louisiana (SLA) and Gulf Coast (GC). Some isolates were not assigned to populations because of mixed ancestry. Population structure was associated with toxin genotype and geographic origin. The NA1, NA2, and SLA populations are differentiated (FST 0.385 - 0.551) but the presence of admixed isolates indicates that the populations are not reproductively isolated. Patterns of linkage disequilibrium (LD) decay suggest frequent recombination within populations. Fusarium graminearum populations from the US have great evolutionary potential given the high recombination rate and a large proportion of admixed isolates. The NA1, NA2, and Southern Louisiana (SLA) populations separated from their common ancestral population roughly at the same time in the past and are evolving with moderate levels of subsequent gene flow between them. Genome-wide selection scans in all three populations revealed outlier regions with the strongest signatures of recent positive natural selection. These outlier regions include many genes with unknown function and some genes with known roles in plant-microbe interaction, fungicide/drug resistance, cellular transport and genes that are related to cellular organelles. Only a very small proportion of outlier regions are shared as outliers among the three populations, suggesting unique host-pathogen interactions and environmental adaptation.

Divergence and Gene Flow Between Fusarium subglutinans and F. temperatum Isolated from Maize in Argentina.

Fusarium subglutinans and F. temperatum are two important fungal pathogens of maize whose distinctness as separate species has been difficult to assess. We isolated strains of these species from commercial and native maize varieties in Argentina and sequenced >28,000 loci to estimate genetic variation in the sample. Our objectives were to measure genetic divergence between the species, infer demographic parameters related to their split, and describe the population structure of the sample. When analyzed together, over 30% of each species' polymorphic sites (>2,500 sites) segregate as polymorphisms in the other. Demographic modeling confirmed the species split predated maize domestication, but subsequent between-species gene flow has occurred, with gene flow from F. subglutinans into F. temperatum greater than gene flow in the reverse direction. In F. subglutinans, little evidence exists for substructure or recent selective sweeps, but there is evidence for limited sexual reproduction. In F. temperatum, there is clear evidence for population substructure and signals of abundant recent selective sweeps, with sexual reproduction probably less common than in F. subglutinans. Both genetic variation and the relative number of polymorphisms shared between species increase near the telomeres of all 12 chromosomes, where genes related to plant-pathogen interactions often are located. Our results suggest that species boundaries between closely related Fusarium species can be semipermeable and merit further study. Such semipermeability could facilitate unanticipated genetic exchange between species and enable quicker permanent responses to changes in the agro-ecosystem, e.g., pathogen-resistant host varieties, new chemical and biological control agents, and agronomic practices.

Fumonisin and Beauvericin Chemotypes and Genotypes of the Sister Species Fusarium subglutinans and Fusarium temperatum.

Fusarium subglutinans and Fusarium temperatum are common maize pathogens that produce mycotoxins and cause plant disease. The ability of these species to produce beauvericin and fumonisin mycotoxins is not settled, as reports of toxin production are not concordant. Our objective was to clarify this situation by determining both the chemotypes and genotypes for strains from both species. We analyzed 25 strains from Argentina, 13 F. subglutinans and 12 F. temperatum strains, for toxin production by ultraperformance liquid chromatography mass spectrometry (UPLC-MS). We used new genome sequences from two strains of F. subglutinans and one strain of F. temperatum, plus genomes of other Fusarium species, to determine the presence of functional gene clusters for the synthesis of these toxins. None of the strains examined from either species produced fumonisins. These strains also lack Fum biosynthetic genes but retain homologs of some genes that flank the Fum cluster in Fusarium verticillioides None of the F. subglutinans strains we examined produced beauvericin although 9 of 12 F. temperatum strains did. A complete beauvericin (Bea) gene cluster was present in all three new genome sequences. The Bea1 gene was presumably functional in F. temperatum but was not functional in F. subglutinans due to a large insertion and multiple mutations that resulted in premature stop codons. The accumulation of only a few mutations expected to disrupt Bea1 suggests that the process of its inactivation is relatively recent. Thus, none of the strains of F. subglutinans or F. temperatum we examined produce fumonisins, and the strains of F. subglutinans examined also cannot produce beauvericin. Variation in the ability of strains of F. temperatum to produce beauvericin requires further study and could reflect the recent shared ancestry of these two species.IMPORTANCEFusarium subglutinans and F. temperatum are sister species and maize pathogens commonly isolated worldwide that can produce several mycotoxins and cause seedling disease, stalk rot, and ear rot. The ability of these species to produce beauvericin and fumonisin mycotoxins is not settled, as reports of toxin production are not concordant at the species level. Our results are consistent with previous reports that strains of F. subglutinans produce neither fumonisins nor beauvericin. The status of toxin production by F. temperatum needs further work. Our strains of F. temperatum did not produce fumonisins, while some strains produced beauvericin and others did not. These results enable more accurate risk assessments of potential mycotoxin contamination if strains of these species are present. The nature of the genetic inactivation of BEA1 is consistent with its relatively recent occurrence and the close phylogenetic relationship of the two sister species.

Molecular, genetic and evolutionary analysis of a paracentric inversion in Arabidopsis thaliana.

Chromosomal inversions can provide windows onto the cytogenetic, molecular, evolutionary and demographic histories of a species. Here we investigate a paracentric 1.17-Mb inversion on chromosome 4 of Arabidopsis thaliana with nucleotide precision of its borders. The inversion is created by Vandal transposon activity, splitting an F-box and relocating a pericentric heterochromatin segment in juxtaposition with euchromatin without affecting the epigenetic landscape. Examination of the RegMap panel and the 1001 Arabidopsis genomes revealed more than 170 inversion accessions in Europe and North America. The SNP patterns revealed historical recombinations from which we infer diverse haplotype patterns, ancient introgression events and phylogenetic relationships. We find a robust association between the inversion and fecundity under drought. We also find linkage disequilibrium between the inverted region and the early flowering Col-FRIGIDA allele. Finally, SNP analysis elucidates the origin of the inversion to South-Eastern Europe approximately 5000 years ago and the FRI-Col allele to North-West Europe, and reveals the spreading of a single haplotype to North America during the 17th to 19th century. The 'American haplotype' was identified from several European localities, potentially due to return migration.

Multiple reference genomes and transcriptomes for Arabidopsis thaliana.

Genetic differences between Arabidopsis thaliana accessions underlie the plant's extensive phenotypic variation, and until now these have been interpreted largely in the context of the annotated reference accession Col-0. Here we report the sequencing, assembly and annotation of the genomes of 18 natural A. thaliana accessions, and their transcriptomes. When assessed on the basis of the reference annotation, one-third of protein-coding genes are predicted to be disrupted in at least one accession. However, re-annotation of each genome revealed that alternative gene models often restore coding potential. Gene expression in seedlings differed for nearly half of expressed genes and was frequently associated with cis variants within 5 kilobases, as were intron retention alternative splicing events. Sequence and expression variation is most pronounced in genes that respond to the biotic environment. Our data further promote evolutionary and functional studies in A. thaliana, especially the MAGIC genetic reference population descended from these accessions.

Recombination and linkage disequilibrium in Arabidopsis thaliana.

Linkage disequilibrium (LD) is a major aspect of the organization of genetic variation in natural populations. Here we describe the genome-wide pattern of LD in a sample of 19 Arabidopsis thaliana accessions using 341,602 non-singleton SNPs. LD decays within 10 kb on average, considerably faster than previously estimated. Tag SNP selection algorithms and 'hide-the-SNP' simulations suggest that genome-wide association mapping will require only 40%-50% of the observed SNPs, a reduction similar to estimates in a sample of African Americans. An Affymetrix genotyping array containing 250,000 SNPs has been designed based on these results; we demonstrate that it should have more than adequate coverage for genome-wide association mapping. The extent of LD is highly variable, and we find clear evidence of recombination hotspots, which seem to occur preferentially in intergenic regions. LD also reflects the action of selection, and it is more extensive between nonsynonymous polymorphisms than between synonymous polymorphisms.

The landscape and predicted roles of structural variants in Fusarium graminearum genomes.

Structural rearrangements, such as inversions, translocations, duplications, and large insertions and deletions, are large-scale genomic variants that can play an important role in shaping phenotypic variation and in genome adaptation and evolution. We used chromosomal-level assemblies from eight Fusarium graminearum isolates to study structural variants and their role in fungal evolution. We generated the assemblies of four of these genomes after Oxford Nanopore sequencing. A total of 87 inversions, 159 translocations, 245 duplications, 58,489 insertions, and 34,102 deletions were detected. Regions of high recombination rate are associated with structural rearrangements, and a significant proportion of inversions, translocations, and duplications overlap with the repeat content of the genome, suggesting recombination and repeat elements are major factors in the origin of structural rearrangements in F. graminearum. Large insertions and deletions introduce presence-absence polymorphisms for many genes, including secondary metabolite biosynthesis cluster genes and predicted effectors genes. Translocation events were found to be shuffling predicted effector-rich regions of the genomes and are likely contributing to the gain and loss of effectors facilitated by recombination. Breakpoints of some structural rearrangements fall within coding sequences and are likely altering the protein products. Structural rearrangements in F. graminearum thus have an important role to play in shaping pathogen-host interactions and broader evolution through genome reorganization, the introduction of presence-absence polymorphisms, and changing protein products and gene regulation.

An Arabidopsis example of association mapping in structured samples.

A potentially serious disadvantage of association mapping is the fact that marker-trait associations may arise from confounding population structure as well as from linkage to causative polymorphisms. Using genome-wide marker data, we have previously demonstrated that the problem can be severe in a global sample of 95 Arabidopsis thaliana accessions, and that established methods for controlling for population structure are generally insufficient. Here, we use the same sample together with a number of flowering-related phenotypes and data-perturbation simulations to evaluate a wider range of methods for controlling for population structure. We find that, in terms of reducing the false-positive rate while maintaining statistical power, a recently introduced mixed-model approach that takes genome-wide differences in relatedness into account via estimated pairwise kinship coefficients generally performs best. By combining the association results with results from linkage mapping in F2 crosses, we identify one previously known true positive and several promising new associations, but also demonstrate the existence of both false positives and false negatives. Our results illustrate the potential of genome-wide association scans as a tool for dissecting the genetics of natural variation, while at the same time highlighting the pitfalls. The importance of study design is clear; our study is severely under-powered both in terms of sample size and marker density. Our results also provide a striking demonstration of confounding by population structure. While statistical methods can be used to ameliorate this problem, they cannot always be effective and are certainly not a substitute for independent evidence, such as that obtained via crosses or transgenic experiments. Ultimately, association mapping is a powerful tool for identifying a list of candidates that is short enough to permit further genetic study.

Low levels of polymorphism in genes that control the activation of defense response in Arabidopsis thaliana.

Plants use signaling pathways involving salicylic acid, jasmonic acid, and ethylene to defend against pathogen and herbivore attack. Many defense response genes involved in these signaling pathways have been characterized, but little is known about the selective pressures they experience. A representative set of 27 defense response genes were resequenced in a worldwide set of 96 Arabidopsis thaliana accessions, and patterns of single nucleotide polymorphisms (SNPs) were evaluated in relation to an empirical distribution of SNPs generated from either 876 fragments or 236 fragments with >400 bp coding sequence (this latter set was selected for comparisons with coding sequences) distributed across the genomes of the same set of accessions. Defense response genes have significantly fewer protein variants, display lower levels of nonsynonymous nucleotide diversity, and have fewer nonsynonymous segregating sites. The majority of defense response genes appear to be experiencing purifying selection, given the dearth of protein variation in this set of genes. Eight genes exhibit some evidence of partial selective sweeps or transient balancing selection. These results therefore provide a strong contrast to the high levels of balancing selection exhibited by genes at the upstream positions in these signaling pathways.

A nonparametric test reveals selection for rapid flowering in the Arabidopsis genome.

The detection of footprints of natural selection in genetic polymorphism data is fundamental to understanding the genetic basis of adaptation, and has important implications for human health. The standard approach has been to reject neutrality in favor of selection if the pattern of variation at a candidate locus was significantly different from the predictions of the standard neutral model. The problem is that the standard neutral model assumes more than just neutrality, and it is almost always possible to explain the data using an alternative neutral model with more complex demography. Today's wealth of genomic polymorphism data, however, makes it possible to dispense with models altogether by simply comparing the pattern observed at a candidate locus to the genomic pattern, and rejecting neutrality if the pattern is extreme. Here, we utilize this approach on a truly genomic scale, comparing a candidate locus to thousands of alleles throughout the Arabidopsis thaliana genome. We demonstrate that selection has acted to increase the frequency of early-flowering alleles at the vernalization requirement locus FRIGIDA. Selection seems to have occurred during the last several thousand years, possibly in response to the spread of agriculture. We introduce a novel test statistic based on haplotype sharing that embraces the problem of population structure, and so should be widely applicable.

The pattern of polymorphism in Arabidopsis thaliana.

We resequenced 876 short fragments in a sample of 96 individuals of Arabidopsis thaliana that included stock center accessions as well as a hierarchical sample from natural populations. Although A. thaliana is a selfing weed, the pattern of polymorphism in general agrees with what is expected for a widely distributed, sexually reproducing species. Linkage disequilibrium decays rapidly, within 50 kb. Variation is shared worldwide, although population structure and isolation by distance are evident. The data fail to fit standard neutral models in several ways. There is a genome-wide excess of rare alleles, at least partially due to selection. There is too much variation between genomic regions in the level of polymorphism. The local level of polymorphism is negatively correlated with gene density and positively correlated with segmental duplications. Because the data do not fit theoretical null distributions, attempts to infer natural selection from polymorphism data will require genome-wide surveys of polymorphism in order to identify anomalous regions. Despite this, our data support the utility of A. thaliana as a model for evolutionary functional genomics.

Genome-wide association mapping in Arabidopsis identifies previously known flowering time and pathogen resistance genes.

There is currently tremendous interest in the possibility of using genome-wide association mapping to identify genes responsible for natural variation, particularly for human disease susceptibility. The model plant Arabidopsis thaliana is in many ways an ideal candidate for such studies, because it is a highly selfing hermaphrodite. As a result, the species largely exists as a collection of naturally occurring inbred lines, or accessions, which can be genotyped once and phenotyped repeatedly. Furthermore, linkage disequilibrium in such a species will be much more extensive than in a comparable outcrossing species. We tested the feasibility of genome-wide association mapping in A. thaliana by searching for associations with flowering time and pathogen resistance in a sample of 95 accessions for which genome-wide polymorphism data were available. In spite of an extremely high rate of false positives due to population structure, we were able to identify known major genes for all phenotypes tested, thus demonstrating the potential of genome-wide association mapping in A. thaliana and other species with similar patterns of variation. The rate of false positives differed strongly between traits, with more clinal traits showing the highest rate. However, the false positive rates were always substantial regardless of the trait, highlighting the necessity of an appropriate genomic control in association studies.

Plant organ evolution revealed by phylotranscriptomics in Arabidopsis thaliana.

The evolution of phenotypes occurs through changes both in protein sequence and gene expression levels. Though much of plant morphological evolution can be explained by changes in gene expression, examining its evolution has challenges. To gain a new perspective on organ evolution in plants, we applied a phylotranscriptomics approach. We combined a phylostratigraphic approach with gene expression based on the strand-specific RNA-seq data from seedling, floral bud, and root of 19 Arabidopsis thaliana accessions to examine the age and sequence divergence of transcriptomes from these organs and how they adapted over time. Our results indicate that, among the sense and antisense transcriptomes of these organs, the sense transcriptomes of seedlings are the evolutionarily oldest across all accessions and are the most conserved in amino acid sequence for most accessions. In contrast, among the sense transcriptomes from these same organs, those from floral bud are evolutionarily youngest and least conserved in sequence for most accessions. Different organs have adaptive peaks at different stages in their evolutionary history; however, all three show a common adaptive signal from the Magnoliophyta to Brassicale stage. Our research highlights how phylotranscriptomic analyses can be used to trace organ evolution in the deep history of plant species.

Unbiased K-mer Analysis Reveals Changes in Copy Number of Highly Repetitive Sequences During Maize Domestication and Improvement.

The major component of complex genomes is repetitive elements, which remain recalcitrant to characterization. Using maize as a model system, we analyzed whole genome shotgun (WGS) sequences for the two maize inbred lines B73 and Mo17 using k-mer analysis to quantify the differences between the two genomes. Significant differences were identified in highly repetitive sequences, including centromere, 45S ribosomal DNA (rDNA), knob, and telomere repeats. Genotype specific 45S rDNA sequences were discovered. The B73 and Mo17 polymorphic k-mers were used to examine allele-specific expression of 45S rDNA in the hybrids. Although Mo17 contains higher copy number than B73, equivalent levels of overall 45S rDNA expression indicates that transcriptional or post-transcriptional regulation mechanisms operate for the 45S rDNA in the hybrids. Using WGS sequences of B73xMo17 doubled haploids, genomic locations showing differential repetitive contents were genetically mapped, which displayed different organization of highly repetitive sequences in the two genomes. In an analysis of WGS sequences of HapMap2 lines, including maize wild progenitor, landraces, and improved lines, decreases and increases in abundance of additional sets of k-mers associated with centromere, 45S rDNA, knob, and retrotransposons were found among groups, revealing global evolutionary trends of genomic repeats during maize domestication and improvement.

Sequence variation and haplotype structure at the human HFE locus.

The HFE locus encodes an HLA class-I-type protein important in iron regulation and segregates replacement mutations that give rise to the most common form of genetic hemochromatosis. The high frequency of one disease-associated mutation, C282Y, and the nature of this disease have led some to suggest a selective advantage for this mutation. To investigate the context in which this mutation arose and gain a better understanding of HFE genetic variation, we surveyed nucleotide variability in 11.2 kb encompassing the HFE locus and experimentally determined haplotypes. We fully resequenced 60 chromosomes of African, Asian, or European ancestry as well as one chimpanzee, revealing 41 variable sites and a nucleotide diversity of 0.08%. This indicates that linkage to the HLA region has not substantially increased the level of HFE variation. Although several haplotypes are shared between populations, one haplotype predominates in Asia but is nearly absent elsewhere, causing higher than average genetic differentiation among the three major populations. Our samples show evidence of intragenic recombination, so the scarcity of recombination events within the C282Y allele class is consistent with selection increasing the frequency of a young allele. Otherwise, the pattern of variability in this region does not clearly indicate the action of positive selection at this or linked loci.

A method for detecting recent selection in the human genome from allele age estimates.

Mutations that have recently increased in frequency by positive natural selection are an important component of naturally occurring variation that affects fitness. To identify such variants, we developed a method to test for recent selection by estimating the age of an allele from the extent of haplotype sharing at linked sites. Neutral coalescent simulations are then used to determine the likelihood of this age given the allele's observed frequency. We applied this method to a common disease allele, the hemochromatosis-associated HFE C282Y mutation. Our results allow us to reject neutral models incorporating plausible human demographic histories for HFE C282Y and one other young but common allele, indicating positive selection at HFE or a linked locus. This method will be useful for scanning the human genome for alleles under selection using the haplotype map now being constructed.

Genome Sequencing of Multiple Isolates Highlights Subtelomeric Genomic Diversity within Fusarium fujikuroi.

Comparisons of draft genome sequences of three geographically distinct isolates of Fusarium fujikuroi with two recently published genome sequences from the same species suggest diverse profiles of secondary metabolite production within F. fujikuroi. Species- and lineage-specific genes, many of which appear to exhibit expression profiles that are consistent with roles in host-pathogen interactions and adaptation to environmental changes, are concentrated in subtelomeric regions. These genomic compartments also exhibit distinct gene densities and compositional characteristics with respect to other genomic partitions, and likely play a role in the generation of molecular diversity. Our data provide additional evidence that gene duplication, divergence, and differential loss play important roles in F. fujikuroi genome evolution and suggest that hundreds of lineage-specific genes might have been acquired through horizontal gene transfer.

A haplotype map of allohexaploid wheat reveals distinct patterns of selection on homoeologous genomes.

BACKGROUND: Bread wheat is an allopolyploid species with a large, highly repetitive genome. To investigate the impact of selection on variants distributed among homoeologous wheat genomes and to build a foundation for understanding genotype-phenotype relationships, we performed population-scale re-sequencing of a diverse panel of wheat lines. RESULTS: A sample of 62 diverse lines was re-sequenced using the whole exome capture and genotyping-by-sequencing approaches. We describe the allele frequency, functional significance, and chromosomal distribution of 1.57 million single nucleotide polymorphisms and 161,719 small indels. Our results suggest that duplicated homoeologous genes are under purifying selection. We find contrasting patterns of variation and inter-variant associations among wheat genomes; this, in addition to demographic factors, could be explained by differences in the effect of directional selection on duplicated homoeologs. Only a small fraction of the homoeologous regions harboring selected variants overlapped among the wheat genomes in any given wheat line. These selected regions are enriched for loci associated with agronomic traits detected in genome-wide association studies. CONCLUSIONS: Evidence suggests that directional selection in allopolyploids rarely acted on multiple parallel advantageous mutations across homoeologous regions, likely indicating that a fitness benefit could be obtained by a mutation at any one of the homoeologs. Additional advantageous variants in other homoelogs probably either contributed little benefit, or were unavailable in populations subjected to directional selection. We hypothesize that allopolyploidy may have increased the likelihood of beneficial allele recovery by broadening the set of possible selection targets.

Parental genetic distance and patterns in nonrandom mating and seed yield in predominately selfing Arabidopsis thaliana.

In this study, we ask two questions: (1) Is reproductive success independent of parental genetic distance in predominately selfing plants? (2) In the absence of early inbreeding depression, is there substantial maternal and/or paternal variation in reproductive success in natural populations? Seed yield in single pollinations and proportion of seeds sired in mixed pollinations were studied in genetically defined accessions of the predominately selfing plant Arabidopsis thaliana by conducting two diallel crosses. The first diallel was a standard, single pollination design that we used to examine variance in seed yield. The second diallel was a mixed pollination design that utilized a standard pollen competitor to examine variance in proportion of seeds sired. We found no correlation between reproductive success and parental genetic distance, and self-pollen does not systematically differ in reproductive success compared to outcross pollen, suggesting that Arabidopsis populations do not experience embryo lethality due to early-acting inbreeding or outbreeding depression. We used these data to partition the contributions to total phenotypic variation from six sources, including maternal contributions, paternal contributions and parental interactions. For seed yield in single pollinations, maternal effects accounted for the most significant source of variance (16.6 %). For proportion of seeds sired in mixed pollinations, the most significant source of variance was paternal effects (17.9 %). Thus, we show that population-level genetic similarities, including selfing, do not correlate with reproductive success, yet there is still significant paternal variance under competition. This suggests two things. First, since these differences are unlikely due to early-acting inbreeding depression or differential pollen viability, this implicates natural variation in pollen germination and tube growth dynamics. Second, this strongly supports a model of fixation of pollen performance genes in populations, offering a focus for future genetic studies in differential reproductive success.

Genome-wide patterns of genetic variation in worldwide Arabidopsis thaliana accessions from the RegMap panel.

Arabidopsis thaliana is native to Eurasia and is naturalized across the world. Its ability to be easily propagated and its high phenotypic variability make it an ideal model system for functional, ecological and evolutionary genetics. To date, analyses of the natural genetic variation of A. thaliana have involved small numbers of individual plants or genetic markers. Here we genotype 1,307 worldwide accessions, including several regional samples, using a 250K SNP chip. This allowed us to produce a high-resolution description of the global pattern of genetic variation. We applied three complementary selection tests and identified new targets of selection. Further, we characterized the pattern of historical recombination in A. thaliana and observed an enrichment of hotspots in its intergenic regions and repetitive DNA, which is consistent with the pattern that is observed for humans but which is strikingly different from that observed in other plant species. We have made the seeds we used to produce this Regional Mapping (RegMap) panel publicly available. This panel comprises one of the largest genomic mapping resources currently available for global natural isolates of a non-human species.