Multigenic phylogeny and analysis of tree incongruences in Triticeae (Poaceae)

Identification and molecular characterization of novel sucrose transporters in the hexaploid wheat (Triticum aestivum L.)

Common wheat (Triticum aestivum) is a major cereal crop grown and consumed globally. Recent advances in sequencing technology have facilitated the exploration of large and repetitive genomes. Plant sucrose transporter (SUT) genes are vital components of energy transport systems that play prominent roles in various plant functions, such as signaling and stress regulation. In this study, we identified and analyzed five novel sucrose transporter genes in wheat. The wheat sucrose transporter genes were divided into five clades based on their phylogenetic relationships. Synteny analysis revealed that synteny in the genome is highly conserved between wheat and rye, barley, and Brachypodium. Furthermore, the cis-element analysis indicated that sucrose transporter genes might be regulated by light and some phytohormone-related transcriptional factors. Overall, plant tissue-specific gene expression revealed enhanced expression of the transporter genes in the root and stem, whereas they were differentially expressed under abiotic stress treatments (cold, heat, NaCl, PEG-6000, and sucrose). These results indicate that each TaSUT gene may play a crucial role in stabilizing plants under stress by actively regulating the energy demands of cells. The findings of this study may provide a basis for further research on sucrose transporters and their significant roles in plant energy metabolism as well as in abiotic stress response, signaling, and regulation.

The Plastome Sequences of Triticum sphaerococcum (ABD) and Triticum turgidum subsp. durum (AB) Exhibit Evolutionary Changes, Structural Characterization, Comparative Analysis, Phylogenomics and Time Divergence

The mechanism and course of Triticum plastome evolution is currently unknown; thus, it remains unclear how Triticum plastomes evolved during recent polyploidization. Here, we report the complete plastomes of two polyploid wheat species, Triticum sphaerococcum (AABBDD) and Triticum turgidum subsp. durum (AABB), and compare them with 19 available and complete Triticum plastomes to create the first map of genomic structural variation. Both T. sphaerococcum and T. turgidum subsp. durum plastomes were found to have a quadripartite structure, with plastome lengths of 134,531 bp and 134,015 bp, respectively. Furthermore, diploid (AA), tetraploid (AB, AG) and hexaploid (ABD, AGA^m) Triticum species plastomes displayed a conserved gene content and commonly harbored an identical set of annotated unique genes. Overall, there was a positive correlation between the number of repeats and plastome size. In all plastomes, the number of tandem repeats was higher than the number of palindromic and forward repeats. We constructed a Triticum phylogeny based on the complete plastomes and 42 shared genes from 71 plastomes. We estimated the divergence of Hordeum vulgare from wheat around 11.04-11.9 million years ago (mya) using a well-resolved plastome tree. Similarly, Sitopsis species diverged 2.8-2.9 mya before Triticum urartu (AA) and Triticum monococcum (AA). Aegilops speltoides was shown to be the maternal donor of polyploid wheat genomes and diverged ~0.2-0.9 mya. The phylogeny and divergence time estimates presented here can act as a reference framework for future studies of Triticum evolution.

When tropical and subtropical congeners met: Multiple ancient hybridization events within Eriobotrya in the Yunnan-Guizhou Plateau, a tropical-subtropical transition area in China

Global climate changes during the Miocene may have created ample opportunities for hybridization between members of tropical and subtropical biomes at the boundary between these zones. Yet, very few studies have explored this possibility. The Yunnan-Guizhou Plateau (YGP) in Southwest China is a biodiversity hotspot for vascular plants, located in a transitional area between the floristic regions of tropical Southeast Asia and subtropical East Asia. The genus Eriobotrya (Rosaceae) comprises both tropical and subtropical taxa, with 12 species recorded in the YGP, making it a suitable basis for testing the hypothesis of between-biome hybridization. Therefore, we surveyed the evolutionary history of Eriobotrya by examining three chloroplast regions and five nuclear genes for 817 individuals (47 populations) of 23 Eriobotrya species (including 19 populations of 12 species in the YGP), plus genome re-sequencing of 33 representative samples. We concluded that: (1) phylogenetic positions for 16 species exhibited strong cytonuclear conflicts, most probably due to ancient hybridization; (2) the YGP is a hotspot for hybridization, with 11 species showing clear evidence of chloroplast capture; and (3) Eriobotrya probably originated in tropical Asia during the Eocene. From the Miocene onwards, the intensification of the Eastern Asia monsoon and global cooling may have shifted the tropical-subtropical boundary and caused secondary contact between species, thus providing ample opportunity for hybridization and diversification of Eriobotrya, especially in the YGP. Our study highlights the significant role that paleoclimate changes probably played in driving hybridization and generating rich species diversity in climate transition zones.

The decreased expression of GW2 homologous genes contributed to the increased grain width and thousand‑grain weight in wheat-Dasypyrum villosum 6VS·6DL translocation lines

This study demonstrated that the aberrant transcription of DvGW2 contributed to the increased grain width and thousand-grain weight in wheat-Dasypyrum villosum T6VS·6DL translocation lines. Due to the high immunity to powdery mildew, Dasypyrum villosum 6VS has been one of the most successful applications of the wild relatives in modern wheat breeding. Along with the desired traits, side-effects could be brought when large alien chromosome fragments are introduced into wheat, but little is known about effects of 6VS on agronomic traits. Here, we found that T6VS·6DL translocation had significantly positive effects on grain weight, plant heightand spike length, and small negative effects on total spikelet number and spikelet compactness using recipient and wheat-D. villosum T6VS·6DL allohexaploid wheats, Wan7107 and Pm97033. Further analysis showed that the 6VS segment might exert direct genetic effect on grain width, then driving the increase of thousand-grain weight. Furthermore, comparative transcriptome analysis identified 2549 and 1282 differentially expressed genes (DEGs) and 2220 and 1496 specifically expressed genes (SEGs) at 6 days after pollination (DAP) grains and 15 DAP endosperms, respectively. Enrichment analysis indicated that the process of cell proliferation category was over-represented in the DEGs. Notably, two homologous genes, TaGW2-D1 and DvGW2, were identified as putative candidate genes associated with grain weight and yield. The expression analysis showed that DvGW2 had an aberrant expression in Pm97033, resulting in significantly lower total expression level of GW2 than Wan7107, which drives the increase of grain weight and width in Pm97033. Collectively, our data indicated that the compromised expression of DvGW2 is critical for increased grain width and weight in T6VS·6DL translocation lines.

Genomic Variation Shaped by Environmental and Geographical Factors in Prairie Cordgrass Natural Populations Collected across Its Native Range in the USA

