Polyploidization as a retraction force in plant genome evolution: sequence rearrangements in triticale

Background

Polyploidization is a major evolutionary process in plants where hybridization and chromosome doubling induce enormous genomic stress and can generate genetic and epigenetic modifications. However, proper evaluation of DNA sequence restructuring events and the precise characterization of sequences involved are still sparse.

Methodology/Principal Findings

Inter Retrotransposons Amplified Polymorphism (IRAP), Retrotransposons Microsatellite Amplified Polymorphism (REMAP) and Inter Simple Sequence Repeat (ISSR) largely confirmed the absence of any intraspecific variation in wheat, rye and triticale. The comparative analysis of banding profiles between wheat and rye inbred lines revealed 34% of monomorphic (common to both parental species) bands for the ten different primer combinations used. The analysis of triticale plants uncovered nearly 51% of rearranged bands in the polyploid, being the majority of these modifications, due to the loss of rye bands (83%). Sequence analysis of rye fragments absent in triticale revealed for instance homology with hydroxyproline-rich glycoproteins (HRGP), a protein that belongs to a major family of inducible defence response proteins. Conversely, a wheat-specific band absent in triticale comprises a nested structure of copia-like retrotransposons elements, namely Claudia and Barbara. Sequencing of a polyploid-specific band (absent in both parents) revealed a microsatellite related sequence. Cytological studies using Fluorescent In Situ Hybridization (FISH) with REMAP products revealed a widespread distribution of retrotransposon and/or microsatellite flanking sequences on rye chromosomes, with a preferential accumulation in heterochromatic sub-telomeric domains.

Conclusions/Significance

Here, we used PCR-based molecular marker techniques involving retrotransposons and microsatellites to uncover polyploidization induced genetic restructuring in triticale. Sequence analysis of rearranged genomic fragments either from rye or wheat origin showed these to be retrotransposon-related as well as coding sequences. Further FISH analysis revealed possible chromosome hotspots for sequence rearrangements. The role of chromatin condensation on the origin of genomic rearrangements mediated by polyploidization in triticale is also discussed.     

Spontaneous and Divergent Hexaploid Triticales Derived from Common Wheat × Rye by Complete Elimination of D-Genome Chromosomes

Background

Hexaploid triticale could be either synthesized by crossing tetraploid wheat with rye, or developed by crossing hexaploid wheat with a hexaploid triticale or an octoploid triticale.

Methodology/Principal Findings

Here two hexaploid triticales with great morphologic divergence derived from common wheat cultivar M8003 (Triticum aestivum L.) × Austrian rye (Secale cereale L.) were reported, exhibiting high resistance for powdery mildew and stripe rust and potential for wheat improvement. Sequential fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH) karyotyping revealed that D-genome chromosomes were completely eliminated and the whole A-genome, B-genome and R-genome chromosomes were retained in both lines. Furthermore, plentiful alterations of wheat chromosomes including 5A and 7B were detected in both triticales and additionally altered 5B, 7A chromosome and restructured chromosome 2A was assayed in N9116H and N9116M, respectively, even after selfing for several decades. Besides, meiotic asynchrony was displayed and a variety of storage protein variations were assayed, especially in the HMW/LMW-GS region and secalins region in both triticales.

Conclusion

This study confirms that whole D-genome chromosomes could be preferentially eliminated in the hybrid of common wheat × rye, “genome shock” was accompanying the allopolyploidization of nascent triticales, and great morphologic divergence might result from the genetic variations. Moreover, new hexaploid triticale lines contributing potential resistance resources for wheat improvement were produced.     

Characterization of the caleosin gene family in the Triticeae

Background

The caleosin genes encode proteins with a single conserved EF hand calcium-binding domain and comprise small gene families found in a wide range of plant species. Some members of the gene family have been shown to be upregulated by environmental stresses including low water availability and high salinity. Caleosin 3 from wheat has been shown to interact with the α-subunit of the heterotrimeric G proteins, and to act as a GTPase activating protein (GAP). This study characterizes the size and diversity of the gene family in wheat and related species and characterizes the differential tissue-specific expression of members of the gene family.

Results

A total of 34 gene family members that belong to eleven paralogous groups of caleosins were identified in the hexaploid bread wheat, T. aestivum. Each group was represented by three homeologous copies of the gene located on corresponding homeologous chromosomes, except the caleosin 10, which has four gene copies. Ten gene family members were identified in diploid barley, Hordeum vulgare, and in rye, Secale cereale, seven in Brachypodium distachyon, and six in rice, Oryza sativa. The analysis of gene expression was assayed in triticale and rye by RNA-Seq analysis of 454 sequence sets and members of the gene family were found to have diverse patterns of gene expression in the different tissues that were sampled in rye and in triticale, the hybrid hexaploid species derived from wheat and rye. Expression of the gene family in wheat and barley was also previously determined by microarray analysis, and changes in expression during development and in response to environmental stresses are presented.

