Sequence composition, organization, and evolution of the core Triticeae genome

Bread wheat satellitome: a complex scenario in a huge genome

In bread wheat (Triticum aestivum L.), chromosome associations during meiosis are extremely regulated and initiate at the telomeres and subtelomeres, which are enriched in satellite DNA (satDNA). We present the study and characterization of the bread wheat satellitome to shed light on the molecular organization of wheat subtelomeres. Our results revealed that the 2.53% of bread wheat genome is composed by satDNA and subtelomeres are particularly enriched in such DNA sequences. Thirty-four satellite DNA (21 for the first time in this work) have been identified, analyzed and cytogenetically validated. Many of the satDNAs were specifically found at particular subtelomeric chromosome regions revealing the asymmetry in subtelomere organisation among the wheat subgenomes, which might play a role in proper homologous recognition and pairing during meiosis. An integrated physical map of the wheat satellitome was also constructed. To the best of our knowledge, our results show that the combination of both cytogenetics and genome research allowed the first comprehensive analysis of the wheat satellitome, shedding light on the complex wheat genome organization, especially on the polymorphic nature of subtelomeres and their putative implication in chromosome recognition and pairing during meiosis.

Telomere and subtelomere high polymorphism might contribute to the specificity of homologous recognition and pairing during meiosis in barley in the context of breeding

Barley (Hordeum vulgare) is one of the most popular cereal crops globally. Although it is a diploid species, (2n = 2x = 14) the study of its genome organization is necessary in the framework of plant breeding since barley is often used in crosses with other cereals like wheat to provide them with advantageous characters. We already have an extensive knowledge on different stages of the meiosis, the cell division to generate the gametes in species with sexual reproduction, such as the formation of the synaptonemal complex, recombination, and chromosome segregation. But meiosis really starts with the identification of homologous chromosomes and pairing initiation, and it is still unclear how chromosomes exactly choose a partner to appropriately pair for additional recombination and segregation. In this work we present an exhaustive molecular analysis of both telomeres and subtelomeres of barley chromosome arms 2H-L, 3H-L and 5H-L. As expected, the analysis of multiple features, including transposable elements, repeats, GC content, predicted CpG islands, recombination hotspots, G4 quadruplexes, genes and targeted sequence motifs for key DNA-binding proteins, revealed a high degree of variability both in telomeres and subtelomeres. The molecular basis for the specificity of homologous recognition and pairing occurring in the early chromosomal interactions at the start of meiosis in barley may be provided by these polymorphisms. A more relevant role of telomeres and most distal part of subtelomeres is suggested.

Domestication of newly evolved hexaploid wheat—A journey of wild grass to cultivated wheat

Domestication of wheat started with the dawn of human civilization. Since then, improvement in various traits including resistance to diseases, insect pests, saline and drought stresses, grain yield, and quality were improved through selections by early farmers and then planned hybridization after the discovery of Mendel’s laws. In the 1950s, genetic variability was created using mutagens followed by the selection of superior mutants. Over the last 3 decades, research was focused on developing superior hybrids, initiating marker-assisted selection and targeted breeding, and developing genetically modified wheat to improve the grain yield, tolerance to drought, salinity, terminal heat and herbicide, and nutritive quality. Acceptability of genetically modified wheat by the end-user remained a major hurdle in releasing into the environment. Since the beginning of the 21st century, changing environmental conditions proved detrimental to achieving sustainability in wheat production particularly in developing countries. It is suggested that high-tech phenotyping assays and genomic procedures together with speed breeding procedures will be instrumental in achieving food security beyond 2050.

Comparative repeatome analysis reveals new evidence on genome evolution in wild diploid Arachis (Fabaceae) species

Opposing changes in the abundance of satellite DNA and long terminal repeat (LTR) retroelements are the main contributors to the variation in genome size and heterochromatin amount in Arachis diploids. The South American genus Arachis (Fabaceae) comprises 83 species organized in nine taxonomic sections. Among them, section Arachis is characterized by species with a wide genome and karyotype diversity. Such diversity is determined mainly by the amount and composition of repetitive DNA. Here we performed computational analysis on low coverage genome sequencing to infer the dynamics of changes in major repeat families that led to the differentiation of genomes in diploid species (x = 10) of genus Arachis, focusing on section Arachis. Estimated repeat content ranged from 62.50 to 71.68% of the genomes. Species with different genome composition tended to have different landscapes of repeated sequences. Athila family retrotransposons were the most abundant and variable lineage among Arachis repeatomes, with peaks of transpositional activity inferred at different times in the evolution of the species. Satellite DNAs (satDNAs) were less abundant, but differentially represented among species. High rates of evolution of an AT-rich superfamily of satDNAs led to the differential accumulation of heterochromatin in Arachis genomes. The relationship between genome size variation and the repetitive content is complex. However, largest genomes presented a higher accumulation of LTR elements and lower contents of satDNAs. In contrast, species with lowest genome sizes tended to accumulate satDNAs in detriment of LTR elements. Phylogenetic analysis based on repetitive DNA supported the genome arrangement of section Arachis. Altogether, our results provide the most comprehensive picture on the repeatome dynamics that led to the genome differentiation of Arachis species.

High genetic diversity in Aegilops tauschii Coss. accessions from North Iran as revealed by IRAP and REMAP markers

BACKGROUND

Aegilops tauschii Coss. as a donor of wheat D genome has an important role in wheat breeding programs. Genetic and phylogeographic diversity of 79 Ae. tauschii accessions collected from north and northwest of Iran were analyzed based on retroelement insertional polymorphisms using inter-retrotransposon amplified polymorphism (IRAP) and retrotransposon-microsatellite amplified polymorphism (REMAP) markers.

RESULTS

In total, 306 and 151 polymorphic bands were amplified in IRAP and REMAP analyses, respectively. As a result, a high level of polymorphism was observed among the studied accessions as revealed by an average of 25.5 bands per primer/primer combination and mean PIC value of 0.47 in IRAP and an average of 25.16 bands per primer combination and mean PIC value of 0.47 in REMAP. Genetic relationships of the accessions were analyzed using distance- and model-based cluster analyses.

CONCLUSION

The result showed that genetic distance did not seem to be related to geographic distribution, and the accessions could be divided into three groups, which was further supported by principal coordinate analysis. These results on genetic diversity and population structure of Ae. tauschii in Iran should provide important knowledge on genetic resources and their applications in wheat breeding programs.

