The {beta}-amylase genes of grasses and a phylogenetic analysis of the Triticeae (Poaceae)

Comprehensive Identification of the β-Amylase (BAM) Gene Family in Response to Cold Stress in White Clover

White clover (Trifolium repens L.) is an allopolyploid plant and an excellent perennial legume forage. However, white clover is subjected to various stresses during its growth, with cold stress being one of the major limiting factors affecting its growth and development. Beta-amylase (BAM) is an important starch-hydrolyzing enzyme that plays a significant role in starch degradation and responses to environmental stress. In this study, 21 members of the BAM gene family were identified in the white clover genome. A phylogenetic analysis using BAMs from Arabidopsis divided TrBAMs into four groups based on sequence similarity. Through analysis of conserved motifs, gene duplication, synteny analysis, and cis-acting elements, a deeper understanding of the structure and evolution of TrBAMs in white clover was gained. Additionally, a gene regulatory network (GRN) containing TrBAMs was constructed; gene ontology (GO) annotation analysis revealed close interactions between TrBAMs and AMY (α-amylase) and DPE (4-alpha-glucanotransferase). To determine the function of TrBAMs under various tissues and stresses, RNA-seq datasets were analyzed, showing that most TrBAMs were significantly upregulated in response to biotic and abiotic stresses and the highest expression in leaves. These results were validated through qRT-PCR experiments, indicating their involvement in multiple gene regulatory pathways responding to cold stress. This study provides new insights into the structure, evolution, and function of the white clover BAM gene family, laying the foundation for further exploration of the functional mechanisms through which TrBAMs respond to cold stress.

Genome-Wide Identification and Preliminary Functional Analysis of BAM (β-Amylase) Gene Family in Upland Cotton

The β-amylase (BAM) gene family encodes important enzymes that catalyze the conversion of starch to maltose in various biological processes of plants and play essential roles in regulating the growth and development of multiple plants. So far, BAMs have been extensively studied in Arabidopsis thaliana (A. thaliana). However, the characteristics of the BAM gene family in the crucial economic crop, cotton, have not been reported. In this study, 27 GhBAM genes in the genome of Gossypium hirsutum L (G. hirsutum) were identified by genome-wide identification, and they were divided into three groups according to sequence similarity and phylogenetic relationship. The gene structure, chromosome distribution, and collinearity of all GhBAM genes identified in the genome of G. hirsutum were analyzed. Further sequence alignment of the core domain of glucosyl hydrolase showed that all GhBAM family genes had the glycosyl hydrolase family 14 domain. We identified the BAM gene GhBAM7 and preliminarily investigated its function by transcriptional sequencing analysis, qRT-PCR, and subcellular localization. These results suggested that the GhBAM7 gene may influence fiber strength during fiber development. This systematic analysis provides new insight into the transcriptional characteristics of BAM genes in G. hirsutum. It may lay the foundation for further study of the function of these genes.

Chromosome-specific painting in Thinopyrum species using bulked oligonucleotides

Chromosome-specific painting probes were developed to identify the individual chromosomes from 1 to 7E in Thinopyrum species and detect alien genetic material of the E genome in a wheat background. The E genome of Thinopyrum is closely related to the ABD genome of wheat (Triticum aestivum L.) and harbors genes conferring beneficial traits to wheat, including high yield, disease resistance, and unique end-use quality. Species of Thinopyrum vary from diploid (2n = 2x = 14) to decaploid (2n = 10x = 70), and chromosome structural variation and differentiation have arisen during polyploidization. To investigate the variation and evolution of the E genome, we developed a complete set of E genome-specific painting probes for identification of the individual chromosomes 1E to 7E based on the genome sequences of Th. elongatum (Host) D. R. Dewey and wheat. By using these new probes in oligonucleotide-based chromosome painting, we showed that Th. bessarabicum (PI 531711, EbEb) has a close genetic relationship with diploid Th. elongatum (EeEe), with five chromosomes (1E, 2E, 3E, 6E, and 7E) maintaining complete synteny in the two species except for a reciprocal translocation between 4 and 5Eb. All 14 pairs of chromosomes of tetraploid Th. elongatum have maintained complete synteny with those of diploid Th. elongatum (Thy14), but the two sets of E genomes have diverged. This study also demonstrated that the E genome-specific painting probes are useful for rapid and effective detection of the alien genetic material of E genome in wheat-Thinopyrum derived lines.

A chromosome-scale genome assembly of Dasypyrum villosum provides insights into its application as a broad-spectrum disease resistance resource for wheat improvement

