Unleashing floret fertility in wheat through the mutation of a homeobox gene.

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.

Genome-wide sequence information reveals recurrent hybridization among diploid wheat wild relatives.

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 were obtained from genotyping-by-sequencing and a target-enrichment experiment that returned 244 low-copy nuclear loci. The data were analyzed using 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 except Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa except 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 for tetraploid 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.

Genomic, karyological and morphological changes of South American garlics (Ipheion) provide insights into mechanisms of speciation in the Pampean region.

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.

Resolving relationships in an exceedingly young Neotropical orchid lineage using Genotyping-by-sequencing data.

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.

Dated tribe-wide whole chloroplast genome phylogeny indicates recurrent hybridizations within Triticeae.

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.

Use of genotyping-by-sequencing to determine the genetic structure in the medicinal plant chamomile, and to identify flowering time and alpha-bisabolol associated SNP-loci by genome-wide association mapping.

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 < 0.01) SNPs for flowering time (9) and alpha-bisabolol (71). One sequence harbouring SNPs associated with flowering time was described to play a role in self-pollination in Arabidopsis thaliana, whereas four sequences harbouring SNPs associated with alpha-bisabolol were identified to be involved in plant biotic and abiotic stress response in various plants species. CONCLUSIONS: The first genomic resource for future applications to enhance breeding in chamomile was created, andanalyses of diversity will facilitate the exploitation of these genetic resources. The GWAS data pave the way for future research towards the genetics underlying important traits in chamomile, the identification of marker-trait associations, and development of reliable markers for practical breeding.

Species-Level Phylogeny and Polyploid Relationships in Hordeum (Poaceae) Inferred by Next-Generation Sequencing and In Silico Cloning of Multiple Nuclear Loci.

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.

Dated Phylogeny of Astragalus Section Stereothrix (Fabaceae) and Allied Taxa in the Hypoglottis Clade.

The Astragalus subgenus Hypoglottis Bunge, which consists of several sections, is one of the taxonomically most complicated groups in the genus. The Astragalus section Stereothrix Bunge belongs to this subgenus and is a significant element of the Irano-Turanian floristic region. A molecular phylogenetic analysis of this section and its closely related taxa using nuclear ribosomal DNA internal transcribed spacers (ITS) and external transcribed spacer (ETS) regions as well as plastid matK sequences were conducted. Parsimony analyses and Bayesian phylogenetic inference revealed that the section is not monophyletic in its current form, as some taxa belonging to closely related sections such as Hypoglottidei DC. and the Malacothrix Bunge group within the sect. Stereothrix render it paraphyletic. Moreover, species groups belonging to sect. Stereothrix are placed in different clades within the phylogenetic tree of subgenus Hypoglottis, which indicates polyphyly, i.e., multiple independent origins of taxa placed in the sect. Stereothrix. Molecular dating of the group estimated an age of 3.62 (1.73-5.62) My for this assemblage with the major diversification events happening during the last 2 My. Many species groups separated only within the last 0.5 to 1 My. Based on morphological and molecular data, we discuss the phylogenetic relationships of the groups and synonymy of species. In addition, the included taxa of sect. Hypoglottidei are not monophyletic and include species belonging to sects. Hololeuce, Koelziana, Malacothrix, Onobrychoideae, and Ornithodpodium group within the sect. Stereothrix taxa. We conclude that only an analysis including all groups and nearly all species of the sections within the Hypoglottis clade can finally result in an new evolutionary-based system for these taxa.

Linkage mapping identifies a non-synonymous mutation in FLOWERING LOCUS T (FT-B1) increasing spikelet number per spike.

Total spikelet number per spike (TSN) is a major component of spike architecture in wheat (Triticum aestivum 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 FLOWERING LOCUS T (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.

Prospects of GWAS and predictive breeding for European winter wheat's grain protein content, grain starch content, and grain hardness.

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.

TaAPO-A1, an ortholog of rice ABERRANT PANICLE ORGANIZATION 1, is associated with total spikelet number per spike in elite European hexaploid winter wheat (Triticum aestivum L.) varieties.

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.

Evolutionarily conserved partial gene duplication in the Triticeae tribe of grasses confers pathogen resistance.

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.

Genetic Architecture of Anther Extrusion in Spring and Winter Wheat.

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.

Molecular phylogeny and divergence times of Astragalus section Hymenostegis: An analysis of a rapidly diversifying species group in Fabaceae.

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.

Progenitor-derivative relationships of Hordeum polyploids (Poaceae, Triticeae) inferred from sequences of TOPO6, a nuclear low-copy gene region.

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.