Prairie cordgrass (Spartina pectinata Link) is a native perennial warm-season (C4) grass common in North American prairies. With its high biomass yield and abiotic stress tolerance, there is a high potential of developing prairie cordgrass for conservation practices and as a dedicated bioenergy crop for sustainable cellulosic biofuel production. However, as with many other undomesticated grass species, little information is known about the genetic diversity or population structure of prairie cordgrass natural populations as compared to their ecotypic and geographic adaptation in North America. In this study, we sampled and characterized a total of 96 prairie cordgrass natural populations with 9315 high quality SNPs from a genotyping-by-sequencing (GBS) approach. The natural populations were collected from putative remnant prairie sites throughout the Midwest and Eastern USA, which are the major habitats for prairie cordgrass. Partitioning of genetic variance using SNP marker data revealed significant variance among and within populations. Two potential gene pools were identified as being associated with ploidy levels, geographical separation, and climatic separation. Geographical factors such as longitude and altitude, and environmental factors such as annual temperature, annual precipitation, temperature of the warmest month, precipitation of the wettest month, precipitation of Spring, and precipitation of the wettest month are important in affecting the intraspecific distribution of prairie cordgrass. The divergence of prairie cordgrass natural populations also provides opportunities to increase breeding value of prairie cordgrass as a bioenergy and conservation crop.

Characterizing chloroplast genomes and inferring maternal divergence of the Triticum-Aegilops complex

The Triticum (wheat)–Aegilops (goatgrass) complex has been extensively studied, but the evolutionary history of polyploid wheats has not been fully elucidated. The chloroplast (cp) with maternal inheritance and homoplasy can simplify the sequence-based evolutionary inferences, but informative inferences would require a complete and accurate cp genome sequence. In this study, 16 cp genomes representing five Aegilops and 11 Triticum species and subspecies were sequenced, assembled and annotated, yielding five novel circular cp genome sequences. Analyzing the assembled cp genomes revealed no marked differences in genome structure and gene arrangement across the assayed species. A polymorphism analysis of 72 published cp genome sequences representing 10 Aegilops and 15 Triticum species and subspecies detected 1183 SNPs and 1881 SSRs. More than 80% SNPs detected resided on the downstream and upstream gene regions and only 2.78% or less SNPs were predicted to be deleterious. The largest nucleotide diversity was observed in the short single-copy genomic region. Relatively weak selection pressure on cp coding genes was detected. Different phylogenetic analyses confirmed that the maternal divergence of the Triticum–Aegilops complex had three deep lineages each representing a diploid species with nuclear A, B, or D genome. Dating the maternal divergence yielded age estimates of divergence that matched well with those reported previously. The divergence between emmer and bread wheats occurred at 8200–11,200 years ago. These findings are useful for further genomic studies, provide insight into cp genome evolvability and allow for better understanding of the maternal divergence of the Triticum–Aegilops complex.

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.

Hybridization preceded radiation in diploid wheats

Evolutionary relationships among the Aegilops-Triticum relatives of cultivated wheats have been difficult to resolve owing to incomplete lineage sorting and reticulate evolution. Recent studies have suggested that the wheat D-genome lineage (progenitor of Ae. tauschii) originated through homoploid hybridization between the A-genome lineage (progenitor of Triticum s.str.) and the B-genome lineage (progenitor of Ae. speltoides). This scenario of reticulation has been debated, calling for adequate phylogenetic analyses based on comprehensive sampling. To reconstruct the evolution of Aegilops-Triticum diploids, we here combined high-throughput sequencing of 38 nuclear low-copy loci of multiple accessions of all 13 species with inferences of the species phylogeny using the full-parameterized MCMC_SEQ method. Phylogenies recovered a monophyletic Aegilops-Triticum lineage that began diversifying ~6.6 Ma ago and gave rise to four sublineages, i.e. the A- (2 species), B- (1 species), D- (9 species) and T- (Ae. mutica) genome lineage. Full-parameterized phylogenies as well as patterns of tree dilation and tree compression supported a hybrid origin of the D-genome lineage from A and B ~3.0-4.0 Ma ago, and did not indicate additional hybridization events. Conflicting ABBA-BABA tests suggestive of further reticulation were shown here to result from ancestral population structure rather than hybridization. This comprehensive and dated phylogeny of wheat relatives indicates that the origin of the hybrid D-genome was followed by intense diversification into the majority of extant diploid as well as allopolyploid wild wheats.

Pervasive hybridizations in the history of wheat relatives

Cultivated wheats are derived from an intricate history of three genomes, A, B, and D, present in both diploid and polyploid species. It was recently proposed that the D genome originated from an ancient hybridization between the A and B lineages. However, this result has been questioned, and a robust phylogeny of wheat relatives is still lacking. Using transcriptome data from all diploid species and a new methodological approach, our comprehensive phylogenomic analysis revealed that more than half of the species descend from an ancient hybridization event but with a more complex scenario involving a different parent than previously thought-Aegilops mutica, an overlooked wild species-instead of the B genome. We also detected other extensive gene flow events that could explain long-standing controversies in the classification of wheat relatives.

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.

Characterization of Chromosomal Rearrangement in New Wheat—Thinopyrum intermedium Addition Lines Carrying Thinopyrum—Specific Grain Hardness Genes

The wild species, Thinopyrum intermedium. (Genome StStJSJSJJ), serves as a valuable germplasm resource providing novel genes for wheat improvement. In the current study, non-denaturing fluorescence in situ hybridization (ND-FISH) with multiple probes and comparative molecular markers were applied to characterize two wheat-Th. intermedium chromosome additions. Sequential ND-FISH with new labeled Th. intermedium specific oligo-probes were used to precisely determine the chromosomal constitution of Th. intermedium, wheat—Th. intermedium partial amphiploids and addition lines Hy36 and Hy37. The ND-FISH results showed that the added JS-St translocated chromosomes in Hy36 had minor Oligo-5S ribosomal DNA (rDNA) signals at the short arm, while a pair of J-St chromosomes in Hy37 had major Oligo-pTa71 and minor Oligo-5S rDNA signals. The 90K SNP array and PCR-based molecular markers that mapped on wheat linkage group 5 and 3 facilitated the identification of Thinopyrum chromosome introgressions in the addition lines, and confirmed that added chromosomes in Hy36 and Hy37 were 5JSS.3StS and 5JS.3StS, respectively. Complete coding sequences at the paralogous puroindoline-a (Pina) loci from Th. intermedium were cloned and localized on the short arm of chromosome 5JS of Hy36. Line Hy36 showed a reduction in the hardness index, which suggested that Th. intermedium-specific Pina gene sequences may be associated with the softness trait in wheat background. The molecular cytogenetic identification of novel wheat—Th. intermedium derivatives indicated that the frequent chromosome rearrangement occurred in the progenies of wheat-Thinopyrum hybridization. The new wheat-Thinopyrum derived lines may increase the genetic diversity for wheat breeding.

Sequence Diversity and Identification of Novel Puroindoline and Grain Softness Protein Alleles in Elymus, Agropyron and Related Species

