Molecular organization and comparative analysis of chromosome 5B of the wild wheat ancestor Triticum dicoccoides

Molecular cytological analysis of alien introgressions in common wheat lines created by crossing of Triticum aestivum with T. dicoccoides and T. dicoccum

Wild and domesticated emmer (ВВАА, 2n = 28) are of significant interest for expanding the genetic diversity of common wheat as sources of a high protein and microelement grain content, resistance to many biotic and abiotic factors. Particular interest in these species is also determined by their close relationship with Triticum aestivum L., which facilitates interspecific hybridization. The objective of this work was to analyze the nature of alien introgressions in hybrid lines from crossing common wheat varieties with T. dicoccoides and T. dicoccum, and to assess the effect of their genome fragments on the cytological stability of introgression lines. A C-banding technique and genotyping with SNP and SSR markers were used to determine localization and length of introgression fragments. Assessment of cytological stability was carried out on the basis of chromosome behavior in microsporogenesis. A molecular cytogenetic analysis of introgression wheat lines indicated that the inclusion of the genetic material of wild and domesticated emmer was carried out mainly in the form of whole arms or large fragments in the chromosomes of the B genome and less extended inserts in the A genome. At the same time, the highest frequency of introgressions of the emmer genome was observed in chromosomes 1A, 1B, 2B, and 3B. The analysis of the final stage of meiosis showed a high level of cytological stability in the vast majority of introgression wheat lines (meiotic index was 83.0-99.0 %), which ensures the formation of functional gametes in an amount sufficient for successful reproduction. These lines are of interest for the selection of promising material with agronomically valuable traits and their subsequent inclusion in the breeding process.

Flow karyotyping of wheat-Aegilops additions facilitate dissecting the genomes of Ae. biuncialis and Ae. geniculata into individual chromosomes

Breeding of wheat adapted to new climatic conditions and resistant to diseases and pests is hindered by a limited gene pool due to domestication and thousands of years of human selection. Annual goatgrasses (Aegilops spp.) with M and U genomes are potential sources of the missing genes and alleles. Development of alien introgression lines of wheat may be facilitated by the knowledge of DNA sequences of Aegilops chromosomes. As the Aegilops genomes are complex, sequencing relevant Aegilops chromosomes purified by flow cytometric sorting offers an attractive route forward. The present study extends the potential of chromosome genomics to allotetraploid Ae. biuncialis and Ae. geniculata by dissecting their M and U genomes into individual chromosomes. Hybridization of FITC-conjugated GAA oligonucleotide probe to chromosomes suspensions of the two species allowed the application of bivariate flow karyotyping and sorting some individual chromosomes. Bivariate flow karyotype FITC vs. DAPI of Ae. biuncialis consisted of nine chromosome-populations, but their chromosome content determined by microscopic analysis of flow sorted chromosomes indicated that only 7M^b and 1U^b could be sorted at high purity. In the case of Ae. geniculata, fourteen chromosome-populations were discriminated, allowing the separation of nine individual chromosomes (1M^g, 3M^g, 5M^g, 6M^g, 7M^g, 1U^g, 3U^g, 6U^g, and 7U^g) out of the 14. To sort the remaining chromosomes, a partial set of wheat-Ae. biuncialis and a whole set of wheat-Ae. geniculata chromosome addition lines were also flow karyotyped, revealing clear separation of the GAA-rich Aegilops chromosomes from the GAA-poor A- and D-genome chromosomes of wheat. All of the alien chromosomes represented by individual addition lines could be isolated at purities ranging from 74.5% to 96.6% and from 87.8% to 97.7%, respectively. Differences in flow karyotypes between Ae. biuncialis and Ae. geniculata were analyzed and discussed. Chromosome-specific genomic resources will facilitate gene cloning and the development of molecular tools to support alien introgression breeding of wheat.

Fine structure and transcription dynamics of bread wheat ribosomal DNA loci deciphered by a multi-omics approach