Wheat omics: Classical breeding to new breeding technologies

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

Heteroalleles in Common Wheat: Multiple Differences between Allelic Variants of the Gli-B1 Locus

The Gli-B1-encoded γ-gliadins and non-coding γ-gliadin DNA sequences for 15 different alleles of common wheat have been compared using seven tests: electrophoretic mobility (EM) and molecular weight (MW) of the encoded major γ-gliadin, restriction fragment length polymorphism patterns (RFLPs) (three different markers), Gli-B1-γ-gliadin-pseudogene known SNP markers (Single nucleotide polymorphisms) and sequencing the pseudogene GAG56B. It was discovered that encoded γ-gliadins, with contrasting EM, had similar MWs. However, seven allelic variants (designated from I to VII) differed among them in the other six tests: I (alleles Gli-B1i, k, m, o), II (Gli-B1n, q, s), III (Gli-B1b), IV (Gli-B1e, f, g), V (Gli-B1h), VI (Gli-B1d) and VII (Gli-B1a). Allele Gli-B1c (variant VIII) was identical to the alleles from group IV in four of the tests. Some tests might show a fine difference between alleles belonging to the same variant. Our results attest in favor of the independent origin of at least seven variants at the Gli-B1 locus that might originate from deeply diverged genotypes of the donor(s) of the B genome in hexaploid wheat and therefore might be called “heteroallelic”. The donor’s particularities at the Gli-B1 locus might be conserved since that time and decisively contribute to the current high genetic diversity of common wheat.

The Evolutionary Dynamics of Repetitive DNA and Its Impact on the Genome Diversification in the Genus Sorghum

Polyploidization is an evolutionary event leading to structural changes of the genome(s), particularly allopolyploidization, which combines different genomes of distinct species. The tetraploid species, Sorghum halepense, is assumed an allopolyploid species formed by hybridization between diploid S. bicolor and S. propinquum. The repeat profiles of S. bicolor, S. halepense, and their relatives were compared to elucidate the repeats' role in shaping their genomes. The repeat frequencies and profiles of the three diploid accessions (S. bicolor, S. bicolor ssp. verticilliflorum, and S. bicolor var. technicum) and two tetraploid accessions (S. halepense) are similar. However, the polymorphic distribution of the subtelomeric satellites preferentially enriched in the tetraploid S. halepense indicates drastic genome rearrangements after the allopolyploidization event. Verified by CENH3 chromatin immunoprecipitation (ChIP)-sequencing and fluorescence in situ hybridization (FISH) analysis the centromeres of S. bicolor are mainly composed of the abundant satellite SorSat137 (CEN38) and diverse CRMs, Athila of Ty3_gypsy and Ty1_copia-SIRE long terminal repeat (LTR) retroelements. A similar centromere composition was found in S. halepense. The potential contribution of S. bicolor in the formation of tetraploid S. halepense is discussed.

Telomeres and Subtelomeres Dynamics in the Context of Early Chromosome Interactions During Meiosis and Their Implications in Plant Breeding

Genomic architecture facilitates chromosome recognition, pairing, and recombination. Telomeres and subtelomeres play an important role at the beginning of meiosis in specific chromosome recognition and pairing, which are critical processes that allow chromosome recombination between homologs (equivalent chromosomes in the same genome) in later stages. In plant polyploids, these terminal regions are even more important in terms of homologous chromosome recognition, due to the presence of homoeologs (equivalent chromosomes from related genomes). Although telomeres interaction seems to assist homologous pairing and consequently, the progression of meiosis, other chromosome regions, such as subtelomeres, need to be considered, because the DNA sequence of telomeres is not chromosome-specific. In addition, recombination operates at subtelomeres and, as it happens in rye and wheat, homologous recognition and pairing is more often correlated with recombining regions than with crossover-poor regions. In a plant breeding context, the knowledge of how homologous chromosomes initiate pairing at the beginning of meiosis can contribute to chromosome manipulation in hybrids or interspecific genetic crosses. Thus, recombination in interspecific chromosome associations could be promoted with the aim of transferring desirable agronomic traits from related genetic donor species into crops. In this review, we summarize the importance of telomeres and subtelomeres on chromatin dynamics during early meiosis stages and their implications in recombination in a plant breeding framework.

Sequence analysis of wheat subtelomeres reveals a high polymorphism among homoeologous chromosomes

Bread wheat, Triticum aestivum L., is one of the most important crops in the world. Understanding its genome organization (allohexaploid; AABBDD; 2n = 6x = 42) is essential for geneticists and plant breeders. Particularly, the knowledge of how homologous chromosomes (equivalent chromosomes from the same genome) specifically recognize each other to pair at the beginning of meiosis, the cellular process to generate gametes in sexually reproducing organisms, is fundamental for plant breeding and has a big influence on the fertility of wheat plants. Initial homologous chromosome interactions contribute to specific recognition and pairing between homologues at the onset of meiosis. Understanding the molecular basis of these critical processes can help to develop genetic tools in a breeding context to promote interspecific chromosome associations in hybrids or interspecific genetic crosses to facilitate the transfer of desirable agronomic traits from related species into a crop like wheat. The terminal regions of chromosomes, which include telomeres and subtelomeres, participate in chromosome recognition and pairing. We present a detailed molecular analysis of subtelomeres of wheat chromosome arms 1AS, 4AS, 7AS, 7BS and 7DS. Results showed a high polymorphism in the subtelomeric region among homoeologues (equivalent chromosomes from related genomes) for all the features analyzed, including genes, transposable elements, repeats, GC content, predicted CpG islands, recombination hotspots and targeted sequence motifs for relevant DNA-binding proteins. These polymorphisms might be the molecular basis for the specificity of homologous recognition and pairing in initial chromosome interactions at the beginning of meiosis in wheat.

Development of genome wide transposable elements based repeat junction markers in Jatropha (Jatropha curcas L.)

