Rare occurrence of cryptic 5' splice sites by downstream 3' splice site/exon boundary mutations in a heavy-ion-induced egy1-4 allele of Arabidopsis thaliana.

Pre-mRNA splicing is a fundamental process in eukaryotic gene expression, and the mechanism of intron definition, involving the recognition of the canonical GU (5'-splice site) and AG (3'-splice site) dinucleotides by splicing factors, has been postulated for most cases of splicing initiation in plants. Splice site mutations have played crucial roles in unraveling the mechanism of pre-mRNA splicing in planta. Typically, splice site mutations abolish splicing events or activate one or more cryptic splice sites surrounding the mutated region. In this report, we investigated the splicing pattern of the EGY1 gene in an Ar-ion-induced egy1-4 allele of Arabidopsis thaliana. egy1-4 has an AG-to-AC mutation in the 3'-end of intron 3, along with 4-bp substitutions and a 5-bp deletion in adjacent exon 4. RT-PCR, cDNA cloning, and amplicon sequencing analyses of EGY1 revealed that while most wild-type EGY1 mRNAs had a single splicing pattern, egy1-4 mRNAs had multiple splicing defects. Almost half of EGY1 transcripts showed 'intron retention' at intron 3, while the other half exhibited activation of 3' cryptic splice sites either upstream or downstream of the original 3'-splice site. Unexpectedly, around 8% of EGY1 transcripts in egy1-4 exhibited activation of cryptic 5'-splice sites positioned upstream of the authentic 5'-splice site of intron 3. Whole genome resequencing of egy1-4 indicated that it has no other known impactful mutations. These results may provide a rare, but real case of activation of cryptic 5'-splice sites by downstream 3'-splice site/exon mutations in planta.

Transcriptome profiles of rice roots under simulated microgravity conditions and following gravistimulation.

Root system architecture affects the efficient uptake of water and nutrients in plants. The root growth angle, which is a critical component in determining root system architecture, is affected by root gravitropism; however, the mechanism of root gravitropism in rice remains largely unknown. In this study, we conducted a time-course transcriptome analysis of rice roots under conditions of simulated microgravity using a three-dimensional clinostat and following gravistimulation to detect candidate genes associated with the gravitropic response. We found that HEAT SHOCK PROTEIN (HSP) genes, which are involved in the regulation of auxin transport, were preferentially up-regulated during simulated microgravity conditions and rapidly down-regulated by gravistimulation. We also found that the transcription factor HEAT STRESS TRANSCRIPTION FACTOR A2s (HSFA2s) and HSFB2s, showed the similar expression patterns with the HSPs. A co-expression network analysis and an in silico motif search within the upstream regions of the co-expressed genes revealed possible transcriptional control of HSPs by HSFs. Because HSFA2s are transcriptional activators, whereas HSFB2s are transcriptional repressors, the results suggest that the gene regulatory networks governed by HSFs modulate the gravitropic response through transcriptional control of HSPs in rice roots.

iPOTs: Internet of Things-based pot system controlling optional treatment of soil water condition for plant phenotyping under drought stress.

A cultivation facility that can assist users in controlling the soil water condition is needed for accurately phenotyping plants under drought stress in an artificial environment. Here we report the Internet of Things-based pot system controlling optional treatment of soil water condition (iPOTs), an automatic irrigation system that mimics the drought condition in a growth chamber. The Wi-Fi-enabled iPOTs system allows water supply from the bottom of the pot, based on the soil water level set by the user, and automatically controls the soil water level at a desired depth. The iPOTs also allows users to monitor environmental parameters, such as soil temperature, air temperature, humidity, and light intensity, in each pot. To verify whether the iPOTs mimics the drought condition, we conducted a drought stress test on rice (Oryza sativa L.) varieties and near-isogenic lines, with diverse root system architecture, using the iPOTs system installed in a growth chamber. Similar to the results of a previous drought stress field trial, the growth of shallow-rooted rice accessions was severely affected by drought stress compared with that of deep-rooted accessions. The microclimate data obtained using the iPOTs system increased the accuracy of plant growth evaluation. Transcriptome analysis revealed that pot positions in the growth chamber had little impact on plant growth. Together, these results suggest that the iPOTs system is a reliable platform for phenotyping plants under drought stress.

