RNA-seq analysis reveals considerable genetic diversity and provides genetic markers saturating all chromosomes in the diploid wild wheat relative Aegilops umbellulata

Discrepancy of flowering time between genetically close sublineages of Aegilops umbellulata Zhuk

Aegilops umbellulata Zhuk., a wild diploid wheat-related species, has been used as a genetic resource for several important agronomic traits. However, its genetic variations have not been comprehensively studied. We sequenced RNA from 114 accessions of Ae. umbellulata to evaluate DNA polymorphisms and phenotypic variations. Bayesian clustering and phylogenetic analysis based on SNPs detected by RNA sequencing revealed two divergent lineages, UmbL1 and UmbL2. The main differences between them were in the sizes of spikes and spikelets, and culm diameter. UmbL1 is divided into two sublineages, UmbL1e and UmbL1w. These genetic differences corresponded to geographic distributions. UmbL1e, UmbL1w, and UmbL2 are found in Turkey, Iran/Iraq, and Greece, respectively. Although UmbL1e and UmbL1w were genetically similar, flowering time and other morphological traits were more distinct between these sublineages than those between the lineages. This discrepancy can be explained by the latitudinal and longitudinal differences in habitats. Specifically, latitudinal clines of flowering time were clearly observed in Ae. umbellulata, strongly correlated with solar radiation in the winter season. This observation implies that latitudinal differences are a factor in differences in the flowering times of Ae. umbellulata. Differences in flowering time could influence other morphological differences and promote genetic divergence between sublineages.

Chromosome-scale assembly of the wild wheat relative Aegilops umbellulata

Wild wheat relatives have been explored in plant breeding to increase the genetic diversity of bread wheat, one of the most important food crops. Aegilops umbellulata is a diploid U genome-containing grass species that serves as a genetic reservoir for wheat improvement. In this study, we report the construction of a chromosome-scale reference assembly of Ae. umbellulata accession TA1851 based on corrected PacBio HiFi reads and chromosome conformation capture. The total assembly size was 4.25 Gb with a contig N50 of 17.7 Mb. In total, 36,268 gene models were predicted. We benchmarked the performance of hifiasm and LJA, two of the most widely used assemblers using standard and corrected HiFi reads, revealing a positive effect of corrected input reads. Comparative genome analysis confirmed substantial chromosome rearrangements in Ae. umbellulata compared to bread wheat. In summary, the Ae. umbellulata assembly provides a resource for comparative genomics in Triticeae and for the discovery of agriculturally important genes.

Comparative genomic analysis of 5Mg chromosome of Aegilops geniculata and 5Uu chromosome of Aegilops umbellulata reveal genic diversity in the tertiary gene pool

Wheat is one of the most important cereal crops for the global food security. Due to its narrow genetic base, modern bread wheat cultivars face challenges from increasing abiotic and biotic stresses. Since genetic improvement is the most sustainable approach, finding novel genes and alleles is critical for enhancing the genetic diversity of wheat. The tertiary gene pool of wheat is considered a gold mine for genetic diversity as novel genes and alleles can be identified and transferred to wheat cultivars. Aegilops geniculata and Ae. umbellulata are the key members of the tertiary gene pool of wheat and harbor important genes against abiotic and biotic stresses. Homoeologous-group five chromosomes (5U^u and 5M^g) have been extensively studied from Ae. geniculata and Ae. umbellulata as they harbor several important genes including Lr57, Lr76, Yr40, Yr70, Sr53 and chromosomal pairing loci. In the present study, using chromosome DNA sequencing and RNAseq datasets, we performed comparative analysis to study homoeologous gene evolution in 5M^g, 5U^u, and group 5 wheat chromosomes. Our findings highlight the diversity of transcription factors and resistance genes, resulting from the differential expansion of the gene families. Both the chromosomes were found to be enriched with the "response to stimulus" category of genes providing resistance against biotic and abiotic stress. Phylogenetic study positioned the M genome closer to the D genome, with higher proximity to the A genome than the B genome. Over 4000 genes were impacted by SNPs on 5D, with 4-5% of those genes displaying non-disruptive variations that affect gene function.

RNA-seq analysis revealed considerable genetic diversity and enabled the development of specific KASP markers for Psathyrostachys huashanica