The puroindoline proteins, PINA and PINB, which are encoded by the Pina and Pinb genes located at the Ha locus on chromosome 5D of bread wheat, are considered to be the most important determinants of grain hardness. However, the recent identification of Pinb-2 genes on group 7 chromosomes has stressed the importance of considering the effects of related genes and proteins. Several species related to wheat (two diploid Agropyron spp., four tetraploid Elymus spp. and five hexaploid Elymus and Agropyron spp.) were therefore analyzed to identify novel variation in Pina, Pinb and Pinb-2 genes which could be exploited for the improvement of cultivated wheat. A novel sequence for the Pina gene was detected in Elymus burchan-buddae, Elymus dahuricus subsp. excelsus and Elymus nutans and novel PINB sequences in Elymus burchan-buddae, Elymus dahuricus subsp. excelsus, and Elymus nutans. A novel PINB-2 variant was also detected in Agropyron repens and Elymus repens. The encoded proteins detected all showed changes in the tryptophan-rich domain as well as changes in and/or deletions of basic and hydrophobic residues. In addition, two new AGP sequences were identified in Elymus nutans and Elymus wawawaiensis. The data presented therefore highlight the sequence diversity in this important gene family and the potential to exploit this diversity to modify grain texture and end-use quality in wheat.

A 250 plastome phylogeny of the grass family (Poaceae): topological support under different data partitions

The systematics of grasses has advanced through applications of plastome phylogenomics, although studies have been largely limited to subfamilies or other subgroups of Poaceae. Here we present a plastome phylogenomic analysis of 250 complete plastomes (179 genera) sampled from 44 of the 52 tribes of Poaceae. Plastome sequences were determined from high throughput sequencing libraries and the assemblies represent over 28.7 Mbases of sequence data. Phylogenetic signal was characterized in 14 partitions, including (1) complete plastomes; (2) protein coding regions; (3) noncoding regions; and (4) three loci commonly used in single and multi-gene studies of grasses. Each of the four main partitions was further refined, alternatively including or excluding positively selected codons and also the gaps introduced by the alignment. All 76 protein coding plastome loci were found to be predominantly under purifying selection, but specific codons were found to be under positive selection in 65 loci. The loci that have been widely used in multi-gene phylogenetic studies had among the highest proportions of positively selected codons, suggesting caution in the interpretation of these earlier results. Plastome phylogenomic analyses confirmed the backbone topology for Poaceae with maximum bootstrap support (BP). Among the 14 analyses, 82 clades out of 309 resolved were maximally supported in all trees. Analyses of newly sequenced plastomes were in agreement with current classifications. Five of seven partitions in which alignment gaps were removed retrieved Panicoideae as sister to the remaining PACMAD subfamilies. Alternative topologies were recovered in trees from partitions that included alignment gaps. This suggests that ambiguities in aligning these uncertain regions might introduce a false signal. Resolution of these and other critical branch points in the phylogeny of Poaceae will help to better understand the selective forces that drove the radiation of the BOP and PACMAD clades comprising more than 99.9% of grass diversity.

Conserved molecular structure of the centromeric histone CENH3 in Secale and its phylogenetic relationships

It has been repeatedly demonstrated that the centromere-specific histone H3 (CENH3), a key component of the centromere, shows considerable variability between species within taxa. We determined the molecular structure and phylogenetic relationships of CENH3 in 11 Secale species and subspecies that possess distinct pollination systems and are adapted to a wide range of abiotic and biotic stresses. The rye (Secale cereale) genome encodes two paralogous CENH3 genes, which differ in intron-exon structure and are transcribed into two main forms of the protein, αCENH3 and βCENH3. These two forms differ in size and amino acid substitutions. In contrast to the reported differences in CENH3 structure between species within other taxa, the main forms of this protein in Secale species and subspecies have a nearly identical structure except some nonsynonymous substitutions. The CENH3 proteins are strictly controlled by genetic factors responsible for purifying selection. A comparison between Hordeum, Secale and Triticum species demonstrates that the structure of CENH3 in the subtribes Hordeinae and Triticinae evolved at different rates. The assumption that reticulate evolution served as a factor stabilizing the structure and evolutionary rate of CENH3 and that this factor was more powerful within Secale and Triticum than in Hordeum, is discussed.

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross-species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next-generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross-genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter- and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.

Transcriptional activity of PIF and Pong-like Class II transposable elements in Triticeae

Background

Transposable elements are major contributors to genome size and variability, accounting for approximately 70–80% of the maize, barley, and wheat genomes. PIF and Pong-like elements belong to two closely-related element families within the PIF/Harbinger superfamily of Class II (DNA) transposons. Both elements contain two open reading frames; one encodes a transposase (ORF2) that catalyzes transposition of the functional elements and their related non-autonomous elements, while the function of the second is still debated. In this work, we surveyed for PIF- and Pong-related transcriptional activity in 13 diploid Triticeae species, all of which have been previously shown to harbor extensive within-genome diversity of both groups of elements.

Results

The results revealed that PIF elements have considerable transcriptional activity in Triticeae, suggesting that they can escape the initial levels of plant cell control and are regulated at the post-transcriptional level. Phylogenetic analysis of 156 PIF cDNA transposase fragments along with 240 genomic partial transposase sequences showed that most, if not all, PIF clades are transcriptionally competent, and that multiple transposases coexisting within a single genome have the potential to act simultaneously. In contrast, we did not detect any transcriptional activity of Pong elements in any sample.

Conclusions

The lack of Pong element transcription shows that even closely related transposon families can exhibit wide variation in their transposase transcriptional activity within the same genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1028-6) contains supplementary material, which is available to authorized users.

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.

Contrasting evolutionary patterns of multiple loci uncover new aspects in the genome origin and evolutionary history of Leymus (Triticeae; Poaceae)

Leymus Hochst. (Triticeae: Poaceae), a group of allopolyploid species with the NsXm genomes, is a perennial genus with diversity in morphology, cytology, ecology, and distribution in the Triticeae. To investigate the genome origin and evolutionary history of Leymus, three unlinked low-copy nuclear genes (Acc1, Pgk1, and GBSSI) and three chloroplast regions (trnL-F, matK, and rbcL) of 32 Leymus species were analyzed with those of 36 diploid species representing 18 basic genomes in the Triticeae. The phylogenetic relationships were reconstructed using Bayesian inference, Maximum parsimony, and NeighborNet methods. A time-calibrated phylogeny was generated to estimate the evolutionary history of Leymus. The results suggest that reticulate evolution has occurred in Leymus species, with several distinct progenitors contributing to the Leymus. The molecular data in resolution of the Xm-genome lineage resulted in two apparently contradictory results, with one placing the Xm-genome lineage as closely related to the P/F genome and the other splitting the Xm-genome lineage as sister to the Ns-genome donor. Our results suggested that (1) the Ns genome of Leymus was donated by Psathyrostachys, and additional Ns-containing alleles may be introgressed into some Leymus polyploids by recurrent hybridization; (2) The phylogenetic incongruence regarding the resolution of the Xm-genome lineage suggested that the Xm genome of Leymus was closely related to the P genome of Agropyron; (3) Both Ns- and Xm-genome lineages served as the maternal donor during the speciation of Leymus species; (4) The Pseudoroegneria, Lophopyrum and Australopyrum genomes contributed to some Leymus species.

Molecular Characterization, Gene Evolution, and Expression Analysis of the Fructose-1, 6-bisphosphate Aldolase (FBA) Gene Family in Wheat (Triticum aestivum L.)