Dasypyrum villosum is one of the most valuable gene resources in wheat improvement, especially for disease resistance. The mining of favorable genes from D. villosum is frustrated by the lack of a whole genome sequence. In this study, we generated a doubled-haploid line, 91C43^DH, using microspore culture and obtained a 4.05-GB high-quality, chromosome-scale genome assembly for D. villosum. The assembly contains39 727 high-confidence genes, and 85.31% of the sequences are repetitive. Two reciprocal translocation events were detected, and 7VS-4VL is a unique translocation in D. villosum. The prolamin seed storage protein-coding genes were found to be duplicated; in particular, the genes encoding low-molecular-weight glutenin at the Glu-V3 locus were significantly expanded. RNA sequencing (RNA-seq) analysis indicated that, after Blumeria graminearum f.sp tritici (Bgt) inoculation, there were more upregulated genes involved in the pattern-triggered immunity and effector-triggered immunity defense pathways in D. villosum than in Triticum urartu. MNase hypersensitive sequencing (MH-seq) identified two Bgt-inducible MH sites (MHSs), one in the promoter and one in the 3' terminal region of the powdery mildew resistance (Pm) gene NLR1-V. Each site had two subpeaks and they were termed MHS1 (MHS1.1/1.2) and MHS2 (MHS2.1/2.2). Bgt-inducible MHS2.2 was uniquely present in D. villosum, and MHS1.1 was more inducible in D. villosum than in wheat, suggesting that MHSs may be critical for regulation of NLR1-V expression and plant defense. In summary, this study provides a valuable genome resource for functional genomics studies and wheat-D. villosum introgression breeding. The identified regulatory mechanisms may also be exploited to develop new strategies for enhancing Pm resistance by optimizing gene expression in wheat.

Comparative transcriptome analysis of resistant and susceptible Kentucky bluegrass varieties in response to powdery mildew infection

BACKGROUND

Poa pratensis is one of the most common cold-season turfgrasses used for urban turf building, and it is also widely used in ecological environment management worldwide. Powdery mildew is a common disease of P. pratensis. To scientifically and ecologically control lawn powdery mildew, the molecular mechanism underlying the response of P. pratensis to powdery mildew infection must better understood.

RESULTS

To explore molecular mechanism underlying the response of P. pratensis to powdery mildew infection, this study compared physiological changes and transcriptomic level differences between the highly resistant variety 'BlackJack' and the extremely susceptible variety 'EverGlade' under powdery mildew infection conditions. We analyzed DEGs using reference canonical pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and the results showed that "starch and sucrose metabolism", "photosynthesis" and "fatty acid metabolism"pathways were only enriched in 'BlackJack', and the expression of DEGs such as HXK, INV, GS, SS, AGpase and β-amylase in "starch and sucrose metabolism" pathway of 'BlackJack' were closely related to powdery mildew resistance. Meanwhile, compared with 'EverGlade', powdery mildew infection promoted synthesis of sucrose, expression of photosynthesis parameters and photosynthesis-related enzymes in leaves of 'BlackJack' and decreased accumulation of monosaccharides such as glucose and fructose.

CONCLUSIONS

This study identified the key metabolic pathways of a P. pratensis variety with high resistance to powdery mildew infection and explored the differences in physiological characteristics and key genes related to sugar metabolism pathways under powdery mildew stress. These findings provide important insights for studying underlying molecular response mechanism.

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 Identification and Expression Analysis of the β-Amylase Gene Family in Chenopodium quinoa

β-Amylase (BAM) is an important starch hydrolase, playing a role in a variety of plant growth and development processes. In this study, 22 BAM gene family members (GFMs) were identified in quinoa (Chenopodium quinoa), an ancient crop gaining modern consumer acceptance because of its nutritional qualities. The genetic structure, phylogenetic and evolutionary relationships, and expression patterns of CqBAM GFMs in different tissues, were analyzed. Phylogenetic analyses assigned the CqBAMs, AtBAMs, and OsBAMs into four clades. The CqBAM gene family had expanded due to segmental duplication. RNA-seq analysis revealed expression of the duplicated pairs to be similar, with the expression of CqBAM GFM pairs showing a degree of tissue specificity that was confirmed by reverse transcription quantitative PCR (RT-qPCR). Several CqBAM GFMs were also responsive to abiotic stresses in shoots and/or roots. In conclusion, the BAM gene family in quinoa was identified and systematically analyzed using bioinformatics and experimental methods. These results will help to elucidate the evolutionary relationship and biological functions of the BAM gene family in quinoa.

Origins and chromosome differentiation of Thinopyrum elongatum revealed by PepC and Pgk1 genes and ND-FISH

