Search

Speciation proceeds through mechanisms that promote reproductive isolation and shape the extent of genetic variation in natural populations, and thus its study is essential to understand the evolutionary processes leading to increased biodiversity. Chromosomal rearrangements are known to facilitate reproductive isolation by hybrid sterility and favour speciation events. The genus Ipheion (Amaryllidaceae, Allioideae) is unique as its species exhibit a remarkable karyological variability but lack population-level genetic data. To unveil the diversification processes acting upon the formation of new lineages within Ipheion in the Pampas of South America, we combined morphology and karyology approaches with genotyping-by-sequencing. Our phylogenomic and population genomics results supported the taxonomic division of Ipheion into three morphological and genetically well-differentiated groups. The origin of Ipheion uniflorum was traced back to its current southern distribution area in the southern Pampean region (in Argentina), from where it had expanded to the north reaching Uruguay. Our results further suggested that chromosome rearrangements and ploidy shifts had triggered speciation events, first during the origin of I. uniflorum and later during its subsequent diversification into I. recurvifolium and I. tweedieanum, in both cases reinforced by extrinsic factors and biogeographical settings. The current study illustrates the analytical power of multidisciplinary approaches integrating phylo- and population genomics with classic analyses to reveal evolutionary processes in plants. Fil: Sassone, Agostina Belén. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Botánica Darwinion. Academia Nacional de Ciencias Exactas, Físicas y Naturales. Instituto de Botánica Darwinion; Argentina. Leibniz Institute of Plant Genetics and Crop Plant Research; Alemania Fil: Hojsgaard, Diego H.. Leibniz Institute of Plant Genetics and Crop Plant Research; Alemania Fil: Giussani, Liliana Mónica. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Botánica Darwinion. Academia Nacional de Ciencias Exactas, Físicas y Naturales. Instituto de Botánica Darwinion; Argentina Fil: Brassac, Jonathan. Leibniz Institute of Plant Genetics and Crop Plant Research; Alemania Fil: Blattner, Frank R.. Leibniz Institute of Plant Genetics and Crop Plant Research; Alemania

Significance Grain number is a key determinant of cereal grain yield, but its underlying genetic basis in wheat remains undefined. This study demonstrates a direct association between increased floret fertility, higher grain number per spike, and higher plot yields of field-grown wheat. The GNI1 gene, encoding an HD-Zip I transcription factor, was identified as responsible for increased floret fertility. The WT allele acts specifically during rachilla development, with its product serving to lower grain yield potential; in contrast, the reduced-function variant increased both floret and grain number. GNI1 evolved through gene duplication in the Triticeae, and its mutations have been under parallel selection in both wheat and barley over the course of domestication., Floret fertility is a key determinant of the number of grains per inflorescence in cereals. During the evolution of wheat (Triticum sp.), floret fertility has increased, such that current bread wheat (Triticum aestivum) cultivars set three to five grains per spikelet. However, little is known regarding the genetic basis of floret fertility. The locus Grain Number Increase 1 (GNI1) is shown here to be an important contributor to floret fertility. GNI1 evolved in the Triticeae through gene duplication. The gene, which encodes a homeodomain leucine zipper class I (HD-Zip I) transcription factor, was expressed most abundantly in the most apical floret primordia and in parts of the rachilla, suggesting that it acts to inhibit rachilla growth and development. The level of GNI1 expression has decreased over the course of wheat evolution under domestication, leading to the production of spikes bearing more fertile florets and setting more grains per spikelet. Genetic analysis has revealed that the reduced-function allele GNI-A1 contributes to the increased number of fertile florets per spikelet. The RNAi-based knockdown of GNI1 led to an increase in the number of both fertile florets and grains in hexaploid wheat. Mutants carrying an impaired GNI-A1 allele out-yielded WT allele carriers under field conditions. The data show that gene duplication generated evolutionary novelty affecting floret fertility while mutations favoring increased grain production have been under selection during wheat evolution under domestication.