Psathyrostachys huashanica, which grows exclusively in Huashan, China, is an important wild relative of common wheat that has many desirable traits relevant for wheat breeding. However, the poorly characterized interspecific phylogeny and genomic variations and the relative lack of species-specific molecular markers have limited the utility of P. huashanica as a genetic resource for enhancing wheat germplasm. In this study, we sequenced the P. huashanica transcriptome, resulting in 50,337,570 clean reads that were assembled into 65,617 unigenes, of which 38,428 (58.56%) matched at least one sequence in public databases. The phylogenetic analysis of P. huashanica, Triticeae species, and Poaceae species was conducted using 68 putative orthologous gene clusters. The data revealed the distant evolutionary relationship between P. huashanica and common wheat as well as the substantial diversity between the P. huashanica genome and the wheat D genome. By comparing the transcriptomes of P. huashanica and Chinese Spring, 750,759 candidate SNPs between P. huashanica Ns genes and their common wheat orthologs were identified. Among the 90 SNPs in the exon regions with different functional annotations, 58 (64.4%) were validated as Ns genome-specific SNPs in the common wheat background by KASP genotyping assays. Marker validation analyses indicated that six specific markers can discriminate between P. huashanica and the other wheat-related species. In addition, five markers are unique to P. huashanica, P. juncea, and Leymus species, which carry the Ns genome. The Ns genome-specific markers in a wheat background were also validated regarding their specificity and stability for detecting P. huashanica chromosomes in four wheat-P. huashanica addition lines. Four and eight SNP markers were detected in wheat-P. huashanica 2Ns and 7Ns addition lines, respectively, and one marker was specific to both wheat-P. huashanica 3Ns, 4Ns, and 7Ns addition lines. These markers developed using transcriptome data may be used to elucidate the genetic relationships among Psathyrostachys, Leymus, and other closely-related species. They may also facilitate precise introgressions and the high-throughput monitoring of P. huashanica exogenous chromosomes or segments in future crop breeding programs.

Evolution of weedy giant ragweed (Ambrosia trifida): Multiple origins and gene expression variability facilitates weediness

Agricultural weeds may originate from wild populations, but the origination patterns and genetics underlying this transition remain largely unknown. Analysis of weedy-wild paired populations from independent locations may provide evidence to identify key genetic variation contributing to this adaptive shift. We performed genetic variation and expression analyses on transcriptome data from 67 giant ragweed samples collected from different locations in Ohio, Iowa, and Minnesota and found geographically separated weedy populations likely originated independently from their adjacent wild populations, but subsequent spreading of weedy populations also occurred locally. By using eight closely related weedy-wild paired populations, we identified thousands of unique transcripts in weedy populations that reflect shared or specific functions corresponding, respectively, to both convergently evolved and population-specific weediness processes. In addition, differential expression of specific groups of genes was detected between weedy and wild giant ragweed populations using gene expression diversity and gene co-expression network analyses. Our study suggests an integrated route of weedy giant ragweed origination, consisting of independent origination combined with the subsequent spreading of certain weedy populations, and provides several lines of evidence to support the hypothesis that gene expression variability plays a key role in the evolution of weedy species.

Molecular cytogenetics for a wheat-Aegilops geniculata 3Mg alien addition line with resistance to stripe rust and powdery mildew

BACKGROUND

Aegilops geniculata Roth is closely related to common wheat (Triticum aestivum L.) and is a valuable genetic resource for improvement of wheat.

RESULTS

In this study, the W19513 line was derived from the BC₁F₁₀ progeny of a cross between wheat 'Chinese Spring' and Ae. geniculata SY159. Cytological examination showed that W19513 contained 44 chromosomes. Twenty-two bivalents were formed at the first meiotic metaphase I in the pollen mother cellsand the chromosomes were evenly distributed to opposite poles at meiotic anaphase I. Genomic in situ hybridization demonstrated that W19513 carried a pair of alien chromosomes from the M genome. Fluorescence in situ hybridization confirmed detection of variation in chromosomes 4A and 6B. Functional molecular marker analysis using expressed sequence tag-sequence-tagged site and PCR-based landmark unique gene primers revealed that the alien gene belonged to the third homologous group. The marker analysis confirmed that the alien chromosome pair was 3M^g. In addition, to further explore the molecular marker specificity of chromosome 3M^g, based on the specific locus amplified fragment sequencing technique, molecular markers specific for W19513 were developed with efficiencies of up to 47.66%. The W19513 line was inoculated with the physiological race E09 of powdery mildew (Blumeria graminis f. sp. tritici) at the seedling stage and showed moderate resistance. Field inoculation with a mixture of the races CYR31, CYR32, CYR33, and CYR34 of the stripe rust fungus (Puccinia striiformis f. sp. triticii) revealed that the line W19513 showed strong resistance.