Thinopyrum elongatum is an important gene pool for wheat genetic improvement. However, the origins of the Thinopyrum genomes and the nature of the genus' intraspecific relationships are still controversial. In this study, we used single-copy nuclear genes and non-denaturing fluorescence in situ hybridization (ND-FISH) to characterize genome constitution and chromosome differentiation in Th. elongatum. According to phylogenetic analyses based on PepC and Pgk1 genes, there was an E genome with three versions (E^e, E^b, E^x) and St genomes in the polyploid Th. elongatum. The ND-FISH results of pSc119.2 and pAs1 revealed that the karyotypes of diploid Th. elongatum and Th. bessarabicum were different, and the chromosome differentiation occurred among accessions of the diploid Th. elongatum. In addition, the tetraploid Th. elongatum has two groups of ND-FISH karyotype, indicating that the tetraploid Th. elongatum might be a segmental allotetraploid. In summary, our results suggested that the diploid Th. elongatum, Th. Bessarabicum, and Pseudoroegneria were the donors of the E^e, E^b, and St genomes to the polyploid Th. elongatum, respectively.

Resequencing of 145 Landmark Cultivars Reveals Asymmetric Sub-genome Selection and Strong Founder Genotype Effects on Wheat Breeding in China

Controlled pedigrees and the multi-decade timescale of national crop plant breeding programs offer a unique experimental context for examining how selection affects plant genomes. More than 3000 wheat cultivars have been registered, released, and documented since 1949 in China. In this study, a set of 145 elite cultivars selected from historical points of wheat breeding in China were re-sequenced. A total of 43.75 Tb of sequence data were generated with an average read depth of 17.94× for each cultivar, and more than 60.92 million SNPs and 2.54 million InDels were captured, based on the Chinese Spring RefSeq genome v1.0. Seventy years of breeder-driven selection led to dramatic changes in grain yield and related phenotypes, with distinct genomic regions and phenotypes targeted by different breeders across the decades. There are very clear instances illustrating how introduced Italian and other foreign germplasm was integrated into Chinese wheat programs and reshaped the genomic landscape of local modern cultivars. Importantly, the resequencing data also highlighted significant asymmetric breeding selection among the three sub-genomes: this was evident in both the collinear blocks for homeologous chromosomes and among sets of three homeologous genes. Accumulation of more newly assembled genes in newer cultivars implied the potential value of these genes in breeding. Conserved and extended sharing of linkage disequilibrium (LD) blocks was highlighted among pedigree-related cultivars, in which fewer haplotype differences were detected. Fixation or replacement of haplotypes from founder genotypes after generations of breeding was related to their breeding value. Based on the haplotype frequency changes in LD blocks of pedigree-related cultivars, we propose a strategy for evaluating the breeding value of any given line on the basis of the accumulation (pyramiding) of beneficial haplotypes. Collectively, our study demonstrates the influence of "founder genotypes" on the output of breeding efforts over many decades and also suggests that founder genotype perspectives are in fact more dynamic when applied in the context of modern genomics-informed breeding.

Transcriptomic and evolutionary analyses of white pear (Pyrus bretschneideri) β-amylase genes reveals their importance for cold and drought stress responses

β-amylase (BAM) genes play essential roles in plant abiotic stress responses. Although the genome of Chinese white pear (Pyrus bretschneideri) has recently been made available, knowledge regarding the BAM family in pear, including gene function, evolutionary history and patterns of gene expression remains limited. In this study, we identified 17 PbBAMs in the pear genome. Of these, 12 PbBAM members were mapped onto 9 chromosomes and 5 PbBAM genes were located on scaffold contigs. Based on gene structure, protein motif analysis, and the topology of the phylogenetic tree of the PbBAM family, we classified member genes into 4 groups. All PbBAM genes were found to contain typical glycosyl hydrolysis 14 domain motifs. Interfamilial comparisons revealed that the phylogenetic relationships of BAM genes in other Rosaceae species were similar those found in pear. We also found that whole-genome duplication (WGD)/segmental duplication events played critical roles in the expansion of the BAM family. Next, we used transcriptomic data to study gene expression during the response of drought and low temperate responses, and found that genes in Group B were related to drought and cold stress. We identified four PbBAM genes associated with abiotic stress in Pear. Finally, by analyzing co-expression networks and co-regulatory genes, we found that PbBAM1a and PbBAM1b were associated with the pear abiotic stress response.

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.

Phylogeny and maternal donors of Elytrigia Desv. sensu lato (Triticeae; Poaceae) inferred from nuclear internal-transcribed spacer and trnL-F sequences