Total spikelet number per spike (TSN) is a major component of spike architecture in wheat (Triticumaestivum L.). A major and consistent quantitative trait locus (QTL) was discovered for TSN in a doubled haploid spring wheat population grown in the field over 4 years. The QTL on chromosome 7B explained up to 20.5% of phenotypic variance. In its physical interval (7B: 6.37–21.67 Mb), the gene FLOWERINGLOCUST (FT-B1) emerged as candidate for the observed effect. In one of the parental lines, FT-B1 carried a non-synonymous substitution on position 19 of the coding sequence. This mutation modifying an aspartic acid (D) into a histidine (H) occurred in a highly conserved position. The mutation was observed with a frequency of ca. 68% in a set of 135 hexaploid wheat varieties and landraces, while it was not found in other plant species. FT-B1 only showed a minor effect on heading and flowering time (FT) which were dominated by a major QTL on chromosome 5A caused by segregation of the vernalization gene VRN-A1. Individuals carrying the FT-B1 allele with amino acid histidine had, on average, a higher number of spikelets (15.1) than individuals with the aspartic acid allele (14.3) independent of their VRN-A1 allele. We show that the effect of TSN is not mainly related to flowering time; however, the duration of pre-anthesis phases may play a major role.

Grain quality traits determine the classification of registered wheat (Triticum aestivum L.) varieties. Although environmental factors and crop management practices exert a considerable influence on wheat quality traits, a significant proportion of the variance is attributed to the genetic factors. To identify the underlying genetic factors of wheat quality parameters viz., grain protein content (GPC), grain starch content (GSC), and grain hardness (GH), we evaluated 372 diverse European wheat varieties in replicated field trials in up to eight environments. We observed that all of the investigated traits hold a wide and significant genetic variation, and a significant negative correlation exists between GPC and GSC plus grain yield. Our association analyses based on 26,694 high-quality single nucleotide polymorphic markers revealed a strong quantitative genetic nature of GPC and GSC with associations on groups 2, 3, and 6 chromosomes. The identification of known Puroindoline-b gene for GH provided a positive analytic proof for our studies. We report that a locus QGpc.ipk-6A controls both GPC and GSC with opposite allelic effects. Based on wheat's reference and pan-genome sequences, the physical characterization of two loci viz., QGpc.ipk-2B and QGpc.ipk-6A facilitated the identification of the candidate genes for GPC. Furthermore, by exploiting additive and epistatic interactions of loci, we evaluated the prospects of predictive breeding for the investigated traits that suggested its efficient use in the breeding programs.

Poor morphological and molecular differentiation in recently diversified lineages is a widespread phenomenon in plants. Phylogenetic relationships within such species complexes are often difficult to resolve because of the low variability in traditional molecular loci. Furthermore, biological phenomena responsible for topological incongruence such as Incomplete Lineage Sorting (ILS) and hybridisation complicate the resolution of phylogenetic relationships among closely related taxa. In this study, we employ a Genotyping-by-sequencing (GBS) approach to disentangle evolutionary relationships within a species complex belonging to the Neotropical orchid genus Cycnoches. This complex includes seven taxa distributed through Central America and the Colombian Chocó, and is nested within a clade estimated to have first diversified in the early Quaternary. Previous phylogenies inferred from few loci failed to provide support for internal relationships within the complex. Our Neighbour-net and coalescent-based analyses inferred from ca. 13,000 GBS loci obtained from 31 individuals belonging to six of the seven traditionally accepted Cycnoches taxa provided a robust phylogeny for this group. The genus Cycnoches includes three main clades that are further supported by morphological traits and geographic distributions. Similarly, a topology reconstructed through maximum likelihood (ML) inference of concatenated GBS loci produced results that are comparable with those reconstructed through coalescence and network-based methods. Our comparative phylogenetic informativeness analyses suggest that the low support evident in the ML phylogeny might be attributed to the abundance of uninformative GBS loci, which can account for up to 50% of the total number of loci recovered. The phylogenomic framework provided here, as well as morphological evidence and geographical patterns, suggest that the six entities previously thought to be different species or subspecies might actually represent only three distinct segregates. We further discuss the limited phylogenetic informativeness found in our GBS approach and its utility to disentangle relationships within recent and rapidly evolving species complexes. Our study is the first to demonstrate the utility of GBS data to reconstruct relationships within young (~2 Ma) Neotropical plant clades, opening new avenues for studies of species complexes that populate the species-rich orchid family UCR::Vicerrectoría de Investigación::Unidades de Investigación::Ciencias Agroalimentarias::Jardín Botánico Lankester (JBL)