The transcriptomic landscapes of rice cultivars with diverse root system architectures grown in upland field conditions.

Root system architecture affects plant drought resistance and other key agronomic traits such as lodging. However, although phenotypic and genomic variation has been extensively analyzed, few field studies have integrated phenotypic and transcriptomic information, particularly for below-ground traits such as root system architecture. Here, we report the phenotypic and transcriptomic landscape of 61 rice (Oryza sativa) accessions with highly diverse below-ground traits grown in an upland field. We found that four principal components explained the phenotypic variation and that accessions could be classified into four subpopulations (indica, aus, japonica and admixed) based on their tiller numbers and crown root diameters. Transcriptome analysis revealed that differentially expressed genes associated with specific subpopulations were enriched with stress response-related genes, suggesting that subpopulations have distinct stress response mechanisms. Root growth was negatively correlated with auxin-inducible genes, suggesting an association between auxin signaling and upland field conditions. A negative correlation between crown root diameter and stress response-related genes suggested that thicker crown root diameter is associated with resistance to mild drought stress. Finally, co-expression network analysis implemented with DNA affinity purification followed by sequencing analysis identified phytohormone signaling networks and key transcription factors negatively regulating crown root diameter. Our datasets provide a useful resource for understanding the genomic and transcriptomic basis of phenotypic variation under upland field conditions.

GRAS-Di system facilitates high-density genetic map construction and QTL identification in recombinant inbred lines of the wheat progenitor Aegilops tauschii.

Due to large and complex genomes of Triticeae species, skim sequencing approaches have cost and analytical advantages for detecting genetic markers and building linkage maps. Here, we develop a high-density linkage map and identify quantitative trait loci (QTLs) for recombinant inbred lines of Aegilops tauschii, a D-genome donor of bread wheat, using the recently developed genotyping by Random Amplicon Sequencing-Direct (GRAS-Di) system, which facilitates skimming of the large and complicated genome and generates a large number of genetic markers. The deduced linkage groups based on the GRAS-Di genetic markers corresponded to the chromosome number of Ae. tauschii. We successfully identified stable QTLs for flowering time and spikelet shape-related traits. Genotype differences of RILs at the QTL-linked markers were significantly associated with the trait variations. In particular, one of the QTL-linked markers for flowering time was mapped close to VRN3 (also known as FLOWERING LOCUS T), which controls flowering. The GRAS-Di system is, therefore, an efficient and useful application for genotyping and linkage mapping in species with large and complex genomes, such as Triticeae species.

Low-cost RNA extraction method for highly scalable transcriptome studies.

RNA extraction has been improved by integration of a variety of materials in the protocol, such as phenol, guanidine thiocyanate, and silica, according to the case-specific demands. However, few methods have been designed for high-throughput RNA preparation for large-scale transcriptome studies. In this study, we established a high-throughput guanidinium thiocyanate and isopropyl alcohol based RNA extraction method (HighGI). HighGI is based on simple and phenol-free homemade buffers and the cost is substantially lower than a column-based commercial kit. We demonstrated that the quality and quantity of RNA extracted with HighGI were comparable to those extracted with a conventional phenol/chloroform-based method and a column-based commercial kit. HighGI retained small RNAs less than 200 bp, which are lost with a commercial column-based kit. We also demonstrated that HighGI is readily applicable to semi-automated RNA extraction. HighGI enables high-throughput RNA extraction for large-scale RNA preparation with high yield and quality.

De novo Genome Assembly of the indica Rice Variety IR64 Using Linked-Read Sequencing and Nanopore Sequencing.