Fructose-1, 6-bisphosphate aldolase (FBA) is a key plant enzyme that is involved in glycolysis, gluconeogenesis, and the Calvin cycle. It plays significant roles in biotic and abiotic stress responses, as well as in regulating growth and development processes. In the present paper, 21 genes encoding TaFBA isoenzymes were identified, characterized, and categorized into three groups: class I chloroplast/plastid FBA (CpFBA), class I cytosol FBA (cFBA), and class II chloroplast/plastid FBA. By using a prediction online database and genomic PCR analysis of Chinese Spring nulli-tetrasomic lines, we have confirmed the chromosomal location of these genes in 12 chromosomes of four homologous groups. Sequence and genomic structure analysis revealed the high identity of the allelic TaFBA genes and the origin of different TaFBA genes. Numerous putative environment stimulus-responsive cis-elements have been identified in 1,500-bp regions of TaFBA gene promoters, of which the most abundant are the light-regulated elements (LREs). Phylogenetic reconstruction using the deduced protein sequence of 245 FBA genes indicated an independent evolutionary pathway for the class I and class II groups. Although, earlier studies have indicated that class II FBA only occurs in prokaryote and fungi, our results have demonstrated that a few class II CpFBAs exist in wheat and other closely related species. Class I TaFBA was predicted to be tetramers and class II to be dimers. Gene expression analysis based on microarray and transcriptome databases suggested the distinct role of TaFBAs in different tissues and developmental stages. The TaFBA 4-9 genes were highly expressed in leaves and might play important roles in wheat development. The differential expression patterns of the TaFBA genes in light/dark and a few abiotic stress conditions were also analyzed. The results suggested that LRE cis-elements of TaFBA gene promoters were not directly related to light responses. Most TaFBA genes had higher expression levels in the roots than in the shoots when under various stresses. Class I cytosol TaFBA genes, particularly TaFBA10/12/18 and TaFBA13/16, and three class II TaFBA genes are involved in responses to various abiotic stresses. Class I CpFBA genes in wheat are apparently sensitive to different stress conditions.

Population and phylogenomic decomposition via genotyping-by-sequencing in Australian Pelargonium

Species delimitation has seen a paradigm shift as increasing accessibility of genomic-scale data enables separation of lineages with convergent morphological traits and the merging of recently diverged ecotypes that have distinguishing characteristics. We inferred the process of lineage formation among Australian species in the widespread and highly variable genus Pelargonium by combining phylogenomic and population genomic analyses along with breeding system studies and character analysis. Phylogenomic analysis and population genetic clustering supported seven of the eight currently described species but provided little evidence for differences in genetic structure within the most widely distributed group that containing P. australe. In contrast, morphometric analysis detected three deep lineages within Australian Pelargonium; with P. australe consisting of five previously unrecognized entities occupying separate geographic ranges. The genomic approach enabled elucidation of parallel evolution in some traits formerly used to delineate species, as well as identification of ecotypic morphological differentiation within recognized species. Highly variable morphology and trait convergence each contribute to the discordance between phylogenomic relationships and morphological taxonomy. Data suggest that genetic divergence among species within the Australian Pelargonium may result from allopatric speciation while morphological differentiation within and among species may be more strongly driven by environmental differences.

Characterization of wheat-Secale africanum chromosome 5R(a) derivatives carrying Secale specific genes for grain hardness

New wheat- Secale africanum chromosome 5R (a) substitution and translocation lines were developed and identified by fluorescence in situ hybridization and molecular markers, and chromosome 5R (a) specific genes responsible for grain hardness were isolated. The wild species, Secale africanum Stapf. (genome R(a)R(a)), serves as a valuable germplasm resource for increasing the diversity of cultivated rye (S. cereale L., genome RR) and providing novel genes for wheat improvement. In the current study, fluorescence in situ hybridization (FISH) and molecular markers were applied to characterize new wheat-S. africanum chromosome 5R(a) derivatives. Labeled rye genomic DNA (GISH) and the Oligo-probes pSc119.2 and pTa535 (FISH) were used to study a wheat-S. africanum amphiploid and a disomic 5R(a) (5D) substitution, and to identify a T5DL.5R(a)S translocation line and 5R(a)S and 5R(a)L isotelosome lines. Twenty-one molecular markers were mapped to chromosome 5R(a) arms which will facilitate future rapid identification of 5R(a) introgressions in wheat backgrounds. Comparative analysis of the molecular markers mapped on 5R(a) with homoeologous regions in wheat confirmed a deletion on the chromosome T5DL.5R(a)S, which suggests that the wheat-S. africanum Robertsonian translocation involving homologous group 5 may not be fully compensating. Complete coding sequences at the paralogous puroindoline-a (Pina) and grain softness protein gene (Gsp-1) loci from S. africanum were cloned and localized onto the short arm of chromosome 5R(a). The S. africanum chromosome 5R(a) substitution and translocation lines showed a reduction in the hardness index, which may be associated with the S. africanum- specific Pina and Gsp-1 gene sequences. The present study reports the production of novel wheat-S. africanum chromosome 5R(a) stripe rust resistant derivatives and new rye-specific molecular markers, which may find application in future use of wild Secale genome resources for grain quality studies and disease resistance breeding.

Wheat miRNA ancestors: evident by transcriptome analysis of A, B, and D genome donors

MicroRNAs are critical players of post-transcriptional gene regulation with profound effects on the fundamental processes of cellular life. Their identification and characterization, together with their targets, hold great significance in exploring and exploiting their roles on a functional context, providing valuable clues into the regulation of important biological processes, such as stress tolerance or environmental adaptation. Wheat is a hardy crop, extensively harvested in temperate regions, and is a major component of the human diet. With the advent of the next generation sequencing technologies considerably decreasing sequencing costs per base-pair, genomic, and transcriptomic data from several wheat species, including the progenitors and wild relatives have become available. In this study, we performed in silico identification and comparative analysis of microRNA repertoires of bread wheat (Triticum aestivum L.) and its diploid progenitors and relatives, Aegilops sharonensis, Aegilops speltoides, Aegilops tauschii, Triticum monococcum, and Triticum urartu through the utilization of publicly available transcriptomic data. Over 200 miRNA families were identified, majority of which have not previously been reported. Ancestral relationships expanded our understanding of wheat miRNA evolution, while T. monococcum miRNAs delivered important clues on the effects of domestication on miRNA expression. Comparative analyses on wild Ae. sharonensis accessions highlighted candidate miRNAs that can be linked to stress tolerance. The miRNA repertoires of bread wheat and its diploid progenitors and relatives provide important insight into the diversification and distribution of miRNA genes, which should contribute to the elucidation of miRNA evolution of Poaceae family. A thorough understanding of the convergent and divergent expression profiles of miRNAs in different genetic backgrounds can provide unique opportunities to modulation of gene regulation for better crop performance.

The Role of Vertical and Horizontal Transfer in the Evolutionary Dynamics of PIF-Like Transposable Elements in Triticeae