Jatropha curcas is a potential biodiesel crop and a highly adaptable species to various agro-climatic conditions. In this study, we have utilized transposable elements' (TE) repeat junctions (RJs) which are an important constituent of the genome, used to form a genome-wide molecular marker platform owing to its use in genomic studies of plants. We screened our previously generated Jatropha hybrid genome assembly of size 265 Mbp using RJPrimers pipeline software and identified a total of 1274 TE junctions. For the predicted RJs, we designed 2868 polymerase chain reaction (PCR) based RJ markers (RJMs) flanking the junction regions. In addition to marker design, the identified RJs were utilized to detect 225,517 TEs across the genome. The different types of transposable repeat elements mainly were scattered into Retro, LTR, Copia and Gypsy categories. The efficacy of the designed markers was tested by utilizing a subset of RJMs selected randomly. We have validated 96 randomly selected RJ primers in a group of 32 J. curcas genotypes and more than 90% of the markers effectively intensified as amplicons. Of these, 10 primers were shown to be polymorphic in estimating genetic diversity among the 32 Jatropha lines. UPGMA cluster analysis revealed the formation of two clusters such as A and B exhibiting 85.5% and 87% similarity coefficient respectively. The various RJMs identified in this study could be utilized as a significant asset in Jatropha functional genomics including genome determination, mapping and marker-assisted selection.

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

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

CRISPR/Cas9 Gene Editing of Gluten in Wheat to Reduce Gluten Content and Exposure-Reviewing Methods to Screen for Coeliac Safety

Ingestion of gluten proteins (gliadins and glutenins) from wheat, barley and rye can cause coeliac disease (CD) in genetically predisposed individuals. The only remedy is a strict and lifelong gluten-free diet. There is a growing desire for coeliac-safe, whole-grain wheat-based products, as consumption of whole-grain foods reduces the risk of chronic diseases. However, due to the large number of gluten genes and the complexity of the wheat genome, wheat that is coeliac-safe but retains baking quality cannot be produced by conventional breeding alone. CD is triggered by immunogenic epitopes, notably those present in α-, γ-, and ω-gliadins. RNA interference (RNAi) silencing has been used to down-regulate gliadin families. Recently, targeted gene editing using CRISPR/Cas9 has been applied to gliadins. These methods produce offspring with silenced, deleted, and/or edited gliadins, that overall may reduce the exposure of patients to CD epitopes. Here we review methods to efficiently screen and select the lines from gliadin gene editing programs for CD epitopes at the DNA and protein level, for baking quality, and ultimately in clinical trials. The application of gene editing for the production of coeliac-safe wheat is further considered within the context of food production and in view of current national and international regulatory frameworks.

Copy Number Variation of Transposable Elements in Thinopyrum intermedium and Its Diploid Relative Species

Diploid and polyploid wild species of Triticeae have complex relationships, and the understanding of their evolution and speciation could help to increase the usability of them in wheat breeding as a source of genetic diversity. The diploid species Pseudoroegneria spicata (St), Thinopyrum bessarabicum (Jb), Dasypyrum villosum (V) derived from a hypothetical common ancestor are considered to be possible subgenome donors in hexaploid species Th. intermedium (JrJvsSt, where indices r, v, and s stand for the partial relation to the genomes of Secale, Dasypyrum, and Pseudoroegneria, respectively). We quantified 10 families of transposable elements (TEs) in P. spicata, Th. bessarabicum, D. villosum (per one genome), and Th. intermedium (per one average subgenome) using the quantitative real time PCR assay and compared their abundance within the studied genomes as well as between them. Sabrina was the most abundant among all studied elements in P. spicata, D. villosum, and Th. intermedium, and among Ty3/Gypsy elements in all studied species. Among Ty1/Copia elements, Angela-A and WIS-A showed the highest and close abundance with the exception of D. villosum, and comprised the majority of all studied elements in Th. bessarabicum. Sabrina, BAGY2, and Angela-A showed similar abundance among diploids and in Th. intermedium hexaploid; Latidu and Barbara demonstrated sharp differences between diploid genomes. The relationships between genomes of Triticeae species based on the studied TE abundance and the role of TEs in speciation and polyploidization in the light of the current phylogenetic models is discussed.

Low coverage sequencing for repetitive DNA analysis in Passiflora edulis Sims: citogenomic characterization of transposable elements and satellite DNA

BACKGROUND

The cytogenomic study of repetitive regions is fundamental for the understanding of morphofunctional mechanisms and genome evolution. Passiflora edulis a species of relevant agronomic value, this work had its genome sequenced by next generation sequencing and bioinformatics analysis performed by RepeatExplorer pipeline. The clusters allowed the identification and characterization of repetitive elements (predominant contributors to most plant genomes). The aim of this study was to identify, characterize and map the repetitive DNA of P. edulis, providing important cytogenomic markers, especially sequences associated with the centromere.

RESULTS

Three clusters of satellite DNAs (69, 118 and 207) and seven clusters of Long Terminal Repeat (LTR) retrotransposons of the superfamilies Ty1/Copy and Ty3/Gypsy and families Angela, Athila, Chromovirus and Maximus-Sire (6, 11, 36, 43, 86, 94 and 135) were characterized and analyzed. The chromosome mapping of satellite DNAs showed two hybridization sites co-located in the 5S rDNA region (PeSat_1), subterminal hybridizations (PeSat_3) and hybridization in four sites, co-located in the 45S rDNA region (PeSat_2). Most of the retroelements hybridizations showed signals scattered in the chromosomes, diverging in abundance, and only the cluster 6 presented pericentromeric regions marking. No satellite DNAs and retroelement associated with centromere was observed.

CONCLUSION

P. edulis has a highly repetitive genome, with the predominance of Ty3/Gypsy LTR retrotransposon. The satellite DNAs and LTR retrotransposon characterized are promising markers for investigation of the evolutionary patterns and genetic distinction of species and hybrids of Passiflora.

The first near-complete assembly of the hexaploid bread wheat genome, Triticum aestivum

Common bread wheat, Triticum aestivum, has one of the most complex genomes known to science, with 6 copies of each chromosome, enormous numbers of near-identical sequences scattered throughout, and an overall haploid size of more than 15 billion bases. Multiple past attempts to assemble the genome have produced assemblies that were well short of the estimated genome size. Here we report the first near-complete assembly of T. aestivum, using deep sequencing coverage from a combination of short Illumina reads and very long Pacific Biosciences reads. The final assembly contains 15 344 693 583 bases and has a weighted average (N50) contig size of 232 659 bases. This represents by far the most complete and contiguous assembly of the wheat genome to date, providing a strong foundation for future genetic studies of this important food crop. We also report how we used the recently published genome of Aegilops tauschii, the diploid ancestor of the wheat D genome, to identify 4 179 762 575 bp of T. aestivum that correspond to its D genome components.

Characterization of leaf transcriptome, development and utilization of unigenes-derived microsatellite markers in sugarcane (Saccharum sp. hybrid)

