Reconstructing the origins and the biogeography of species’ genomes in the highly reticulate allopolyploid-rich model grass genus Brachypodium using minimum evolution, coalescence and maximum likelihood approaches

• Published in: Molecular phylogenetics and evolution 2018-10, Vol.127, p.256-271
• Main Authors: Díaz-Pérez, Antonio, López-Álvarez, Diana, Sancho, Rubén, Catalán, Pilar
• Format: Article
• Language: English
• Subjects: Alleles; Bayes Theorem; Biological Evolution; Brachypodium - genetics; Diploidy; Divergence times (lineages) and coalescence ages (allopolyploids); Genome, Plant; Genomic biogeography; Grafted species tree and species network; Homeologous genomes; Likelihood Functions; Model Brachypodium grasses; Models, Biological; Phylogeny; Phylogeography; Polyploidy; Reticulate historical scenarios; Species Specificity; Time Factors
• ISSN: 1055-7903; 1095-9513
• DOI: 10.1016/j.ympev.2018.06.003

[Display omitted] •A comprehensive 5-gene phylogeny of the model grass genus Brachypodium was built.•Minimum-Evolution and Species-Network approaches identified the concurring genomes.•Detected homeologous genomes matched the expected ploidy levels of allopolyploids.•Splittings and mergings of genomes occurred in different spatio-temporal scenarios.•Our biogeographical study infers dispersals of diploids but not of allopolyploids. The identification of homeologous genomes and the biogeographical analyses of highly reticulate allopolyploid-rich groups face the challenge of incorrectly inferring the genomic origins and the biogeographical patterns of the polyploids from unreliable strictly bifurcating trees. Here we reconstruct a plausible evolutionary scenario of the diverging and merging genomes inherited by the diploid and allopolyploid species and cytotypes of the model grass genus Brachypodium. We have identified the ancestral Brachypodium genomes and inferred the paleogeographical ranges for potential hybridization events that originated its allopolyploid taxa. We also constructed a comprehensive phylogeny of Brachypodium from five nuclear and plastid genes using Species Tree Minimum Evolution allele grafting and Species Network analysis. The divergence ages of the lineages were estimated from a consensus maximum clade credibility tree using fossil calibrations, whereas ages of origin of the diploid and allopolyploid species were inferred from coalescence Bayesian methods. The biogeographical events of the genomes were reconstructed using a stratified Dispersal-Extinction-Colonization model with three temporal windows. Our combined Minimum Evolution-coalescence-Bayesian approach allowed us to infer the origins and the identities of the homeologous genomes of the Brachypodium allopolyploids, matching the expected ploidy levels of the hybrids. To date, the current extant progenitor genomes (species) are only known for B. hybridum. Putative ancestral homeologous genome have been inherited by B. mexicanum, ancestral and recent genomes by B. boissieri, and only recently evolved genomes by B. retusum and the core perennial clade allopolyploids (B. phoenicoides, B. pinnatum 4x, B. rupestre 4x). We dissected the complex spatio-temporal evolution of ancestral and recent genomes and have detected successive splitting, dispersal and merging events for dysploid homeologous genomes in diverse geographical scenarios that have led to the current extant taxa. Our data