Separating homeologs by phasing in the tetraploid wheat transcriptome

Plant Tolerance to Drought Stress with Emphasis on Wheat

Environmental stresses, such as drought, have negative effects on crop yield. Drought is a stress whose impact tends to increase in some critical regions. However, the worldwide population is continuously increasing and climate change may affect its food supply in the upcoming years. Therefore, there is an ongoing effort to understand the molecular processes that may contribute to improving drought tolerance of strategic crops. These investigations should contribute to delivering drought-tolerant cultivars by selective breeding. For this reason, it is worthwhile to review regularly the literature concerning the molecular mechanisms and technologies that could facilitate gene pyramiding for drought tolerance. This review summarizes achievements obtained using QTL mapping, genomics, synteny, epigenetics, and transgenics for the selective breeding of drought-tolerant wheat cultivars. Synthetic apomixis combined with the msh1 mutation opens the way to induce and stabilize epigenomes in crops, which offers the potential of accelerating selective breeding for drought tolerance in arid and semi-arid regions.

A Robust Methodology for Assessing Homoeolog-Specific Expression

Angiosperm evolution is marked by numerous, recurring polyploidization events. While hybridization and polyploidization have greatly increased the degree of genetic and phenotypic diversity in plants, the mechanisms underlying changes in the genotype-to-phenotype relationships remain unclear. As the field of natural sciences continues to expand during the post-genomic era, large datasets are becoming increasingly common. However, the development of tools and workflows available to robustly assess these changes have lagged behind data production. A robust homoeolog-specific expression analysis strongly depends upon proper homoeolog calling, the ability to account for reference sequence biases, flexible and accurate methods for dealing with residual bias, and a reproducible workflow. To that end, this chapter aims to provide a detailed description of the potential pitfalls encountered while estimating homoeolog-specific expression as well as provide a workflow that allows for robust inferences based on precise estimates of expression changes.

Molecular evidence for adaptive evolution of drought tolerance in wild cereals

The considerable drought tolerance of wild cereal crop progenitors has diminished during domestication in the pursuit of higher productivity. Regaining this trait in cereal crops is essential for global food security but requires novel genetic insight. Here, we assessed the molecular evidence for natural variation of drought tolerance in wild barley (Hordeum spontaneum), wild emmer wheat (Triticum dicoccoides), and Brachypodium species collected from dry and moist habitats at Evolution Canyon, Israel (ECI). We report that prevailing moist vs dry conditions have differentially shaped the stomatal and photosynthetic traits of these wild cereals in their respective habitats. We present the genomic and transcriptomic evidence accounting for differences, including co-expression gene modules, correlated with physiological traits, and selective sweeps, driven by the xeric site conditions on the African Slope (AS) at ECI. Co-expression gene module 'circadian rhythm' was linked to significant drought-induced delay in flowering time in Brachypodium stacei genotypes. African Slope-specific differentially expressed genes are important in barley drought tolerance, verified by silencing Disease-Related Nonspecific Lipid Transfer 1 (DRN1), Nonphotochemical Quenching 4 (NPQ4), and Brassinosteroid-Responsive Ring-H1 (BRH1). Our results provide new genetic information for the breeding of resilient wheat and barley in a changing global climate with increasingly frequent drought events.

Equilibrated evolution of the mixed auto-/allopolyploid haplotype-resolved genome of the invasive hexaploid Prussian carp

Understanding genome evolution of polyploids requires dissection of their often highly similar subgenomes and haplotypes. Polyploid animal genome assemblies so far restricted homologous chromosomes to a ‘collapsed’ representation. Here, we sequenced the genome of the asexual Prussian carp, which is a close relative of the goldfish, and present a haplotype-resolved chromosome-scale assembly of a hexaploid animal. Genome-wide comparisons of the 150 chromosomes with those of two ancestral diploid cyprinids and the allotetraploid goldfish and common carp revealed the genomic structure, phylogeny and genome duplication history of its genome. It consists of 25 syntenic, homeologous chromosome groups and evolved by a recent autoploid addition to an allotetraploid ancestor. We show that de-polyploidization of the alloploid subgenomes on the individual gene level occurred in an equilibrated fashion. Analysis of the highly conserved actinopterygian gene set uncovered a subgenome dominance in duplicate gene loss of one ancestral chromosome set.

The haplotype-resolved assembly of the asexual invasive Prussian carp shows six genome copies (AAABBB), evolved from two ancestral species by a recent self-addition (AB) to its hybrid-tetraploid (AABB) goldfish ancestor. Equilibrated gene loss led to subgenome dominance.

Expression Level Dominance and Homeolog Expression Bias Upon Cold Stress in the F1 Hybrid Between the Invasive Sphagneticola trilobata and the Native S. calendulacea in South China, and Implications for Its Invasiveness

The role of hybridization is significant in biological invasion, and thermotolerance is a trait critical to range expansions. The South American Sphagneticola trilobata is now widespread in South China, threatening the native S. calendulacea by competition and hybridization. Furthermore, upon formation, their F1 hybrid can quickly replace both parents. In this study, the three taxa were used as a model to investigate the consequences of hybridization on cold tolerance, particularly the effect of subgenome dominance in the hybrid. Upon chilling treatments, physiological responses and transcriptome profiles were compared across different temperature points to understand their differential responses to cold. While both parents showed divergent responses, the hybrid’s responses showed an overall resemblance to S. calendulacea, but the contribution of homeolog expression bias to cold stress was not readily evident in the F1 hybrid possibly due to inherent bias that comes with the sampling location. Our findings provided insights into the role of gene expression in differential cold tolerance, and further contribute to predicting the invasive potential of other hybrids between S. trilobata and its congeners around the world.

Rubiscosome gene expression is balanced across the hexaploid wheat genome

Functional and active Rubisco is essential for CO₂ fixation and is a primary target for engineering approaches to increasing crop yields. However, the assembly and maintenance of active Rubisco are dependent on the coordinated biosynthesis of at least 11 nuclear-encoded proteins, termed the 'Rubiscosome'. Using publicly available gene expression data for wheat (Triticum aestivum L.), we show that the expression of Rubiscosome genes is balanced across the three closely related subgenomes that form the allohexaploid genome. Each subgenome contains a near complete set of homoeologous genes and contributes equally to overall expression, both under optimal and under heat stress conditions. The expression of the wheat thermo-tolerant Rubisco activase isoform 1β increases under heat stress and remains balanced across the subgenomes, albeit with a slight shift towards greater contribution from the D subgenome. The findings show that the gene copies in all three subgenomes need to be accounted for when designing strategies for crop improvement.

Mapping and Characterization of a Wheat Stem Rust Resistance Gene in Durum Wheat "Kronos"

Wheat stem (or black) rust is one of the most devastating fungal diseases, threatening global wheat production. Identification, mapping, and deployment of effective resistance genes are critical to addressing this challenge. In this study, we mapped and characterized one stem rust resistance (Sr) gene from the tetraploid durum wheat variety Kronos (temporary designation SrKN). This gene was mapped on the long arm of chromosome 2B and confers resistance to multiple virulent Pgt races, such as TRTTF and BCCBC. Using a large mapping population (3,366 gametes), we mapped SrKN within a 0.29 cM region flanked by the sequenced-based markers pku4856F2R2 and pku4917F3R3, which corresponds to 5.6- and 7.2-Mb regions in the Svevo and Chinese Spring reference genomes, respectively. Both regions include a cluster of nucleotide binding leucine-repeat (NLR) genes that likely includes the candidate gene. An allelism test failed to detect recombination between SrKN and the previously mapped Sr9e gene. This result, together with the similar seedling resistance responses and resistance profiles, suggested that SrKN and Sr9e may represent the same gene. We introgressed SrKN into common wheat and developed completely linked markers to accelerate its deployment in the wheat breeding programs. SrKN can be a valuable component of transgenic cassettes or gene pyramids that includes multiple resistance genes to control this devastating disease.

Relaxed selection on male mitochondrial genes in DUI bivalves eases the need for mitonuclear coevolution