PIF-like transposable elements are members of the PIF/Harbinger superfamily of DNA transposons found in the genomes of many plants, animals, and fungi. The evolution of the gene that encodes the transposase responsible for mobilizing PIF-like elements has been studied in both plants and animals, but the elements' history in flowering plants remains poorly known. In this work, we describe the phylogenetic distribution and evolution of PIF-like elements in the genomes of 21 diploid species from the wheat tribe, Triticeae, and we present the first convincing evidence of horizontal transfer of PIF elements in plant genomes. A phylogenetic analysis of 240 PIF sequences based on the conserved region of the transposase domain revealed at least four main transposase lineages. Their complex evolutionary history can be best explained by a combination of vertical transmission with differential evolutionary success among lineages, and occasional horizontal transfer between phylogenetically distant Triticeae genera. In addition, we identified 127 potentially functional transposase sequences indicating possible recent activity of PIF.

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.

Diversity, abundance, and evolutionary dynamics of Pong-like transposable elements in Triticeae

Pong-like elements are members of the PIF/Harbinger superfamily of DNA transposons that has been described in many plants, animals, and fungi. Most Pong elements contain two open reading frames (ORFs). One encodes a transposase (ORF2) that catalyzes transposition of Pong and related non-autonomous elements, while the function of the second is unknown. Little is known about the evolutionary history of Pong elements in flowering plants. In this work, we present the first comprehensive analysis of the diversity, abundance, and evolution of the Pong-like transposase gene in the genomes of 21 diploid species from the wheat tribe, Triticeae, and we present the first convincing evidence of horizontal transfer of nuclear-encoded Pong elements in any organism. A phylogenetic analysis of nearly 300 Pong sequences based on a conserved region of the transposase domain revealed a complex evolutionary history of Pong elements that can be best explained by ancestral polymorphism, followed by differential evolutionary success of some transposase lineages, and by occasional horizontal transfer between phylogenetically distant genera. In addition, we used transposon display to estimate the abundance of the transposase gene within Triticeae genomes, and our results revealed varying levels of Pong proliferation, with numbers of transposase copies ranging from 22 to 92. Comparisons of Pong transposase abundance to flow cytometry estimates of genome size revealed that larger Triticeae genome size was not correlated with transposase abundance.

Plastid phylogenomics of the cool-season grass subfamily: clarification of relationships among early-diverging tribes

Whole plastid genomes are being sequenced rapidly from across the green plant tree of life, and phylogenetic analyses of these are increasing resolution and support for relationships that have varied among or been unresolved in earlier single- and multi-gene studies. Pooideae, the cool-season grass lineage, is the largest of the 12 grass subfamilies and includes important temperate cereals, turf grasses and forage species. Although numerous studies of the phylogeny of the subfamily have been undertaken, relationships among some 'early-diverging' tribes conflict among studies, and some relationships among subtribes of Poeae have not yet been resolved. To address these issues, we newly sequenced 25 whole plastomes, which showed rearrangements typical of Poaceae. These plastomes represent 9 tribes and 11 subtribes of Pooideae, and were analysed with 20 existing plastomes for the subfamily. Maximum likelihood (ML), maximum parsimony (MP) and Bayesian inference (BI) robustly resolve most deep relationships in the subfamily. Complete plastome data provide increased nodal support compared with protein-coding data alone at nodes that are not maximally supported. Following the divergence of Brachyelytrum, Phaenospermateae, Brylkinieae-Meliceae and Ampelodesmeae-Stipeae are the successive sister groups of the rest of the subfamily. Ampelodesmeae are nested within Stipeae in the plastome trees, consistent with its hybrid origin between a phaenospermatoid and a stipoid grass (the maternal parent). The core Pooideae are strongly supported and include Brachypodieae, a Bromeae-Triticeae clade and Poeae. Within Poeae, a novel sister group relationship between Phalaridinae and Torreyochloinae is found, and the relative branching order of this clade and Aveninae, with respect to an Agrostidinae-Brizinae clade, are discordant between MP and ML/BI trees. Maximum likelihood and Bayesian analyses strongly support Airinae and Holcinae as the successive sister groups of a Dactylidinae-Loliinae clade.

From gene trees to a dated allopolyploid network: insights from the angiosperm genus Viola (Violaceae)

Allopolyploidization accounts for a significant fraction of speciation events in many eukaryotic lineages. However, existing phylogenetic and dating methods require tree-like topologies and are unable to handle the network-like phylogenetic relationships of lineages containing allopolyploids. No explicit framework has so far been established for evaluating competing network topologies, and few attempts have been made to date phylogenetic networks. We used a four-step approach to generate a dated polyploid species network for the cosmopolitan angiosperm genus Viola L. (Violaceae Batch.). The genus contains ca 600 species and both recent (neo-) and more ancient (meso-) polyploid lineages distributed over 16 sections. First, we obtained DNA sequences of three low-copy nuclear genes and one chloroplast region, from 42 species representing all 16 sections. Second, we obtained fossil-calibrated chronograms for each nuclear gene marker. Third, we determined the most parsimonious multilabeled genome tree and its corresponding network, resolved at the section (not the species) level. Reconstructing the "correct" network for a set of polyploids depends on recovering all homoeologs, i.e., all subgenomes, in these polyploids. Assuming the presence of Viola subgenome lineages that were not detected by the nuclear gene phylogenies ("ghost subgenome lineages") significantly reduced the number of inferred polyploidization events. We identified the most parsimonious network topology from a set of five competing scenarios differing in the interpretation of homoeolog extinctions and lineage sorting, based on (i) fewest possible ghost subgenome lineages, (ii) fewest possible polyploidization events, and (iii) least possible deviation from expected ploidy as inferred from available chromosome counts of the involved polyploid taxa. Finally, we estimated the homoploid and polyploid speciation times of the most parsimonious network. Homoploid speciation times were estimated by coalescent analysis of gene tree node ages. Polyploid speciation times were estimated by comparing branch lengths and speciation rates of lineages with and without ploidy shifts. Our analyses recognize Viola as an old genus (crown age 31 Ma) whose evolutionary history has been profoundly affected by allopolyploidy. Between 16 and 21 allopolyploidizations are necessary to explain the diversification of the 16 major lineages (sections) of Viola, suggesting that allopolyploidy has accounted for a high percentage-between 67% and 88%-of the speciation events at this level. The theoretical and methodological approaches presented here for (i) constructing networks and (ii) dating speciation events within a network, have general applicability for phylogenetic studies of groups where allopolyploidization has occurred. They make explicit use of a hitherto underexplored source of ploidy information from chromosome counts to help resolve phylogenetic cases where incomplete sequence data hampers network inference. Importantly, the coalescent-based method used herein circumvents the assumption of tree-like evolution required by most techniques for dating speciation events.

A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome

An ordered draft sequence of the 17-gigabase hexaploid bread wheat (Triticum aestivum) genome has been produced by sequencing isolated chromosome arms. We have annotated 124,201 gene loci distributed nearly evenly across the homeologous chromosomes and subgenomes. Comparative gene analysis of wheat subgenomes and extant diploid and tetraploid wheat relatives showed that high sequence similarity and structural conservation are retained, with limited gene loss, after polyploidization. However, across the genomes there was evidence of dynamic gene gain, loss, and duplication since the divergence of the wheat lineages. A high degree of transcriptional autonomy and no global dominance was found for the subgenomes. These insights into the genome biology of a polyploid crop provide a springboard for faster gene isolation, rapid genetic marker development, and precise breeding to meet the needs of increasing food demand worldwide.