Background Chamomile (Matricaria recutita L.) has a long history of use in herbal medicine with various applications, and the flower heads contain numerous secondary metabolites which are medicinally active. In the major crop plants, next generation sequencing (NGS) approaches are intensely applied to exploit genetic resources, to develop genomic resources and to enhance breeding. Here, genotyping-by-sequencing (GBS) has been used in the non-model medicinal plant chamomile to evaluate the genetic structure of the cultivated varieties/populations, and to perform genome wide association study (GWAS) focusing on genes with large effect on flowering time and the medicinally important alpha-bisabolol content. Results GBS analysis allowed the identification of 6495 high-quality SNP-markers in our panel of 91 M. recutita plants from 33 origins (2–4 genotypes each) and 4 M. discoidea plants as outgroup, grown in the greenhouse in Gatersleben, Germany. M. recutita proved to be clearly distinct from the outgroup, as was demonstrated by different cluster and principal coordinate analyses using the SNP-markers. Chamomile genotypes from the same origin were mostly genetically similar. Model-based cluster analysis revealed one large group of tetraploid genotypes with low genetic differentiation including 39 plants from 14 origins. Tetraploids tended to display lower genetic diversity than diploids, probably reflecting their origin by artificial polyploidisation from only a limited set of genetic backgrounds. Analyses of flowering time demonstrated that diploids generally flowered earlier than tetraploids, and the analysis of alpha-bisabolol identified several tetraploid genotypes with a high content. GWAS identified highly significant (P

We dissected the genetic basis of total spikelet number (TSN) along with other traits, viz. spike length (SL) and flowering time (FT) in a panel of 518 elite European winter wheat varieties. Genome-wide association studies (GWAS) based on 39,908 SNP markers revealed highly significant quantitative trait loci (QTL) for TSN on chromosomes 2D, 7A, and 7B, for SL on 5A, and FT on 2D, with 2D-QTL being the functional marker for the gene Ppd-D1. The physical region of the 7A-QTL for TSN revealed the presence of a wheat ortholog (TaAPO-A1) to APO1–a rice gene that positively controls the spikelet number on the panicles. Interspecific analyses of the TaAPO-A1 orthologs showed that it is a highly conserved gene important for floral development and present in a wide range of terrestrial plants. Intraspecific studies of the TaAPO-A1 across wheat genotypes revealed a polymorphism in the conserved F-box domain, defining two haplotypes. A KASP marker developed on the polymorphic site showed a highly significant association of TaAPO-A1 with TSN, explaining 23.2% of the total genotypic variance. Also, the TaAPO-A1 alleles showed weak but significant differences for SL and grain yield. Our results demonstrate the importance of wheat sequence resources to identify candidate genes for important traits based on genetic analyses.

Many conflicting hypotheses regarding the relationships among crops and wild species closely related to wheat (the genera Aegilops, Amblyopyrum, and Triticum) have been postulated. The contribution of hybridization to the evolution of these taxa is intensely discussed. To determine possible causes for this, and provide a phylogeny of the diploid taxa based on genome-wide sequence information, independent data was obtained from genotyping-by-sequencing and a target-enrichment experiment that returned 244 low-copy nuclear loci. The data were analyzed with Bayesian, likelihood and coalescent-based methods. D statistics were used to test if incomplete lineage sorting alone or together with hybridization is the source for incongruent gene trees. Here we present the phylogeny of all diploid species of the wheat wild relatives. We hypothesize that most of the wheat-group species were shaped by a primordial homoploid hybrid speciation event involving the ancestral Triticum and Am. muticum lineages to form all other species but Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa but Triticum. Mostly progenitors of the extant species were involved in these processes, while recent interspecific gene flow seems insignificant. The composite nature of many genomes of wheat group taxa results in complicated patterns of diploid contributions when these lineages are involved in polyploid formation, which is, for example, the case in the tetra-and hexaploid wheats. Our analysis provides phylogenetic relationships and a testable hypothesis for the genome compositions in the basic evolutionary units within the wheat group of Triticeae.

We dissected the genetic basis of total spikelet number (TSN) along with other traits, namely spike length (SL) and flowering time (FT) in a panel of 518 elite European winter wheat varieties. Genome-wide association studies based on 39,908 SNP markers revealed highly significant quantitative trait loci (QTL) for TSN on chromosomes 2D, 7A, and 7B, for SL on 5A, and FT on 2D, with 2D-QTL being the functional marker for the genePpd-D1. The physical region of the 7A-QTL for TSN revealed the presence of an ortholog toAPO1– a rice gene that positively controls spikelet number on panicles. Interspecific analyses ofTaAPO-A1orthologs showed that it is a highly conserved gene important for floral development, and present in a wide range of terrestrial plants. Intraspecific studies of the wheat orthologTaAPO-A1across wheat genotypes revealed a polymorphism in the highly conserved F-box domain, defining two haplotypes. A KASP maker developed on the polymorphic site showed a highly significant association ofTaAPO-A1with TSN, explaining 23.2% of the genotypic variance. Also, theTaAPO-A1alleles showed weak but significant differences for SL and grain yield. Our results demonstrate the importance of wheat sequence resources to identify candidate genes for important traits based on genetic analyses.

Poor morphological and molecular differentiation in recently diversified lineages is a widespread phenomenon in plants. Phylogenetic relationships within such species complexes are often difficult to resolve because of the low variability in traditional molecular loci, as well as various other biological phenomena responsible for topological incongruence such as ILS and hybridization. In this study, we employ a Genotyping-by-sequencing (GBS) approach to disentangle evolutionary relationships within a species complex belonging to the Neotropical orchid genus Cycnoches. The complex includes seven taxa distributed in Central America and the adjacent Chocó biogeographic region, nested within a clade estimated to have first diversified in the early Quaternary. Previous phylogenies inferred from a handful of loci failed to provide support for internal relationships within the complex. Our Neighbor-net and coalescent-based analyses inferred from ca. 13,000 GBS loci obtained from 31 individuals belonging to six of the seven traditionally accepted Cycnoches species provided a robustly supported network. The resulting three main clades are corroborated by morphological traits and geographical distributions. Similarly, Maximum Likelihood (ML) inferences of concatenated GBS-loci produced results comparable with those derived from coalescence and network-based methods, albeit always with poor statistical support. The low support evident in the ML phylogeny might be attributed to the abundance of uninformative GBS loci, which can account for up to 50% of the total number of loci recovered. The phylogenomic framework provided here, as well as morphological evidence and geographical patterns, suggest that the six entities previously thought to be different species might actually represent only three distinct segregates. Our study is the first to demonstrate the utility of GBS data in phylogenomic research of a very young Neotropical plant clade (~2 Ma), and it paves the way for the study of the many other species complexes that populate the species-rich orchid family.

Background The large and highly repetitive genomes of the cultivated species Hordeum vulgare (barley), Triticum aestivum (wheat), and Secale cereale (rye) belonging to the Triticeae tribe of grasses appear to be particularly rich in gene-like sequences including partial duplicates. Most of them have been classified as putative pseudogenes. In this study we employ transient and stable gene silencing- and over-expression systems in barley to study the function of HvARM1 (for H. vulgare Armadillo 1), a partial gene duplicate of the U-box/armadillo-repeat E3 ligase HvPUB15 (for H. vulgare Plant U-Box 15). Results The partial ARM1 gene is derived from a gene-duplication event in a common ancestor of the Triticeae and contributes to quantitative host as well as nonhost resistance to the biotrophic powdery mildew fungus Blumeria graminis. In barley, allelic variants of HvARM1 but not of HvPUB15 are significantly associated with levels of powdery mildew infection. Both HvPUB15 and HvARM1 proteins interact in yeast and plant cells with the susceptibility-related, plastid-localized barley homologs of THF1 (for Thylakoid formation 1) and of ClpS1 (for Clp-protease adaptor S1) of Arabidopsis thaliana. A genome-wide scan for partial gene duplicates reveals further events in barley resulting in stress-regulated, potentially neo-functionalized, genes. Conclusion The results suggest neo-functionalization of the partial gene copy HvARM1 increases resistance against powdery mildew infection. It further links plastid function with susceptibility to biotrophic pathogen attack. These findings shed new light on a novel mechanism to employ partial duplication of protein-protein interaction domains to facilitate the expansion of immune signaling networks. Electronic supplementary material The online version of this article (10.1186/s13059-018-1472-7) contains supplementary material, which is available to authorized users.

The taxa of Astragalus section Hymenostegis are an important element of mountainous and steppe habitats in Southwest Asia. A phylogenetic hypothesis of sect. Hymenostegis has been obtained from nuclear ribosomal DNA internal transcribed spacer (ITS) and plastid ycf1 sequences of up to 303 individuals from 106 species, including all 89 taxa currently assigned to sect. Hymenostegis, 14 species of other Astragalus sections, and two species of Oxytropis and one Biserrula designated as outgroups. Bayesian phylogenetic inference and parsimony analyses reveal that three species from two other closely related sections group within sect. Hymenostegis, making the section paraphyletic. DNA sequence diversity is generally very low among Hymenostegis taxa, which is consistent with recent diversification of the section. We estimate that diversification in sect. Hymenostegis occurred in the middle to late Pleistocene, with many species arising only during the last one million years, when environmental conditions in the mountain regions of Southwest and Central Asia cycled repeatedly between dry and more humid conditions.

Hybrid wheat breeding is gaining prominence worldwide because it ensures higher and more static yield than conventionally bred varieties. The cleistogamous floral architecture of wheat (Triticum aestivum L.) impedes anthers inside the floret, making it largely an inbreeder. For hybrid seed production, high anther extrusion is needed to promote cross pollination and to ensure a high level of pollen availability for the seed plant. This study, therefore, aimed at the genetic dissection of anther extrusion (AE) in panels of spring (SP), and winter wheat (WP) accessions by genome wide association studies (GWAS). We performed GWAS to identify the SNP markers potentially linked with AE in each panel separately. Phenotypic data were collected for 3 years for each panel. The average levels of Pearson's correlation (r) among all years and their best linear unbiased estimates (BLUEs) within both panels were high (r(SP) = 0.75, P < 0.0001;r(WP) = 0.72, P < 0.0001). Genotypic data (with minimum of 0.05 minor allele frequency applied) included 12,066 and 12,191 SNP markers for SP and WP, respectively. Both genotypes and environment influenced the magnitude of AE. In total, 23 significant (|log10(P)| > 3.0) marker trait associations (MTAs) were detected (SP = 11; WP = 12). Anther extrusion behaved as a complex trait with significant markers having either favorable or unfavorable additive effects and imparting minor to moderate levels of phenotypic variance (R2(SP) = 9.75−14.24%; R2 (WP) = 9.44−16.98%). All mapped significant markers as well as the markers within their significant linkage disequilibrium (r2 ≥ 0.30) regions were blasted against wheat genome assembly (IWGSC1+popseq) to find the corresponding genes and their high confidence descriptions were retrieved. These genes and their orthologs in Hordeum vulgare, Brachypodium distachyon, Oryza sativa, and Sorghum bicolor revealed syntenic genomic regions potentially involved in flowering-related traits. Moreover, the expression data of these genes suggested potential candidates for AE. Our results suggest that the use of significant markers can help to introduce AE in high yielding varieties to increase cross fertilization rates and improve hybrid-seed production in wheat.

Background Triticeae, the tribe of wheat grasses, harbours the cereals barley, rye and wheat and their wild relatives. Although economically important, relationships within the tribe are still not understood. We analysed the phylogeny of chloroplast lineages among nearly all monogenomic Triticeae taxa and polyploid wheat species aiming at a deeper understanding of the tribe’s evolution. We used on- and off-target reads of a target-enrichment experiment followed by Illumina sequencing. Results The read data was used to assemble the plastid locus ndhF for 194 individuals and the whole chloroplast genome for 183 individuals, representing 53 Triticeae species and 15 genera. We conducted Bayesian and multispecies coalescent analyses to infer relationships and estimate divergence times of the taxa. We present the most comprehensive dated Triticeae chloroplast phylogeny and review previous hypotheses in the framework of our results. Monophyly of Triticeae chloroplasts could not be confirmed, as either Bromus or Psathyrostachys captured a chloroplast from a lineage closely related to a Bromus-Triticeae ancestor. The most recent common ancestor of Triticeae occurred approximately between ten and 19 million years ago. Conclusions The comparison of the chloroplast phylogeny with available nuclear data in several cases revealed incongruences indicating past hybridizations. Recent events of chloroplast capture were detected as individuals grouped apart from con-specific accessions in otherwise monopyhletic groups. Electronic supplementary material The online version of this article (doi:10.1186/s12862-017-0989-9) contains supplementary material, which is available to authorized users.

Polyploidization is an important speciation mechanism in the barley genus Hordeum. To analyze evolutionary changes after allopolyploidization, knowledge of parental relationships is essential. One chloroplast and 12 nuclear single-copy loci were amplified by polymerase chain reaction (PCR) in all Hordeum plus six out-group species. Amplicons from each of 96 individuals were pooled, sheared, labeled with individual-specific barcodes and sequenced in a single run on a 454 platform. Reference sequences were obtained by cloning and Sanger sequencing of all loci for nine supplementary individuals. The 454 reads were assembled into contigs representing the 13 loci and, for polyploids, also homoeologues. Phylogenetic analyses were conducted for all loci separately and for a concatenated data matrix of all loci. For diploid taxa, a Bayesian concordance analysis and a coalescent-based dated species tree was inferred from all gene trees. Chloroplast matK was used to determine the maternal parent in allopolyploid taxa. The relative performance of different multilocus analyses in the presence of incomplete lineage sorting and hybridization was also assessed. The resulting multilocus phylogeny reveals for the first time species phylogeny and progenitor-derivative relationships of all di- and polyploid Hordeum taxa within a single analysis. Our study proves that it is possible to obtain a multilocus species-level phylogeny for di- and polyploid taxa by combining PCR with next-generation sequencing, without cloning and without creating a heavy load of sequence data.

Polyploidization is a major mechanism of speciation in plants. Within the barley genus Hordeum, approximately half of the taxa are polyploids. While for diploid species a good hypothesis of phylogenetic relationships exists, there is little information available for the polyploids (4×, 6×) of Hordeum. Relationships among all 33 diploid and polyploid Hordeum species were analyzed with the low-copy nuclear marker region TOPO6 for 341 Hordeum individuals and eight outgroup species. PCR products were either directly sequenced or cloned and on average 12 clones per individual were included in phylogenetic analyses. In most diploid Hordeum species TOPO6 is probably a single-copy locus. Most sequences found in polyploid individuals phylogenetically cluster together with sequences derived from diploid species and thus allow the identification of parental taxa of polyploids. Four groups of sequences occurring only in polyploid taxa are interpreted as footprints of extinct diploid taxa, which contributed to allopolyploid evolution. Our analysis identifies three key species involved in the evolution of the American polyploids of the genus. (i) All but one of the American tetraploids have a TOPO6 copy originating from the Central Asian diploid H. roshevitzii, the second copy clustering with different American diploid species. (ii) All hexaploid species from the New World have a copy of an extinct close relative of H. californicum and (iii) possess the TOPO6 sequence pattern of tetraploid H. jubatum, each with an additional copy derived from different American diploids. Tetraploid H. bulbosum is an autopolyploid, while the assumed autopolyploid H. brevisubulatum (4×, 6×) was identified as allopolyploid throughout most of its distribution area. The use of a proof-reading DNA polymerase in PCR reduced the proportion of chimerical sequences in polyploids in comparison to Taq polymerase.