On computing relative effective population size estimates and parameters from an equilibrium cycle of hermaphrodite frequency fluctuation due to mixed reproductive modes in filamentous fungi

Many filamentous ascomycete fungi reproduce primarily asexually, with only occasional sexual generations. This can lead to a departure from the 1:1 mating type ratio that is expected in obligate sexual populations. The relaxed selection on sexual traits also can lead to a decrease in the frequency of female fertile strains in field populations, while male fertility does not similarly decrease since male gametes also can serve as asexual spores. Both changes ultimately impact the strength of genetic drift in populations. The frequency of female sterility likely increases with the time since the last generation of sexual reproduction, such that it can be used to estimate the relative frequency of sexual reproduction. Here I provide additional details relevant to Leslie and Klein’s (1996) model of mixed sexual reproduction and vegetative propagation as related to the frequency of female sterility. This includes new or modified equations that allow for simpler calculations of i) two estimates of relative Ne, ii) the expected range of hermaphrodite frequencies during the cycles of mixed reproduction, and iii) the relative frequency of sex. These equations also are included in spreadsheet templates into which researchers can directly enter frequencies computed from their population data to estimate these parameters for their own populations. These resources will make the results of Leslie and Klein (1996) more accessible and should increase the use of this model in evaluating the frequency of sexual reproduction of filamentous fungi

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.Fil: Fumero, María Verónica. Universidad Nacional de Río Cuarto. Facultad de Ciencias Exactas Fisicoquímicas y Naturales. Instituto de Investigación en Micología y Micotoxicología. - Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigación en Micología y Micotoxicología; ArgentinaFil: Villani, Alessandra. Institute of Sciences of Food Production; ItaliaFil: Susca, Antonia. Institute of Sciences of Food Production; ItaliaFil: Haidukowski, Miriam. Institute of Sciences of Food Production; ItaliaFil: Cimmarusti, María. Institute of Sciences of Food Production; ItaliaFil: Toomajian, Christopher. Kansas State University; Estados UnidosFil: Leslie, John Franklin. Kansas State University; Estados UnidosFil: Chulze, Sofia Noemi. Universidad Nacional de Río Cuarto. Facultad de Ciencias Exactas Fisicoquímicas y Naturales. Instituto de Investigación en Micología y Micotoxicología. - Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigación en Micología y Micotoxicología; ArgentinaFil: Moretti, Antonio. Institute of Sciences of Food Production; Itali

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

Citation: Liu, S. Z., Zheng, J., Migeon, P., Ren, J., Hu, Y., He, C., . . . Wang, G. Y. (2017). Unbiased K-mer Analysis Reveals Changes in Copy Number of Highly Repetitive Sequences During Maize Domestication and Improvement. Scientific Reports, 7, 15. https://doi.org/10.1038/srep42444The 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 allelespecific 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

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

Genome-Wide Association Mapping in <i>Arabidopsis</i> 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.</p

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

Citation: Jordan, K. W., Wang, S., Lun, Y., Gardiner, L. J., MacLachlan, R., Hucl, P., . . . Akhunov, E. (2015). A haplotype map of allohexaploid wheat reveals distinct patterns of selection on homoeologous genomes. Genome Biology, 16(1). https://doi.org/10.1186/s13059-015-0606-4Background: 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. © 2015 Jordan et al.; licensee BioMed Central.Additional Authors: Talbert, L.;Bansal, U. K.;Bariana, H. S.;Hayden, M. J.;Pozniak, C.;Jeddeloh, J. A.;Hall, A

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

The Pattern of Polymorphism in <i>Arabidopsis thaliana</i>

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

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

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