Genetic diversity and relationships among accessions of five crested wheatgrass species (Poaceae: Agropyron) based on gliadin analysis

Complete mitochondrial genome of Agropyron cristatum reveals gene transfer and RNA editing events

BACKGROUND

As an important forage in arid and semi-arid regions, Agropyron cristatum provides livestock with exceptionally high nutritional value. Additionally, A. cristatum exhibits outstanding genetic characteristics to endure drought and disease. Therefore, rich genetic diversity serves as a cornerstone for the improvement of major food crops. The purposes of this study were to systematically describe mitogenome of A.cristatum and preliminarily analyze its internal variations.

RESULT

The A. cristatum mitogenome was a single-ring molecular structure of 381,065 bp that comprised 52 genes, including 35 protein-coding, 3 rRNA and 14 tRNA genes. Among these, two pseudoprotein-coding genes and multiple copies of tRNA genes were observed. A total of 320 repetitive sequences was found to cover more than 10% of the mitogenome (105 simple sequences, 185 dispersed and 30 tandem repeats), which led to a large number of fragment rearrangements in the mitogenome of A. cristatum. Leucine was the most frequent amino acid (n = 1087,10.8%) in the protein-coding genes of A. cristatum mitogenome, and the highest usage codon was ATG (initiation codon). The number of A/T changes at the third base of the codon was much higher than that of G/C. Among 23 PCGs, the range of Pi values is from 0.0021 to 0.0539, with an average of 0.013. Additionally, 81 RNA editing sites were predicted, which were considerably fewer than those reported in other plant mitogenomes. Most of the RNA editing site base positions were concentrated at the first and second codon bases, which were C to T transitions. Moreover, we identified 95 sequence fragments (total length of 34, 343 bp) that were transferred from the chloroplast to mitochondria genes, introns, and intergenic regions. The stability of the tRNA genes was maintained during this process. Selection pressure analysis of 23 protein-coding genes shared by 15 Poaceae plants, showed that most genes were subjected to purifying selection during evolution, whereas rps4, cob, mttB, and ccmB underwent positive selection in different plants. Finally, a phylogenetic tree was constructed based on 22 plant mitogenomes, which showed that Agropyron plants have a high degree of independent heritability in Triticeae.

CONCLUSION

The findings of this study provide new data for a better understanding of A. cristatum genes, and demonstrate that mitogenomes are suitable for the study of plant classifications, such as those of Agropyron. Moreover, it provides a reference for further exploration of the phylogenetic relationships within Agropyron species, and establishes a theoretical basis for the subsequent development and utilization of A. cristatum plant germplasm resources.

Genotyping-by-Sequencing Enhances Genetic Diversity Analysis of Crested Wheatgrass [Agropyron cristatum (L.) Gaertn.]

Molecular characterization of unsequenced plant species with complex genomes is now possible by genotyping-by-sequencing (GBS) using recent next generation sequencing technologies. This study represents the first use of GBS application to sample genome-wide variants of crested wheatgrass [Agropyron cristatum (L.) Gaertn.] and assess the genetic diversity present in 192 genotypes from 12 tetraploid lines. Bioinformatic analysis identified 45,507 single nucleotide polymorphism (SNP) markers in this outcrossing grass species. The model-based Bayesian analysis revealed four major clusters of the samples assayed. The diversity analysis revealed 15.8% of SNP variation residing among the 12 lines, and 12.1% SNP variation present among four genetic clusters identified by the Bayesian analysis. The principal coordinates analysis and dendrogram were able to distinguish four lines of Asian origin from Canadian cultivars and breeding lines. These results serve as a valuable resource for understanding genetic variability, and will aid in the genetic improvement of this outcrossing polyploid grass species for forage production. These findings illustrate the potential of GBS application in the characterization of non-model polyploid plants with complex genomes.

RNA-Seq Analysis Provides the First Insights into the Phylogenetic Relationship and Interspecific Variation between Agropyron cristatum and Wheat