Three out of four RNA components of ribosomes are encoded by 45S ribosomal DNA (rDNA) loci, which are organized as long head-to-tail tandem arrays of nearly identical units, spanning several megabases of sequence. Due to this structure, the rDNA loci are the major sources of gaps in genome assemblies, and gene copy number, sequence composition, and expression status of particular arrays remain elusive, especially in complex genomes harboring multiple loci. Here we conducted a multi-omics study to decipher the 45S rDNA loci in hexaploid bread wheat. Coupling chromosomal genomics with optical mapping, we reconstructed individual rDNA arrays, enabling locus-specific analyses of transcription activity and methylation status from RNA- and bisulfite-sequencing data. We estimated a total of 6,650 rDNA units in the bread wheat genome, with approximately 2,321, 3,910, 253, and 50 gene copies located in short arms of chromosomes 1B, 6B, 5D, and 1A, respectively. Only 1B and 6B loci contributed substantially to rRNA transcription at a roughly 2:1 ratio. The ratio varied among five tissues analyzed (embryo, coleoptile, root tip, primary leaf, mature leaf), being the highest (2.64:1) in mature leaf and lowest (1.72:1) in coleoptile. Cytosine methylation was considerably higher in CHG context in the silenced 5D locus as compared with the active 1B and 6B loci. In conclusion, a fine genomic organization and tissue-specific expression of rDNA loci were deciphered, for the first time, in a complex polyploid species. The results are discussed in the context of wheat evolution and transcription regulation.

Pan-Genome miRNomics in Brachypodium

Pan-genomes are efficient tools for the identification of conserved and varying genomic sequences within lineages of a species. Investigating genetic variations might lead to the discovery of genes present in a subset of lineages, which might contribute into beneficial agronomic traits such as stress resistance or yield. The content of varying genomic regions in the pan-genome could include protein-coding genes as well as microRNA(miRNAs), small non-coding RNAs playing key roles in the regulation of gene expression. In this study, we performed in silico miRNA identification from the genomic sequences of 54 lineages of Brachypodium distachyon, aiming to explore varying miRNA contents and their functional interactions. A total of 115 miRNA families were identified in 54 lineages, 56 of which were found to be present in all lineages. The miRNA families were classified based on their conservation among lineages and potential mRNA targets were identified. Obtaining information about regulatory mechanisms stemming from these miRNAs offers strong potential to provide a better insight into the complex traits that were potentially present in some lineages. Future work could lead us to introduce these traits to different lineages or other economically important plant species in order to promote their survival in different environmental conditions.

Chromosome analysis and sorting

Flow cytometric analysis and sorting of plant mitotic chromosomes has been mastered by only a few laboratories worldwide. Yet, it has been contributing significantly to progress in plant genetics, including the production of genome assemblies and the cloning of important genes. The dissection of complex genomes by flow sorting into the individual chromosomes that represent small parts of the genome reduces DNA sample complexity and streamlines projects relying on molecular and genomic techniques. Whereas flow cytometric analysis, that is, chromosome classification according to fluorescence and light scatter properties, is an integral part of any chromosome sorting project, it has rarely been used on its own due to lower resolution and sensitivity as compared to other cytogenetic methods. To perform chromosome analysis and sorting, commercially available electrostatic droplet sorters are suitable. However, in order to resolve and purify chromosomes of interest the instrument must offer high resolution of optical signals as well as stability during long runs. The challenge is thus not the instrumentation, but the adequate sample preparation. The sample must be a suspension of intact mitotic metaphase chromosomes and the protocol, which includes the induction of cell cycle synchrony, accumulation of dividing cells at metaphase, and release of undamaged chromosomes, is time consuming and laborious and needs to be performed very carefully. Moreover, in addition to fluorescent staining chromosomal DNA, the protocol may include specific labelling of DNA repeats to facilitate discrimination of particular chromosomes. This review introduces the applications of chromosome sorting in plants, and discusses in detail sample preparation, chromosome analysis and sorting to achieve the highest purity in flow-sorted fractions, and their suitability for downstream applications.

Capture of complete ciliate chromosomes in single sequencing reads reveals widespread chromosome isoforms

BACKGROUND