Mitonuclear coevolution is an important prerequisite for efficient energy production in eukaryotes. However, many bivalve taxa experience doubly uniparental inheritance (DUI) and have sex-specific mitochondrial (mt) genomes, providing a challenge for mitonuclear coevolution. We examined possible mechanisms to reconcile mitonuclear coevolution with DUI. No nuclear-encoded, sex-specific OXPHOS paralogs were found in the DUI clam Ruditapes philippinarum, refuting OXPHOS paralogy as a solution in this species. It is also unlikely that mt changes causing disruption of nuclear interactions are strongly selected against because sex-specific mt-residues or those under positive selection in M mt genes were not depleted for contacting nuclear-encoded residues. However, M genomes showed consistently higher d_N /d_S ratios compared to putatively ancestral F genomes in all mt OXPHOS genes and across all DUI species. Further analyses indicated that this was consistently due to relaxed, not positive selection on M vs. F mt OXPHOS genes. Similarly, selection was relaxed on the F genome of DUI species compared to species with strict maternal inheritance. Coupled with recent physiological and molecular evolution studies, we suggest that relaxed selection on M mt function limits the need to maintain mitonuclear interactions in M genomes compared to F genomes. We discuss our findings with regard to OXPHOS function and the origin of DUI.

De Novo Transcriptome Assembly, Functional Annotation, and Transcriptome Dynamics Analyses Reveal Stress Tolerance Genes in Mangrove Tree (Bruguiera gymnorhiza)

Their high adaptability to difficult coastal conditions makes mangrove trees a valuable resource and an interesting model system for understanding the molecular mechanisms underlying stress tolerance and adaptation of plants to the stressful environmental conditions. In this study, we used RNA sequencing (RNA-Seq) for de novo assembling and characterizing the Bruguiera gymnorhiza (L.) Lamk leaf transcriptome. B. gymnorhiza is one of the most widely distributed mangrove species from the biggest family of mangroves; Rhizophoraceae. The de novo assembly was followed by functional annotations and identification of individual transcripts and gene families that are involved in abiotic stress response. We then compared the genome-wide expression profiles between two populations of B. gymnorhiza, growing under different levels of stress, in their natural habitats. One population living in high salinity environment, in the shore of the Pacific Ocean- Japan, and the other population living about one kilometre farther from the ocean, and next to the estuary of a river; in less saline and more brackish condition. Many genes involved in response to salt and osmotic stress, showed elevated expression levels in trees growing next to the ocean in high salinity condition. Validation of these genes may contribute to future salt-resistance research in mangroves and other woody plants. Furthermore, the sequences and transcriptome data provided in this study are valuable scientific resources for future comparative transcriptome research in plants growing under stressful conditions.

Relationship between the Phenylpropanoid Pathway and Dwarfism of Paspalum seashore Based on RNA-Seq and iTRAQ

Seashore paspalum is a major warm-season turfgrass requiring frequent mowing. The use of dwarf cultivars with slow growth is a promising method to decrease mowing frequency. The present study was conducted to provide an in-depth understanding of the molecular mechanism of T51 dwarfing in the phenylpropane pathway and to screen the key genes related to dwarfing. For this purpose, we obtained transcriptomic information based on RNA-Seq and proteomic information based on iTRAQ for the dwarf mutant T51 of seashore paspalum. The combined results of transcriptomic and proteomic analysis were used to identify the differential expression pattern of genes at the translational and transcriptional levels. A total of 8311 DEGs were detected at the transcription level, of which 2540 were upregulated and 5771 were downregulated. Based on the transcripts, 2910 proteins were identified using iTRAQ, of which 392 (155 upregulated and 237 downregulated) were DEPs. The phenylpropane pathway was found to be significantly enriched at both the transcriptional and translational levels. Combined with the decrease in lignin content and the increase in flavonoid content in T51, we found that the dwarf phenotype of T51 is closely related to the abnormal synthesis of lignin and flavonoids in the phenylpropane pathway. CCR and HCT may be the key genes for T51 dwarf. This study provides the basis for further study on the dwarfing mechanism of seashore paspalum. The screening of key genes lays a foundation for further studies on the molecular mechanism of seashore paspalum dwarfing.

Mitonuclear Coevolution, but not Nuclear Compensation, Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.

Mutational signatures associated with exposure to carcinogenic microplastic compounds bisphenol A and styrene oxide

Microplastic pollutants in oceans and food chains are concerning to public health. Common plasticizing compounds Bisphenol-A (BPA) and Styrene-7,8-Oxide (SO) are now labeled as carcinogens. We show that BPA and SO cause deoxyribonucleic acid damage and mutagenesis in human cells, and analyze the genome-wide point mutation and genomic rearrangement patterns associated with BPA and SO exposure. A subset of the single- and doublet base substitutions shows mutagenesis near or at guanine, consistent with these compounds' preferences to form guanosine adducts. Presence of other mutational signatures suggest additional mutagenesis probably due to complex effects of BPA and SO on diverse cellular processes. Analyzing data for 19 cancer cohorts, we find that tumors of digestive and urinary organs show relatively high similarity in mutational profiles, and the burden of such mutations increases with age. Even within the same cancer type, proportions of corresponding mutational patterns vary among the cohorts from different countries, as does the amount of microplastic waste in ocean waters. BPA and SO are relatively mild mutagens, and other environmental agents can also potentially generate similar, complex mutational patterns in cancer genomes. Nonetheless, our findings call for systematic evaluation of public health consequences of microplastic exposure worldwide.

Multi-Omics Analysis of Small RNA, Transcriptome, and Degradome in T. turgidum-Regulatory Networks of Grain Development and Abiotic Stress Response

Crop reproduction is highly sensitive to water deficit and heat stress. The molecular networks of stress adaptation and grain development in tetraploid wheat (Triticum turgidum durum) are not well understood. Small RNAs (sRNAs) are important epigenetic regulators connecting the transcriptional and post-transcriptional regulatory networks. This study presents the first multi-omics analysis of the sRNAome, transcriptome, and degradome in T. turgidum developing grains, under single and combined water deficit and heat stress. We identified 690 microRNAs (miRNAs), with 84 being novel, from 118 sRNA libraries. Complete profiles of differentially expressed miRNAs (DEMs) specific to genotypes, stress types, and different reproductive time-points are provided. The first degradome sequencing report for developing durum grains discovered a significant number of new target genes regulated by miRNAs post-transcriptionally. Transcriptome sequencing profiled 53,146 T. turgidum genes, swith differentially expressed genes (DEGs) enriched in functional categories such as nutrient metabolism, cellular differentiation, transport, reproductive development, and hormone transduction pathways. miRNA-mRNA networks that affect grain characteristics such as starch synthesis and protein metabolism were constructed on the basis of integrated analysis of the three omics. This study provides a substantial amount of novel information on the post-transcriptional networks in T. turgidum grains, which will facilitate innovations for breeding programs aiming to improve crop resilience and grain quality.

Evolution of NLR resistance genes with noncanonical N-terminal domains in wild tomato species

Nucleotide-binding and leucine-rich repeat immune receptors (NLRs) provide resistance against diverse pathogens. To create comparative NLR resources, we conducted resistance gene enrichment sequencing (RenSeq) with single-molecule real-time sequencing of PacBio for 18 accessions in Solanaceae, including 15 accessions of five wild tomato species. We investigated the evolution of a class of NLRs, CNLs with extended N-terminal sequences previously named Solanaceae Domain. Through comparative genomic analysis, we revealed that the extended CNLs (exCNLs) anciently emerged in the most recent common ancestor between Asterids and Amaranthaceae, far predating the Solanaceae family. In tomatoes, the exCNLs display exceptional modes of evolution in a clade-specific manner. In the clade G3, exCNLs have substantially elongated their N-termini through tandem duplications of exon segments. In the clade G1, exCNLs have evolved through recent proliferation and sequence diversification. In the clade G6, an ancestral exCNL has lost its N-terminal domains in the course of evolution. Our study provides high-quality NLR gene models for close relatives of domesticated tomatoes that can serve as a useful resource for breeding and molecular engineering for disease resistance. Our findings regarding the exCNLs offer unique backgrounds and insights for future functional studies of the NLRs.

Epistatic interactions between PHOTOPERIOD1, CONSTANS1 and CONSTANS2 modulate the photoperiodic response in wheat