BACKGROUND

Elytrigia Desv. is a genus with a varied array of morphology, cytology, ecology, and distribution in Triticeae. Classification and systematic position of Elytrigia remain controversial. We used nuclear internal-transcribed spacer (nrITS) sequences and chloroplast trnL-F region to study the relationships of phylogenetic and maternal genome donor of Elytrigia Desv. sensu lato.

RESULTS

(1) E, F, P, St, and W genomes bear close relationship with one another and are distant from H and Ns genomes. E^e and E^b are homoeologous. (2) In ESt genome species, E genome is the origin of diploid Elytrigia species with E genome, St genome is the origin of Pseudoroegneria. (3) Diploid species Et. elongata were differentiated. (4) Et. stipifolia and Et. varnensis sequences are diverse based on nrITS data. (5) Et. lolioides contains St and H genomes and belongs to Elymus s. l. (6) E genome diploid species in Elytrigia serve as maternal donors of E genome for Et. nodosa (PI547344), Et. farcta, Et. pontica, Et. pycnantha, Et. scirpea, and Et. scythica. At least two species act as maternal donor of allopolyploids (ESt and EStP genomes).

CONCLUSIONS

Our results suggested that Elytrigia s. l. species contain different genomes, which should be divided into different genera. However, the genomes of Elytrigia species had close relationships with one another. Diploid species were differentiated, because of introgression and different geographical sources. The results also suggested that the same species and the same genomes of different species have different maternal donor. Further study of molecular biology and cytology could facilitate the evaluation of our results of phylogenetic in a more specific and accurate way.

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.

Molecular cytogenetic (FISH) and genome analysis of diploid wheatgrasses and their phylogenetic relationship

This paper reports detailed FISH-based karyotypes for three diploid wheatgrass species Agropyron cristatum (L.) Beauv., Thinopyrum bessarabicum (Savul.&Rayss) A. Löve, Pseudoroegneria spicata (Pursh) A. Löve, the supposed ancestors of hexaploid Thinopyrum intermedium (Host) Barkworth & D.R.Dewey, compiled using DNA repeats and comparative genome analysis based on COS markers. Fluorescence in situ hybridization (FISH) with repetitive DNA probes proved suitable for the identification of individual chromosomes in the diploid JJ, StSt and PP genomes. Of the seven microsatellite markers tested only the (GAA)_n trinucleotide sequence was appropriate for use as a single chromosome marker for the P. spicata A^S chromosome. Based on COS marker analysis, the phylogenetic relationship between diploid wheatgrasses and the hexaploid bread wheat genomes was established. These findings confirmed that the J and E genomes are in neighbouring clusters.

Beta-amylase gene variability in introgressive wheat lines

Variability of the beta-amylase gene in bread wheat, artificial amphidiploids, and derived introgression wheat lines was analyzed. Variation in homeologous beta-amylase sequences caused by the presence of MITE (Miniature Inverted-Repeat Transposable Element) and its footprint has been identified in bread wheat. The previously unknown location of MITE in Triticum urartu and T. aestivum L. beta-amylase gene has been found. These species have a MITE sequence in the third intron of beta-amylase, as opposed to Aegilops comosa and a number of other Triticeae species, which have it in the fourth intron. These two MITEs from Ae. comosa and T. aestivum were shown to have low identity scores. Miosa, an artificial amphidiploid, which has the M genome from Ae. comosa was shown to lose the MITE sequences. This loss might be caused by genomic shock due to allopolyploidization.

Characterization and phylogenetic analysis of α-gliadin gene sequences reveals significant genomic divergence in Triticeae species

Although the unique properties of wheat α-gliadin gene family are well characterized, little is known about the evolution and genomic divergence of α-gliadin gene family within the Triticeae. We isolated a total of 203 α-gliadin gene sequences from 11 representative diploid and polyploid Triticeae species, and found 108 sequences putatively functional. Our results indicate that α-gliadin genes may have possibly originated from wild Secale species, where the sequences contain the shortest repetitive domains and display minimum variation. A miniature inverted-repeat transposable element insertion is reported for the first time in α-gliadin gene sequence of Thinopyrum intermedium in this study, indicating that the transposable element might have contributed to the diversification of α-gliadin genes family among Triticeae genomes. The phylogenetic analyses revealed that the α-gliadin gene sequences of Dasypyrum, Australopyrum, Lophopyrum, Eremopyrum and Pseudoroengeria species have amplified several times. A search for four typical toxic epitopes for celiac disease within the Triticeae α-gliadin gene sequences showed that the α-gliadins of wild Secale, Australopyrum and Agropyron genomes lack all four epitopes, while other Triticeae species have accumulated these epitopes, suggesting that the evolution of these toxic epitopes sequences occurred during the course of speciation, domestication or polyploidization of Triticeae.

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.

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.