Whole-genome shotgun sequencing, which stitches together millions of short sequencing reads into a single genome, ushered in the era of modern genomics and led to a rapid expansion of the number of genome sequences available. Nevertheless, assembly of short reads remains difficult, resulting in fragmented genome sequences. Ultimately, only a sequencing technology capable of capturing complete chromosomes in a single run could resolve all ambiguities. Even "third generation" sequencing technologies produce reads far shorter than most eukaryotic chromosomes. However, the ciliate Oxytricha trifallax has a somatic genome with thousands of chromosomes averaging only 3.2 kbp, making it an ideal candidate for exploring the benefits of sequencing whole chromosomes without assembly.

RESULTS

We used single-molecule real-time sequencing to capture thousands of complete chromosomes in single reads and to update the published Oxytricha trifallax JRB310 genome assembly. In this version, over 50% of the completed chromosomes with two telomeres derive from single reads. The improved assembly includes over 12,000 new chromosome isoforms, and demonstrates that somatic chromosomes derive from variable rearrangements between somatic segments encoded up to 191,000 base pairs away. However, while long reads reduce the need for assembly, a hybrid approach that supplements long-read sequencing with short reads for error correction produced the most complete and accurate assembly, overall.

CONCLUSIONS

This assembly provides the first example of complete eukaryotic chromosomes captured by single sequencing reads and demonstrates that traditional approaches to genome assembly can mask considerable structural variation.

Dissecting the Complex Genome of Crested Wheatgrass by Chromosome Flow Sorting

Wheatgrass (Agropyron sp.) is a potential source of beneficial traits for wheat improvement. Among them, crested wheatgrass [A. cristatum (L.) Gaertn.] comprises a complex of diploid, tetraploid, and hexaploid forms with the basic genome P, with some accessions carrying supernumerary B chromosomes (Bs). In this work, we applied flow cytometry to dissect the complex genome of crested wheatgrass into individual chromosomes to facilitate its analysis. Flow karyotypes obtained after the analysis of 4',6-diamidino-2-phenylindole (DAPI)-stained mitotic chromosomes of diploid and tetraploid accessions consisted of three peaks, each corresponding to a group of two or three chromosomes. To improve the resolution, bivariate flow karyotyping after fluorescent labeling of chromosomes with fluorescein isothiocyanate (FITC)-conjugated probe (GAA) microsatellite was applied and allowed discrimination and sorting of P genome chromosomes from wheat-crested wheatgrass addition lines. Chromosomes 1P-6P and seven telomeric chromosomes could be sorted at purities ranging from 81.7 to 98.2% in disomics and from 44.8 to 87.3% in telosomics. Chromosome 7P was sorted at purities reaching 50.0 and 39.5% in diploid and tetraploid crested wheatgrass, respectively. In addition to the whole complement chromosomes (A), Bs could be easily discriminated and sorted from a diploid accession at 95.4% purity. The sorted chromosomes will streamline genome analysis of crested wheatgrass, facilitating gene cloning and development of molecular tools to support alien introgression into wheat.

Differential Expression of MicroRNAs During Root Formation in Taxus Chinensis Var. mairei Cultivars

MicroRNAs (miRNAs) have been shown to play key roles in the regulation of plant growth and development by modifying the expression of their target genes. However, the influence of miRNAs on root formation and development in woody plants, such as Taxus chinensis, remains largely unknown. In the current study, we explored the phytohormone-response and nutrition-response miRNA expression profiles during T. chinensis rooting by quantitative real-time PCR (qPCR). We identified six phytohormone-response miRNAs, namely, miR164a, miR165, miR167a, miR171b, miR319, and miR391, and eight nutrition-response miRNAs, namely, miR169b, miR395a, miR399c, miR408, miR826, miR827, miR857, and miR2111a, that were differentially expressed at different rooting phases of T. chinensis. Using northern blot analysis of the putative target genes of these miRNAs, we detected the relative gene expression changes of the target genes. Taken together, our results suggest that miRNAs are involved in root formation of T. chinensis and that miRNAs may play important regulatory roles in primary root, crown root, and root hair formation by targeting phytohormone and/or nutrition response genes in T. chinensis. For the first time, these results expand our understanding of the molecular mechanisms of plant root formation and development in a conifer species.

Salinity-associated microRNAs and their potential roles in mediating salt tolerance in rice colonized by the endophytic root fungus Piriformospora indica

