A molecular-cytogenetic method for locating genes to pericentromeric regions facilitates a genomewide comparison of synteny between the centromeric regions of wheat and rice

Comparison of Oryza sativa and Oryza brachyantha Genomes Reveals Selection-Driven Gene Escape from the Centromeric Regions

Centromeres are dynamic chromosomal regions, and the genetic and epigenetic environment of the centromere is often regarded as oppressive to protein-coding genes. Here, we used comparative genomic and phylogenomic approaches to study the evolution of centromeres and centromere-linked genes in the genus Oryza We report a 12.4-Mb high-quality BAC-based pericentromeric assembly for Oryza brachyantha, which diverged from cultivated rice (Oryza sativa) ∼15 million years ago. The synteny analyses reveal seven medium (>50 kb) pericentric inversions in O. sativa and 10 in O. brachyantha Of these inversions, three resulted in centromere movement (Chr1, Chr7, and Chr9). Additionally, we identified a potential centromere-repositioning event, in which the ancestral centromere on chromosome 12 in O. brachyantha jumped ∼400 kb away, possibly mediated by a duplicated transposition event (>28 kb). More strikingly, we observed an excess of syntenic gene loss at and near the centromeric regions (P < 2.2 × 10^-16). Most (33/47) of the missing genes moved to other genomic regions; therefore such excess could be explained by the selective loss of the copy in or near centromeric regions after gene duplication. The pattern of gene loss immediately adjacent to centromeric regions suggests centromere chromatin dynamics (e.g., spreading or microrepositioning) may drive such gene loss.

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.

A New Map Location of Gene Stb3 for Resistance to Septoria Tritici Blotch in Wheat

Septoria tritici blotch (STB), caused by Mycosphaerella graminicola (synonym: Zymoseptoria tritici; asexual stage: Septoria tritici), is an important disease of wheat worldwide. Management of the disease usually is by host resistance or fungicides. However, M. graminicola has developed insensitivity to most commonly applied fungicides so there is a continuing need for well-characterized sources of host resistance to accelerate the development of improved wheat cultivars. Gene Stb3 has been a useful source of major resistance, but its mapping location has not been well characterized. Based on linkage to a single marker, a previous study assigned Stb3 to a location on the short arm of chromosome 6D. However, the results from the present study show that this reported location is incorrect. Instead, linkage analysis revealed that Stb3 is located on the short arm of wheat chromosome 7A, completely linked to microsatellite (SSR) locus Xwmc83 and flanked by loci Xcfa2028 (12.4 cM distal) and Xbarc222 (2.1 cM proximal). Linkage between Stb3 and Xwmc83 was validated in BC₁F₃ progeny of other crosses, and analyses of the flanking markers with deletion stocks showed that the gene is located on 7AS between fraction lengths 0.73 and 0.83. This revised location of Stb3 is different from those for other STB resistance genes previously mapped in hexaploid wheat but is approximately 20 cM proximal to an STB resistance gene mapped on the short arm of chromosome 7A^m in Triticum monococcum. The markers described in this study are useful for accelerating the deployment of Stb3 in wheat breeding programs.

Next-generation survey sequencing and the molecular organization of wheat chromosome 6B

Common wheat (Triticum aestivum L.) is one of the most important cereals in the world. To improve wheat quality and productivity, the genomic sequence of wheat must be determined. The large genome size (∼17 Gb/1 C) and the hexaploid status of wheat have hampered the genome sequencing of wheat. However, flow sorting of individual chromosomes has allowed us to purify and separately shotgun-sequence a pair of telocentric chromosomes. Here, we describe a result from the survey sequencing of wheat chromosome 6B (914 Mb/1 C) using massively parallel 454 pyrosequencing. From the 4.94 and 5.51 Gb shotgun sequence data from the two chromosome arms of 6BS and 6BL, 235 and 273 Mb sequences were assembled to cover ∼55.6 and 54.9% of the total genomic regions, respectively. Repetitive sequences composed 77 and 86% of the assembled sequences on 6BS and 6BL, respectively. Within the assembled sequences, we predicted a total of 4798 non-repetitive gene loci with the evidence of expression from the wheat transcriptome data. The numbers and chromosomal distribution patterns of the genes for tRNAs and microRNAs in wheat 6B were investigated, and the results suggested a significant involvement of DNA transposon diffusion in the evolution of these non-protein-coding RNA genes. A comparative analysis of the genomic sequences of wheat 6B and monocot plants clearly indicated the evolutionary conservation of gene contents.

Sequence organization and evolutionary dynamics of Brachypodium-specific centromere retrotransposons

Brachypodium distachyon is a wild annual grass belonging to the Pooideae, more closely related to wheat, barley, and forage grasses than rice and maize. As an experimental model, the completed genome sequence of B. distachyon provides a unique opportunity to study centromere evolution during the speciation of grasses. Centromeric satellite sequences have been identified in B. distachyon, but little is known about centromeric retrotransposons in this species. In the present study, bacterial artificial chromosome (BAC)-fluorescence in situ hybridization was conducted in maize, rice, barley, wheat, and rye using B. distachyon (Bd) centromere-specific BAC clones. Eight Bd centromeric BAC clones gave no detectable fluorescence in situ hybridization (FISH) signals on the chromosomes of rice and maize, and three of them also did not yield any FISH signals in barley, wheat, and rye. In addition, four of five Triticeae centromeric BAC clones did not hybridize to the B. distachyon centromeres, implying certain unique features of Brachypodium centromeres. Analysis of Brachypodium centromeric BAC sequences identified a long terminal repeat (LTR)-centromere retrotransposon of B. distachyon (CRBd1). This element was found in high copy number accounting for 1.6 % of the B. distachyon genome, and is enriched in Brachypodium centromeric regions. CRBd1 accumulated in active centromeres, but was lost from inactive ones. The LTR of CRBd1 appears to be specific to B. distachyon centromeres. These results reveal different evolutionary events of this retrotransposon family across grass species.