Phylogeny of a genomically diverse group of elymus (poaceae) allopolyploids reveals multiple levels of reticulation

The grass tribe Triticeae (=Hordeeae) comprises only about 300 species, but it is well known for the economically important crop plants wheat, barley, and rye. The group is also recognized as a fascinating example of evolutionary complexity, with a history shaped by numerous events of auto- and allopolyploidy and apparent introgression involving diploids and polyploids. The genus Elymus comprises a heterogeneous collection of allopolyploid genome combinations, all of which include at least one set of homoeologs, designated St, derived from Pseudoroegneria. The current analysis includes a geographically and genomically diverse collection of 21 tetraploid Elymus species, and a single hexaploid species. Diploid and polyploid relationships were estimated using four molecular data sets, including one that combines two regions of the chloroplast genome, and three from unlinked nuclear genes: phosphoenolpyruvate carboxylase, β-amylase, and granule-bound starch synthase I. Four gene trees were generated using maximum likelihood, and the phylogenetic placement of the polyploid sequences reveals extensive reticulation beyond allopolyploidy alone. The trees were interpreted with reference to numerous phenomena known to complicate allopolyploid phylogenies, and introgression was identified as a major factor in their history. The work illustrates the interpretation of complicated phylogenetic results through the sequential consideration of numerous possible explanations, and the results highlight the value of careful inspection of multiple independent molecular phylogenetic estimates, with particular focus on the differences among them.

What does the 5S rRNA multigene family tell us about the origin of the annual Triticeae (Poaceae)?

We have investigated the complex relationships among the annual genera within the tribe Triticeae through phylogenetic analyses of the 5S rRNA multigene family. Cloned sequences were assigned to groups of orthologous sequences, called unit classes, that were subjected to several analyses including BLAST (Basic Local Alignment Search Tool) searches to assess possible ancestral relationships with perennial genera; phylogenetic analyses using parsimony (Pars), maximum likelihood (ML), and Bayesian methods; and minimum reticulation networks from the Pars, ML, and Bayesian trees. In this study, we included genera with both annual and perennial species, such as Dasypyrum, Hordeum, and Secale. BLAST pointed to Pseudoroegneria (carrier of the St genome) and possibly Thinopyrum (carrier of the J genome) as the potential next of kin. However, Thinopyrum and Pseudoroegneria have never fallen together on the individual trees with the former generally associated with Crithopsis, Aegilops, Triticum, and Dasypyrum, while the latter is usually associated with the rest of the genera within Triticeae. The "long" unit classes placed Dasypyrum breviaristatum together with Dasypyrum villosum, whereas the "short" unit classes put them far apart on the trees. None of the gene trees alone was able to summarize the complex relationships among the genera, in line with previous results in the Triticeae. However, the application of tools designed to display phylogenetic networks was able to depict the complex links among the genera based on the short and the long gene trees, including the close link between Thinopyrum and Pseudoroegneria suggested by the phylogenetic analyses. In addition, our analyses provide support for the hypothesis that at least some annual Triticeae taxa are derived from their perennial relatives.

Reprint of: using nuclear gene data for plant phylogenetics: progress and prospects

The paper reviews the current state of low and single copy nuclear markers that have been applied successfully in plant phylogenetics to date, and discusses case studies highlighting the potential of massively parallel high throughput or next-generation sequencing (NGS) approaches for molecular phylogenetic and evolutionary investigations. The current state, prospects and challenges of specific single- or low-copy plant nuclear markers as well as phylogenomic case studies are presented and evaluated.

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

Background

Introgressive events (e.g., hybridization, gene flow, horizontal gene transfer) and incomplete lineage sorting of ancestral polymorphisms are a challenge for phylogenetic analyses since different genes may exhibit conflicting genealogical histories. Grasses of the Triticeae tribe provide a particularly striking example of incongruence among gene trees. Previous phylogenies, mostly inferred with one gene, are in conflict for several taxon positions. Therefore, obtaining a resolved picture of relationships among genera and species of this tribe has been a challenging task. Here, we obtain the most comprehensive molecular dataset to date in Triticeae, including one chloroplastic and 26 nuclear genes. We aim to test whether it is possible to infer phylogenetic relationships in the face of (potentially) large-scale introgressive events and/or incomplete lineage sorting; to identify parts of the evolutionary history that have not evolved in a tree-like manner; and to decipher the biological causes of gene-tree conflicts in this tribe.