Piriformospora indica (P. indica), an endophytic root fungus, supports the growth and enhanced tolerance of plants to biotic and abiotic stresses. Several recent studies showed the significant role of small RNA (sRNA) molecules including microRNAs (miRNAs) in plant adaption to environmental stress, but little is known concerning the symbiosis-mediated salt stress tolerance regulated at miRNAs level. The overarching goal of this research is to elucidate the impact of miRNAs in regulating the P. indica-mediated salt tolerance in rice. Applying sRNA-seq analysis led to identify a set of 547 differentially abundant miRNAs in response to P. indica inoculation and salt stress. These included 206 rice-specific and 341 previously known miRNAs from other plant species. In silico analysis of miRNAs predictions of the differentially abundant miRNAs led to identifying of 193 putatively target genes, most of which were encoded either genes or transcription factors involved in nutrient uptake, sodium ion transporters, growth regulators, and auxin- responsive proteins. The rice-specific miRNAs targeted the transcription factors involved in the import of potassium ions into the root cells, the export of sodium ions, and plant growth and development. Interestingly, P. indica affected the differential abundance of miRNAs regulated genes and transcription factors linked to salt stress tolerance. Our data helps to understand the molecular basis of salt stress tolerance mediated by symbionts in plant and the potential impact of miRNAs for genetic improvement of rice varieties for tolerance to salt stress.

Chromosome-based survey sequencing reveals the genome organization of wild wheat progenitor Triticum dicoccoides

Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is the progenitor of wheat. We performed chromosome-based survey sequencing of the 14 chromosomes, examining repetitive sequences, protein-coding genes, miRNA/target pairs and tRNA genes, as well as syntenic relationships with related grasses. We found considerable differences in the content and distribution of repetitive sequences between the A and B subgenomes. The gene contents of individual chromosomes varied widely, not necessarily correlating with chromosome size. We catalogued candidate agronomically important loci, along with new alleles and flanking sequences that can be used to design exome sequencing. Syntenic relationships and virtual gene orders revealed several small-scale evolutionary rearrangements, in addition to providing evidence for the 4AL-5AL-7BS translocation in wild emmer wheat. Chromosome-based sequence assemblies contained five novel miRNA families, among 59 families putatively encoded in the entire genome which provide insight into the domestication of wheat and an overview of the genome content and organization.

Genetic Contribution of Emmer Wheat (Triticum dicoccon Schrank) to Heat Tolerance of Bread Wheat

Rising global temperatures cause substantial yield losses in many wheat growing environments. Emmer wheat (Triticum dicoccon Schrank), one of the first wheat species domesticated, carries significant variation for tolerance to abiotic stresses. This study identified new genetic variability for high-temperature tolerance in hexaploid progeny derived from crosses with emmer wheat. Eight hexaploid and 11 tetraploid parents were recombined in 43 backcross combinations using the hexaploid as the recurrent parent. A total of 537 emmer-based hexaploid lines were developed by producing approximately 10 doubled haploids on hexaploid like BC1F1 progeny and subsequent selection for hexaploid morphology. These materials and 17 commercial cultivars and hexaploid recurrent parents were evaluated under two times of sowing in the field, in 2014-2016. The materials were genotyped using a 90K SNP platform and these data were used to estimate the contribution of emmer wheat to the progeny. Significant phenotypic and genetic variation for key agronomical traits including grain yield, TKW and screenings was observed. Many of the emmer derived lines showed improved performance under heat stress (delayed sowing) compared with parents and commercial cultivars. Emmer derived lines were the highest yielding material in both sowing dates. The emmer wheat parent contributed between 1 and 44% of the genome of the derived lines. Emmer derived lines with superior kernel weight and yield generally had a greater genetic contribution from the emmer parent compared to those with lower trait values. The study showed that new genetic variation for key traits such as yield, kernel weight and screenings can be introduced to hexaploid wheat from emmer wheat. These genetic resources should be explored more systematically to stabilize grain yield and quality in a changing climate.

Hotspots in the genomic architecture of field drought responses in wheat as breeding targets