Conclusions

The caleosin gene family had a greater degree of expansion in the Triticeae than in the other monocot species, Brachypodium and rice. The prior implication of one member of the gene family in the stress response and heterotrimeric G protein signaling, points to the potential importance of the caleosin gene family. The complexity of the family and differential expression in various tissues and under conditions of abiotic stress suggests the possibility that caleosin family members may play diverse roles in signaling and development that warrants further investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-239) contains supplementary material, which is available to authorized users.     

A high density consensus map of rye (Secale cereale L.) based on DArT markers

Background

Rye (Secale cereale L.) is an economically important crop, exhibiting unique features such as outstanding resistance to biotic and abiotic stresses and high nutrient use efficiency. This species presents a challenge to geneticists and breeders due to its large genome containing a high proportion of repetitive sequences, self incompatibility, severe inbreeding depression and tissue culture recalcitrance. The genomic resources currently available for rye are underdeveloped in comparison with other crops of similar economic importance. The aim of this study was to create a highly saturated, multilocus linkage map of rye via consensus mapping, based on Diversity Arrays Technology (DArT) markers.

Methodology/Principal Findings

Recombinant inbred lines (RILs) from 5 populations (564 in total) were genotyped using DArT markers and subjected to linkage analysis using Join Map 4.0 and Multipoint Consensus 2.2 software. A consensus map was constructed using a total of 9703 segregating markers. The average chromosome map length ranged from 199.9 cM (2R) to 251.4 cM (4R) and the average map density was 1.1 cM. The integrated map comprised 4048 loci with the number of markers per chromosome ranging from 454 for 7R to 805 for 4R. In comparison with previously published studies on rye, this represents an eight-fold increase in the number of loci placed on a consensus map and a more than two-fold increase in the number of genetically mapped DArT markers.

Conclusions/Significance

Through the careful choice of marker type, mapping populations and the use of software packages implementing powerful algorithms for map order optimization, we produced a valuable resource for rye and triticale genomics and breeding, which provides an excellent starting point for more in-depth studies on rye genome organization.     

Next generation sequencing provides rapid access to the genome of Puccinia striiformis f. sp. tritici, the causal agent of wheat stripe rust

Background

The wheat stripe rust fungus (Puccinia striiformis f. sp. tritici, PST) is responsible for significant yield losses in wheat production worldwide. In spite of its economic importance, the PST genomic sequence is not currently available. Fortunately Next Generation Sequencing (NGS) has radically improved sequencing speed and efficiency with a great reduction in costs compared to traditional sequencing technologies. We used Illumina sequencing to rapidly access the genomic sequence of the highly virulent PST race 130 (PST-130).

Methodology/Principal Findings

We obtained nearly 80 million high quality paired-end reads (>50x coverage) that were assembled into 29,178 contigs (64.8 Mb), which provide an estimated coverage of at least 88% of the PST genes and are available through GenBank. Extensive micro-synteny with the Puccinia graminis f. sp. tritici (PGTG) genome and high sequence similarity with annotated PGTG genes support the quality of the PST-130 contigs. We characterized the transposable elements present in the PST-130 contigs and using an ab initio gene prediction program we identified and tentatively annotated 22,815 putative coding sequences. We provide examples on the use of comparative approaches to improve gene annotation for both PST and PGTG and to identify candidate effectors. Finally, the assembled contigs provided an inventory of PST repetitive elements, which were annotated and deposited in Repbase.

Conclusions/Significance

The assembly of the PST-130 genome and the predicted proteins provide useful resources to rapidly identify and clone PST genes and their regulatory regions. Although the automatic gene prediction has limitations, we show that a comparative genomics approach using multiple rust species can greatly improve the quality of gene annotation in these species. The PST-130 sequence will also be useful for comparative studies within PST as more races are sequenced. This study illustrates the power of NGS for rapid and efficient access to genomic sequence in non-model organisms.     

The physical map of wheat chromosome 5DS revealed gene duplications and small rearrangements

Background

The substantially large bread wheat genome, organized into highly similar three sub-genomes, renders genomic research challenging. The construction of BAC-based physical maps of individual chromosomes reduces the complexity of this allohexaploid genome, enables elucidation of gene space and evolutionary relationships, provides tools for map-based cloning, and serves as a framework for reference sequencing efforts. In this study, we constructed the first comprehensive physical map of wheat chromosome arm 5DS, thereby exploring its gene space organization and evolution.