In Arabidopsis, CONSTANS (CO) integrates light and circadian clock signals to promote flowering under long days (LD). In the grasses, a duplication generated two paralogs designated as CONSTANS1 (CO1) and CONSTANS2 (CO2). Here we show that in tetraploid wheat plants grown under LD, combined loss-of-function mutations in the A and B-genome homeologs of CO1 and CO2 (co1 co2) result in a small (3 d) but significant (P<0.0001) acceleration of heading time both in PHOTOPERIOD1 (PPD1) sensitive (Ppd-A1b, functional ancestral allele) and insensitive (Ppd-A1a, functional dominant allele) backgrounds. Under short days (SD), co1 co2 mutants headed 13 d earlier than the wild type (P<0.0001) in the presence of Ppd-A1a. However, in the presence of Ppd-A1b, spikes from both genotypes failed to emerge by 180 d. These results indicate that CO1 and CO2 operate mainly as weak heading time repressors in both LD and SD. By contrast, in ppd1 mutants with loss-of-function mutations in both PPD1 homeologs, the wild type Co1 allele accelerated heading time >60 d relative to the co1 mutant allele under LD. We detected significant genetic interactions among CO1, CO2 and PPD1 genes on heading time, which were reflected in complex interactions at the transcriptional and protein levels. Loss-of-function mutations in PPD1 delayed heading more than combined co1 co2 mutations and, more importantly, PPD1 was able to perceive and respond to differences in photoperiod in the absence of functional CO1 and CO2 genes. Similarly, CO1 was able to accelerate heading time in response to LD in the absence of a functional PPD1. Taken together, these results indicate that PPD1 and CO1 are able to respond to photoperiod in the absence of each other, and that interactions between these two photoperiod pathways at the transcriptional and protein levels are important to fine-tune the flowering response in wheat.

phylostratr: a framework for phylostratigraphy

MOTIVATION

The goal of phylostratigraphy is to infer the evolutionary origin of each gene in an organism. This is done by searching for homologs within increasingly broad clades. The deepest clade that contains a homolog of the protein(s) encoded by a gene is that gene's phylostratum.

RESULTS

We have created a general R-based framework, phylostratr, to estimate the phylostratum of every gene in a species. The program fully automates analysis: selecting species for balanced representation, retrieving sequences, building databases, inferring phylostrata and returning diagnostics. Key diagnostics include: detection of genes with inferred homologs in old clades, but not intermediate ones; proteome quality assessments; false-positive diagnostics, and checks for missing organellar genomes. phylostratr allows extensive customization and systematic comparisons of the influence of analysis parameters or genomes on phylostrata inference. A user may: modify the automatically generated clade tree or use their own tree; provide custom sequences in place of those automatically retrieved from UniProt; replace BLAST with an alternative algorithm; or tailor the method and sensitivity of the homology inference classifier. We show the utility of phylostratr through case studies in Arabidopsis thaliana and Saccharomyces cerevisiae.

AVAILABILITY AND IMPLEMENTATION

Source code available at https://github.com/arendsee/phylostratr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Apomixis Technology: Separating the Wheat from the Chaff

Projections indicate that current plant breeding approaches will be unable to incorporate the global crop yields needed to deliver global food security. Apomixis is a disruptive innovation by which a plant produces clonal seeds capturing heterosis and gene combinations of elite phenotypes. Introducing apomixis into hybrid cultivars is a game-changing development in the current plant breeding paradigm that will accelerate the generation of high-yield cultivars. However, apomixis is a developmentally complex and genetically multifaceted trait. The central problem behind current constraints to apomixis breeding is that the genomic configuration and molecular mechanism that initiate apomixis and guide the formation of a clonal seed are still unknown. Today, not a single explanation about the origin of apomixis offer full empirical coverage, and synthesizing apomixis by manipulating individual genes has failed or produced little success. Overall evidence suggests apomixis arise from a still unknown single event molecular mechanism with multigenic effects. Disentangling the genomic basis and complex genetics behind the emergence of apomixis in plants will require the use of novel experimental approaches benefiting from Next Generation Sequencing technologies and targeting not only reproductive genes, but also the epigenetic and genomic configurations associated with reproductive phenotypes in homoploid sexual and apomictic carriers. A comprehensive picture of most regulatory changes guiding apomixis emergence will be central for successfully installing apomixis into the target species by exploiting genetic modification techniques.

Dissecting the Genetic Complexity of Fusarium Crown Rot Resistance in Wheat

Fusarium crown rot (FCR) is one of the most important diseases of wheat (Triticum aestivum L.). FCR is mainly caused by the fungal pathogens Fusarium culmorum and F. pseudograminearum. In order to identify new sources of resistance to FCR and to dissect the complexity of FCR resistance, a panel of 161 wheat accessions was phenotyped under growth room (GR) and greenhouse conditions (GH). Analysis of variance showed significant differences in crown rot development among wheat accessions and high heritability of genotype-environment interactions for GR (0.96) and GH (0.91). Mixed linear model analysis revealed seven novel quantitative trait loci (QTLs) linked to F. culmorum on chromosomes 2AL, 3AS, 4BS, 5BS, 5DS, 5DL and 6DS for GR and eight QTLs on chromosomes on 3AS, 3BS, 3DL, 4BS (2), 5BS, 6BS and 6BL for GH. Total phenotypic variances (R²) explained by the QTLs linked to GR and GH were 48% and 59%, respectively. In addition, five favorable epistasis interactions among the QTLs were detected for both GR and GH with and without main effects. Epistatic interaction contributed additional variation up to 21% under GR and 7% under GH indicating strong effects of environment on the expression of QTLs. Our results revealed FCR resistance responses in wheat to be complex and controlled by multiple QTLs.

Signatures of the Evolution of Parthenogenesis and Cryptobiosis in the Genomes of Panagrolaimid Nematodes

Most animal species reproduce sexually and fully parthenogenetic lineages are usually short lived in evolution. Still, parthenogenesis may be advantageous as it avoids the cost of sex and permits colonization by single individuals. Panagrolaimid nematodes have colonized environments ranging from arid deserts to Arctic and Antarctic biomes. Many are obligatory meiotic parthenogens, and most have cryptobiotic abilities, being able to survive repeated cycles of complete desiccation and freezing. To identify systems that may contribute to these striking abilities, we sequenced and compared the genomes and transcriptomes of parthenogenetic and outcrossing panagrolaimid species, including cryptobionts and non-cryptobionts. The parthenogens are triploids, most likely originating through hybridization. Adaptation to cryptobiosis shaped the genomes of panagrolaimid nematodes and is associated with the expansion of gene families and signatures of selection on genes involved in cryptobiosis. All panagrolaimids have acquired genes through horizontal gene transfer, some of which are likely to contribute to cryptobiosis.

Comprehensive Stress-Based De Novo Transcriptome Assembly and Annotation of Guar (Cyamopsis tetragonoloba (L.) Taub.): An Important Industrial and Forage Crop

The forage crop Guar (Cyamopsis tetragonoloba (L.) Taub.) has the ability to endure heat, drought, and mild salinity. A complete image on its genic architecture will promote our understanding about gene expression networks and different tolerance mechanisms at the molecular level. Therefore, whole mRNA sequence approach on the Guar plant was conducted to provide a snapshot of the mRNA information in the cell under salinity, heat, and drought stresses to be integrated with previous transcriptomic studies. RNA-Seq technology was employed to perform a 2 × 100 paired-end sequencing using an Illumina HiSeq 2500 platform for the transcriptome of leaves of C. tetragonoloba under normal, heat, drought, and salinity conditions. Trinity was used to achieve a de novo assembly followed by gene annotation, functional classification, metabolic pathway analysis, and identification of SSR markers. A total of 218.2 million paired-end raw reads (~44 Gbp) were generated. Of those, 193.5M paired-end reads of high quality were used to reconstruct a total of 161,058 transcripts (~266 Mbp) with N50 of 2552 bp and 61,508 putative genes. There were 6463 proteins having >90% full-length coverage against the Swiss-Prot database and 94% complete orthologs against Embryophyta. Approximately, 62.87% of transcripts were blasted, 50.46% mapped, and 43.50% annotated. A total of 4715 InterProScan families, 3441 domains, 74 repeats, and 490 sites were detected. Biological processes, molecular functions, and cellular components comprised 64.12%, 25.42%, and 10.4%, respectively. The transcriptome was associated with 985 enzymes and 156 KEGG pathways. A total of 27,066 SSRs were gained with an average frequency of one SSR/9.825 kb in the assembled transcripts. This resulting data will be helpful for the advanced analysis of Guar to multi-stress tolerance.