Wheat can adapt to most agricultural conditions across temperate regions. This success is the result of phenotypic plasticity conferred by a large and complex genome composed of three homoeologous genomes (A, B, and D). Although drought is a major cause of yield and quality loss in wheat, the adaptive mechanisms and gene networks underlying drought responses in the field remain largely unknown. Here, we addressed this by utilizing an interdisciplinary approach involving field water status phenotyping, sampling, and gene expression analyses. Overall, changes at the transcriptional level were reflected in plant spectral traits amenable to field-level physiological measurements, although changes in photosynthesis-related pathways were found likely to be under more complex post-transcriptional control. Examining homoeologous genes with a 1:1:1 relationship across the A, B, and D genomes (triads), we revealed a complex genomic architecture for drought responses under field conditions, involving gene homoeolog specialization, multiple gene clusters, gene families, miRNAs, and transcription factors coordinating these responses. Our results provide a new focus for genomics-assisted breeding of drought-tolerant wheat cultivars.

Complete chloroplast genome sequence and comparative analysis of loblolly pine (Pinus taeda L.) with related species

Pinaceae, the largest family of conifers, has a diversified organization of chloroplast (cp) genomes with two typical highly reduced inverted repeats (IRs). In the current study, we determined the complete sequence of the cp genome of an economically and ecologically important conifer tree, the loblolly pine (Pinus taeda L.), using Illumina paired-end sequencing and compared the sequence with those of other pine species. The results revealed a genome size of 121,531 base pairs (bp) containing a pair of 830-bp IR regions, distinguished by a small single copy (42,258 bp) and large single copy (77,614 bp) region. The chloroplast genome of P. taeda encodes 120 genes, comprising 81 protein-coding genes, four ribosomal RNA genes, and 35 tRNA genes, with 151 randomly distributed microsatellites. Approximately 6 palindromic, 34 forward, and 22 tandem repeats were found in the P. taeda cp genome. Whole cp genome comparison with those of other Pinus species exhibited an overall high degree of sequence similarity, with some divergence in intergenic spacers. Higher and lower numbers of indels and single-nucleotide polymorphism substitutions were observed relative to P. contorta and P. monophylla, respectively. Phylogenomic analyses based on the complete genome sequence revealed that 60 shared genes generated trees with the same topologies, and P. taeda was closely related to P. contorta in the subgenus Pinus. Thus, the complete P. taeda genome provided valuable resources for population and evolutionary studies of gymnosperms and can be used to identify related species.

A large-scale multiomics analysis of wheat stem solidness and the wheat stem sawfly feeding response, and syntenic associations in barley, Brachypodium, and rice

The wheat stem sawfly (WSS), Cephus cinctus Norton (Hymenoptera: Cephidae), is an important pest of wheat and other cereals, threatening the quality and quantity of grain production. WSS larvae feed and develop inside the stem where they are protected from the external environment; therefore, pest management strategies primarily rely on host plant resistance. A major locus on the long arm of wheat chromosome 3B underlies most of the variation in stem solidness; however, the impact of stem solidness on WSS feeding has not been completely characterized. Here, we used a multiomics approach to examine the response to WSS in both solid- and semi-solid-stemmed wheat varieties. The combined transcriptomic, proteomic, and metabolomic data revealed that two important molecular pathways, phenylpropanoid and phosphate pentose, are involved in plant defense against WSS. We also detected a general downregulation of several key defense transcripts, including those encoding secondary metabolites such as DIMBOA, tricetin, and lignin, which suggested that the WSS larva might interfere with plant defense. We comparatively analyzed the stem solidness genomic region known to be associated with WSS tolerance in wild emmer, durum, and bread wheats, and described syntenic regions in the close relatives barley, Brachypodium, and rice. Additionally, microRNAs identified from the same genomic region revealed potential regulatory pathways associated with the WSS response. We propose a model outlining the molecular responses of the WSS–wheat interactions. These findings provide insight into the link between stem solidness and WSS feeding at the molecular level.

Mapping QTLs conferring salt tolerance and micronutrient concentrations at seedling stagein wheat