Results

The physical map of 5DS was comprised of 164 contigs, of which 45 were organized into 21 supercontigs, covering 176 Mb with an N50 value of 2,173 kb. Fifty-eight of the contigs were larger than 1 Mb, with the largest contig spanning 6,649 kb. A total of 1,864 molecular markers were assigned to the map at a density of 10.5 markers/Mb, anchoring 100 of the 120 contigs (>5 clones) that constitute ~95 % of the cumulative length of the map. Ordering of 80 contigs along the deletion bins of chromosome arm 5DS revealed small-scale breaks in syntenic blocks. Analysis of the gene space of 5DS suggested an increasing gradient of genes organized in islands towards the telomere, with the highest gene density of 5.17 genes/Mb in the 0.67-0.78 deletion bin, 1.4 to 1.6 times that of all other bins.

Conclusions

Here, we provide a chromosome-specific view into the organization and evolution of the D genome of bread wheat, in comparison to one of its ancestors, revealing recent genome rearrangements. The high-quality physical map constructed in this study paves the way for the assembly of a reference sequence, from which breeding efforts will greatly benefit.

Electronic supplementary material

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

Characterization of Ferredoxin-Dependent Glutamine-Oxoglutarate Amidotransferase (Fd-GOGAT) Genes and Their Relationship with Grain Protein Content QTL in Wheat

Background

In higher plants, inorganic nitrogen is assimilated via the glutamate synthase cycle or GS-GOGAT pathway. GOGAT enzyme occurs in two distinct forms that use NADH (NADH-GOGAT) or Fd (Fd-GOGAT) as electron carriers. The goal of the present study was to characterize wheat Fd-GOGAT genes and to assess the linkage with grain protein content (GPC), an important quantitative trait controlled by multiple genes.

Results

We report the complete genomic sequences of the three homoeologous A, B and D Fd-GOGAT genes from hexaploid wheat (Triticum aestivum) and their localization and characterization. The gene is comprised of 33 exons and 32 introns for all the three homoeologues genes. The three genes show the same exon/intron number and size, with the only exception of a series of indels in intronic regions. The partial sequence of the Fd-GOGAT gene located on A genome was determined in two durum wheat (Triticum turgidum ssp. durum) cvs Ciccio and Svevo, characterized by different grain protein content. Genomic differences allowed the gene mapping in the centromeric region of chromosome 2A. QTL analysis was conducted in the Svevo×Ciccio RIL mapping population, previously evaluated in 5 different environments. The study co-localized the Fd-GOGAT-A gene with the marker GWM-339, identifying a significant major QTL for GPC.

Conclusions

The wheat Fd-GOGAT genes are highly conserved; both among the three homoeologous hexaploid wheat genes and in comparison with other plants. In durum wheat, an association was shown between the Fd-GOGAT allele of cv Svevo with increasing GPC - potentially useful in breeding programs.     

A high-resolution physical map integrating an anchored chromosome with the BAC physical maps of wheat chromosome 6B

Background

A complete genome sequence is an essential tool for the genetic improvement of wheat. Because the wheat genome is large, highly repetitive and complex due to its allohexaploid nature, the International Wheat Genome Sequencing Consortium (IWGSC) chose a strategy that involves constructing bacterial artificial chromosome (BAC)-based physical maps of individual chromosomes and performing BAC-by-BAC sequencing. Here, we report the construction of a physical map of chromosome 6B with the goal of revealing the structural features of the third largest chromosome in wheat.

Results

We assembled 689 informative BAC contigs (hereafter reffered to as contigs) representing 91 % of the entire physical length of wheat chromosome 6B. The contigs were integrated into a radiation hybrid (RH) map of chromosome 6B, with one linkage group consisting of 448 loci with 653 markers. The order and direction of 480 contigs, corresponding to 87 % of the total length of 6B, were determined. We also characterized the contigs that contained a part of the nucleolus organizer region or centromere based on their positions on the RH map and the assembled BAC clone sequences. Analysis of the virtual gene order along 6B using the information collected for the integrated map revealed the presence of several chromosomal rearrangements, indicating evolutionary events that occurred on chromosome 6B.

Conclusions

We constructed a reliable physical map of chromosome 6B, enabling us to analyze its genomic structure and evolutionary progression. More importantly, the physical map should provide a high-quality and map-based reference sequence that will serve as a resource for wheat chromosome 6B.

Electronic supplementary material

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