Results

We obtain resolved phylogenetic hypotheses using the supermatrix and Bayesian Concordance Factors (BCF) approaches despite numerous incongruences among gene trees. These phylogenies suggest the existence of 4-5 major clades within Triticeae, with Psathyrostachys and Hordeum being the deepest genera. In addition, we construct a multigenic network that highlights parts of the Triticeae history that have not evolved in a tree-like manner. Dasypyrum, Heteranthelium and genera of clade V, grouping Secale, Taeniatherum, Triticum and Aegilops, have evolved in a reticulated manner. Their relationships are thus better represented by the multigenic network than by the supermatrix or BCF trees. Noteworthy, we demonstrate that gene-tree incongruences increase with genetic distance and are greater in telomeric than centromeric genes. Together, our results suggest that recombination is the main factor decoupling gene trees from multigenic trees.

Conclusions

Our study is the first to propose a comprehensive, multigenic phylogeny of Triticeae. It clarifies several aspects of the relationships among genera and species of this tribe, and pinpoints biological groups with likely reticulate evolution. Importantly, this study extends previous results obtained in Drosophila by demonstrating that recombination can exacerbate gene-tree conflicts in phylogenetic reconstructions.

An integrative test of the dead-end hypothesis of selfing evolution in Triticeae (Poaceae)

Self-fertilization is hypothesized to be an evolutionary dead end because reversion to outcrossing can rarely happen, and selfing lineages are thought to rapidly become extinct because of limited potential for adaptation and/or accumulation of deleterious mutations. We tested these two assumptions by combining morphological characters and molecular-evolution analyses in a tribe of hermaphroditic grasses (Triticeae). First, we determined the mating system of the 19 studied species. Then, we sequenced 27 protein-coding loci and compared base composition and substitution patterns between selfers and outcrossers. We found that the evolution of the mating system is best described by a model including outcrossing-to-selfing transitions only. At the molecular level, we showed that regions of low recombination exhibit signatures of relaxed selection. However, we did not detect any evidence of accumulation of nonsynonymous substitutions in selfers compared to outcrossers. Additionally, we tested for the potential deleterious effects of GC-biased gene conversion in outcrossing species. We found that recombination and not the mating system affected substitution patterns and base composition. We suggest that, in Triticeae, although recombination patterns have remained stable, selfing lineages are of recent origin and inbreeding may have persisted for insufficient time for differences between the two mating systems to evolve.

Reticulate evolutionary history of a complex group of grasses: phylogeny of Elymus StStHH allotetraploids based on three nuclear genes

BACKGROUND

Elymus (Poaceae) is a large genus of polyploid species in the wheat tribe Triticeae. It is polyphyletic, exhibiting many distinct allopolyploid genome combinations, and its history might be further complicated by introgression and lineage sorting. We focus on a subset of Elymus species with a tetraploid genome complement derived from Pseudoroegneria (genome St) and Hordeum (H). We confirm the species' allopolyploidy, identify possible genome donors, and pinpoint instances of apparent introgression or incomplete lineage sorting.

METHODOLOGY/PRINCIPAL FINDINGS

We sequenced portions of three unlinked nuclear genes-phosphoenolpyruvate carboxylase, beta-amylase, and granule-bound starch synthase I-from 27 individuals, representing 14 Eurasian and North American StStHHElymus species. Elymus sequences were combined with existing data from monogenomic representatives of the tribe, and gene trees were estimated separately for each data set using maximum likelihood. Trees were examined for evidence of allopolyploidy and additional reticulate patterns. All trees confirm the StStHH genome configuration of the Elymus species. They suggest that the StStHH group originated in North America, and do not support separate North American and European origins. Our results point to North American Pseudoroegneria and Hordeum species as potential genome donors to Elymus. Diploid P. spicata is a prospective St-genome donor, though conflict among trees involving P. spicata and the Eurasian P. strigosa suggests either introgression of GBSSI sequences from P. strigosa into North American Elymus and Pseudoroegneria, or incomplete lineage sorting of ancestral GBSSI polymorphism. Diploid H. californicum and/or allotetraploid H. jubatum are possible H-genome donors; direct involvement of an allotetraploid Hordeum species would simultaneously introduce two distinct H genomes to Elymus, consistent with some of the relationships among H-genome sequences in Hordeum and Elymus.

CONCLUSIONS/SIGNIFICANCE