CONCLUSIONS

This study provides a foundation for use of the line W19513 in future genetic research and wheat improvement.

RNA-Seq-based DNA marker analysis of the genetics and molecular evolution of Triticeae species

The release of high-quality chromosome-level genome sequences of members of the Triticeae tribe has greatly facilitated genetic and genomic analyses of important crops such as wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to the large diploid genome size of Triticeae plants (ca. 5 Gbp), transcript analysis is an important method for identifying genetic and genomic differences among Triticeae species. In this review, we summarize our results of RNA-Seq analyses of diploid wheat accessions belonging to the genera Aegilops and Triticum. We also describe studies of the molecular relationships among these accessions and provide insight into the evolution of common hexaploid wheat. DNA markers based on polymorphisms within species can be used to map loci of interest. Even though the genome sequence of diploid Aegilops tauschii, the D-genome donor of common wheat, has been released, the diploid barley genome continues to provide key information about the physical structures of diploid wheat genomes. We describe how a series of RNA-Seq analyses of wheat relatives has helped uncover the structural and evolutionary features of genomic and genetic systems in wild and cultivated Triticeae species.

Transcriptome analysis of the impact of exogenous methyl jasmonate on the opening of sorghum florets

Background

Methyl Jasmonate (MeJA) could promote the opening of sorghum florets, but the molecular mechanism remains unclear.

Objective

We aimed to investigate the molecular mechanism of exogenous MeJA in promoting the opening of sorghum florets.

Methods

Hybrid sorghum Aikang-8 was selected as the test material in this study. Sorghum plants of uniform growth with approximately 20%-25% florets open were selected and treated with 0, 0.5 and 2.0 mmol/L of MeJA. Totally there were 27 samples with lodicules removed were obtained at different time points and used for the transcriptome analysis using the BGISEQ_500RS platform.

Results

The results showed the sorghum florets opened earlier than the control after the treatment with exogenous MeJA, and the promotive effect increased along with the increase of exogenous MeJA concentration. The number of differentially expressed genes (DEGs) in plasma cells increased with the increase of MeJA concentration, whether up- or down-regulated, after the exogenous MeJA treatment. Besides, the number of metabolic pathways was also positively correlated with the concentration of MeJA. GO and KEGG analysis suggested the DEGs were mainly enriched in starch and sucrose metabolism-related pathways (i.e., LOC8063704, LOC8083539 and LOC8056206), plant hormone signal transduction pathways (i.e., LOC8084842, LOC8072010, and LOC8057408), energy metabolic pathway (i.e., LOC8076139) and the α-linolenic acid metabolic pathway (i.e., LOC8055636, LOC8057399, LOC8063048 and LOC110430730). Functional analysis of target genes showed that two genes named LOC-1 (LOC8063704) and LOC-2 (LOC8076139) could induce the earlier flowering of Arabidopsis thaliana.

Conclusion

The results of this study suggest that exogenous MeJA treatments could induce the up- or down- regulation of genes related to starch and sucrose metabolism, -linolenic acid metabolism and plant hormone signal transduction pathways in the plasma cells of sorghum florets, thereby promoting the opening of sorghum florets.

Wheat omics: Classical breeding to new breeding technologies