IR64 is a rice variety with high-yield that has been widely cultivated around the world. IR64 has been replaced by modern varieties in most growing areas. Given that modern varieties are mostly progenies or relatives of IR64, genetic analysis of IR64 is valuable for rice functional genomics. However, chromosome-level genome sequences of IR64 have not been available previously. Here, we sequenced the IR64 genome using synthetic long reads obtained by linked-read sequencing and ultra-long reads obtained by nanopore sequencing. We integrated these data and generated the de novo assembly of the IR64 genome of 367 Mb, equivalent to 99% of the estimated size. Continuity of the IR64 genome assembly was improved compared with that of a publicly available IR64 genome assembly generated by short reads only. We annotated 41,458 protein-coding genes, including 657 IR64-specific genes, that are missing in other high-quality rice genome assemblies IRGSP-1.0 of japonica cultivar Nipponbare or R498 of indica cultivar Shuhui498. The IR64 genome assembly will serve as a genome resource for rice functional genomics as well as genomics-driven and/or molecular breeding.

Backhoe-assisted monolith method for plant root phenotyping under upland conditions.

Root system architecture (RSA) is one of the most important traits determining water and nutrient availability for plants. Modification of RSA is known to be a useful approach for improving root performance of crops. However, for conducting root phenotyping, there are few alternatives for the rapid collection of root samples from a constant soil volume. In this report, we propose a rapid root-sampling method, which uses a steel cylinder known as round monolith and backhoes to reduce the physical effort. The monolith was set on the ground surrounding individual rice plants and vertically driven back by a backhoe. Soil samples with 20 cm width and 25 cm depth were excavated by the monolith, from which root samples were then isolated. This backhoe-assisted monolith method requires at most five minutes to collect root samples from one plant. Using this method, we quantified the root traits of three rice lines, reported to form different types of root system such as shallow-, intermediate-, and deep-roots, using a root image analysis software. The data obtained through this method, which showed the same trend as previously reported, clearly demonstrated that this method is useful for quantitative evaluation of roots in the soil.

RNA Sequencing-Based Bulked Segregant Analysis Facilitates Efficient D-genome Marker Development for a Specific Chromosomal Region of Synthetic Hexaploid Wheat.

Common wheat originated from interspecific hybridization between cultivated tetraploid wheat and its wild diploid relative Aegilops tauschii followed by amphidiploidization. This evolutionary process can be reproduced artificially, resulting in synthetic hexaploid wheat lines. Here we performed RNA sequencing (RNA-seq)-based bulked segregant analysis (BSA) using a bi-parental mapping population of two synthetic hexaploid wheat lines that shared identical A and B genomes but included with D-genomes of distinct origins. This analysis permitted identification of D-genome-specific polymorphisms around the Net2 gene, a causative locus to hybrid necrosis. The resulting single nucleotide polymorphisms (SNPs) were classified into homoeologous polymorphisms and D-genome allelic variations, based on the RNA-seq results of a parental tetraploid and two Ae. tauschii accessions. The difference in allele frequency at the D-genome-specific SNP sites between the contrasting bulks (ΔSNP-index) was higher on the target chromosome than on the other chromosomes. Several SNPs with the highest ΔSNP-indices were converted into molecular markers and assigned to the Net2 chromosomal region. These results indicated that RNA-seq-based BSA can be applied efficiently to a synthetic hexaploid wheat population to permit molecular marker development in a specific chromosomal region of the D genome.

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