Comparisons among molecular phylogenetic trees confirm allopolyploidy, identify potential genome donors, and highlight cases of apparent introgression or incomplete lineage sorting. The complicated history of this group emphasizes an inherent problem with interpreting conflicts among bifurcating trees-identifying introgression and determining its direction depend on which tree is chosen as a starting point of comparison. In spite of difficulties with interpretation, differences among gene trees allow us to identify reticulate species and develop hypotheses about underlying evolutionary processes.

Phylogeny and evolutionary history of Leymus (Triticeae; Poaceae) based on a single-copy nuclear gene encoding plastid acetyl-CoA carboxylase

Background

Single- and low- copy genes are less likely subject to concerted evolution, thus making themselves ideal tools for studying the origin and evolution of polyploid taxa. Leymus is a polyploid genus with a diverse array of morphology, ecology and distribution in Triticeae. The genomic constitution of Leymus was assigned as NsXm, where Ns was presumed to be originated from Psathyrostachys, while Xm represented a genome of unknown origin. In addition, little is known about the evolutionary history of Leymus. Here, we investigate the phylogenetic relationship, genome donor, and evolutionary history of Leymus based on a single-copy nuclear Acc1 gene.

Results

Two homoeologues of the Acc1 gene were isolated from nearly all the sampled Leymus species using allele-specific primer and were analyzed with those from 35 diploid taxa representing 18 basic genomes in Triticeae. Sequence diversity patterns and genealogical analysis suggested that (1) Leymus is closely related to Psathyrostachys, Agropyron, and Eremopyrum; (2) Psathyrostachys juncea is an ancestral Ns-genome donor of Leymus species; (3) the Xm genome in Leymus may be originated from an ancestral lineage of Agropyron and Eremopyrum triticeum; (4) the Acc1 sequences of Leymus species from the Qinghai-Tibetan plateau are evolutionarily distinct; (5) North America Leymus species might originate from colonization via the Bering land bridge; (6) Leymus originated about 11-12MYA in Eurasia, and adaptive radiation might have occurred in Leymus during the period of 3.7-4.3 MYA and 1.7-2.1 MYA.

Conclusion

Leymus species have allopolyploid origin. It is hypothesized that the adaptive radiation of Leymus species might have been triggered by the recent upliftings of the Qinghai-Tibetan plateau and subsequent climatic oscillations. Adaptive radiation may have promoted the rapid speciation, as well as the fixation of unique morphological characters in Leymus. Our results shed new light on our understanding of the origin of Xm genome, the polyploidization events and evolutionary history of Leymus that could account for the rich diversity and ecological adaptation of Leymus species.

Catalytically-inactive beta-amylase BAM4 required for starch breakdown in Arabidopsis leaves is a starch-binding-protein

Of the four chloroplast beta-amylase (BAM) proteins identified in Arabidopsis, BAM3 and BAM4 were previously shown to play the major roles in leaf starch breakdown, although BAM4 apparently lacks key active site residues and beta-amylase activity. Here we tested multiple BAM4 proteins with different N-terminal sequences with a range of glucan substrates and assay methods, but detected no alpha-1,4-glucan hydrolase activity. BAM4 did not affect BAM1, BAM2 or BAM3 activity even when added in 10-fold excess, nor the BAM3-catalysed release of maltose from isolated starch granules in the presence of glucan water dikinase. However, BAM4 binds to amylopectin and to amylose-Sepharose whereas BAM2 has very low beta-amylase activity and poor glucan binding. The low activity of BAM2 may be explained by poor glucan binding but absence of BAM4 activity is not. These results suggest that BAM4 facilitates starch breakdown by a mechanism involving direct interaction with starch or other alpha-1,4-glucan.

Allohexaploidy, introgression, and the complex phylogenetic history of Elymus repens (Poaceae)

The phylogenetic position of hexaploid Elymus repens within the tribe Triticeae (Poaceae) was examined using cloned sequences from the low-copy nuclear genes encoding phosphoenolpyruvate carboxylase (pepC) and beta-amylase. A previous analysis of E. repens using data from the nuclear granule-bound starch synthase I (GBSSI) gene had yielded five phylogenetically distinct gene copies, two more than expected from hexaploidy alone. The three gene trees share three distinct E. repens clades, suggesting that E. repens contains three phylogenetically divergent genomes, contributed by Hordeum, Pseudoroegneria, and an unknown donor. The two additional GBSSI sequences, including one that was apparently derived from outside of the tribe, appear to reflect past introgression of GBSSI sequences into the E. repens genome. On all three trees, the Hordeum-like E. repens sequences are polyphyletic within Hordeum, and the trees are in conflict with regard to the placement of these sequences within Hordeum, highlighting multiple contributions from Hordeum to E. repens.