Wheat is an important cereal crop, and its significance is more due to compete for dietary products in the world. Many constraints facing by the wheat crop due to environmental hazardous, biotic, abiotic stress and heavy matters factors, as a result, decrease the yield. Understanding the molecular mechanism related to these factors is significant to figure out genes regulate under specific conditions. Classical breeding using hybridization has been used to increase the yield but not prospered at the desired level. With the development of newly emerging technologies in biological sciences i.e., marker assisted breeding (MAB), QTLs mapping, mutation breeding, proteomics, metabolomics, next-generation sequencing (NGS), RNA_sequencing, transcriptomics, differential expression genes (DEGs), computational resources and genome editing techniques i.e. (CRISPR cas9; Cas13) advances in the field of omics. Application of new breeding technologies develops huge data; considerable development is needed in bioinformatics science to interpret the data. However, combined omics application to address physiological questions linked with genetics is still a challenge. Moreover, viroid discovery opens the new direction for research, economics, and target specification. Comparative genomics important to figure gene of interest processes are further discussed about considering the identification of genes, genomic loci, and biochemical pathways linked with stress resilience in wheat. Furthermore, this review extensively discussed the omics approaches and their effective use. Integrated plant omics technologies have been used viroid genomes associated with CRISPR and CRISPR-associated Cas13a proteins system used for engineering of viroid interference along with high-performance multidimensional phenotyping as a significant limiting factor for increasing stress resistance in wheat.

Analysis of Structural Genomic Diversity in Aegilops umbellulata, Ae. markgrafii, Ae. comosa, and Ae. uniaristata by Fluorescence In Situ Hybridization Karyotyping

Fluorescence in situ hybridization karyotypes have been widely used for evolutionary analysis on chromosome organization and genetic/genomic diversity in the wheat alliance (tribe Triticeae of Poaceae). The karyotpic diversity of Aegilops umbellulata, Ae. markgrafii, Ae. comosa subsp. comosa and subsp. subventricosa, and Ae. uniaristata was evaluated by the fluorescence in situ hybridization (FISH) probes oligo-pSc119.2 and pTa71 in combination with (AAC)5, (ACT)7, and (CTT)12, respectively. Abundant intra- and interspecific genetic variation was discovered in Ae. umbellulata, Ae. markgrafii, and Ae. comosa, but not Ae. uniaristata. Chromosome 7 of Ae. umbellulata had more variants (six variants) than the other six U chromosomes (2-3 variants) as revealed by probes oligo-pSc119.2 and (AAC)5. Intraspecific variation in Ae. markgrafii and Ae. comosa was revealed by oligo-pSc119.2 in combination with (ACT)7 and (CTT)12, respectively. At least five variants were found in every chromosome of Ae. markgrafii and Ae. comosa, and up to 18, 10, and 15 variants were identified for chromosomes 2 of Ae. markgrafii, 4 of Ae. comosa subsp. comosa, and 6 of Ae. comosa subsp. subventricosa. The six Ae. uniaristata accessions showed identical FISH signal patterns. A large number of intra-specific polymorphic FISH signals were observed between the homologous chromosomes of Ae. markgrafii and Ae. comosa, especially chromosomes 1, 2, 4, and 7 of Ae. markgrafii, chromosome 4 of Ae. comosa subsp. comosa, and chromosome 6 of Ae. comosa subsp. subventricosa. Twelve Ae. comosa and 24 Ae. markgrafii accessions showed heteromorphism between homologous chromosomes. Additionally, a translocation between the short arms of chromosomes 1 and 7 of Ae. comosa PI 551038 was identified. The FISH karyotypes can be used to clearly identify the chromosome variations of each chromosome in these Aegilops species and also provide valuable information for understanding the evolutionary relationships and structural genomic variation among Aegilops species.

Appraisal of wheat genomics for gene discovery and breeding applications: a special emphasis on advances in Asia

We discussed the most recent efforts in wheat functional genomics to discover new genes and their deployment in breeding with special emphasis on advances in Asian countries. Wheat research community is making significant progress to bridge genotype-to-phenotype gap and then applying this knowledge in genetic improvement. The advances in genomics and phenomics have intrigued wheat researchers in Asia to make best use of this knowledge in gene and trait discovery. These advancements include, but not limited to, map-based gene cloning, translational genomics, gene mapping, association genetics, gene editing and genomic selection. We reviewed more than 57 homeologous genes discovered underpinning important traits and multiple strategies used for their discovery. Further, the complementary advancements in wheat phenomics and analytical approaches to understand the genetics of wheat adaptability, resilience to climate extremes and resistance to pest and diseases were discussed. The challenge to build a gold standard reference genome sequence of bread wheat is now achieved and several de novo reference sequences from the cultivars representing different gene pools will be available soon. New pan-genome sequencing resources of wheat will strengthen the foundation required for accelerated gene discovery and provide more opportunities to practice the knowledge-based breeding.