Soil salinization and degradation is one of the consequences of climate change. Identification of major salt tolerance genes and marker assisted selection (MAS) can accelerate wheat breeding for this trait. We genotyped 154 wheat F₂ lines derived from a cross between salt tolerant and susceptible cultivars using the Axiom Wheat Breeder's Genotyping Array. A high-density linkage map of 988 single nucleotide polymorphisms (SNPs) was constructed and utilized for quantitative trait loci (QTL) mapping for salt tolerance traits and mineral concentrations under salinity. Of 49 mapped QTLs, six were for Na⁺ exclusion (NAX) and two QTLs (qSNAX.2 A.1, qSNAX.2 A.2) on chromosome 2 A coincided with a reported major NAX QTL (Nax1 or HKT1;4). Two other major NAX QTLs were mapped on 7 A, which contributed 11.23 and 18.79% of the salt tolerance respectively. In addition to Ca⁺² and Mg⁺² QTLs, twenty-seven QTLs for tissue Phosphorus, Zinc, Iron, Manganese, Copper, Sulphur and Boron concentrations under salinity were also mapped. The 1293 segregating SNPs were annotated/located within genes for various ion channels, signalling pathways, transcription factors (TFs), metabolic pathways and 258 of them showed differential expression in silico under salinity. These findings will create new opportunities for salt tolerance breeding programs.

Mapping QTLs conferring salt tolerance and micronutrient concentrations at seedling stagein wheat

Soil salinization and degradation is one of the consequences of climate change. Identification of major salt tolerance genes and marker assisted selection (MAS) can accelerate wheat breeding for this trait. We genotyped 154 wheat F₂ lines derived from a cross between salt tolerant and susceptible cultivars using the Axiom Wheat Breeder’s Genotyping Array. A high-density linkage map of 988 single nucleotide polymorphisms (SNPs) was constructed and utilized for quantitative trait loci (QTL) mapping for salt tolerance traits and mineral concentrations under salinity. Of 49 mapped QTLs, six were for Na⁺ exclusion (NAX) and two QTLs (qSNAX.2 A.1, qSNAX.2 A.2) on chromosome 2 A coincided with a reported major NAX QTL (Nax1 or HKT1;4). Two other major NAX QTLs were mapped on 7 A, which contributed 11.23 and 18.79% of the salt tolerance respectively. In addition to Ca⁺² and Mg⁺² QTLs, twenty-seven QTLs for tissue Phosphorus, Zinc, Iron, Manganese, Copper, Sulphur and Boron concentrations under salinity were also mapped. The 1293 segregating SNPs were annotated/located within genes for various ion channels, signalling pathways, transcription factors (TFs), metabolic pathways and 258 of them showed differential expression in silico under salinity. These findings will create new opportunities for salt tolerance breeding programs.

High-throughput SNP genotyping of modern and wild emmer wheat for yield and root morphology using a combined association and linkage analysis

Durum wheat (Triticum turgidum var. durum Desf.) is a major world crop that is grown primarily in areas of the world that experience periodic drought, and therefore, breeding climate-resilient durum wheat is a priority. High-throughput single nucleotide polymorphism (SNP) genotyping techniques have greatly increased the power of linkage and association mapping analyses for bread wheat, but as yet there is no durum wheat-specific platform available. In this study, we evaluate the new 384HT Wheat Breeders Array for its usefulness in tetraploid wheat breeding by genotyping a breeding population of F₆ hybrids, derived from multiple crosses between T. durum cultivars and wild and cultivated emmer wheat accessions. Using a combined linkage and association mapping approach, we generated a genetic map including 1345 SNP markers, and identified markers linked to 6 QTLs for coleoptile length (2), heading date (1), anthocyanin accumulation (1) and osmotic stress tolerance (2). We also developed a straightforward approach for combining genetic data from multiple families of limited size that will be useful for evaluating and mapping pre-existing breeding material.

DNApod: DNA polymorphism annotation database from next-generation sequence read archives