Ancient hybridizations among the ancestral genomes of bread wheat

The allohexaploid bread wheat genome consists of three closely related subgenomes (A, B, and D), but a clear understanding of their phylogenetic history has been lacking. We used genome assemblies of bread wheat and five diploid relatives to analyze genome-wide samples of gene trees, as well as to estimate evolutionary relatedness and divergence times. We show that the A and B genomes diverged from a common ancestor ~7 million years ago and that these genomes gave rise to the D genome through homoploid hybrid speciation 1 to 2 million years later. Our findings imply that the present-day bread wheat genome is a product of multiple rounds of hybrid speciation (homoploid and polyploid) and lay the foundation for a new framework for understanding the wheat genome as a multilevel phylogenetic mosaic.

Comparison of Acetyl-CoA carboxylase 1 (Acc-1) gene diversity among different Triticeae genomes

It has widely been documented that life form and mating system have significant influences on genetic diversity. In the tribe Triticeae, several genera contain both annual and perennial species, whereas other genera comprise strictly annual or perennial species. It was suggested that Triticeae annuals have originated from Triticeae perennials. The present study aims to analyze nucleotide diversity of Acc-1 gene among different Triticeae genomes, and attempts to link effects of life history (annuals and perennials) and mating systems. The nucleotide diversity of 364 Acc-1 sequences in Triticeae species was characterized. The highest estimates of nucleotide diversity values (π=0.01919, θ=0.03515) were found for the Ns genome among the genomes analyzed. Nucleotide diversities in the D genome and Ns genome of polyploids are higher than those in respective genomes of diploids, while in the St genome of polyploids, it is lower than that in the St genome of diploids. The averaged π value (0.013705) in the genomes of perennials is more than twice of the value (0.00508) in the genomes of annuals. The averaged π value (0.01323) in the genomes of outcrossing species is two-fold of the value (0.005664) in the genomes of selfer. Our results suggested that the evolutionary history and mating system may play an important role in determining nucleotide diversity of Acc-1 gene in each genome.

Limited genetic divergence among Australian alpine Poa tussock grasses coupled with regional structuring points to ongoing gene flow and taxonomic challenges

BACKGROUND AND AIMS

While molecular approaches can often accurately reconstruct species relationships, taxa that are incompletely differentiated pose a challenge even with extensive data. Such taxa are functionally differentiated, but may be genetically differentiated only at small and/or patchy regions of the genome. This issue is considered here in Poa tussock grass species that dominate grassland and herbfields in the Australian alpine zone.

METHODS

Previously reported tetraploidy was confirmed in all species by sequencing seven nuclear regions and five microsatellite markers. A Bayesian approach was used to co-estimate nuclear and chloroplast gene trees with an overall dated species tree. The resulting species tree was used to examine species structure and recent hybridization, and intertaxon fertility was tested by experimental crosses.

KEY RESULTS

Species tree estimation revealed Poa gunnii, a Tasmanian endemic species, as sister to the rest of the Australian alpine Poa. The taxa have radiated in the last 0·5-1·2 million years and the non-gunnii taxa are not supported as genetically distinct. Recent hybridization following past species divergence was also not supported. Ongoing gene flow is suggested, with some broad-scale geographic structure within the group.

CONCLUSIONS

The Australian alpine Poa species are not genetically distinct despite being distinguishable phenotypically, suggesting recent adaptive divergence with ongoing intertaxon gene flow. This highlights challenges in using conventional molecular taxonomy to infer species relationships in recent, rapid radiations.

Evolutionary pattern of rDNA following polyploidy in Leymus (Triticeae: Poaceae)

Ribosomal ITS polymorphism and its ancestral genome origin of polyploid Leymus were examined to infer the evolutionary outcome of rDNA gene following allopolyploid speciation and to elucidate the geographic pattern of ITS variation. The results demonstrated that different polyploids have experienced varying fates, including maintenance or homogenization of divergent arrays, occurrence of chimeric repeats and potential pseudogenes. Our data suggested that (1) the Ns, P/F, and St genomic types in Leymus were originated from Psathyrostachys, Agropyron/Eremopyrum, and Pseudoroegneria, respectively; (2) the occurrence of a higher proportion of Leymus species with predominant uniparental rDNA type might associate with the segmental allopolyploid origin, nucleolar dominance of alloploids, and rapid radiation of Leymus; (3) maintenance of multiple parental ITS types in allopolyploid might result from long generation times associated to vegetative multiplication, number and chromosomal location of ribosomal loci and/or recurrent hybridization; (4) the rDNA genealogical structure of Leymus species might associate with the geographic origins; and (5) ITS sequence clade shared by Leymus species from Central Asia, North America, and Nordic might be an outcome of ancestral ITS homogenization. Our results shed new light on understanding evolutionary outcomes of rDNA following allopolyploid speciation and geographic isolation.

Sequencing of chloroplast genomes from wheat, barley, rye and their relatives provides a detailed insight into the evolution of the Triticeae tribe

Using Roche/454 technology, we sequenced the chloroplast genomes of 12 Triticeae species, including bread wheat, barley and rye, as well as the diploid progenitors and relatives of bread wheat Triticum urartu, Aegilops speltoides and Ae. tauschii. Two wild tetraploid taxa, Ae. cylindrica and Ae. geniculata, were also included. Additionally, we incorporated wild Einkorn wheat Triticum boeoticum and its domesticated form T. monococcum and two Hordeum spontaneum (wild barley) genotypes. Chloroplast genomes were used for overall sequence comparison, phylogenetic analysis and dating of divergence times. We estimate that barley diverged from rye and wheat approximately 8–9 million years ago (MYA). The genome donors of hexaploid wheat diverged between 2.1–2.9 MYA, while rye diverged from Triticum aestivum approximately 3–4 MYA, more recently than previously estimated. Interestingly, the A genome taxa T. boeoticum and T. urartu were estimated to have diverged approximately 570,000 years ago. As these two have a reproductive barrier, the divergence time estimate also provides an upper limit for the time required for the formation of a species boundary between the two. Furthermore, we conclusively show that the chloroplast genome of hexaploid wheat was contributed by the B genome donor and that this unknown species diverged from Ae. speltoides about 980,000 years ago. Additionally, sequence alignments identified a translocation of a chloroplast segment to the nuclear genome which is specific to the rye/wheat lineage. We propose the presented phylogeny and divergence time estimates as a reference framework for future studies on Triticeae.

Development of a wheat single gene FISH map for analyzing homoeologous relationship and chromosomal rearrangements within the Triticeae