Characterization of a new 4BS.7HL wheat–barley translocation line using GISH, FISH, and SSR markers and its effect on the β-glucan content of wheat

A spontaneous interspecific Robertsonian translocation was revealed by genomic in situ hybridization (GISH) in the progenies of a monosomic 7H addition line originating from a new wheat ‘Asakaze komugi’ × barley ‘Manas’ hybrid. Fluorescence in situ hybridization (FISH) with repetitive DNA sequences (Afa family, pSc119.2, and pTa71) allowed identification of all wheat chromosomes, including wheat chromosome arm 4BS involved in the translocation. FISH using barley telomere- and centromere-specific repetitive DNA probes (HvT01 and (AGGGAG)n) confirmed that one of the arms of barley chromosome 7H was involved in the translocation. Simple sequence repeat (SSR) markers specific to the long (L) and short (S) arms of barley chromosome 7H identified the translocated chromosome segment as 7HL. Further analysis of the translocation chromosome clarified the physical position of genetically mapped SSRs within 7H, with a special focus on its centromeric region. The presence of the HvCslF6 gene, responsible for (1,3;1,4)-β-d-glucan production, was revealed in the centromeric region of 7HL. An increased (1,3;1,4)-β-d-glucan level was also detected in the translocation line, demonstrating that the HvCslF6 gene is of potential relevance for the manipulation of wheat (1,3;1,4)-β-d-glucan levels.

Asymmetric distribution of gene expression in the centromeric region of rice chromosome 5

There is controversy as to whether gene expression is silenced in the functional centromere. The complete genomic sequences of the centromeric regions in higher eukaryotes have not been fully elucidated, because the presence of highly repetitive sequences complicates many aspects of genomic sequencing. We performed resequencing, assembly, and sequence finishing of two P1-derived artificial chromosome clones in the centromeric region of rice (Oryza sativa L.) chromosome 5 (Cen5). The pericentromeric region, where meiotic recombination is silenced, is located at the center of chromosome 5 and is 2.14 Mb long; a total of six restriction-fragment-length polymorphism markers (R448, C1388, S20487S, E3103S, C53260S, and R2059) genetically mapped at 54.6 cM were located in this region. In the pericentromeric region, 28 genes were annotated on the short arm and 45 genes on the long arm. To quantify all transcripts in this region, we performed massive parallel sequencing of mRNA. Transcriptional density (total length of transcribed region/length of the genomic region) and expression level (number of uniquely mapped reads/length of transcribed region) were calculated on the basis of the mapped reads on the rice genome. Transcriptional density and expression level were significantly lower in Cen5 than in the average of the other chromosomal regions. Moreover, transcriptional density in Cen5 was significantly lower on the short arm than on the long arm; the distribution of transcriptional density was asymmetric. The genomic sequence of Cen5 has been integrated into the most updated reference rice genome sequence constructed by the International Rice Genome Sequencing Project.

Molecular mapping of the suppressor gene Igc1 to the gametocidal gene Gc3-C1 in common wheat

Several species of the genus Aegilops, wild relatives of wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) carry gametocidal (Gc) genes. Gc genes kill the gametes without themselves by causing chromosomal breakage during post-meiotic cell divisions, and therefore are strong segregation distorters. The Gc gene Gc3-C1 derived from chromosome 3C of Ae. triuncialis (2n = 4x = 28, CCUU) induces chromosomal breakage in wheat cultivar 'Chinese Spring' (CS) but not in cultivar 'Norin 26' (N26). This cultivar-specific inhibition of Gc function is caused by a suppressor gene Igc1 located on chromosome 3B of N26. Igc1 is presumed to be a modified Gc gene without breakage function because of its homoeology to Gc3-C1. Here we report the results of linkage and physical mapping of Igc1 to help elucidate the molecular mechanisms underlying Gc action. Segregation analysis of the phenotypic data in BC(1)F(1) mapping population of the cross between (CSxN26)F(1) and CS + 3C" showed a 1:1 segregation ratio indicating that Igc1 is a dominant gene. In the linkage analysis, three molecular marker loci Xgwm285, Xgwm376, and Xcfp1886 cosegregated with the Igc1 locus. Bin mapping assigned the loci Xgwm285 and Xcfp1886 to bin C-3BS1-0.33 and Xgwm376 to bin C-3BL2-0.22. Physical mapping using Gc-induced chromosomal deletion lines of chromosome 3B of N26 revealed that the Igc1 locus resides in 52.0% or 2.1% of bins C-3BS1-0.33 and C-3BL2-0.22, respectively. Pericentromeric localization of Igc1 in chromosome 3B of N26 may have a positive effect to keep the two-component system of the Gc action. Map-based cloning approach to isolate the Igc1 may be difficult because recombination is depleted in the pericentromeric region. As is shown in this study, the combination of genetic and physical mapping offers high efficiency to identify the regions where genes are located especially in regions with low levels of recombination.