Sugarcane (Saccharum species hybrid) is the major source of sugar (> 80% sugar) in the world and is cultivated in more than 115 countries. It has recently gained attention as a source of biofuel (ethanol). Due to genomic complexity, the development of new genomic resources is imperative in understanding the gene regulation and function, and to fine tune the genetic improvement of sugarcane. In this study, a cDNA library was constructed from mature leaves so as to develop ESTs resources which were further compared with nucleotide and protein databases to explore the functional identity of sugarcane genes. The non-redundant ESTs (unigenes) were categorized into 18 metabolic functions. The major categories were bioenergetics and photosynthesis (4%), cell metabolism (5%), development related protein (3%), membrane-related, mobile genetic elements (5%), signal transduction (2%), DNA (1%), RNA (1%) and protein (2%) metabolism, other metabolic processes (3%), transcription factors (1%), transport (4%) and proteins related to stress/defense (4%). From 540 unique ESTs, 212 simple sequence repeats were identified, of which 206 were from 463 singlets and six were mined from 77 contig sequences. A total of 540 unique EST sequences were used for SSR search of which 97 (17.9%) contained specified SSR motifs, generating 212 unique SSRs. The genes characterized in this study and the EST-derived microsatellite markers identified from the cDNA library will enrich genomic resources for association- and linkage-mapping studies in sugarcane.

Fluorescence in situ hybridization karyotyping reveals the presence of two distinct genomes in the taxon Aegilops tauschii

Background

Aegilops tauschii is the donor of the bread wheat D genome. Based on spike morphology, the taxon has conventionally been subdivided into ssp. tauschii and ssp. strangulata. The present study was intended to address the poor match between this whole plant morphology-based subdivision and genetic relationships inferred from genotyping by fluorescence in situ hybridization karyotyping a set of 31 Ae. tauschii accessions.

Results

The distribution of sites hybridizing to the two probes oligo-pTa-535 and (CTT)₁₀ split the Ae. tauschii accessions into two clades, designated D^t and D^s, which corresponded perfectly with a previously assembled phylogeny based on marker genotype. The D^t cluster was populated exclusively by ssp. tauschii accessions, while the D^s cluster harbored both ssp. strangulata and morphologically intermediate accessions. As a result, it is proposed that Ae. tauschii ssp. tauschii is restricted to carriers of the D^t karyotype: their spikelets are regularly spaced along the rachis, at least in the central portion of their spike. Accessions classified as Ae. tauschii ssp. strangulata carry the D^s karyotype; their spikelets are irregularly spaced. Based on this criterion, forms formerly classified as ssp. tauschii var. meyeri have been re-designated ssp. strangulata var. meyeri.

Conclusions

According to the reworking of the taxon, the bread wheat D genome was most probably donated by ssp. strangulata var. meyeri. Chromosomal differentiation reveals intra-species taxon of Ae. tauschii. Ae. tauschii ssp. tauschii has more distant relationship with breed wheat than ssp. strangulata and can be used for breeding improving effectively.

PlaNC-TE: a comprehensive knowledgebase of non-coding RNAs and transposable elements in plants

Transposable elements (TEs) play an essential role in the genetic variability of eukaryotic species. In plants, they may comprise up to 90% of the total genome. Non-coding RNAs (ncRNAs) are known to control gene expression and regulation. Although the relationship between ncRNAs and TEs is known, obtaining the organized data for sequenced genomes is not straightforward. In this study, we describe the PlaNC-TE (http://planc-te.cp.utfpr.edu.br), a user-friendly portal harboring a knowledgebase created by integrating and analysing plant ncRNA-TE data. We identified a total of 14 350 overlaps between ncRNAs and TEs in 40 plant genomes. The database allows users to browse, search and download all ncRNA and TE data analysed. Overall, PlaNC-TE not only organizes data and provides insights about the relationship between ncRNA and TEs in plants but also helps improve genome annotation strategies. Moreover, this is the first database to provide resources to broadly investigate functions and mechanisms involving TEs and ncRNAs in plants.

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross-species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next-generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross-genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter- and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.

Translational genomics of grain size regulation in wheat

Identifying and mapping grain size candidate genes in the wheat genome greatly empowers reverse genetics approaches to improve grain yield potential of wheat. Grain size (GS) or grain weight is believed to be a major driving force for further improvement of wheat yield. Although the large, polyploid genome of wheat poses an obstacle to identifying GS determinants using map-based cloning, a translational genomics approach using GS regulators identified in the model plants rice and Arabidopsis as candidate genes appears to be effective and supports a hypothesis that a conserved genetic network regulates GS in rice and wheat. In this review, we summarize the progress in the studies on GS in the model plants and wheat and identify 45 GS candidate loci in the wheat genome. In silico mapping of these GS loci in the diploid wheat and barley genomes showed (1) several gene families amplified in the wheat lineage, (2) a significant number of the GS genes located in the proximal regions surrounding the centromeres, and (3) more than half of candidate genes to be negative regulators, or their expression negatively related by microRNAs. Identifying and mapping the wheat GS gene homologs will not only facilitate candidate gene analysis, but also open the door to improving wheat yield using reverse genetics approaches by mining desired alleles in landraces and wild ancestors and to developing novel germplasm by TILLING and genome editing technologies.

Diversity, distribution and dynamics of full-length Copia and Gypsy LTR retroelements in Solanum lycopersicum

Transposable elements are the most abundant components of plant genomes and can dramatically induce genetic changes and impact genome evolution. In the recently sequenced genome of tomato (Solanum lycopersicum), the estimated fraction of elements corresponding to retrotransposons is nearly 62%. Given that tomato is one of the most important vegetable crop cultivated and consumed worldwide, understanding retrotransposon dynamics can provide insight into its evolution and domestication processes. In this study, we performed a genome-wide in silico search of full-length LTR retroelements in the tomato nuclear genome and annotated 736 full-length Gypsy and Copia retroelements. The dispersion level across the 12 chromosomes, the diversity and tissue-specific expression of those elements were estimated. Phylogenetic analysis based on the retrotranscriptase region revealed the presence of 12 major lineages of LTR retroelements in the tomato genome. We identified 97 families, of which 77 and 20 belong to the superfamilies Copia and Gypsy, respectively. Each retroelement family was characterized according to their element size, relative frequencies and insertion time. These analyses represent a valuable resource for comparative genomics within the Solanaceae, transposon-tagging and for the design of cultivar-specific molecular markers in tomato.

Transcriptional activity of PIF and Pong-like Class II transposable elements in Triticeae

Background

Transposable elements are major contributors to genome size and variability, accounting for approximately 70–80% of the maize, barley, and wheat genomes. PIF and Pong-like elements belong to two closely-related element families within the PIF/Harbinger superfamily of Class II (DNA) transposons. Both elements contain two open reading frames; one encodes a transposase (ORF2) that catalyzes transposition of the functional elements and their related non-autonomous elements, while the function of the second is still debated. In this work, we surveyed for PIF- and Pong-related transcriptional activity in 13 diploid Triticeae species, all of which have been previously shown to harbor extensive within-genome diversity of both groups of elements.

Results

The results revealed that PIF elements have considerable transcriptional activity in Triticeae, suggesting that they can escape the initial levels of plant cell control and are regulated at the post-transcriptional level. Phylogenetic analysis of 156 PIF cDNA transposase fragments along with 240 genomic partial transposase sequences showed that most, if not all, PIF clades are transcriptionally competent, and that multiple transposases coexisting within a single genome have the potential to act simultaneously. In contrast, we did not detect any transcriptional activity of Pong elements in any sample.

Conclusions

The lack of Pong element transcription shows that even closely related transposon families can exhibit wide variation in their transposase transcriptional activity within the same genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1028-6) contains supplementary material, which is available to authorized users.