A cytogenetic map of wheat was constructed using FISH with cDNA probes. FISH markers detected homoeology and chromosomal rearrangements of wild relatives, an important source of genes for wheat improvement. To transfer agronomically important genes from wild relatives to bread wheat (Triticum aestivum L., 2n = 6 x = 42, AABBDD) by induced homoeologous recombination, it is important to know the chromosomal relationships of the species involved. Fluorescence in situ hybridization (FISH) can be used to study chromosome structure. The genomes of allohexaploid bread wheat and other species from the Triticeae tribe are colinear to some extent, i.e., composed of homoeoloci at similar positions along the chromosomes, and with genic regions being highly conserved. To develop cytogenetic markers specific for genic regions of wheat homoeologs, we selected more than 60 full-length wheat cDNAs using BLAST against mapped expressed sequence tags and used them as FISH probes. Most probes produced signals on all three homoeologous chromosomes at the expected positions. We developed a wheat physical map with several cDNA markers located on each of the 14 homoeologous chromosome arms. The FISH markers confirmed chromosome rearrangements within wheat genomes and were successfully used to study chromosome structure and homoeology in wild Triticeae species. FISH analysis detected 1 U-6 U chromosome translocation in the genome of Aegilops umbellulata, showed colinearity between chromosome A of Ae. caudata and group-1 wheat chromosomes, and between chromosome arm 7S#3 L of Thinopyrum intermedium and the long arm of the group-7 wheat chromosomes.

Development of a wheat single gene FISH map for analyzing homoeologous relationship and chromosomal rearrangements within the Triticeae

A cytogenetic map of wheat was constructed using FISH with cDNA probes. FISH markers detected homoeology and chromosomal rearrangements of wild relatives, an important source of genes for wheat improvement. To transfer agronomically important genes from wild relatives to bread wheat (Triticum aestivum L., 2n = 6 x = 42, AABBDD) by induced homoeologous recombination, it is important to know the chromosomal relationships of the species involved. Fluorescence in situ hybridization (FISH) can be used to study chromosome structure. The genomes of allohexaploid bread wheat and other species from the Triticeae tribe are colinear to some extent, i.e., composed of homoeoloci at similar positions along the chromosomes, and with genic regions being highly conserved. To develop cytogenetic markers specific for genic regions of wheat homoeologs, we selected more than 60 full-length wheat cDNAs using BLAST against mapped expressed sequence tags and used them as FISH probes. Most probes produced signals on all three homoeologous chromosomes at the expected positions. We developed a wheat physical map with several cDNA markers located on each of the 14 homoeologous chromosome arms. The FISH markers confirmed chromosome rearrangements within wheat genomes and were successfully used to study chromosome structure and homoeology in wild Triticeae species. FISH analysis detected 1 U-6 U chromosome translocation in the genome of Aegilops umbellulata, showed colinearity between chromosome A of Ae. caudata and group-1 wheat chromosomes, and between chromosome arm 7S#3 L of Thinopyrum intermedium and the long arm of the group-7 wheat chromosomes.

A review of the occurrence of Grain softness protein-1 genes in wheat (Triticum aestivum L.)

Grain softness protein-1 (Gsp-1) is a small, 495-bp intronless gene found throughout the Triticeae tribe at the distal end of group 5 chromosomes. With the Puroindolines, it constitutes a key component of the Hardness locus. Gsp-1 likely plays little role in grain hardness, but has direct interest due to its utility in phylogeny and its role in arabinogalactan peptides. Further role(s) remain to be identified. In the polyploid wheats, Triticum aestivum and T. turgidum, the gene is present in a homoeologous series. Since its discovery, there have been conflicting reports and data as to the number of Gsp-1 genes and the level of sequence polymorphism. Little is known about allelic variation within a species. In the simplest model, a single Gsp-1 gene is present in each wheat and Aegilops tauschii genome. The present review critically re-examines the published and some unpublished data (sequence available in the NCBI nucleotide and MIPS Wheat Genome Databases). A number of testable hypotheses are identified, and include the level of polymorphism that may represent (and define) different Gsp-1 alleles, the existence of a fourth Gsp-1 gene, and the apparent, at times, high level of naturally-occurring or artifactual gene chimeras. In summary, the present data provide firm evidence for at most, three Gsp-1 genes in wheat, although there are numerous data that suggest a more complex model.

Sina and Sinb genes in triticale do not determine grain hardness contrary to their orthologs Pina and Pinb in wheat

BACKGROUND

Secaloindoline a (Sina) and secaloindoline b (Sinb) genes of hexaploid triticale (x Triticosecale Wittmack) are orthologs of puroindoline a (Pina) and puroindoline b (Pinb) in hexaploid wheat (Triticum aestivum L.). It has already been proven that RNA interference (RNAi)-based silencing of Pina and Pinb genes significantly decreased the puroindoline a and puroindoline b proteins in wheat and essentially increased grain hardness (J Exp Bot 62:4025-4036, 2011). The function of Sina and Sinb in triticale was tested by means of RNAi silencing and compared to wheat.

RESULTS

Novel Sina and Sinb alleles in wild-type plants of cv. Wanad were identified and their expression profiles characterized. Alignment with wheat Pina-D1a and Pinb-D1a alleles showed 95% and 93.3% homology with Sina and Sinb coding sequences. Twenty transgenic lines transformed with two hpRNA silencing cassettes directed to silence Sina or Sinb were obtained by the Agrobacterium-mediated method. A significant decrease of expression of both Sin genes in segregating progeny of tested T1 lines was observed independent of the silencing cassette used. The silencing was transmitted to the T4 kernel generation. The relative transcript level was reduced by up to 99% in T3 progeny with the mean for the sublines being around 90%. Silencing of the Sin genes resulted in a substantial decrease of secaloindoline a and secaloindoline b content. The identity of SIN peptides was confirmed by mass spectrometry. The hardness index, measured by the SKCS (Single Kernel Characterization System) method, ranged from 22 to 56 in silent lines and from 37 to 49 in the control, and the mean values were insignificantly lower in the silent ones, proving increased softness. Additionally, the mean total seed protein content of silenced lines was about 6% lower compared with control lines. Correlation coefficients between hardness and transcript level were weakly positive.

CONCLUSIONS

We documented that RNAi-based silencing of Sin genes resulted in significant decrease of their transcripts and the level of both secaloindoline proteins, however did not affect grain hardness. The unexpected, functional differences of Sin genes from triticale compared with their orthologs, Pin of wheat, are discussed.

Reticulate evolution of the rye genome

Rye (Secale cereale) is closely related to wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to its large genome (~8 Gb) and its regional importance, genome analysis of rye has lagged behind other cereals. Here, we established a virtual linear gene order model (genome zipper) comprising 22,426 or 72% of the detected set of 31,008 rye genes. This was achieved by high-throughput transcript mapping, chromosome survey sequencing, and integration of conserved synteny information of three sequenced model grass genomes (Brachypodium distachyon, rice [Oryza sativa], and sorghum [Sorghum bicolor]). This enabled a genome-wide high-density comparative analysis of rye/barley/model grass genome synteny. Seventeen conserved syntenic linkage blocks making up the rye and barley genomes were defined in comparison to model grass genomes. Six major translocations shaped the modern rye genome in comparison to a putative Triticeae ancestral genome. Strikingly dissimilar conserved syntenic gene content, gene sequence diversity signatures, and phylogenetic networks were found for individual rye syntenic blocks. This indicates that introgressive hybridizations (diploid or polyploidy hybrid speciation) and/or a series of whole-genome or chromosome duplications played a role in rye speciation and genome evolution.