Campylobacter group II phage CP21 is the prototype of a new subgroup revealing a distinct modular genome organization and host specificity

Background

The application of phages is a promising tool to reduce the number of Campylobacter along the food chain. Besides the efficacy against a broad range of strains, phages have to be safe in terms of their genomes. Thus far, no genes with pathogenic potential (e.g., genes encoding virulence factors) have been detected in Campylobacter phages. However, preliminary studies suggested that the genomes of group II phages may be diverse and prone to genomic rearrangements.

Results

We determined and analysed the genomic sequence (182,761 bp) of group II phage CP21 that is closely related to the already characterized group II phages CP220 and CPt10. The genomes of these phages are comprised of four modules separated by very similar repeat regions, some of which harbouring open reading frames (ORFs). Though, the arrangement of the modules and the location of some ORFs on the genomes are different in CP21 and in CP220/CPt10. In this work, a PCR system was established to study the modular genome organization of other group II phages demonstrating that they belong to different subgroups of the CP220-like virus genus, the prototypes of which are CP21 and CP220. The subgroups revealed different restriction patterns and, interestingly enough, also distinct host specificities, tail fiber proteins and tRNA genes. We additionally analysed the genome of group II phage vB_CcoM-IBB_35 (IBB_35) for which to date only five individual contigs could be determined. We show that the contigs represent modules linked by long repeat regions enclosing some yet not identified ORFs (e.g., for a head completion protein). The data suggest that IBB_35 is a member of the CP220 subgroup.

Conclusion

Campylobacter group II phages are diverse regarding their genome organization. Since all hitherto characterized group II phages contain numerous genes for transposases and homing endonucleases as well as similar repeat regions, it cannot be excluded that these phages are genetically unstable. To answer this question, further experiments and sequencing of more group II phages should be performed.

Electronic supplementary material

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

Alterations and Abnormal Mitosis of Wheat Chromosomes Induced by Wheat-Rye Monosomic Addition Lines

Background

Wheat-rye addition lines are an old topic. However, the alterations and abnormal mitotic behaviours of wheat chromosomes caused by wheat-rye monosomic addition lines are seldom reported.

Methodology/Principal Findings

Octoploid triticale was derived from common wheat T. aestivum L. ‘Mianyang11’×rye S. cereale L. ‘Kustro’ and some progeny were obtained by the controlled backcrossing of triticale with ‘Mianyang11’ followed by self-fertilization. Genomic in situ hybridization (GISH) using rye genomic DNA and fluorescence in situ hybridization (FISH) using repetitive sequences pAs1 and pSc119.2 as probes were used to analyze the mitotic chromosomes of these progeny. Strong pSc119.2 FISH signals could be observed at the telomeric regions of 3DS arms in ‘Mianyang11’. However, the pSc119.2 FISH signals were disappeared from the selfed progeny of 4R monosomic addition line and the changed 3D chromosomes could be transmitted to next generation stably. In one of the selfed progeny of 7R monosomic addition line, one 2D chromosome was broken and three 4A chromosomes were observed. In the selfed progeny of 6R monosomic addition line, structural variation and abnormal mitotic behaviour of 3D chromosome were detected. Additionally, 1A and 4B chromosomes were eliminated from some of the progeny of 6R monosomic addition line.

Conclusions/Significance

These results indicated that single rye chromosome added to wheat might cause alterations and abnormal mitotic behaviours of wheat chromosomes and it is possible that the stress caused by single alien chromosome might be one of the factors that induced karyotype alteration of wheat.     

Development of a D genome specific marker resource for diploid and hexaploid wheat

Background

Mapping and map-based cloning of genes that control agriculturally and economically important traits remain great challenges for plants with complex highly repetitive genomes such as those within the grass tribe, Triticeae. Mapping limitations in the Triticeae are primarily due to low frequencies of polymorphic gene markers and poor genetic recombination in certain genetic regions. Although the abundance of repetitive sequence may pose common problems in genome analysis and sequence assembly of large and complex genomes, they provide repeat junction markers with random and unbiased distribution throughout chromosomes. Hence, development of a high-throughput mapping technology that combine both gene-based and repeat junction-based markers is needed to generate maps that have better coverage of the entire genome.

Results