Useful parasites: the evolutionary biology and biotechnology applications of transposable elements

Transposable elements usually comprise the most abundant nongenic fraction of eukaryotic genomes. Because of their capacity to selfreplicate and to induce a wide range of mutations, transposable elements have long been considered as 'parasitic' or 'selfish'. Today, we recognize that the findings about genomic changes affected by transposable elements have considerably altered our view of the ways in which genomes evolve and work. Numerous studies have provided evidences that mobile elements have the potential to act as agents of evolution by increasing, rearranging and diversifying the genetic repertoire of their hosts. With large-scale sequencing becoming increasingly available, more and more scientists come across transposable element sequences in their data. I will provide examples that transposable elements, although having signatures of 'selfish' DNA, play a significant biological role in the maintainance of genome integrity and providing novel regulatoty networks. These features, along with the transpositional and mutagenic capacity to produce a raw genetic diversity, make the genome mobile fraction, a key player in species adaptation and microevolution. The last but not least, transposable elements stand as informative DNA markers that may complement other conventional DNA markers. Altogether, transposable elements represent a promising, but still largely unexplored research niche and deserve to be included into the agenda of molecular ecologists, evolutionary geneticists, conservation biologists and plant breeders.

Pseudogenes and Their Genome-Wide Prediction in Plants

Pseudogenes are paralogs generated from ancestral functional genes (parents) during genome evolution, which contain critical defects in their sequences, such as lacking a promoter, having a premature stop codon or frameshift mutations. Generally, pseudogenes are functionless, but recent evidence demonstrates that some of them have potential roles in regulation. The majority of pseudogenes are generated from functional progenitor genes either by gene duplication (duplicated pseudogenes) or retro-transposition (processed pseudogenes). Pseudogenes are primarily identified by comparison to their parent genes. Bioinformatics tools for pseudogene prediction have been developed, among which PseudoPipe, PSF and Shiu's pipeline are publicly available. We compared these three tools using the well-annotated Arabidopsis thaliana genome and its known 924 pseudogenes as a test data set. PseudoPipe and Shiu's pipeline identified ~80% of A. thaliana pseudogenes, of which 94% were shared, while PSF failed to generate adequate results. A need for improvement of the bioinformatics tools for pseudogene prediction accuracy in plant genomes was thus identified, with the ultimate goal of improving the quality of genome annotation in plants.

New molecular markers and cytogenetic probes enable chromosome identification of wheat-Thinopyrum intermedium introgression lines for improving protein and gluten contents

New molecular markers were developed for targeting Thinopyrum intermedium 1St#2 chromosome, and novel FISH probe representing the terminal repeats was produced for identification of Thinopyrum chromosomes. Thinopyrum intermedium has been used as a valuable resource for improving the disease resistance and yield potential of wheat. A wheat-Th. intermedium ssp. trichophorum chromosome 1St#2 substitution and translocation has displayed superior grain protein and wet gluten content. With the aim to develop a number of chromosome 1St#2 specific molecular and cytogenetic markers, a high throughput, low-cost specific-locus amplified fragment sequencing (SLAF-seq) technology was used to compare the sequences between a wheat-Thinopyrum 1St#2 (1D) substitution and the related species Pseudoroegneria spicata (St genome, 2n = 14). A total of 5142 polymorphic fragments were analyzed and 359 different SLAF markers for 1St#2 were predicted. Thirty-seven specific molecular markers were validated by PCR from 50 randomly selected SLAFs. Meanwhile, the distribution of transposable elements (TEs) at the family level between wheat and St genomes was compared using the SLAFs. A new oligo-nucleotide probe named Oligo-pSt122 from high SLAF reads was produced for fluorescence in situ hybridization (FISH), and was observed to hybridize to the terminal region of 1St#L and also onto the terminal heterochromatic region of Th. intermedium genomes. The genome-wide markers and repetitive based probe Oligo-pSt122 will be valuable for identifying Thinopyrum chromosome segments in wheat backgrounds.

detectMITE: A novel approach to detect miniature inverted repeat transposable elements in genomes