BACKGROUND: Aegilops umbellulata Zhuk. (2n = 14), a wild diploid wheat relative, has been the source of trait improvement in wheat breeding. Intraspecific genetic variation of Ae. umbellulata, however, has not been well studied and the genomic information in this species is limited. RESULTS: To develop novel genetic markers distributed over all chromosomes of Ae. umbellulata and to evaluate its genetic diversity, we performed RNA sequencing of 12 representative accessions and reconstructed transcripts by de novo assembly of reads for each accession. A large number of single nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were obtained and anchored to the pseudomolecules of Ae. tauschii and barley (Hordeum vulgare L.), which were regarded as virtual chromosomes of Ae. umbellulata. Interestingly, genetic diversity in Ae. umbellulata was higher than in Ae. tauschii, despite the narrow habitat of Ae. umbellulata. Comparative analyses of nucleotide polymorphisms between Ae. umbellulata and Ae. tauschii revealed no clear lineage differentiation and existence of alleles with rarer frequencies predominantly in Ae. umbellulata, with patterns clearly distinct from those in Ae. tauschii. CONCLUSIONS: The anchored SNPs, covering all chromosomes, provide sufficient genetic markers between Ae. umbellulata accessions. The alleles with rarer frequencies might be the main source of the high genetic diversity in Ae. umbellulata.

Genetic mapping of a novel recessive allele for non-glaucousness in wild diploid wheat Aegilops tauschii: implications for the evolution of common wheat.

Cuticular wax on the aerial surface of plants has a protective function against many environmental stresses. The bluish-whitish appearance of wheat leaves and stems is called glaucousness. Most modern cultivars of polyploid wheat species exhibit the glaucous phenotype, while in a wild wheat progenitor, Ae. tauschii, both glaucous and non-glaucous accessions exist. Iw2, a wax inhibitor locus on the short arm of chromosome 2D, is the main contributor to this phenotypic variation in Ae. tauschii, and the glaucous/non-glaucous phenotype of Ae. tauschii is usually inherited by synthetic hexaploid wheat. However, a few synthetic lines show the glaucous phenotype although the parental Ae. tauschii accessions are non-glaucous. Molecular marker genotypes indicate that the exceptional non-glaucous Ae. tauschii accessions share the same genotype in the Iw2 chromosomal region as glaucous accessions, suggesting that these accessions have a different causal locus for their phenotype. This locus was assigned to the long arm of chromosome 3D using an F2 mapping population and designated W4, a novel glaucous locus in Ae. tauschii. The dominant W4 allele confers glaucousness, consistent with phenotypic observation of Ae. tauschii accessions and the derived synthetic lines. These results implied that glaucous accessions of Ae. tauschii with the W2W2iw2iw2W4W4 genotype could have been the D-genome donor of common wheat.

Genetic mapping reveals a dominant awn-inhibiting gene related to differentiation of the variety anathera in the wild diploid wheat Aegilops tauschii.

Aegilops tauschii, a wild wheat relative, is the D-genome donor of common wheat. Subspecies and varieties of Ae. tauschii are traditionally classified based on differences in their inflorescence architecture. However, the genetic information for their diversification has been quite limited in the wild wheat relatives. The variety anathera has no awn on the lemma, but the genetic basis for this diagnostic character is unknown. Wide variations in awn length traits at the top and middle spikes were found in the Ae. tauschii core collection, and the awn length at the middle spike was significantly smaller in the eastward-dispersed sublineage than in those in other sublineages. To clarify loci controlling the awnless phenotype of var. anathera, we measured awn length of an intervariety F2 mapping population, and found that the F2 individuals could be divided into two groups mainly based on the awn length at the middle of spike, namely short and long awn groups, significantly fitting a 3:1 segregation ratio, which indicated that a single locus controls the awnless phenotype. The awnless locus, Anathera (Antr), was assigned to the distal region of the short arm of chromosome 5D. Quantitative trait locus analysis using the awn length data of each F2 individual showed that only one major locus was at the same chromosomal position as Antr. These results suggest that a single dominant allele determines the awnless diagnostic character in the variety anathera. The Antr dominant allele is a novel gene inhibiting awn elongation in wheat and its relatives.

Quantitative trait locus analysis for spikelet shape-related traits in wild wheat progenitor Aegilops tauschii: Implications for intraspecific diversification and subspecies differentiation.