In this study, the available genomics resource of the diploid Aegilop tauschii, the D genome donor of bread wheat, were used to develop genome specific markers that can be applied for mapping in modern hexaploid wheat. A NimbleGen array containing both gene-based and repeat junction probe sequences derived from Ae. tauschii was developed and used to map the Chinese Spring nullisomic-tetrasomic lines and deletion bin lines of the D genome chromosomes. Based on these mapping data, we have now anchored 5,171 repeat junction probes and 10,892 gene probes, corresponding to 5,070 gene markers, to the delineated deletion bins of the D genome. The order of the gene-based markers within the deletion bins of the Chinese Spring can be inferred based on their positions on the Ae. tauschii genetic map. Analysis of the probe sequences against the Chinese Spring chromosome sequence assembly database facilitated mapping of the NimbleGen probes to the sequence contigs and allowed assignment or ordering of these sequence contigs within the deletion bins. The accumulated length of anchored sequence contigs is about 155 Mb, representing ~ 3.2 % of the D genome. A specific database was developed to allow user to search or BLAST against the probe sequence information and to directly download PCR primers for mapping specific genetic loci.

Conclusions

In bread wheat, aneuploid stocks have been extensively used to assign markers linked with genes/traits to chromosomes, chromosome arms, and their specific bins. Through this study, we added thousands of markers to the existing wheat chromosome bin map, representing a significant step forward in providing a resource to navigate the wheat genome. The database website (http://probes.pw.usda.gov/ATRJM/) provides easy access and efficient utilization of the data. The resources developed herein can aid map-based cloning of traits of interest and the sequencing of the D genome of hexaploid wheat.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1852-2) contains supplementary material, which is available to authorized users.Keyword: Wheat deletion bins, Molecular markers, Repeat junction markers, NimbleGen array, Recombination, Genetic map     

A novel model of double replications and random loss accounts for rearrangements in the Mitogenome of Samariscus latus (Teleostei: Pleuronectiformes)

Background

Although more than one thousand complete mitochondrial DNA (mtDNA) sequences have been determined in teleostean fishes, only a few gene rearrangements have been observed, and genome-scale rearrangements are even rarer. However, flatfishes (Pleuronectiformes) have been identified as having diverse types of mitochondrial gene rearrangements. It has been reported that tongue soles and the blue flounder mitogenomes exhibit different types of large-scale gene rearrangements.

Results

In the present study, the complete mitochondrial genome of another flatfish, Samariscus latus, was sequenced, and genome-scale rearrangements were observed. The genomic features of this flounder are different from those of any other studied vertebrates, including flatfish species too. The mitogenome of S. latus is characterized by the duplication and translocation of the control region (CR). The genes located between the two CRs are divided into two clusters in which their relative orders are maintained.

Conclusions

We propose a “Double Replications and Random Loss” model to explain the rearrangement events in S. latus mitogenome. This model consists of the following steps. First, the CR was duplicated and translocated. Subsequently, double replications of the mitogenome were successively initiated from the two CRs, leading to the duplication of the genes between the two CRs. Finally, one of each pair of duplicated genes was lost in a random event.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-352) contains supplementary material, which is available to authorized users.     

New Types of Wheat Chromosomal Structural Variations in Derivatives of Wheat-Rye Hybrids

Background

Chromosomal rearrangements induced by wheat-rye hybridization is a very well investigated research topic. However, the structural alterations of wheat chromosomes in wheat-rye hybrids are seldom reported.

Methodology/Principal Findings

Octoploid triticale lines were derived from common wheat Triticum. aestivum L. ‘Mianyang11’×rye Secale cereale L. ‘Kustro’. Some progeny were obtained by the controlled backcrossing of triticale with ‘Mianyang11’ and common wheat T. aestivum L. ‘Chuannong27’ followed by self-fertilization. Fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH) using Oligo-pSc119.2-1, Oligo-pTa535-1 and rye genomic DNA as probes were used to analyze the mitotic chromosomes of these progeny. Alterations of wheat chromosomes including 5A, 6A, 1B, 2B, 6B, 7B, 1D, 3D and 7D were observed. 5AL arm carrying intercalary Oligo-pSc119.2-1, Oligo-pTa535-1 or both Oligo-pSc119.2-1 and Oligo-pTa535-1 signals, 6AS, 1BS and 1DL arms containing terminal Oligo-pSc119.2-1 signal, 6BS and 3DS arms without terminal Oligo-pSc119.2-1 signal, 7BS without subtelomeric Oligo-pSc119.2-1 signal and 7DL with intercalary Oligo-pSc119.2-1 signal have been observed. However, these changed wheat chromosomes have not been detected in ‘Mianyang11’ and Chuannong 27. The altered 5A, 6A, 7B and 7D chromosomes in this study have not been reported and represent several new karyotype structures of common wheat chromosomes.

Conclusions/Significance

These rearranged wheat chromosomes in the present study afford some new genetic variations for wheat breeding program and are valuable materials for studying the biological function of tandem repetitive DNA sequences.