Miniature inverted repeat transposable elements (MITEs) are prevalent in eukaryotic genomes, including plants and animals. Classified as a type of non-autonomous DNA transposable elements, they play important roles in genome organization and evolution. Comprehensive and accurate genome-wide detection of MITEs in various eukaryotic genomes can improve our understanding of their origins, transposition processes, regulatory mechanisms, and biological relevance with regard to gene structures, expression, and regulation. In this paper, we present a new MATLAB-based program called detectMITE that employs a novel numeric calculation algorithm to replace conventional string matching algorithms in MITE detection, adopts the Lempel-Ziv complexity algorithm to filter out MITE candidates with low complexity, and utilizes the powerful clustering program CD-HIT to cluster similar MITEs into MITE families. Using the rice genome as test data, we found that detectMITE can more accurately, comprehensively, and efficiently detect MITEs on a genome-wide scale than other popular MITE detection tools. Through comparison with the potential MITEs annotated in Repbase, the widely used eukaryotic repeat database, detectMITE has been shown to find known and novel MITEs with a complete structure and full-length copies in the genome. detectMITE is an open source tool (https://sourceforge.net/projects/detectmite).

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.

The Role of Vertical and Horizontal Transfer in the Evolutionary Dynamics of PIF-Like Transposable Elements in Triticeae

PIF-like transposable elements are members of the PIF/Harbinger superfamily of DNA transposons found in the genomes of many plants, animals, and fungi. The evolution of the gene that encodes the transposase responsible for mobilizing PIF-like elements has been studied in both plants and animals, but the elements' history in flowering plants remains poorly known. In this work, we describe the phylogenetic distribution and evolution of PIF-like elements in the genomes of 21 diploid species from the wheat tribe, Triticeae, and we present the first convincing evidence of horizontal transfer of PIF elements in plant genomes. A phylogenetic analysis of 240 PIF sequences based on the conserved region of the transposase domain revealed at least four main transposase lineages. Their complex evolutionary history can be best explained by a combination of vertical transmission with differential evolutionary success among lineages, and occasional horizontal transfer between phylogenetically distant Triticeae genera. In addition, we identified 127 potentially functional transposase sequences indicating possible recent activity of PIF.

Sequencing chromosome 5D of Aegilops tauschii and comparison with its allopolyploid descendant bread wheat (Triticum aestivum)

Flow cytometric sorting of individual chromosomes and chromosome-based sequencing reduces the complexity of large, repetitive Triticeae genomes. We flow-sorted chromosome 5D of Aegilops tauschii, the D genome donor of bread wheat and sequenced it by Roche 454 GS FLX platform to approximately 2.2x coverage. Repetitive sequences represent 81.09% of the survey sequences of this chromosome, and Class I retroelements are the prominent type, with a particular abundance of LTR/Gypsy superfamily. Nonrepetitive sequences were assembled to cover 17.76% of the total chromosome regions. Up to 6188 nonrepetitive gene loci were predicted to be encoded by the 5D chromosome. The numbers and chromosomal distribution patterns of tRNA genes suggest abundance in tRNA(L) (ys) and tRNA(M) (et) species, while the nonrepetitive assembly reveals tRNA(A) (la) species as the most abundant type. A comparative analysis of the genomic sequences of bread wheat and Aegilops chromosome 5D indicates conservation of gene content. Orthologous unique genes, matching Aegilops 5D sequences, numbered 3730 in barley, 5063 in Brachypodium, 4872 in sorghum and 4209 in rice. In this study, we provide a chromosome-specific view into the structure and organization of the 5D chromosome of Ae. tauschii, the D genome ancestor of bread wheat. This study contributes to our understanding of the chromosome-level evolution of the wheat genome and presents a valuable resource in wheat genomics due to the recent hybridization of Ae. tauschii genome with its tetraploid ancestor.

Diversity, abundance, and evolutionary dynamics of Pong-like transposable elements in Triticeae

Pong-like elements are members of the PIF/Harbinger superfamily of DNA transposons that has been described in many plants, animals, and fungi. Most Pong elements contain two open reading frames (ORFs). One encodes a transposase (ORF2) that catalyzes transposition of Pong and related non-autonomous elements, while the function of the second is unknown. Little is known about the evolutionary history of Pong elements in flowering plants. In this work, we present the first comprehensive analysis of the diversity, abundance, and evolution of the Pong-like transposase gene in the genomes of 21 diploid species from the wheat tribe, Triticeae, and we present the first convincing evidence of horizontal transfer of nuclear-encoded Pong elements in any organism. A phylogenetic analysis of nearly 300 Pong sequences based on a conserved region of the transposase domain revealed a complex evolutionary history of Pong elements that can be best explained by ancestral polymorphism, followed by differential evolutionary success of some transposase lineages, and by occasional horizontal transfer between phylogenetically distant genera. In addition, we used transposon display to estimate the abundance of the transposase gene within Triticeae genomes, and our results revealed varying levels of Pong proliferation, with numbers of transposase copies ranging from 22 to 92. Comparisons of Pong transposase abundance to flow cytometry estimates of genome size revealed that larger Triticeae genome size was not correlated with transposase abundance.

Genome size stability despite high chromosome number variation in Carex gr. laevigata

PREMISE OF THE STUDY

In organisms with holocentric chromosomes like Carex species, chromosome number evolution has been hypothesized to be a result of fission, fusion, and/or translocation events. Negative, positive, or the absence of correlations have been found between chromosome number and genome size in Carex.

METHODS

Using the inferred diploid chromosome number and 80 genome size measurements from 26 individuals and 20 populations of Carex gr. laevigata, we tested the null hypothesis of chromosome number evolution by duplication and deletion of whole chromosomes.

KEY RESULTS

Our results show a significant positive correlation between genome size and chromosome number, but the slope of such correlation supports the hypothesis of proliferation and removal of repetitive DNA fragments to explain genome size variation rather than duplication and deletion of whole chromosomes.

CONCLUSIONS

Our results refine the theory of the holokinetic drive: this mechanism is proposed to facilitate repetitive DNA removal (or any segmental deletion) when smaller homologous chromosomes are preferentially inherited, or repetitive DNA proliferation (or any segmental duplication) when larger homologs are preferred. This study sheds light on how karyotype evolution plays an important role in the diversification of the species of the genus Carex.

Pivoting from Arabidopsis to wheat to understand how agricultural plants integrate responses to biotic stress

In this review, we argue for a research initiative on wheat's responses to biotic stress. One goal is to begin a conversation between the disparate communities of plant pathology and entomology. Another is to understand how responses to a variety of agents of biotic stress are integrated in an important crop. We propose gene-for-gene interactions as the focus of the research initiative. On the parasite's side is an Avirulence (Avr) gene that encodes one of the many effector proteins the parasite applies to the plant to assist with colonization. On the plant's side is a Resistance (R) gene that mediates a surveillance system that detects the Avr protein directly or indirectly and triggers effector-triggered plant immunity. Even though arthropods are responsible for a significant proportion of plant biotic stress, they have not been integrated into important models of plant immunity that come from plant pathology. A roadblock has been the absence of molecular evidence for arthropod Avr effectors. Thirty years after this evidence was discovered in a plant pathogen, there is now evidence for arthropods with the cloning of the Hessian fly's vH13 Avr gene. After reviewing the two models of plant immunity, we discuss how arthropods could be incorporated. We end by showing features that make wheat an interesting system for plant immunity, including 479 resistance genes known from agriculture that target viruses, bacteria, fungi, nematodes, insects, and mites. It is not likely that humans will be subsisting on Arabidopsis in the year 2050. It is time to start understanding how agricultural plants integrate responses to biotic stress.

Contrasting evolutionary trajectories of multiple retrotransposons following independent allopolyploidy in wild wheats

Transposable elements (TEs) are expectedly central to genome evolution. To assess the impact of TEs in driving genome turnover, we used allopolyploid genomes, showing considerable deviation from the predicted additivity of their diploid progenitors and thus having undergone major restructuring. Genome survey sequencing was used to select 17 putatively active families of long terminal repeat retrotransposons. Genome-wide TE insertions were genotyped with sequence-specific amplified polymorphism (SSAP) in diploid progenitors and their derived polyploids, and compared with changes in random sequences to assess restructuring of four independent Aegilops allotetraploid genomes. Generally, TEs with different evolutionary trajectories from those of random sequences were identified. Thus, TEs presented family-specific and species-specific dynamics following polyploidy, as illustrated by Sabine showing proliferation in particular polyploids, but massive elimination in others. Contrasting with that, only a few families (BARE1 and Romani) showed proliferation in all polyploids. Overall, TE divergence between progenitors was strongly correlated with the degree of restructuring in polyploid TE fractions. TE families present evolutionary trajectories that are decoupled from genome-wide changes after allopolyploidy and have a pervasive impact on their restructuring.

Variation, evolution, and correlation analysis of C+G content and genome or chromosome size in different kingdoms and phyla

C+G content (GC content or G+C content) is known to be correlated with genome/chromosome size in bacteria but the relationship for other kingdoms remains unclear. This study analyzed genome size, chromosome size, and base composition in most of the available sequenced genomes in various kingdoms. Genome size tends to increase during evolution in plants and animals, and the same is likely true for bacteria. The genomic C+G contents were found to vary greatly in microorganisms but were quite similar within each animal or plant subkingdom. In animals and plants, the C+G contents are ranked as follows: monocot plants>mammals>non-mammalian animals>dicot plants. The variation in C+G content between chromosomes within species is greater in animals than in plants. The correlation between average chromosome C+G content and chromosome length was found to be positive in Proteobacteria, Actinobacteria (but not in other analyzed bacterial phyla), Ascomycota fungi, and likely also in some plants; negative in some animals, insignificant in two protist phyla, and likely very weak in Archaea. Clearly, correlations between C+G content and chromosome size can be positive, negative, or not significant depending on the kingdoms/groups or species. Different phyla or species exhibit different patterns of correlation between chromosome-size and C+G content. Most chromosomes within a species have a similar pattern of variation in C+G content but outliers are common. The data presented in this study suggest that the C+G content is under genetic control by both trans- and cis- factors and that the correlation between C+G content and chromosome length can be positive, negative, or not significant in different phyla.

Differentially expressed genes in response to gamma-irradiation during the vegetative stage in Arabidopsis thaliana

Biochemical and physiological processes in plants are affected by gamma-irradiation, which causes significant changes in gene transcripts and expression. To identify the differentially expressed Arabidopsis genes in response to gamma-irradiation, we performed a microarray analysis with rosette leaves during the vegetative stage. Arabidopsis plants were exposed to a wide spectrum doses of gamma ray (100, 200, 300, 400, 800, 1200, 1600 or 2000 Gy) for 24 h. At the dose range from 100 to 400 Gy, irradiated plants were found to be shorter than controls after 8 days of irradiation, while doses over 800 Gy caused severe growth retardation. Therefore, 100 and 800 Gy were selected as adequate doses for microarray analysis to identify differentially expressed genes. Among the 20,993 genes used as microarray probes, a total number of 496 and 1,042 genes were up-regulated and down-regulated by gamma-irradiation, respectively (P < 0.05). We identified the characteristics of the genes that were up-and down-regulated fourfold higher genes by gamma irradiation according to The arabidopsis information resource gene ontology. To confirm the microarray results, we performed a northern blot and quantitative real-time PCR with several selected genes that had a large difference in expression after irradiation. In particular, genes associated with lipid transfer proteins, histones and transposons were down-regulated by 100 and/or 800 Gy of gamma irradiation. The expression patterns of selected genes were generally in agreement with the microarray results, although there were quantitative differences in the expression levels.

Organization and evolution of transposable elements along the bread wheat chromosome 3B

BACKGROUND

The 17 Gb bread wheat genome has massively expanded through the proliferation of transposable elements (TEs) and two recent rounds of polyploidization. The assembly of a 774 Mb reference sequence of wheat chromosome 3B provided us with the opportunity to explore the impact of TEs on the complex wheat genome structure and evolution at a resolution and scale not reached so far.

RESULTS

We develop an automated workflow, CLARI-TE, for TE modeling in complex genomes. We delineate precisely 56,488 intact and 196,391 fragmented TEs along the 3B pseudomolecule, accounting for 85% of the sequence, and reconstruct 30,199 nested insertions. TEs have been mostly silent for the last one million years, and the 3B chromosome has been shaped by a succession of bursts that occurred between 1 to 3 million years ago. Accelerated TE elimination in the high-recombination distal regions is a driving force towards chromosome partitioning. CACTAs overrepresented in the high-recombination distal regions are significantly associated with recently duplicated genes. In addition, we identify 140 CACTA-mediated gene capture events with 17 genes potentially created by exon shuffling and show that 19 captured genes are transcribed and under selection pressure, suggesting the important role of CACTAs in the recent wheat adaptation.

CONCLUSION

Accurate TE modeling uncovers the dynamics of TEs in a highly complex and polyploid genome. It provides novel insights into chromosome partitioning and highlights the role of CACTA transposons in the high level of gene duplication in wheat.

Roles of dehydrin genes in wheat tolerance to drought stress

Physiological parameters and expression levels of drought related genes were analyzed in early vegetative stage of two bread wheat cultivars (Sids and Gmiza) differ in drought tolerance capacity. Both cultivars were imposed to gradual water depletion started on day 17 till day 32 after sowing. Sids, the more tolerant cultivar to drought showed higher fresh and dry weights than the drought sensitive genotype, Gmiza. Under water stress, Sids had higher membrane stability index (MSI), lower accumulated H₂O₂ and higher activity of the antioxidant enzymes; catalase (CAT), guaiacol peroxidase (GPX), ascorbate peroxidase (APX) and superoxide dismutase (SOD) than Gmiza. On the other hand, the differential expression patterns of the genes dhn, wcor and dreb were observed due to water deficit intensity according to cultivar’s tolerance to drought. The DNA sequence alignment of dun showed high similarity of about 80–92% identities with other related plants. The most striking overall observed trend was the highly induction in the expression of dun, wcor and dreb in leaves of the tolerant genotype, Sids under severe water stress.

Unique and conserved microRNAs in wheat chromosome 5D revealed by next-generation sequencing

MicroRNAs are a class of short, non-coding, single-stranded RNAs that act as post-transcriptional regulators in gene expression. miRNA analysis of Triticum aestivum chromosome 5D was performed on 454 GS FLX Titanium sequences of flow-sorted chromosome 5D with a total of 3,208,630 good quality reads representing 1.34x and 1.61x coverage of the short (5DS) and long (5DL) arms of the chromosome respectively. In silico and structural analyses revealed a total of 55 miRNAs; 48 and 42 miRNAs were found to be present on 5DL and 5DS respectively, of which 35 were common to both chromosome arms, while 13 miRNAs were specific to 5DL and 7 miRNAs were specific to 5DS. In total, 14 of the predicted miRNAs were identified in wheat for the first time. Representation (the copy number of each miRNA) was also found to be higher in 5DL (1,949) compared to 5DS (1,191). Targets were predicted for each miRNA, while expression analysis gave evidence of expression for 6 out of 55 miRNAs. Occurrences of the same miRNAs were also found in Brachypodium distachyon and Oryza sativa genome sequences to identify syntenic miRNA coding sequences. Based on this analysis, two other miRNAs: miR1133 and miR167 were detected in B. distachyon syntenic region of wheat 5DS. Five of the predicted miRNA coding regions (miR6220, miR5070, miR169, miR5085, miR2118) were experimentally verified to be located to the 5D chromosome and three of them : miR2118, miR169 and miR5085, were shown to be 5D specific. Furthermore miR2118 was shown to be expressed in Chinese Spring adult leaves. miRNA genes identified in this study will expand our understanding of gene regulation in bread wheat.

Retrotransposons represent the most labile fraction for genomic rearrangements in polyploid plant species

Understanding how increased genome size and diversity within polyploid genomes impacts plant evolution and breeding continues to be challenging. Although historical studies by McClintock suggested the importance of transposable elements mediated by polyploidisation on genomic changes, data from plant crosses remain scarce. Despite the absence of a conclusive proof regarding autonomous retrotransposon movement in synthetic allopolyploids, the transposition of retrotransposons and their ubiquitous dispersion in all plant species might explain the positive correlation between the genome size of plants and the prevalence of retrotransposons. Here, we address polyploidisation-mediated rearrangements of retrotransposon-associated sequences and discuss a tendency for a preferential restructuring of large ancestral genomes after polyploidisation. A comparative analysis of the frequency of modifications of retrotransposon-associated sequences in synthetic polyploids with marked differences in genome sizes is presented. Such analyses suggest the absence of a significant difference in the rates of rearrangements despite vast dissimilarities in the retrotransposon copy number between species, which emphasises the high plasticity of this genomic feature. See also the sister article focusing on animals by Arkhipova and Rodriguez in this themed issue.

A 4-gigabase physical map unlocks the structure and evolution of the complex genome of Aegilops tauschii, the wheat D-genome progenitor

The current limitations in genome sequencing technology require the construction of physical maps for high-quality draft sequences of large plant genomes, such as that of Aegilops tauschii, the wheat D-genome progenitor. To construct a physical map of the Ae. tauschii genome, we fingerprinted 461,706 bacterial artificial chromosome clones, assembled contigs, designed a 10K Ae. tauschii Infinium SNP array, constructed a 7,185-marker genetic map, and anchored on the map contigs totaling 4.03 Gb. Using whole genome shotgun reads, we extended the SNP marker sequences and found 17,093 genes and gene fragments. We showed that collinearity of the Ae. tauschii genes with Brachypodium distachyon, rice, and sorghum decreased with phylogenetic distance and that structural genome evolution rates have been high across all investigated lineages in subfamily Pooideae, including that of Brachypodieae. We obtained additional information about the evolution of the seven Triticeae chromosomes from 12 ancestral chromosomes and uncovered a pattern of centromere inactivation accompanying nested chromosome insertions in grasses. We showed that the density of noncollinear genes along the Ae. tauschii chromosomes positively correlates with recombination rates, suggested a cause, and showed that new genes, exemplified by disease resistance genes, are preferentially located in high-recombination chromosome regions.

Integration of mate pair sequences to improve shotgun assemblies of flow-sorted chromosome arms of hexaploid wheat

BACKGROUND

The assembly of the bread wheat genome sequence is challenging due to allohexaploidy and extreme repeat content (>80%). Isolation of single chromosome arms by flow sorting can be used to overcome the polyploidy problem, but the repeat content cause extreme assembly fragmentation even at a single chromosome level. Long jump paired sequencing data (mate pairs) can help reduce assembly fragmentation by joining multiple contigs into single scaffolds. The aim of this work was to assess how mate pair data generated from multiple displacement amplified DNA of flow-sorted chromosomes affect assembly fragmentation of shotgun assemblies of the wheat chromosomes.

RESULTS

Three mate pair (MP) libraries (2 Kb, 3 Kb, and 5 Kb) were sequenced to a total coverage of 89x and 64x for the short and long arm of chromosome 7B, respectively. Scaffolding using SSPACE improved the 7B assembly contiguity and decreased gene space fragmentation, but the degree of improvement was greatly affected by scaffolding stringency applied. At the lowest stringency the assembly N50 increased by ~7 fold, while at the highest stringency N50 was only increased by ~1.5 fold. Furthermore, a strong positive correlation between estimated scaffold reliability and scaffold assembly stringency was observed. A 7BS scaffold assembly with reduced MP coverage proved that assembly contiguity was affected only to a small degree down to ~50% of the original coverage.

CONCLUSION

The effect of MP data integration into pair end shotgun assemblies of wheat chromosome was moderate; possibly due to poor contig assembly contiguity, the extreme repeat content of wheat, and the use of amplified chromosomal DNA for MP library construction.