Genomic tools for durum wheat breeding: de novo assembly of Svevo transcriptome and SNP discovery in elite germplasm

BACKGROUND

The tetraploid durum wheat (Triticum turgidum L. ssp. durum Desf. Husnot) is an important crop which provides the raw material for pasta production and a valuable source of genetic diversity for breeding hexaploid wheat (Triticum aestivum L.). Future breeding efforts to enhance yield potential and climate resilience will increasingly rely on genomics-based approaches to identify and select beneficial alleles. A deeper characterisation of the molecular and functional diversity of the durum wheat transcriptome will be instrumental to more effectively harness its genetic diversity.

RESULTS

We report on the de novo transcriptome assembly of durum wheat cultivar 'Svevo'. The transcriptome of four tissues/organs (shoots and roots at the seedling stage, reproductive organs and developing grains) was assembled de novo, yielding 180,108 contigs, with a N50 length of 1121 bp and mean contig length of 883 bp. Alignment against the transcriptome of nine plant species identified 43% of transcripts with homology to at least one reference transcriptome. The functional annotation was completed by means of a combination of complementary software. The presence of differential expression between the A- and B-homoeolog copies of the durum wheat tetraploid genome was ascertained by phase reconstruction of polymorphic sites based on the T. urartu transcripts and inferring homoeolog-specific sequences. We observed greater expression divergence between A and B homoeologs in grains rather than in leaves and roots. The transcriptomes of 13 durum wheat cultivars spanning the breeding period from 1969 to 2005 were analysed for SNP diversity, leading to 95,358 non-rare, hemi-SNPs shared among two or more cultivars and 33,747 locus-specific (diploid inheritance) SNPs.

CONCLUSIONS

Our study updates and expands the de novo transcriptome reference assembly available for durum wheat. Out of 180,108 assembled transcripts, 13,636 were specific to the Svevo cultivar as compared to the only other reference transcriptome available for durum, thus contributing to the identification of the tetraploid wheat pan-transcriptome. Additionally, the analysis of 13 historically relevant hallmark varieties produced a SNP dataset that could successfully validate the genotyping in tetraploid wheat and provide a valuable resource for genomics-assisted breeding of both tetraploid and hexaploid wheats.

Exploiting next-generation sequencing to solve the haplotyping puzzle in polyploids: a simulation study

Haplotypes are the units of inheritance in an organism, and many genetic analyses depend on their precise determination. Methods for haplotyping single individuals use the phasing information available in next-generation sequencing reads, by matching overlapping single-nucleotide polymorphisms while penalizing post hoc nucleotide corrections made. Haplotyping diploids is relatively easy, but the complexity of the problem increases drastically for polyploid genomes, which are found in both model organisms and in economically relevant plant and animal species. Although a number of tools are available for haplotyping polyploids, the effects of the genomic makeup and the sequencing strategy followed on the accuracy of these methods have hitherto not been thoroughly evaluated.We developed the simulation pipeline haplosim to evaluate the performance of three haplotype estimation algorithms for polyploids: HapCompass, HapTree and SDhaP, in settings varying in sequencing approach, ploidy levels and genomic diversity, using tetraploid potato as the model. Our results show that sequencing depth is the major determinant of haplotype estimation quality, that 1 kb PacBio circular consensus sequencing reads and Illumina reads with large insert-sizes are competitive and that all methods fail to produce good haplotypes when ploidy levels increase. Comparing the three methods, HapTree produces the most accurate estimates, but also consumes the most resources. There is clearly room for improvement in polyploid haplotyping algorithms.

Resistance gene cloning from a wild crop relative by sequence capture and association genetics

Disease resistance (R) genes from wild relatives could be used to engineer broad-spectrum resistance in domesticated crops. We combined association genetics with R gene enrichment sequencing (AgRenSeq) to exploit pan-genome variation in wild diploid wheat and rapidly clone four stem rust resistance genes. AgRenSeq enables R gene cloning in any crop that has a diverse germplasm panel.

Atlantic salmon (Salmo salar L.) genetics in the 21st century: taking leaps forward in aquaculture and biological understanding

Atlantic salmon (Salmo salar L.) is among the most iconic and economically important fish species and was the first member of Salmonidae to have a high-quality reference genome assembly published. Advances in genomics have become increasingly central to the genetic improvement of farmed Atlantic salmon as well as conservation of wild salmon stocks. The salmon genome has also been pivotal in shaping our understanding of the evolutionary and functional consequences arising from an ancestral whole-genome duplication event characterising all Salmonidae members. Here, we provide a review of the current status of Atlantic salmon genetics and genomics, focussed on progress made from genome-wide research aimed at improving aquaculture production and enhancing understanding of salmonid ecology, physiology and evolution. We present our views on the future direction of salmon genomics, including the role of emerging technologies (e.g. genome editing) in elucidating genetic features that underpin functional variation in traits of commercial and evolutionary importance.

De novo assembly of wheat root transcriptomes and transcriptional signature of longitudinal differentiation

Hidden underground, root systems constitute an important part of the plant for its development, nourishment and sensing the soil environment around it, but we know very little about its genetic regulation in crop plants like wheat. In the present study, we de novo assembled the root transcriptomes in reference cultivar Chinese Spring from RNA-seq reads generated by the 454-GS-FLX and HiSeq platforms. The FLX reads were assembled into 24,986 transcripts with completeness of 54.84%, and the HiSeq reads were assembled into 91,543 high-confidence protein-coding transcripts, 2,404 low-confidence protein-coding transcripts, and 13,181 non-coding transcripts with the completeness of >90%. Combining the FLX and HiSeq assemblies, we assembled a root transcriptome of 92,335 ORF-containing transcripts. Approximately 7% of the coding transcripts and ~2% non-coding transcripts are not present in the current wheat genome assembly. Functional annotation of both assemblies showed similar gene ontology patterns and that ~7% coding and >5% non-coding transcripts are root-specific. Transcription quantification identified 1,728 differentially expressed transcripts between root tips and maturation zone, and functional annotation of these transcripts captured a transcriptional signature of longitudinal development of wheat root. With the transcriptomic resources developed, this study provided the first view of wheat root transcriptome under different developmental zones and laid a foundation for molecular studies of wheat root development and growth using a reverse genetic approach.

Comparative Transcriptomics in Two Bivalve Species Offers Different Perspectives on the Evolution of Sex-Biased Genes

Comparative genomics has become a central tool for evolutionary biology, and a better knowledge of understudied taxa represents the foundation for future work. In this study, we characterized the transcriptome of male and female mature gonads in the European clam Ruditapes decussatus, compared with that in the Manila clam Ruditapes philippinarum providing, for the first time in bivalves, information about transcription dynamics and sequence evolution of sex-biased genes. In both the species, we found a relatively low number of sex-biased genes (1,284, corresponding to 41.3% of the orthologous genes between the two species), probably due to the absence of sexual dimorphism, and the transcriptional bias is maintained in only 33% of the orthologs. The dN/dS is generally low, indicating purifying selection, with genes where the female-biased transcription is maintained between the two species showing a significantly higher dN/dS. Genes involved in embryo development, cell proliferation, and maintenance of genome stability show a faster sequence evolution. Finally, we report a lack of clear correlation between transcription level and evolutionary rate in these species, in contrast with studies that reported a negative correlation. We discuss such discrepancy and call into question some methodological approaches and rationales generally used in this type of comparative studies.

A Genome-wide View of Transcriptome Dynamics During Early Spike Development in Bread Wheat