Agropyron cristatum, which is a wild grass of the tribe Triticeae, grows widely in harsh environments and provides many desirable genetic resources for wheat improvement. However, unclear interspecific phylogeny and genome-wide variation has limited the utilization of A. cristatum in the production of superior wheat varieties. In this study, by sequencing the transcriptome of the representative tetraploid A. cristatum Z559 and the common wheat variety Fukuhokomugi (Fukuho), which are often used as parents in a wide cross, their phylogenetic relationship and interspecific variation were dissected. First, 214,854 transcript sequences were assembled, and 3,457 orthologous genes related to traits of interest were identified in A. cristatum. Second, a total of 72 putative orthologous gene clusters were used to construct phylogenetic relationships among A. cristatum, Triticeae and other genomes. A clear division between A. cristatum and the other Triticeae species was revealed. Third, the sequence similarity of most genes related to traits of interest is greater than 95% between A. cristatum and wheat. Therefore, using the 5% mismatch parameter for A. cristatum, we mapped the transcriptome sequencing data to wheat reference sequences to discover the variations between A. cristatum and wheat and 862,340 high-quality variants were identified. Additionally, compared with the wheat A and B genomes, the P and D genomes displayed an obviously larger variant density and a longer evolutionary distance, suggesting that A. cristatum is more distantly related to the wheat D genome. Finally, by using Kompetitive Allele Specific PCR array (KASPar) technology, 37 of 53 (69.8%) SNPs were shown to be genuine in Z559, Fukuho, and additional lines with seven different P chromosomes, and function of the genes in which these SNPs are located were also determined. This study provides not only the first insights into the phylogenetic relationships between the P genome and Triticeae but also genetic resources for gene discovery and specific marker development in A. cristatum, and this information will be vital for future wheat-breeding efforts. The sequence data have been deposited in the Sequence Read Archive (SRA) database at the NCBI under accession number SRP090613.

A high-density genetic map for P genome of Agropyron Gaertn. based on specific-locus amplified fragment sequencing (SLAF-seq)

This genetic map for Agropyron Gaertn. contained 1023 markers on seven linkage groups, with a total of 907.8 cM and an average distance of 1.5 cM between adjacent loci. Many wheat- Agropyron cristatum derivative lines exhibit superior agronomic traits, and part of them are valuable for future wheat breeding. To date, no high-density genetic map for Agropyron Gaertn. has been published. Specific-locus amplified fragment sequencing (SLAF-seq), a recently developed strategy for large scale de novo discovery and genotyping of single nucleotide polymorphisms (SNPs), was employed in this study to develop sufficient markers for a segregating Agropyron F1 population derived from an interspecific cross between two cross-pollinated diploid collections A. cristatum (L.) Beauv. 'Z1842' and A. mongolicum Keng 'Z2098'. In total, we obtained raw data consisting of 128,932,358 pair-end reads of ~80 bp long after sequencing. Then 69,325 high-quality SLAFs were detected, of which 26,248 SLAFs were polymorphic and 1752 of the polymorphic markers were used for the genetic map construction. The final map contained 1023 markers on the seven linkage groups (LGs), which spanned a total of 907.8 cM with an average number of 146 markers and 89 loci per LG and an average distance of 1.5 cM between adjacent loci. To our knowledge, this map is the densest genetic linkage map for Agropyron so far. Through BLAT alignment of Agropyron SLAF marker sequences with the draft genome assemblies of wheat and barley, the Agropyron LGs were assigned as LG1-7 according to their corresponding homoeologous chromosomal groups of wheat. Results of this study will not only provide a platform for gene/QTL fine mapping, but also serve as a reference to assist the assembling of the P genome sequence in future.

De novo transcriptome sequencing of Agropyron cristatum to identify available gene resources for the enhancement of wheat

Agropyron cristatum is a wild grass of the tribe Triticeae that is widely grown in harsh environments. As a wild relative of wheat, A. cristatum carries many resistance genes that could be used to broaden the genetic diversity of wheat. Here, we report the transcriptome sequencing of the flag leaf and young spike tissues of a representative tetraploid A. cristatum. More than 90 million reads from the two tissues were assembled into 73,664 unigenes. All unigenes were functionally annotated against the KEGG, COG, and Gene Ontology databases and predicted long non-coding RNAs. Pfam prediction demonstrates that A. cristatum carries an abundance of stress resistance genes. The extent of specific genes and rare alleles make A. cristatum a vital genetic reservoir for the improvement of wheat. Altogether, the available gene resources in A. cristatum facilitate efforts to harness the genetic diversity of wild relatives to enhance wheat.