With the rapid advances in next-generation sequencing (NGS), datasets for DNA polymorphisms among various species and strains have been produced, stored, and distributed. However, reliability varies among these datasets because the experimental and analytical conditions used differ among assays. Furthermore, such datasets have been frequently distributed from the websites of individual sequencing projects. It is desirable to integrate DNA polymorphism data into one database featuring uniform quality control that is distributed from a single platform at a single place. DNA polymorphism annotation database (DNApod; http://tga.nig.ac.jp/dnapod/) is an integrated database that stores genome-wide DNA polymorphism datasets acquired under uniform analytical conditions, and this includes uniformity in the quality of the raw data, the reference genome version, and evaluation algorithms. DNApod genotypic data are re-analyzed whole-genome shotgun datasets extracted from sequence read archives, and DNApod distributes genome-wide DNA polymorphism datasets and known-gene annotations for each DNA polymorphism. This new database was developed for storing genome-wide DNA polymorphism datasets of plants, with crops being the first priority. Here, we describe our analyzed data for 679, 404, and 66 strains of rice, maize, and sorghum, respectively. The analytical methods are available as a DNApod workflow in an NGS annotation system of the DNA Data Bank of Japan and a virtual machine image. Furthermore, DNApod provides tables of links of identifiers between DNApod genotypic data and public phenotypic data. To advance the sharing of organism knowledge, DNApod offers basic and ubiquitous functions for multiple alignment and phylogenetic tree construction by using orthologous gene information.

Transcriptome-based SNP discovery by GBS and the construction of a genetic map for olive

Molecular markers located in the genic regions of plants are valuable tools for the identification of candidate genes of economically important traits and consequent use in marker-assisted selection (MAS). In the past, simple sequence repeat markers (SSRs) and single-nucleotide polymorphisms (SNPs) located in expressed sequence tags (ESTs) were developed by sequencing RNA derived from different plant tissues, which involves laborious RNA extraction, mRNA isolation, and cDNA synthesis. In order to develop SNP markers located in olive transcriptomes, we used the recently developed genotyping-by-sequencing (GBS) technique. An analysis was done for 125 olive DNA samples (123 DNA samples from a cross-pollinated F₁ mapping population, and two samples from parents). From 45 to 66% of Illumina reads from GBS analysis were aligned to the olive transcriptome. A total of 22,033 transcriptome-based SNP markers were identified, and 3384 of these were mapped in the olive genome. The genetic linkage map constructed in this study consists of 1 cleaved amplified polymorphic sequence (CAPS), 19 SSR, and 3384 transcriptome-based SNP markers. The map covers 3340.8 cM of the olive genome in 23 linkage groups, with the length of the linkage groups ranging from 55.6 to 248.7 cM. Average map distance between flanking markers was 0.98 cM. This genetic linkage map is a saturated genetic map and will be a useful tool for the localization of quantitative trait loci (QTLs) and gene(s) of interest and for the identification of candidate genes for economically important traits.

A large-scale chromosome-specific SNP discovery guideline

Single-nucleotide polymorphisms (SNPs) are the most prevalent type of variation in genomes that are increasingly being used as molecular markers in diversity analyses, mapping and cloning of genes, and germplasm characterization. However, only a few studies reported large-scale SNP discovery in Aegilops tauschii, restricting their potential use as markers for the low-polymorphic D genome. Here, we report 68,592 SNPs found on the gene-related sequences of the 5D chromosome of Ae. tauschii genotype MvGB589 using genomic and transcriptomic sequences from seven Ae. tauschii accessions, including AL8/78, the only genotype for which a draft genome sequence is available at present. We also suggest a workflow to compare SNP positions in homologous regions on the 5D chromosome of Triticum aestivum, bread wheat, to mark single nucleotide variations between these closely related species. Overall, the identified SNPs define a density of 4.49 SNPs per kilobyte, among the highest reported for the genic regions of Ae. tauschii so far. To our knowledge, this study also presents the first chromosome-specific SNP catalog in Ae. tauschii that should facilitate the association of these SNPs with morphological traits on chromosome 5D to be ultimately targeted for wheat improvement.

Flow Analysis and Sorting of Plant Chromosomes

Analysis and sorting of plant chromosomes (plant flow cytogenetics) is a special application of flow cytometry in plant genomics and its success depends critically on sample quality. This unit describes the methodology in a stepwise manner, starting with the induction of cell cycle synchrony and accumulation of dividing cells in mitotic metaphase, and continues with the preparation of suspensions of intact mitotic chromosomes, flow analysis and sorting of chromosomes, and finally processing of the sorted chromosomes. Each step of the protocol is described in detail as some procedures have not been used widely. Supporting histograms are presented as well as hints on dealing with plant material; the utility of sorted chromosomes for plant genomics is also discussed. © 2016 by John Wiley & Sons, Inc.

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

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

Dissecting miRNAs in Wheat D Genome Progenitor, Aegilops tauschii

As the post-transcriptional regulators of gene expression, microRNAs or miRNAs comprise an integral part of understanding how genomes function. Although miRNAs have been a major focus of recent efforts, miRNA research is still in its infancy in most plant species. Aegilops tauschii, the D genome progenitor of bread wheat, is a wild diploid grass exhibiting remarkable population diversity. Due to the direct ancestry and the diverse gene pool, A. tauschii is a promising source for bread wheat improvement. In this study, a total of 87 Aegilops miRNA families, including 51 previously unknown, were computationally identified both at the subgenomic level, using flow-sorted A. tauschii 5D chromosome, and at the whole genome level. Predictions at the genomic and subgenomic levels suggested A. tauschii 5D chromosome as rich in pre-miRNAs that are highly associated with Class II DNA transposons. In order to gain insights into miRNA evolution, putative 5D chromosome miRNAs were compared to its modern ortholog, Triticum aestivum 5D chromosome, revealing that 48 of the 58 A. tauschii 5D miRNAs were conserved in orthologous T. aestivum 5D chromosome. The expression profiles of selected miRNAs (miR167, miR5205, miR5175, miR5523) provided the first experimental evidence for miR5175, miR5205 and miR5523, and revealed differential expressional changes in response to drought in different genetic backgrounds for miR167 and miR5175. Interestingly, while miR5523 coding regions were present and expressed as pre-miR5523 in both T. aestivum and A. tauschii, the expression of mature miR5523 was observed only in A. tauschii under normal conditions, pointing out to an interference at the downstream processing of pre-miR5523 in T. aestivum. Overall, this study expands our knowledge on the miRNA catalog of A. tauschii, locating a subset specifically to the 5D chromosome, with ample functional and comparative insight which should contribute to and complement efforts to develop drought tolerant wheat varieties.

Transcriptome profiling of wheat glumes in wild emmer, hulled landraces and modern cultivars

Background

Wheat domestication is considered as one of the most important events in the development of human civilization. Wheat spikelets have undergone significant changes during evolution under domestication, resulting in soft glumes and larger kernels that are released easily upon threshing. Our main goal was to explore changes in transcriptome expression in glumes that accompanied wheat evolution under domestication.

Methods

A total of six tetraploid wheat accessions were selected for transcriptome profiling based on their rachis brittleness and glumes toughness. RNA pools from glumes of the central spikelet at heading time were used to construct cDNA libraries for sequencing. The trimmed reads from each library were separately aligned to the reference sub-genomes A and B, which were extracted from wheat survey sequence. Differentially expression analysis and functional annotation were performed between wild and domesticated wheat, to identity candidate genes associated with evolution under domestication. Selected candidate genes were validated using real time PCR.

Results

Transcriptome profiles of wild emmer wheat, wheat landraces, and wheat cultivars were compared using next generation sequencing (RNA-seq). We have found a total of 194,893 transcripts, of which 73,150 were shared between wild, landraces, and cultivars. From 781 differentially expressed genes (DEGs), 336 were down-regulated and 445 were up-regulated in the domesticated compared to wild wheat genotypes. Gene Ontology (GO) annotation assigned 293 DEGs (37.5 %) to GO term groups, of which 134 (17.1 %) were down-regulated and 159 (20.4 %) up-regulated in the domesticated wheat. Some of the down-regulated DEGs in domesticated wheat are related to the biosynthetic pathways that eventually define the mechanical strength of the glumes, such as cell wall, lignin, pectin and wax biosynthesis. The reduction in gene expression of such genes, may explain the softness of the glumes in the domesticated forms. In addition, we have identified genes involved in nutrient remobilization that may affect grain size and other agronomic traits evolved under domestication.

Conclusions

The comparison of RNA-seq profiles between glumes of wheat groups differing in glumes toughness and rachis brittleness revealed a few DEGs that may be involved in glumes toughness and nutrient remobilization. These genes may be involved in processes of wheat improvement under domestication.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1996-0) contains supplementary material, which is available to authorized users.