Wheat spike development is a coordinated process of cell proliferation and differentiation with distinctive phases and architecture changes. However, the dynamic alteration of gene expression in this process remains enigmatic. Here, we characterized and dissected bread wheat spike into six developmental stages, and used genome-wide gene expression profiling, to investigate the underlying regulatory mechanisms. High gene expression correlations between any two given stages indicated that wheat early spike development is controlled by a small subset of genes. Throughout, auxin signaling increased, while cytokinin signaling decreased. Besides, many genes associated with stress responses highly expressed during the double ridge stage. Among the differentially expressed genes (DEGs), were identified 375 transcription factor (TF) genes, of which some homologs in rice or Arabidopsis are proposed to function in meristem maintenance, flowering time, meristem initiation or transition, floral organ development or response to stress. Gene expression profiling demonstrated that these genes had either similar or distinct expression pattern in wheat. Several genes regulating spike development were expressed in the early spike, of which Earliness per se 3 (Eps-3) was found might function in the initiation of spikelet meristem. Our study helps uncover important genes associated with apical meristem morphology and development in wheat.

Genome-wide exploration of C2H2 zinc finger family in durum wheat (Triticum turgidum ssp. Durum): insights into the roles in biological processes especially stress response

C2H2 type of zinc finger transcription factors (C2H2-ZFP TFs) play crucial roles in plant developments and stress response. Regarding its importance, genome-wide study of C2H2-ZFs were performed in multiple important plant species, but any such investigation was not fulfilled in Triticum turgidum ssp. Durum (durum wheat) as an important nutritional crop. The present study identified 122 C2H2-ZFs in durum wheat and physically mapped them onto the genome. The phylogenetic analysis classified these TFs into six major groups. Genes structure and conserved motifs assay showed TtC2H2-ZF involvement in the important cellular functions. Comparative phylogeny between durum wheat TtC2H2-ZF genes and the orthologs in rice revealed the evolutionary relationships of C2H2-ZF proteins. The gene ontology and promoter cis-element analysis indicated that most of TtC2H2-ZF genes are involved in multiple molecular functions including metal ion-binding and various stimuli responses. Further, the miRNAs targeting TtC2H2-ZF transcripts, homology modeling and proteins interaction network were also demonstrated, suggesting the vital cellular functions of TtC2H2-ZFs during various circumstances. The expression heatmap demonstrated differential and tissue-specific expression patterns of these genes. Expression profiling of this gene family members in response to dehydration and heat stresses showed differential expression pattern of these genes at multiple time points of stresses. This study can prepare a comprehensive overview of the durum wheat C2H2-ZF gene family and may provide a new perspective on the evolution of them, which will form the basis for further investigation of the roles of this family members and future genetic engineering studies in crops.

Effects of different nitrogen fertilizers on two wheat cultivars: An integrated approach

Investigation of cultivated plant physiology grown under low energy input plays an important role to indicate their fitness to the new environmental conditions. The durum-wheat cultivars Creso and Dylan were tested to evaluate the growth, production, and proteomic and transcriptomic profiles of the crop under different synthetic and organic nitrogen fertilization regimes. In this work, a two-dimensional gel electrophoresis (2-DE) approach combined with liquid chromatography-mass spectrometry (LC-MS) was used to investigate the protein changes induced by the use of different nitrogen sources (hydrolysate of proteins 1 and 2, rhizovit, synthesis, leather) on wheat plants. Proteomic studies were integrated with qPCR analysis of genes related to glutamine synthetase/glutamine-2-oxoglutarate aminotransferase (GS-GOGAT) and tricarboxylic acid (TCA) metabolic pathways because most relevant for nitrogen-dependent plants growth. The proteomic analysis lead to the isolation of 23 spots that were able to distinguish the analyzed samples. These spots yielded the identification of 60 proteins involved in photosynthesis, glycolysis, and nitrogen metabolism. As an example, the quinone oxidoreductase-like protein and probable glutathione S-transferase GSTU proteins were identified in two spots that represents the most statistically significant ones in Dylan samples. Transcript analysis indicated that related genes exhibited different expression trends; the heat map also revealed the different behaviors of the hydrolysates of the proteins 1 and 2 nitrogen sources. The effects of nitrogenous fertilizers at the proteomic and agronomic levels revealed that plants fertilized with synthesis or rhizovit gave the best results concerning yield, whereas rhizovit and protein hydrolysates were most effective for proteins content in the grain (% of dry weight). Therefore, all parameters measured in this study indicated that different kinds of nitrogen fertilization used have a relevant impact on plant growth and production.

A Robust Methodology for Assessing Differential Homeolog Contributions to the Transcriptomes of Allopolyploids

Polyploidy has played a pivotal and recurring role in angiosperm evolution. Allotetraploids arise from hybridization between species and possess duplicated gene copies (homeologs) that serve redundant roles immediately after polyploidization. Although polyploidization is a major contributor to plant evolution, it remains poorly understood. We describe an analytical approach for assessing homeolog-specific expression that begins with de novo assembly of parental transcriptomes and effectively (i) reduces redundancy in de novo assemblies, (ii) identifies putative orthologs, (iii) isolates common regions between orthologs, and (iv) assesses homeolog-specific expression using a robust Bayesian Poisson-Gamma model to account for sequence bias when mapping polyploid reads back to parental references. Using this novel methodology, we examine differential homeolog contributions to the transcriptome in the recently formed allopolyploids Tragopogon mirus and T. miscellus (Compositae). Notably, we assess a larger Tragopogon gene set than previous studies of this system. Using carefully identified orthologous regions and filtering biased orthologs, we find in both allopolyploids largely balanced expression with no strong parental bias. These new methods can be used to examine homeolog expression in any tetrapolyploid system without requiring a reference genome.

Assessment of an Organ-Specific de Novo Transcriptome of the Nematode Trap-Crop, Solanum sisymbriifolium

Solanum sisymbriifolium, also known as “Litchi Tomato” or “Sticky Nightshade,” is an undomesticated and poorly researched plant related to potato and tomato. Unlike the latter species, S. sisymbriifolium induces eggs of the cyst nematode, Globodera pallida, to hatch and migrate into its roots, but then arrests further nematode maturation. In order to provide researchers with a partial blueprint of its genetic make-up so that the mechanism of this response might be identified, we used single molecule real time (SMRT) sequencing to compile a high quality de novo transcriptome of 41,189 unigenes drawn from individually sequenced bud, root, stem, and leaf RNA populations. Functional annotation and BUSCO analysis showed that this transcriptome was surprisingly complete, even though it represented genes expressed at a single time point. By sequencing the 4 organ libraries separately, we found we could get a reliable snapshot of transcript distributions in each organ. A divergent site analysis of the merged transcriptome indicated that this species might have undergone a recent genome duplication and re-diploidization. Further analysis indicated that the plant then retained a disproportionate number of genes associated with photosynthesis and amino acid metabolism in comparison to genes with characteristics of R-proteins or involved in secondary metabolism. The former processes may have given S. sisymbriifolium a bigger competitive advantage than the latter did.

Genome-specific histories of divergence and introgression between an allopolyploid unisexual salamander lineage and two ancestral sexual species

Quantifying introgression between sexual species and polyploid lineages traditionally thought to be asexual is an important step in understanding what drives the longevity of putatively asexual groups. Here, we capitalize on three recent innovations-ultraconserved element (UCE) sequencing, bioinformatic techniques for identifying genome-specific variation in polyploids, and model-based methods for evaluating historical gene flow-to measure the extent and tempo of introgression over the evolutionary history of an allopolyploid lineage of all-female salamanders and two ancestral sexual species. Our analyses support a scenario in which the genomes sampled in unisexual salamanders last shared a common ancestor with genomes in their parental species ∼3.4 million years ago, followed by a period of divergence between homologous genomes. Recently, secondary introgression has occurred at different times with each sexual species during the last 500,000 years. Sustained introgression of sexual genomes into the unisexual lineage is the defining characteristic of their reproductive mode, but this study provides the first evidence that unisexual genomes have undergone long periods of divergence without introgression. Unlike other sperm-dependent taxa in which introgression is rare, the alternating periods of divergence and introgression between unisexual salamanders and their sexual relatives could explain why these salamanders are among the oldest described unisexual animals.

Wheat genetic resources in the post-genomics era: promise and challenges

Background

Wheat genetic resources have been used for genetic improvement since 1876, when Stephen Wilson (Transactions and Proceedings of the Botanical Society of Edinburgh 12: 286) consciously made the first wide hybrid involving wheat and rye in Scotland. Wide crossing continued with sporadic attempts in the first half of 19th century and became a sophisticated scientific discipline during the last few decades with considerable impact in farmers' fields. However, a large diversity of untapped genetic resources could contribute in meeting future wheat production challenges.