Wild diploid wheat Aegilops tauschii, the D-genome progenitor of common wheat, carries large genetic variation in spikelet and grain morphology. Two differentiated subspecies of Ae. tauschii, subspecies tauschii and strangulata, have been traditionally defined based on differences in spikelet morphology. Here, we first assessed six spikelet shape-related traits among 199 Ae. tauschii accessions, and found that the accessions belonging to TauL1major lineage produced significantly longer spikes, higher spikelet density, and shorter, narrower spikelets than another major lineage, TauL2, in which the strangulata accessions are included. Next, we performed quantitative trait locus (QTL) analysis of the spikelet and grain shape using three mapping populations derived from interlineage crosses between TauL1 and TauL2 to identify the genetic loci for the morphological variations of the spikelet and grain shape in Ae. tauschii. Three major QTL regions for the examined traits were detected on chromosomes 3D, 4D and 7D. The 3D and 4D QTL regions for several spikelet shape-related traits were conserved in the three mapping populations, which indicated that the 3D and 4D QTLs contribute to divergence of the two major lineages. The 7D QTLs were found only in a mapping population from a cross of the two subspecies, suggesting that these 7D QTLs may be closely related to subspecies differentiation in Ae. tauschii. Thus, QTL analysis for spikelet and grain morphology may provide useful information to elucidate the evolutionary processes of intraspecific differentiation.

Fine mapping and genetic association analysis of Net2, the causative D-genome locus of low temperature-induced hybrid necrosis in interspecific crosses between tetraploid wheat and Aegilops tauschii.

Hybrid necrosis has been observed in many interspecific hybrids from crosses between tetraploid wheat and the wheat D-genome donor Aegilops tauschii. Type II necrosis is a kind of hybrid incompatibility that is specifically characterized by low-temperature induction and growth suppression. Two complementary genes, Net1 on the AB genome and Net2 on the D genome, putatively control type II necrosis in ABD triploids and synthetic hexaploid wheat. Toward map-based cloning of Net2, a fine map around the Net2 region on 2DS was constructed in this study. Using the draft genome sequence of Ae. tauschii and the physical map of the barley genome, the Net2 locus was mapped within a 0.6 cM interval between two closely linked markers. Although local chromosomal rearrangements were observed in the Net2-corresponding region between the barley/Brachypodium and Ae. tauschii genomes, the two closely linked markers were significantly associated with type II necrosis in Ae. tauschii. These results suggest that these markers will aid efficient selection of Net2 non-carrier individuals from the Ae. tauschii population and intraspecific progeny, and could help with introgression of agriculturally important genes from Ae. tauschii to common wheat.

Genome-wide identification of novel genetic markers from RNA sequencing assembly of diverse Aegilops tauschii accessions.

The wild species in the Triticeae tribe are tremendous resources for crop breeding due to their abundant natural variation. However, their huge and highly repetitive genomes have hindered the establishment of physical maps and the completeness of their genome sequences. To develop molecular markers for the efficient utilization of their valuable traits while avoiding their genome complexity, we assembled RNA sequences of ten representative accessions of Aegilops tauschii, the progenitor of the wheat D genome, and estimated single nucleotide polymorphisms (SNPs) and insertions/deletions (indels). The deduced unigenes were anchored to the chromosomes of Ae. tauschii and barley. The SNPs and indels in the anchored unigenes, covering entire chromosomes, were sufficient for linkage map construction, even in combinations between the genetically closest accessions. Interestingly, the resolution of SNP and indel distribution on barley chromosomes was slightly higher than on Ae. tauschii chromosomes. Since barley chromosomes are regarded as virtual chromosomes of Triticeae species, our strategy allows capture of genetic markers arranged on the chromosomes in order based on the conserved synteny. The resolution of these genetic markers will be comparable to that of the Ae. tauschii whose draft genome sequence is available. Our procedure should be applicable to marker development for Triticeae species, which have no draft sequences available.

Fine mapping of Hch1, the causal D-genome gene for hybrid chlorosis in interspecific crosses between tetraploid wheat and Aegilops tauschii.