Developing nuclear DNA phylogenetic markers in the angiosperm genus Leucadendron (Proteaceae): a next-generation sequencing transcriptomic approach

Despite the recent advances in generating molecular data, reconstructing species-level phylogenies for non-models groups remains a challenge. The use of a number of independent genes is required to resolve phylogenetic relationships, especially for groups displaying low polymorphism. In such cases, low-copy nuclear exons and non-coding regions, such as 3' untranslated regions (3'-UTRs) or introns, constitute a potentially interesting source of nuclear DNA variation. Here, we present a methodology meant to identify new nuclear orthologous markers using both public-nucleotide databases and transcriptomic data generated for the group of interest by using next generation sequencing technology. To identify PCR primers for a non-model group, the genus Leucadendron (Proteaceae), we adopted a framework aimed at minimizing the probability of paralogy and maximizing polymorphism. We anchored when possible the right-hand primer into the 3'-UTR and the left-hand primer into the coding region. Seven new nuclear markers emerged from this search strategy, three of those included 3'-UTRs. We further compared the phylogenetic potential between our new markers and the ribosomal internal transcribed spacer region (ITS). The sequenced 3'-UTRs yielded higher polymorphism rates than the ITS region did. We did not find strong incongruences with the phylogenetic signal contained in the ITS region and the seven new designed markers but they strongly improved the phylogeny of the genus Leucadendron. Overall, this methodology is efficient in isolating orthologous loci and is valid for any non-model group given the availability of transcriptomic data.

Less is more in mammalian phylogenomics: AT-rich genes minimize tree conflicts and unravel the root of placental mammals

Despite the rapid increase of size in phylogenomic data sets, a number of important nodes on animal phylogeny are still unresolved. Among these, the rooting of the placental mammal tree is still a controversial issue. One difficulty lies in the pervasive phylogenetic conflicts among genes, with each one telling its own story, which may be reliable or not. Here, we identified a simple criterion, that is, the GC content, which substantially helps in determining which gene trees best reflect the species tree. We assessed the ability of 13,111 coding sequence alignments to correctly reconstruct the placental phylogeny. We found that GC-rich genes induced a higher amount of conflict among gene trees and performed worse than AT-rich genes in retrieving well-supported, consensual nodes on the placental tree. We interpret this GC effect mainly as a consequence of genome-wide variations in recombination rate. Indeed, recombination is known to drive GC-content evolution through GC-biased gene conversion and might be problematic for phylogenetic reconstruction, for instance, in an incomplete lineage sorting context. When we focused on the AT-richest fraction of the data set, the resolution level of the placental phylogeny was greatly increased, and a strong support was obtained in favor of an Afrotheria rooting, that is, Afrotheria as the sister group of all other placentals. We show that in mammals most conflicts among gene trees, which have so far hampered the resolution of the placental tree, are concentrated in the GC-rich regions of the genome. We argue that the GC content-because it is a reliable indicator of the long-term recombination rate-is an informative criterion that could help in identifying the most reliable molecular markers for species tree inference.

Expansion of the gamma-gliadin gene family in Aegilops and Triticum

BACKGROUND

The gamma-gliadins are considered to be the oldest of the gliadin family of storage proteins in Aegilops/Triticum. However, the expansion of this multigene family has not been studied in an evolutionary perspective.

RESULTS

We have cloned 59 gamma-gliadin genes from Aegilops and Triticum species (Aegilops caudata L., Aegilops comosa Sm. in Sibth. & Sm., Aegilops mutica Boiss., Aegilops speltoides Tausch, Aegilops tauschii Coss., Aegilops umbellulata Zhuk., Aegilops uniaristata Vis., and Triticum monococcum L.) representing eight different genomes: Am, B/S, C, D, M, N, T and U. Overall, 15% of the sequences contained internal stop codons resulting in pseudogenes, but this percentage was variable among genomes, up to over 50% in Ae. umbellulata. The most common length of the deduced protein, including the signal peptide, was 302 amino acids, but the length varied from 215 to 362 amino acids, both obtained from Ae. speltoides. Most genes encoded proteins with eight cysteines. However, all Aegilops species had genes that encoded a gamma-gliadin protein of 302 amino acids with an additional cysteine. These conserved nine-cysteine gamma-gliadins may perform a specific function, possibly as chain terminators in gluten network formation in protein bodies during endosperm development. A phylogenetic analysis of gamma-gliadins derived from Aegilops and Triticum species and the related genera Lophopyrum, Crithopsis, and Dasypyrum showed six groups of genes. Most Aegilops species contained gamma-gliadin genes from several of these groups, which also included sequences from the genera Lophopyrum, Crithopsis, and Dasypyrum. Hordein and secalin sequences formed separate groups.

CONCLUSIONS

We present a model for the evolution of the gamma-gliadins from which we deduce that the most recent common ancestor (MRCA) of Aegilops/Triticum-Dasypyrum-Lophopyrum-Crithopsis already had four groups of gamma-gliadin sequences, presumably the result of two rounds of duplication of the locus.

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.

Generic phylogeny, historical biogeography and character evolution of the cosmopolitan aquatic plant family Hydrocharitaceae

Background

Hydrocharitaceae is a fully aquatic monocot family, consists of 18 genera with approximately 120 species. The family includes both fresh and marine aquatics and exhibits great diversity in form and habit including annual and perennial life histories; submersed, partially submersed and floating leaf habits and linear to orbicular leaf shapes. The family has a cosmopolitan distribution and is well represented in the Tertiary fossil record in Europe. At present, the historical biogeography of the family is not well understood and the generic relationships remain controversial. In this study we investigated the phylogeny and biogeography of Hydrocharitaceae by integrating fossils and DNA sequences from eight genes. We also conducted ancestral state reconstruction for three morphological characters.

Results

Phylogenetic analyses produced a phylogeny with most branches strongly supported by bootstrap values greater than 95 and Bayesian posterior probability values of 1.0. Stratiotes is the first diverging lineage with the remaining genera in two clades, one clade consists of Lagarosiphon, Ottelia, Blyxa, Apalanthe, Elodea and Egeria; and the other consists of Hydrocharis-Limnobium, Thalassia, Enhalus, Halophila, Najas, Hydrilla, Vallisneria, Nechamandra and Maidenia. Biogeographic analyses (DIVA, Mesquite) and divergence time estimates (BEAST) resolved the most recent common ancestor of Hydrocharitaceae as being in Asia during the Late Cretaceous and Palaeocene (54.7-72.6 Ma). Dispersals (including long-distance dispersal and migrations through Tethys seaway and land bridges) probably played major roles in the intercontinental distribution of this family. Ancestral state reconstruction suggested that in Hydrocharitaceae evolution of dioecy is bidirectional, viz., from dioecy to hermaphroditism, and from hermaphroditism to dioecy, and that the aerial-submerged leaf habit and short-linear leaf shape are the ancestral states.

Conclusions

Our study has shed light on the previously controversial generic phylogeny of Hydrocharitaceae. The study has resolved the historical biogeography of this family and supported dispersal as the most likely explanation for the intercontinental distribution. We have also provided valuable information for understanding the evolution of breeding system and leaf phenotype in aquatic monocots.