Perspectives and Conclusion

Recently the complete reference genome of hexaploid (Chinese Spring) and tetraploid (Triticum turgidum ssp. dicoccoides) wheat became publicly available coupled with on-going international efforts on wheat pan-genome sequencing. We anticipate that an objective appraisal is required in the post-genomics era to prioritize genetic resources for use in the improvement of wheat production if the goal of doubling yield by 2050 is to be met. Advances in genomics have resulted in the development of high-throughput genotyping arrays, improved and efficient methods of gene discovery, genomics-assisted selection and gene editing using endonucleases. Likewise, ongoing advances in rapid generation turnover, improved phenotyping, envirotyping and analytical methods will significantly accelerate exploitation of exotic genes and increase the rate of genetic gain in breeding. We argue that the integration of these advances will significantly improve the precision and targeted identification of potentially useful variation in the wild relatives of wheat, providing new opportunities to contribute to yield and quality improvement, tolerance to abiotic stresses, resistance to emerging biotic stresses and resilience to weather extremes.

Construction of Agropyron Gaertn. genetic linkage maps using a wheat 660K SNP array reveals a homoeologous relationship with the wheat genome

Agropyron Gaertn. (P genome) is a wild relative of wheat that harbours many genetic variations that could be used to increase the genetic diversity of wheat. To agronomically transfer important genes from the P genome to a wheat chromosome by induced homoeologous pairing and recombination, it is necessary to determine the chromosomal relationships between Agropyron and wheat. Here, we report using the wheat 660K single nucleotide polymorphism (SNP) array to genotype a segregating Agropyron F1 population derived from an interspecific cross between two cross-pollinated diploid collections 'Z1842' [A. cristatum (L.) Beauv.] (male parent) and 'Z2098' [A. mongolicum Keng] (female parent) and 35 wheat-A. cristatum addition/substitution lines. Genetic linkage maps were constructed using 913 SNP markers distributed among seven linkage groups spanning 839.7 cM. The average distance between adjacent markers was 1.8 cM. The maps identified the homoeologous relationship between the P genome and wheat and revealed that the P and wheat genomes are collinear and relatively conserved. In addition, obvious rearrangements and introgression spread were observed throughout the P genome compared with the wheat genome. Combined with genotyping data, the complete set of wheat-A. cristatum addition/substitution lines was characterized according to their homoeologous relationships. In this study, the homoeologous relationship between the P genome and wheat was identified using genetic linkage maps, and the detection mean for wheat-A. cristatum introgressions might significantly accelerate the introgression of genetic variation from Agropyron into wheat for exploitation in wheat improvement programmes.

Amino acid permease 6 modulates host response to cyst nematodes in wheat and Arabidopsis

Cyst nematodes are plant parasites that cause significant crop loss in wheat and other crops. Infective juveniles invade roots and induce syncytial feeding structures as the only source of nutrients throughout their life. A previous genome-wide association study in wheat identified amino acid permease 6 (TaAAP6) to be linked to susceptibility to the cereal cyst nematode Heterodera filipjevi. To characterise the role of AAP6 during nematode parasitism, we analysed the expression of TaAAP6 and the Arabidopsis orthologue AtAAP6. TaAAP6 was found to be highly expressed in nematode-infected roots of susceptible wheat, whereas it was not upregulated in nematode-infected roots of resistant accessions. AtAAP6 was also found to be highly upregulated in nematode-induced syncytia compared with non-infected roots. Infection assays with an AtAAP6 knock-out mutant revealed reduction in developing females, female size, and size of female-associated syncytia, thus indicating the importance of AAP6 in cyst nematode parasitism.

Transcriptome dynamics provide insights into long-term salinity stress tolerance in Triticum aestivum cv. Kharchia Local

Kharchia Local, a wheat (Triticum aestivum) cultivar, is native to the saline-sodic soils of Pali district, Rajasthan, India and well known for its salinity stress tolerance. In the present study, we performed transcriptome sequencing to compare genome wide differential expression pattern between flag leaves of salinity stressed (15 EC) and control plants at anthesis stage. The 63.9 million paired end raw reads were assembled into 74,106 unigenes, of which, 3197 unigenes were found to be differentially expressed. Functional annotation analysis revealed the upregulation of genes associated with various biological processes including signal transduction, phytohormones signaling, osmoregulation, flavonoid biosynthesis, ion transport and ROS homeostasis. Expression pattern of fourteen differentially expressed genes was validated using qRT-PCR and was found to be consistent with the results of the transcriptome sequencing. Present study is the primary report on transcriptome profiling of Kharchia Local flag leaf under long-term salinity stress at anthesis stage. In conclusion, the data generated in this study can improve our knowledge in understanding the molecular mechanism of salinity stress tolerance. It will also serve as a valuable genomic resource in wheat breeding programs.

A wheat NBS-LRR gene TaRGA19 participates in Lr19-mediated resistance to Puccinia triticina

Wheat leaf rust, caused by Puccinia triticina (Pt), is one of the most severe fungal diseases on wheat globally. Rational utilization of wheat leaf rust resistance (Lr) genes is still the best choice for control this disease. Wheat seedlings carrying Lr19 showed a high resistance phenotype to all Pt races in China. So far, all the cloned seedling Lr genes including Lr1, Lr10 and Lr21, encode protein with NBS-LRR domain. In this study, a wheat gene with NBS-LRR domain from previously established Lr19-resistance-related cDNA library was cloned and designated as TaRGA19. Full length of this gene was amplified by rapid amplification of cDNA ends (RACE). By blast against IWGSC wheat genome database, we have noticed that TaRGA19 was located on chromosome 2DS, which was different from Lr19 located on chromosome 7DL. Compared with susceptible Thatcher line, expression level of TaRGA19 was upregulated in wheat isogenic lines carrying Lr19 (TcLr19) after inoculation of Pt race THTS. By particle bombardment, TaRGA19-GFP fused protein was localized on plasma membrane of epidermal cells. Using virus-induced gene silencing (VIGS), TaRGA19-knockdown plants of TcLr19 showed reduced resistance and few sporulation phenotype upon Pt challenge. Further histological observation indicated that Pt hyphal growth at the infection sites was less suppressed in the TaRGA19-knockdown plants. In conclusion, we speculate this TaRGA19 gene was involved in the Lr19-mediated resistance to wheat leaf rust along with other components.

Haplotype Analysis of the Pre-harvest Sprouting Resistance Locus Phs-A1 Reveals a Causal Role of TaMKK3-A in Global Germplasm

Pre-harvest sprouting (PHS) is an important cause of quality loss in many cereal crops and is particularly prevalent and damaging in wheat. Resistance to PHS is therefore a valuable target trait in many breeding programs. The Phs-A1 locus on wheat chromosome arm 4AL has been consistently shown to account for a significant proportion of natural variation to PHS in diverse mapping populations. However, the deployment of sprouting resistance is confounded by the fact that different candidate genes, including the tandem duplicated Plasma Membrane 19 (PM19) genes and the mitogen-activated protein kinase kinase 3 (TaMKK3-A) gene, have been proposed to underlie Phs-A1. To further define the Phs-A1 locus, we constructed a physical map across this interval in hexaploid and tetraploid wheat. We established close proximity of the proposed candidate genes which are located within a 1.2 Mb interval. Genetic characterization of diverse germplasm used in previous genetic mapping studies suggests that TaMKK3-A, and not PM19, is the major gene underlying the Phs-A1 effect in European, North American, Australian and Asian germplasm. We identified the non-dormant TaMKK3-A allele at low frequencies within the A-genome diploid progenitor Triticum urartu genepool, and show an increase in the allele frequency in modern varieties. In United Kingdom varieties, the frequency of the dormant TaMKK3-A allele was significantly higher in bread-making quality varieties compared to feed and biscuit-making cultivars. Analysis of exome capture data from 58 diverse hexaploid wheat accessions identified fourteen haplotypes across the extended Phs-A1 locus and four haplotypes for TaMKK3-A. Analysis of these haplotypes in a collection of United Kingdom and Australian cultivars revealed distinct major dormant and non-dormant Phs-A1 haplotypes in each country, which were either rare or absent in the opposing germplasm set. The diagnostic markers and haplotype information reported in the study will help inform the choice of germplasm and breeding strategies for the deployment of Phs-A1 resistance into breeding germplasm.