Hybrid chlorosis, one of the reproductive barriers between tetraploid wheat and its D-genome progenitor, Aegilops tauschii, inhibits normal growth of synthetic wheat hexaploids. Hybrid chlorosis appears to be due to an epistatic interaction of two loci from the AB and D wheat genomes. Our previous study assigned the causal D-genome gene for hybrid chlorosis, Hch1, to the short arm of chromosome 7D. Here, we constructed a fine map of 7DS near Hch1 using 280 F2 individuals from a cross of two wheat synthetic lines, one showing normal growth and the other showing hybrid chlorosis. The hybrid chlorosis phenotype was controlled by a single dominant allele of the Hch1 locus in the synthetic hexaploids. Hch1 was closely linked to four new markers within 0.2 cM, and may be localized near or within the two Ae. tauschii scaffolds containing the linked markers on 7DS. Comparative analysis of the Hch1 chromosomal region for Ae. tauschii, barley and Brachypodium showed that a local inversion occurred in the region proximal to Hch1 during the divergence between barley and Ae. tauschii, and that the Hch1 region on wheat 7DS is syntenic to Brachypodium chromosome 1. These observations provide useful information for further studies toward map-based cloning of Hch1.

Quantitative trait locus analysis for flowering-related traits using two F2 populations derived from crosses between Japanese common wheat cultivars and synthetic hexaploids.

Flowering time is an important trait for Japanese wheat breeding. Aegilops tauschii, the D-genome donor of hexaploid wheat, is a useful resource to enlarge the D-genome diversity of common wheat. Previously, we identified flowering-related QTLs in F2 populations of synthetic hexaploid wheat lines between the tetraploid wheat cultivar Langdon and Ae. tauschii accessions. Here, to evaluate the usefulness of the early-flowering alleles from Ae. tauschii for Japanese wheat breeding, QTL analyses were conducted in two F2 populations derived from crosses between Japanese wheat cultivars and early-flowering lines of synthetic hexaploid wheat. Only two chromosomal regions controlling flowering-related traits were identified, on chromosomes 2DS and 5AL in the mapping populations, and no previously identified QTLs were found in the synthetic hexaploid lines. The strong effect of the 2DS QTL, putatively corresponding to Ppd-D1, was considered to hide any significant expression of other QTLs with small effects on flowering-related traits. When F2 individuals carrying Ae. tauschii-homozygous alleles around the 2DS QTL region were selected, the Ae. tauschii-derived alleles of the previously identified flowering QTLs partly showed an early-flowering phenotype compared with the Japanese wheat-derived alleles. Thus, some early-flowering alleles from Ae. tauschii may be useful for production of early-flowering Japanese wheat cultivars.

Comparison of gene expression profiles and responses to zinc chloride among inter- and intraspecific hybrids with growth abnormalities in wheat and its relatives.

Hybrid necrosis is a well-known reproductive isolation mechanism in plant species, and an autoimmune response is generally considered to trigger hybrid necrosis through epistatic interaction between disease resistance-related genes in hybrids. In common wheat, the complementary Ne1 and Ne2 genes control hybrid necrosis, defined as type I necrosis. Two other types of hybrid necrosis (type II and type III) have been observed in interspecific hybrids between tetraploid wheat and Aegilops tauschii. Another type of hybrid necrosis, defined here as type IV necrosis, has been reported in F1 hybrids between Triticum urartu and some accessions of Triticum monococcum ssp. aegilopoides. In types I, III and IV, cell death occurs gradually starting in older tissues, whereas type II necrosis symptoms occur only under low temperature. To compare comprehensive gene expression patterns of hybrids showing growth abnormalities, transcriptome analysis of type I and type IV necrosis was performed using a wheat 38k oligo-DNA microarray. Defense-related genes including many WRKY transcription factor genes were dramatically up-regulated in plants showing type I and type IV necrosis, similarly to other known hybrid abnormalities, suggesting an association with an autoimmune response. Reactive oxygen species generation and necrotic cell death were effectively inhibited by ZnCl2 treatment in types I, III and IV necrosis, suggesting a significant association of Ca(2+) influx in upstream signaling of necrotic cell death in wheat hybrid necrosis.