Phenotypic effects of the U-genome variation in nascent synthetic hexaploids derived from interspecific crosses between durum wheat and its diploid relative Aegilops umbellulata

Aegilops umbellulata is a wild diploid wheat species with the UU genome that is an important genetic resource for wheat breeding. To exploit new synthetic allohexaploid lines available as bridges for wheat breeding, a total of 26 synthetic hexaploid lines were generated through crossing between the durum wheat cultivar Langdon and 26 accessions of Ae. umbellulata. In nascent synthetic hexaploids with the AABBUU genome, the presence of the set of seven U-genome chromosomes was confirmed with U-genome chromosome-specific markers developed based on RNA-seq-derived data from Ae. umbellulata. The AABBUU synthetic hexaploids showed large variations in flowering- and morphology-related traits, and these large variations transmitted well from the parental Ae. umbellulata accessions. However, the variation ranges in most traits examined were reduced under the AABBUU hexaploid background compared with under the diploid parents. The AABBUU and AABBDD synthetic hexaploids were clearly discriminated by several morphological traits, and an increase of plant height and in the number of spikes and a decrease of spike length were commonly observed in the AABBUU synthetics. Thus, interspecific differences in several morphological traits between Ae. umbellulata and A. tauschii largely affected the basic plant architecture of the synthetic hexaploids. In conclusion, the AABBUU synthetic hexaploid lines produced in the present study are useful resources for the introgression of desirable genes from Ae. umbellulata to common wheat.

Diploid genome differentiation conferred by RNA sequencing-based survey of genome-wide polymorphisms throughout homoeologous loci in Triticum and Aegilops

Background

Triticum and Aegilops diploid species have morphological and genetic diversity and are crucial genetic resources for wheat breeding. According to the chromosomal pairing-affinity of these species, their genome nomenclatures have been defined. However, evaluations of genome differentiation based on genome-wide nucleotide variations are still limited, especially in the three genomes of the genus Aegilops: Ae. caudata L. (CC genome), Ae. comosa Sibth. et Sm. (MM genome), and Ae. uniaristata Vis. (NN genome). To reveal the genome differentiation of these diploid species, we first performed RNA-seq-based polymorphic analyses for C, M, and N genomes, and then expanded the analysis to include the 12 diploid species of Triticum and Aegilops.

Results

Genetic divergence of the exon regions throughout the entire chromosomes in the M and N genomes was larger than that between A- and A^m-genomes. Ae. caudata had the second highest genetic diversity following Ae. speltoides, the putative B genome donor of common wheat. In the phylogenetic trees derived from the nuclear and chloroplast genome-wide polymorphism data, the C, D, M, N, U, and S genome species were connected with short internal branches, suggesting that these diploid species emerged during a relatively short period in the evolutionary process. The highly consistent nuclear and chloroplast phylogenetic topologies indicated that nuclear and chloroplast genomes of the diploid Triticum and Aegilops species coevolved after their diversification into each genome, accounting for most of the genome differentiation among the diploid species.

Conclusions

RNA-sequencing-based analyses successfully evaluated genome differentiation among the diploid Triticum and Aegilops species and supported the chromosome-pairing-based genome nomenclature system, except for the position of Ae. speltoides. Phylogenomic and epigenetic analyses of intergenic and centromeric regions could be essential for clarifying the mechanisms behind this inconsistency.

Cytogenetic Analysis and Molecular Marker Development for a New Wheat-Thinopyrum ponticum 1Js (1D) Disomic Substitution Line With Resistance to Stripe Rust and Powdery Mildew