Multiple hybrid de novo genome assembly of finger millet, an orphan allotetraploid crop

Finger millet (Eleusine coracana (L.) Gaertn) is an important crop for food security because of its tolerance to drought, which is expected to be exacerbated by global climate changes. Nevertheless, it is often classified as an orphan/underutilized crop because of the paucity of scientific attention. Among several small millets, finger millet is considered as an excellent source of essential nutrient elements, such as iron and zinc; hence, it has potential as an alternate coarse cereal. However, high-quality genome sequence data of finger millet are currently not available. One of the major problems encountered in the genome assembly of this species was its polyploidy, which hampers genome assembly compared with a diploid genome. To overcome this problem, we sequenced its genome using diverse technologies with sufficient coverage and assembled it via a novel multiple hybrid assembly workflow that combines next-generation with single-molecule sequencing, followed by whole-genome optical mapping using the Bionano Irys® system. The total number of scaffolds was 1,897 with an N50 length >2.6 Mb and detection of 96% of the universal single-copy orthologs. The majority of the homeologs were assembled separately. This indicates that the proposed workflow is applicable to the assembly of other allotetraploid genomes.

Increased pericarp cell length underlies a major quantitative trait locus for grain weight in hexaploid wheat

Crop yields must increase to address food insecurity. Grain weight, determined by grain length and width, is an important yield component, but our understanding of the underlying genes and mechanisms is limited. We used genetic mapping and near isogenic lines (NILs) to identify, validate and fine-map a major quantitative trait locus (QTL) on wheat chromosome 5A associated with grain weight. Detailed phenotypic characterisation of developing and mature grains from the NILs was performed. We identified a stable and robust QTL associated with a 6.9% increase in grain weight. The positive interval leads to 4.0% longer grains, with differences first visible 12 d after fertilization. This grain length effect was fine-mapped to a 4.3 cM interval. The locus also has a pleiotropic effect on grain width (1.5%) during late grain development that determines the relative magnitude of the grain weight increase. Positive NILs have increased maternal pericarp cell length, an effect which is independent of absolute grain length. These results provide direct genetic evidence that pericarp cell length affects final grain size and weight in polyploid wheat. We propose that combining genes that control distinct biological mechanisms, such as cell expansion and proliferation, will enhance crop yields.

Transcriptomic Profiles Reveal the Interactions of Cd/Zn in Dwarf Polish Wheat (Triticum polonicum L.) Roots

Different intra- or interspecific wheat show different interactions of Cd/Zn. Normally, Zn has been/being widely utilized to reduce the Cd toxicity. In the present study, the DPW seedlings exhibited strong Cd tolerance. Zn and Cd mutually inhibited their uptake in the roots, showed antagonistic Cd/Zn interactions. However, Zn promoted the Cd transport from the roots to shoots, showed synergistic. In order to discover the interactive molecular responses, a transcriptome, including 123,300 unigenes, was constructed using RNA-Sequencing (RNA-Seq). Compared with CK, the expression of 1,269, 820, and 1,254 unigenes was significantly affected by Cd, Zn, and Cd+Zn, respectively. Only 381 unigenes were co-induced by these three treatments. Several metal transporters, such as cadmium-transporting ATPase and plant cadmium resistance 4, were specifically regulated by Cd+Zn. Other metal-related unigenes, such as ABC transporters, metal chelator, nicotianamine synthase (NAS), vacuolar iron transporters (VIT), metal-nicotianamine transporter YSL (YSL), and nitrate transporter (NRT), were regulated by Cd, but were not regulated by Cd+Zn. These results indicated that these transporters participated in the mutual inhibition of the Cd/Zn uptake in the roots, and also participated in the Cd transport, accumulation and detoxification. Meanwhile, some unigenes involved in other processes, such as oxidation-reduction, auxin metabolism, glutathione (GSH) metabolism nitrate transport, played different and important roles in the detoxification of these heavy metals.

Identification and expression analysis of the NAC transcription factor family in durum wheat (Triticum turgidum L. ssp. durum)

The NAC (NAM, ATAF and CUC) proteins belong to one of the largest plant-specific transcription factor (TF) families and play important roles in plant development processes, response to biotic and abiotic cues and hormone signaling. Our analysis led to the identification of 168 NAC genes in durum wheat, including nine putative membrane-bound TFs and 48 homeologous genes pairs. Phylogenetic analyses of TtNACs along with their Arabidopsis, grape, barley and rice counterparts divided these proteins into 8 phylogenetic groups and allowed the identification of TtNAC-A7, TtNAC-B35, TtNAC-A68, TtNAC-B69 and TtNAC-A43 as homologs of OsNAC1, OsNAC8, OsNTL2, OsNTL5 and ANAC025/NTL14, respectively. In silico expression analysis, using RNA-seq data, revealed tissue-specific and stress responsive TtNAC genes. The expression of ten selected genes was analyzed under salt and drought stresses in two contrasting tolerance cultivars. This analysis is the first report of NAC gene family in durum wheat and will be useful for the identification and selection of candidate genes associated with stress tolerance.

Uncovering hidden variation in polyploid wheat

Comprehensive reverse genetic resources, which have been key to understanding gene function in diploid model organisms, are missing in many polyploid crops. Young polyploid species such as wheat, which was domesticated less than 10,000 y ago, have high levels of sequence identity among subgenomes that mask the effects of recessive alleles. Such redundancy reduces the probability of selection of favorable mutations during natural or human selection, but also allows wheat to tolerate high densities of induced mutations. Here we exploited this property to sequence and catalog more than 10 million mutations in the protein-coding regions of 2,735 mutant lines of tetraploid and hexaploid wheat. We detected, on average, 2,705 and 5,351 mutations per tetraploid and hexaploid line, respectively, which resulted in 35-40 mutations per kb in each population. With these mutation densities, we identified an average of 23-24 missense and truncation alleles per gene, with at least one truncation or deleterious missense mutation in more than 90% of the captured wheat genes per population. This public collection of mutant seed stocks and sequence data enables rapid identification of mutations in the different copies of the wheat genes, which can be combined to uncover previously hidden variation. Polyploidy is a central phenomenon in plant evolution, and many crop species have undergone recent genome duplication events. Therefore, the general strategy and methods developed herein can benefit other polyploid crops.

Enabling Molecular Technologies for Trait Improvement in Wheat

Wheat is the major staple food crop and a source of calories for humans worldwide. A steady increase in the wheat production is essential to meet the demands of an ever-increasing global population and to achieve food security. The large size and structurally intricate genome of polyploid wheat had hindered the genomic analysis. However, with the advent of new genomic technologies such as next generation sequencing has led to genome drafts for bread wheat and its progenitors and has paved the way to design new strategies for crop improvement. Here we provide an overview of the advancements made in wheat genomics together with the available "omics approaches" and bioinformatics resources developed for wheat research. Advances in genomic, transcriptomic, and metabolomic technologies are highlighted as options to circumvent existing bottlenecks in the phenotypic and genomic selection and gene transfer. The contemporary reverse genetics approaches, including the novel genome editing techniques to inform targeted manipulation of a single/multiple genes and strategies for generating marker-free transgenic wheat plants, emphasize potential to revolutionize wheat improvement shortly.

De novo assembly, functional annotation, and analysis of the giant reed (Arundo donax L.) leaf transcriptome provide tools for the development of a biofuel feedstock

BACKGROUND

Arundo donax has attracted renewed interest as a potential candidate energy crop for use in biomass-to-liquid fuel conversion processes and biorefineries. This is due to its high productivity, adaptability to marginal land conditions, and suitability for biofuel and biomaterial production. Despite its importance, the genomic resources currently available for supporting the improvement of this species are still limited.

RESULTS