The cuticular wax inhibitor locus Iw2 in wild diploid wheat Aegilops tauschii: phenotypic survey, genetic analysis, and implications for the evolution of common wheat.

BACKGROUND: Cuticular wax production on plant surfaces confers a glaucous appearance and plays important roles in plant stress tolerance. Most common wheat cultivars, which are hexaploid, and most tetraploid wheat cultivars are glaucous; in contrast, a wild wheat progenitor, Aegilops tauschii, can be glaucous or non-glaucous. A dominant non-glaucous allele, Iw2, resides on the short arm of chromosome 2D, which was inherited from Ae. tauschii through polyploidization. Iw2 is one of the major causal genes related to variation in glaucousness among hexaploid wheat. Detailed genetic and phylogeographic knowledge of the Iw2 locus in Ae. tauschii may provide important information and lead to a better understanding of the evolution of common wheat. RESULTS: Glaucous Ae. tauschii accessions were collected from a broad area ranging from Armenia to the southwestern coastal part of the Caspian Sea. Linkage analyses with five mapping populations showed that the glaucous versus non-glaucous difference was mainly controlled by the Iw2 locus in Ae. tauschii. Comparative genomic analysis of barley and Ae. tauschii was then used to develop molecular markers tightly linked with Ae. tauschii Iw2. Chromosomal synteny around the orthologous Iw2 regions indicated that some chromosomal rearrangement had occurred during the genetic divergence leading to Ae. tauschii, barley, and Brachypodium. Genetic associations between specific Iw2-linked markers and respective glaucous phenotypes in Ae. tauschii indicated that at least two non-glaucous accessions might carry other glaucousness-determining loci outside of the Iw2 locus. CONCLUSION: Allelic differences at the Iw2 locus were the main contributors to the phenotypic difference between the glaucous and non-glaucous accessions of Ae. tauschii. Our results supported the previous assumption that the D-genome donor of common wheat could have been any Ae. tauschii variant that carried the recessive iw2 allele.

Genome-wide marker development for the wheat D genome based on single nucleotide polymorphisms identified from transcripts in the wild wheat progenitor Aegilops tauschii.

KEY MESSAGE: 13,347 high-confidence SNPs were discovered through transcriptome sequencing of Aegilops tauschii, which are useful for genomic analysis and molecular breeding of hexaploid wheat. In organisms with large and complex genomes, such as wheat, RNA-seq analysis is cost-effective for discovery of genome-wide single nucleotide polymorphisms (SNPs). In this study, deep sequencing of the spike transcriptome from two Aegilops tauschii accessions representing two major lineages led to the discovery of 13,347 high-confidence (HC) SNPs in 4,872 contigs. After removing redundant SNPs detected in the leaf transcriptome from the same accessions in an earlier study, 10,589 new SNPs were discovered. In total, 5,642 out of 5,808 contigs with HC SNPs were assigned to the Ae. tauschii draft genome sequence. On average, 732 HC polymorphic contigs were mapped in silico to each Ae. tauschii chromosome. Based on the polymorphic data, we developed markers to target the short arm of chromosome 2D and validated the polymorphisms using 20 Ae. tauschii accessions. Of the 29 polymorphic markers, 28 were successfully mapped to 2DS in the diploid F2 population of Ae. tauschii. Among ten hexaploid wheat lines, which included wheat synthetics and common wheat cultivars, 25 of the 43 markers were polymorphic. In the hexaploid F2 population between a common wheat cultivar and a synthetic wheat line, 23 of the 25 polymorphic markers between the parents were available for genotyping of the F2 plants and 22 markers mapped to chromosome 2DS. These results indicate that molecular markers that developed from polymorphisms between two distinct lineages of Ae. tauschii might be useful for analysis not only of the diploid, but also of the hexaploid wheat genome.