Thinopyrum ponticum (2n = 10x = 70), a member of the tertiary gene pool of wheat (Triticum aestivum L.), harbors many biotic and abiotic stress resistance genes. CH10A5, a novel disomic substitution line from a cross of T. aestivum cv. 7182 and Th. ponticum, was characterized by cytogenetic identification, in situ hybridization, molecular marker analysis, and morphological investigation of agronomic traits and disease resistance. Cytological observations showed that CH10A5 contained 42 chromosomes and formed 21 bivalents at meiotic metaphase I. Genome in situ hybridization (GISH) analysis indicated that two of its chromosomes came from the J^s genome of Th. ponticum, and wheat 15K array mapping and fluorescence in situ hybridization (FISH) revealed that chromosome 1D was absent from CH10A5. Polymorphic analysis of molecular markers indicated that the pair of alien chromosomes belonged to homoeologous group one, designated as 1J^s. Thus, CH10A5 was a wheat-Th. ponticum 1J^s (1D) disomic substitution line. Field disease resistance trials demonstrated that the introduced Th. ponticum chromosome 1J^s was probably responsible for resistance to both stripe rust and powdery mildew at the adult stage. Based on specific-locus amplified fragment sequencing (SLAF-seq), 507 STS molecular markers were developed to distinguish chromosome 1J^s genetic material from that of wheat. Of these, 49 STS markers could be used to specifically identify the genetic material of Th. ponticum. CH10A5 will increase the resistance gene diversity of wheat breeding materials, and the markers developed here will permit further tracing of heterosomal chromosome fragments in the future.

Genome-wide polymorphisms from RNA sequencing assembly of leaf transcripts facilitate phylogenetic analysis and molecular marker development in wild einkorn wheat

A survey of genome-wide polymorphisms between closely related species is required to understand the molecular basis of the evolutionary differentiation of their genomes. Two wild diploid wheat species, namely Triticum monococcum ssp. aegilopoides and T. urartu, are closely related and harbour the A^m and A genomes, respectively. The A-genome donor of tetraploid and common wheat is T. urartu, and T. monococcum ssp. monococcum is the cultivated form derived from the wild einkorn wheat subspecies aegilopoides. Although subspecies aegilopoides has been a useful genetic resource in wheat breeding, genome-wide molecular markers for this subspecies have not been sufficiently developed. Here, we describe the detection of genome-wide polymorphisms such as single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) from RNA sequencing (RNA-seq) data of leaf transcripts in 15 accessions of the two diploid wheat species. The SNPs and indels, detected using the A genome of common wheat as the reference genome, covered the entire chromosomes of these species. The polymorphism information facilitated a comparison of the genetic diversity of einkorn wheat with that of two related diploid Aegilops species, namely, Ae. tauschii and Ae. umbellulata. Cleaved amplified polymorphic sequence (CAPS) markers converted from the SNP data were efficiently developed to confirm the addition of aegilopoides subspecies chromosomes to tetraploid wheat in nascent allohexaploid lines with AABBA^mA^m genomes. In addition, the CAPS markers permitted linkage map construction in mapping populations of aegilopoides subspecies accessions. Therefore, these RNA-seq data provide information for further breeding of closely related species with no reference genome sequence data.

Origin of wheat B-genome chromosomes inferred from RNA sequencing analysis of leaf transcripts from section Sitopsis species of Aegilops

Dramatic changes occasionally occur in intergenic regions leading to genomic alterations during speciation and will consequently obscure the ancestral species that have contributed to the formation of allopolyploid organisms. The S genome of five species of section Sitopsis of genus Aegilops is considered to be an origin of B-genome in cultivated tetraploid and hexaploid wheat species, although its actual donor is still unclear. Here, we attempted to elucidate phylogenetic relationship among Sitopsis species by performing RNA sequencing of the coding regions of each chromosome. Thus, genome-wide polymorphisms were extensively analyzed in 19 accessions of the Sitopsis species in reference to the tetraploid and hexaploid wheat B genome sequences and consequently were efficiently anchored to the B-genome chromosomes. The results of our genome-wide exon sequencing and resultant phylogenetic analysis indicate that Ae. speltoides is likely to be the direct donor of all chromosomes of the wheat B genome. Our results also indicate that the genome differentiation during wheat allopolyploidization from S to B proceeds at different speeds over the chromosomes rather than at constant rate and recombination could be a factor determining the speed. This observation is potentially generalized to genome differentiation during plant allopolyploid evolution.