We used RNA sequencing (RNA-Seq) to de novo assemble and characterize the A. donax leaf transcriptome. The sequencing generated 1249 million clean reads that were assembled using single-k-mer and multi-k-mer approaches into 62,596 unique sequences (unitranscripts) with an N50 of 1134 bp. TransDecoder and Trinotate software suites were used to obtain putative coding sequences and annotate them by mapping to UniProtKB/Swiss-Prot and UniRef90 databases, searching for known transcripts, proteins, protein domains, and signal peptides. Furthermore, the unitranscripts were annotated by mapping them to the NCBI non-redundant, GO and KEGG pathway databases using Blast2GO. The transcriptome was also characterized by BLAST searches to investigate homologous transcripts of key genes involved in important metabolic pathways, such as lignin, cellulose, purine, and thiamine biosynthesis and carbon fixation. Moreover, a set of homologous transcripts of key genes involved in stomatal development and of genes coding for stress-associated proteins (SAPs) were identified. Additionally, 8364 simple sequence repeat (SSR) markers were identified and surveyed. SSRs appeared more abundant in non-coding regions (63.18%) than in coding regions (36.82%). This SSR dataset represents the first marker catalogue of A. donax. 53 SSRs (PolySSRs) were then predicted to be polymorphic between ecotype-specific assemblies, suggesting genetic variability in the studied ecotypes.

CONCLUSIONS

This study provides the first publicly available leaf transcriptome for the A. donax bioenergy crop. The functional annotation and characterization of the transcriptome will be highly useful for providing insight into the molecular mechanisms underlying its extreme adaptability. The identification of homologous transcripts involved in key metabolic pathways offers a platform for directing future efforts in genetic improvement of this species. Finally, the identified SSRs will facilitate the harnessing of untapped genetic diversity. This transcriptome should be of value to ongoing functional genomics and genetic studies in this crop of paramount economic importance.

Transcriptome Sequencing of Diverse Peanut (Arachis) Wild Species and the Cultivated Species Reveals a Wealth of Untapped Genetic Variability

To test the hypothesis that the cultivated peanut species possesses almost no molecular variability, we sequenced a diverse panel of 22 Arachis accessions representing Arachis hypogaea botanical classes, A-, B-, and K- genome diploids, a synthetic amphidiploid, and a tetraploid wild species. RNASeq was performed on pools of three tissues, and de novo assembly was performed. Realignment of individual accession reads to transcripts of the cultivar OLin identified 306,820 biallelic SNPs. Among 10 naturally occurring tetraploid accessions, 40,382 unique homozygous SNPs were identified in 14,719 contigs. In eight diploid accessions, 291,115 unique SNPs were identified in 26,320 contigs. The average SNP rate among the 10 cultivated tetraploids was 0.5, and among eight diploids was 9.2 per 1000 bp. Diversity analysis indicated grouping of diploids according to genome classification, and cultivated tetraploids by subspecies. Cluster analysis of variants indicated that sequences of B genome species were the most similar to the tetraploids, and the next closest diploid accession belonged to the A genome species. A subset of 66 SNPs selected from the dataset was validated; of 782 SNP calls, 636 (81.32%) were confirmed using an allele-specific discrimination assay. We conclude that substantial genetic variability exists among wild species. Additionally, significant but lesser variability at the molecular level occurs among accessions of the cultivated species. This survey is the first to report significant SNP level diversity among transcripts, and may explain some of the phenotypic differences observed in germplasm surveys. Understanding SNP variants in the Arachis accessions will benefit in developing markers for selection.

A Conserved Puccinia striiformis Protein Interacts with Wheat NPR1 and Reduces Induction of Pathogenesis-Related Genes in Response to Pathogens

In Arabidopsis, NPR1 is a key transcriptional coregulator of systemic acquired resistance. Upon pathogen challenge, NPR1 translocates from the cytoplasm to the nucleus, in which it interacts with TGA-bZIP transcription factors to activate the expression of several pathogenesis-related (PR) genes. In a screen of a yeast two-hybrid library from wheat leaves infected with Puccinia striiformis f. sp. tritici, we identified a conserved rust protein that interacts with wheat NPR1 and named it PNPi (for Puccinia NPR1 interactor). PNPi interacts with the NPR1/NIM1-like domain of NPR1 via its C-terminal DPBB_1 domain. Using bimolecular fluorescence complementation assays, we detected the interaction between PNPi and wheat NPR1 in the nucleus of Nicotiana benthamiana protoplasts. A yeast three-hybrid assay showed that PNPi interaction with NPR1 competes with the interaction between wheat NPR1 and TGA2.2. In barley transgenic lines overexpressing PNPi, we observed reduced induction of multiple PR genes in the region adjacent to Pseudomonas syringae pv. tomato DC3000 infection. Based on these results, we hypothesize that PNPi has a role in manipulating wheat defense response via its interactions with NPR1.

Transcriptome Sequencing and Expression Analysis of Cadmium (Cd) Transport and Detoxification Related Genes in Cd-Accumulating Salix integra

Salix integra is a shrub willow native to northeastern China, Japan, Korea, and Primorsky Krai in the far southeast of Russia, and has been identified as cadmium (Cd)-accumulating trees in recent years. Although many physiological studies have been conducted with these plants, little is known about the molecular basis underlying Cd response in this plant, and this is confirmed by the very few number of gene sequences (only 39 nucleotide sequences) available in public databases. Advances in genomics for Salix are promising for future improvement in identification of new candidate genes involved in metal tolerance and accumulation. Thus, high-throughput transcriptome sequencing is essential for generating enormous transcript sequences from S. integra, especially for the purpose of Cd toxicity-responsive genes discovery. Using Illumina paired-end sequencing, approximately 60.05 million high-quality reads were obtained. De novo assembly yielded 80,105 unigenes with an average length of 703 bp, A total of 50,221 (63%) unigenes were further functionally annotated by comparing their sequences to different proteins and functional domain databases. GO annotation reveals 1849 Cd responsive genes involving in Cd binding, transport, and detoxification and cellular Cd homeostasis, and these genes were highly enriched in plant response to Cd ion and Cd ion transport. By searching against the PlantCyc database, 509 unigenes were assigned to 14 PlantCyc pathways related to Cd transport and cellular detoxification, and many of them are genes encoding heavy metal ATPases (HMAs), nature resistance-associated with microphage proteins (NRAMPs), ATP-binding cassette (ABC) transporters, etc., Comprehensive RT-qPCR analysis of these selected genes in different tissues of S. integra under the control and Cd treatment revealed metallothionein-like protein (MT2A and MT2B), Metal tolerance protein (MTP1), ABCB25, NRAMP5, and ZIP1 may be involved in the Cd transport and detoxification in leaves, while NRAMP2, ZIP8, and NRAMP5 may be related to Cd transport in roots. Our study will enrich the sequence information of S. integra in public database, and would provide some new understanding of the molecular mechanisms of heavy metal tolerance and detoxification in willows.

Generation of an integrated Hieracium genomic and transcriptomic resource enables exploration of small RNA pathways during apomixis initiation

Background

Application of apomixis, or asexual seed formation, in crop breeding would allow rapid fixation of complex traits, economizing improved crop delivery. Identification of apomixis genes is confounded by the polyploid nature, high genome complexity and lack of genomic sequence integration with reproductive tissue transcriptomes in most apomicts.

Results

A genomic and transcriptomic resource was developed for Hieracium subgenus Pilosella (Asteraceae) which incorporates characterized sexual, apomictic and mutant apomict plants exhibiting reversion to sexual reproduction. Apomicts develop additional female gametogenic cells that suppress the sexual pathway in ovules. Disrupting small RNA pathways in sexual Arabidopsis also induces extra female gametogenic cells; therefore, the resource was used to examine if changes in small RNA pathways correlate with apomixis initiation. An initial characterization of small RNA pathway genes within Hieracium was undertaken, and ovary-expressed ARGONAUTE genes were identified and cloned. Comparisons of whole ovary transcriptomes from mutant apomicts, relative to the parental apomict, revealed that differentially expressed genes were enriched for processes involved in small RNA biogenesis and chromatin silencing. Small RNA profiles within mutant ovaries did not reveal large-scale alterations in composition or length distributions; however, a small number of differentially expressed, putative small RNA targets were identified.

Conclusions

The established Hieracium resource represents a substantial contribution towards the investigation of early sexual and apomictic female gamete development, and the generation of new candidate genes and markers. Observed changes in small RNA targets and biogenesis pathways within sexual and apomictic ovaries will underlie future functional research into apomixis initiation in Hieracium.

Electronic supplementary material

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