Multiple homoplasious insertions and deletions of a Triticeae (Poaceae) DNA transposon: a phylogenetic perspective

Background

Stowaway elements are short, non-autonomous DNA transposons categorized as miniature inverted-repeat transposable elements (MITEs). The high MITE copy number in grass genomes suggests an active history of amplification and insertion, but ongoing MITE activity has only rarely been seen, and ongoing Stowaway activity has never been observed. Thus, a phylogenetic perspective on presence vs. absence of elements in an aligned data set can provide valuable historical insights into the dynamics of MITE acquisition and loss.

Results

A Stowaway-like element resides within the fourth intron of a β-amylase gene in representatives of five genera in the wheat tribe, Triticeae. Its presence vs. absence was examined with reference to the β-amylase gene tree topology, and in light of sequence comparisons of the β-amylase elements to Triticeae Stowaway elements in the Entrez nucleotide database. Among the sequences lacking the element, there are five distinct putative excision footprints (one widespread and four restricted to unrelated lineages) and two flanking deletions. The sequences that do contain elements are polyphyletic on the β-amylase tree, and their elements are divergent at the sequence level. The β-amylase elements do not form a monophyletic group relative to other Stowaway elements in Entrez; most are more similar to elements from other loci in other Triticeae genomes than they are to one another.

Conclusion

Combined, the phylogenetic distribution, sequence variation, and Entrez database comparisons indicate that a Stowaway-like element has undergone multiple deletions from and insertions into the same site in β-amylase intron 4 during the history of the tribe. The elements currently at the site represent multiple, distinct lineages that transcend generic boundaries. While patterns of Stowaway polymorphism across a phylogenetic data set do not allow evolutionary mechanisms to be inferred with certainty, they do provide insights into the dynamics of element evolution over an extended time scale. The historical perspective provided by a phylogenetic approach is complementary to the few studies in which ongoing MITE activity has been documented.

Widespread Genealogical Nonmonophyly in Species of Pinus Subgenus Strobus

Phylogenetic relationships among Pinus species from subgenus Strobus remain unresolved despite combined efforts based on nrITS and cpDNA. To provide greater resolution among these taxa, a 900-bp intron from a late embryogenesis abundant (LEA)-like gene (IFG8612)was sequenced from 39 pine species, with two or more alleles representing 33 species. Nineteen of 33 species exhibited allelic nonmonphyly in the strict consensus tree, and 10 deviated significantly from allelic monophyly based on topology incongruence tests. Intraspecific nucleotide diversity ranged from 0.0 to 0.0211, and analysis of variance shows that nucleotide diversity was strongly associated (P < 0.0001)with the degree of species monophyly. Although species nonmonophyly complicates phylogenetic interpretations, this nuclear locus offers greater topological support than previously observed for cpDNA or nrITS. Lacking evidence for hybridization, recombination, or imperfect taxonomy, we feel that incomplete lineage sorting remains the best explanation for the polymorphisms shared among species. Depending on the species, coalescent expectations indicate that reciprocal monophyly will be more likely than paraphyly in 1.71 to 24.0 x 10(6) years, and that complete genome-wide coalescence in these species may require up to 76.3 x 10(6) years. The absence of allelic coalescence is a severe constraint in the application of phylogenetic methods in Pinus, and taxa sharing similar life history traits with Pinus are likely to show species nonmonophyly using nuclear markers.

Phylogenetic relationships of Triticum and Aegilops and evidence for the origin of the A, B, and D genomes of common wheat (Triticum aestivum)

Common wheat (Triticum aestivum) has for decades been a textbook example of the evolution of a major crop species by allopolyploidization. Using a sophisticated extension of the PCR technique, we have successfully isolated two single-copy nuclear genes, DMC1 and EF-G, from each of the three genomes found in hexaploid wheat (BA(u)D) and from the two genomes of the tetraploid progenitor Triticum turgidum (BA(u)). By subjecting these sequences to phylogenetic analysis together with sequences from representatives of all the diploid Triticeae genera we are able for the first time to provide simultaneous and strongly supported evidence for the D genome being derived from Aegilops tauschii, the A(u) genome being derived from Triticum urartu, and the hitherto enigmatic B genome being derived from Aegilops speltoides. Previous problems of identifying the B genome donor may be associated with a higher diversification rate of the B genome compared to the A(u) genome in the polyploid wheats. The phylogenetic hypothesis further suggests that neither Triticum, Aegilops, nor Triticum plus Aegilops are monophyletic.