Targeted analysis of orthologous phytochrome A regions of the sorghum,maize, and rice genomes using comparative gene-island sequencing

A "gene-island" sequencing strategy has been developed that expedites the targeted acquisition of orthologous gene sequences from related species for comparative genome analysis. A 152-kb bacterial artificial chromosome (BAC) clone from sorghum (Sorghum bicolor) encoding phytochrome A (PHYA) was fully sequenced, revealing 16 open reading frames with a gene density similar to many regions of the rice (Oryza sativa) genome. The sequences of genes in the orthologous region of the maize (Zea mays) and rice genomes were obtained using the gene-island sequencing method. BAC clones containing the orthologous maize and rice PHYA genes were identified, sheared, subcloned, and probed with the sorghum PHYA-containing BAC DNA. Sequence analysis revealed that approximately 75% of the cross-hybridizing subclones contained sequences orthologous to those within the sorghum PHYA BAC and less than 25% contained repetitive and/or BAC vector DNA sequences. The complete sequence of four genes, including up to 1 kb of their promoter regions, was identified in the maize PHYA BAC. Nine orthologous gene sequences were identified in the rice PHYA BAC. Sequence comparison of the orthologous sorghum and maize genes aided in the identification of exons and conserved regulatory sequences flanking each open reading frame. Within genomic regions where micro-colinearity of genes is absolutely conserved, gene-island sequencing is a particularly useful tool for comparative analysis of genomes between related species.     

Transfer of grain softness from 5U-5A wheat-Aegilops triuncialis substitution line to bread wheat through induced homeologous pairing

Journal of Plant Biochemistry and Biotechnology - Bread wheat sustains genes for grain softness on “Ha” locus of short arm of 5D chromosome. Pina-D1 and Pinb-D1 alleles of...     

Rapid genome divergence at orthologous low molecular weight glutenin loci of the A and Am genomes of wheat

To study genome evolution in wheat, we have sequenced and compared two large physical contigs of 285 and 142 kb covering orthologous low molecular weight (LMW) glutenin loci on chromosome 1AS of a diploid wheat species (Triticum monococcum subsp monococcum) and a tetraploid wheat species (Triticum turgidum subsp durum). Sequence conservation between the two species was restricted to small regions containing the orthologous LMW glutenin genes, whereas >90% of the compared sequences were not conserved. Dramatic sequence rearrangements occurred in the regions rich in repetitive elements. Dating of long terminal repeat retrotransposon insertions revealed different insertion events occurring during the last 5.5 million years in both species. These insertions are partially responsible for the lack of homology between the intergenic regions. In addition, the gene space was conserved only partially, because different predicted genes were identified on both contigs. Duplications and deletions of large fragments that might be attributable to illegitimate recombination also have contributed to the differentiation of this region in both species. The striking differences in the intergenic landscape between the A and A(m) genomes that diverged 1 to 3 million years ago provide evidence for a dynamic and rapid genome evolution in wheat species.     

A bacterial artificial chromosome contig spanning the major domestication locus Q in wheat and identification of a candidate gene

The Q locus played a major role in the domestication of wheat because it confers the free-threshing character and influences many other agronomically important traits. We constructed a physical contig spanning the Q locus using a Triticum monococcum BAC library. Three chromosome walking steps were performed by complete sequencing of BACs and identification of low-copy markers through similarity searches of database sequences. The BAC contig spans a physical distance of approximately 300 kb corresponding to a genetic distance of 0.9 cM. The physical map of T. monococcum had perfect colinearity with the genetic map of wheat chromosome arm 5AL. Recombination data in conjunction with analysis of fast neutron deletions confirmed that the contig spanned the Q locus. The Q gene was narrowed to a 100-kb segment, which contains an APETALA2 (AP2)-like gene that cosegregates with Q. AP2 is known to play a major role in controlling floral homeotic gene expression and thus is an excellent candidate for Q.     

In silico comparative analysis reveals a mosaic conservation of genes within a novel colinear region in wheat chromosome 1AS and rice chromosome 5S

Comparative RFLP mapping has revealed extensive conservation of marker order in different grass genomes. However, microcolinearity studies at the sequence level have shown rapid genome evolution and many exceptions to colinearity. Most of these studies have focused on a limited size of genomic fragment and the extent of microcolinearity over large distances or across entire genomes remains poorly characterized in grasses. Here, we have investigated the microcolinearity between the rice genome and a total of 1,500 kb from physical BAC contigs on wheat chromosome 1AS. Using ESTs mapped in wheat chromosome bins as an additional source of physical data, we have identified 27 conserved orthologous sequences between wheat chromosome 1AS and a region of 1,210 kb located on rice chromosome 5S. Our results extend the orthology described earlier between wheat chromosome group 1S and rice chromosome 5S. Microcolinearity was found to be frequently disrupted by rearrangements which must have occurred after the divergence of wheat and rice. At the Lr10 orthologous loci, microrearrangements were due to the insertion of mobile elements, but also originated from gene movement, amplification, deletion and inversion. These mechanisms of genome evolution are at the origin of the mosaic conservation observed between the orthologous regions. Finally, in silico mapping of wheat genes identified an intragenomic colinearity between fragments from rice chromosome 1L and 5S, suggesting an ancestral segmental duplication in rice.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Comparative mapping of wheat chromosome 1AS which contains the tiller inhibition gene (<Emphasis Type="Italic">tin</Emphasis>) with rice chromosome 5S

The capacity to tiller is a key factor that determines plant architecture. Using molecular markers, a single major gene reducing tiller number, formally named the tiller inhibition gene (tin), was mapped to the short arm of chromosome 1A in wheat. We identified a tightly linked microsatellite marker (Xgwm136) that may be useful in future marker-assisted selection. The tin gene was mapped to the distal deletion bin of chromosome 1AS (FLM value 0.86) and wheat ESTs which were previously mapped to the same deletion bin were used to identify 18 closely related sequences in the syntenic region of rice chromosome 5. For a subset of wheat ESTs that detected flanking markers for tin, we identified closely related sequences within the most distal 300 kb of rice chromosome 5S. The synteny between the distal chromosome ends of wheat 1AS and rice 5S appeared to be disrupted at the hairy glume locus and seed storage protein loci. We compared map position of tin with other reduced tillering mutants characterised in other cereals to identify possible orthologous genes.     

Molecular and biochemical characterization of puroindoline a and b alleles in Chinese landraces and historical cultivars

Kernel hardness that is conditioned by puroindoline genes has a profound effect on milling, baking and end-use quality of bread wheat. In this study, 219 landraces and 166 historical cultivars from China and 12 introduced wheats were investigated for their kernel hardness and puroindoline alleles, using molecular and biochemical markers. The results indicated that frequencies of soft, mixed and hard genotypes were 42.7, 24.3, and 33.0%, respectively, in Chinese landraces and 45.2, 13.9, and 40.9% in historical cultivars. The frequencies of PINA null, Pinb-D1b and Pinb-D1p genotypes were 43.8, 12.3, and 39.7%, respectively, in hard wheat of landraces, while 48.5, 36.8, and 14.7%, respectively, in historical hard wheats. A new Pinb-D1 allele, designated Pinb-D1t, was identified in two landraces, Guangtouxianmai and Hongmai from the Guizhou province, with the characterization of a glycine to arginine substitution at position 47 in the coding region of Pinb gene. Surprisingly, a new Pina-D1 allele, designated Pina-D1m, was detected in the landrace Hongheshang, from the Jiangsu province, with the characterization of a proline to serine substitution at position 35 in the coding region of Pina gene; it was the first novel mutation found in bread wheat, resulting in a hard endosperm with PINA expression. Among the PINA null genotypes, an allele designed as Pina-D1l, was detected in five landraces with a cytosine deletion at position 265 in Pina locus; while another novel Pina-D1 allele, designed as Pina-D1n, was identified in six landraces, with the characterization of an amino acid change from tryptophan-43 to a ‘stop’ codon in the coding region of Pina gene. The study of puroindoline polymorphism in Chinese wheat germplasm could provide useful information for the further understanding of the molecular basis of kernel hardness in bread wheat.     

The wheat D-genome HMW-glutenin locus: BAC sequencing,gene distribution,and retrotransposon clusters

A bacterial-artificial-chromosome (BAC) clone from the genome of Triticum tauschii, the D-genome ancestor of hexaploid bread wheat, was sequenced and the presence of the two paralogous x- and y-type high-molecular-weight (HMW) glutenin genes of the Glu-D1 locus was confirmed. These two genes occur in the same orientation, are 51,893 bp apart, and the separating DNA includes a 31,000-bp cluster of retrotransposons. A second retrotransposon cluster of 32,000 bp follows the x-type HMW-glutenin gene region. Each HMW-glutenin gene is found within a region of mainly unique DNA sequence which includes multiple additional genes including an active endosperm globulin gene not previously reported in the Triticeae family, a leucine-rich-repeat (LRR) type gene truncated at the 5′ end of the BAC, a kinase gene of unknown activity, remnants of a paralogous second globulin gene, and genes similar to two hypothetical rice genes. The newly identified globulin genes are assigned to a locus designated Glo-2. Comparison to available orthologous regions of the wheat A and B genomes show rapid sequence divergences flanking the HMW-glutenin genes, and the absence of two hypothetical and unknown genes found 5′ to the B-genome x-type ortholog. The region surrounding the Glu-D1 locus is similar to other reported Triticeae BAC sequences; i.e. small gene islands separated by retrotransposon clusters. Electronic Publication     

High-resolution mapping of the leaf rust disease resistance gene <Emphasis Type="Italic">Lr1</Emphasis> in wheat and characterization of BAC clones from the <Emphasis Type="Italic">Lr1</Emphasis> locus

Leaf rust is the most common disease in wheat production. There are more than 45 specific resistance genes described and used in wheat breeding to control epidemics of leaf rust, but none of them has been cloned. The leaf rust disease resistance gene 1 ( Lr1) is a good model gene for isolation by map-based cloning because it is a single, dominant gene which is located in the distal region of chromosome 5DL of wheat. As the first step towards the isolation of this gene we constructed a high-resolution genetic map in the region of the Lr1 locus by saturation mapping of two large segregating F(2) populations (Thatcher Lr1 x Thatcher, Thatcher Lr1 x Frisal). The resistance gene Lr1 was delimited in a 0.16-cM region between the RFLP markers ABC718 and PSR567 (0.12 cM from ABC718 and 0.04 cM from PSR567). A genomic BAC library of Aegilops tauschii (D genome) was screened using the RFLP markers ABC718 and PSR567. Five positive BAC clones were identified by ABC718 and four clones by PSR567. Two NBS-LRR type of resistance gene analogs, which encode proteins highly homologous to the bacterial blight disease resistance protein Xa1 of rice, were identified on BAC clones isolated with PSR567. Polymorphic BAC end probes were isolated from both ends of a 105-kb large BAC clone identified by ABC718. The end probes were mapped at the same locus as ABC718, and no recombination event was found within 105 kb around ABC718 in our analysis of more than 4,000 gametes.     

Molecular characterization of a diagnostic DNA marker for domesticated tetraploid wheat provides evidence for gene flow from wild tetraploid wheat to hexaploid wheat

All forms of domesticated tetraploid wheat (Triticum turgidum, genomes AABB) are nearly monomorphic for restriction fragment length polymorphism (RFLP) haplotype a at the Xpsr920 locus on chromosome 4A (Xpsr920-A1a), and wild tetraploid wheat is monomorphic for haplotype b. The Xpsr920-A1a/b dimorphism provides a molecular marker for domesticated and wild tetraploid wheat, respectively. Hexaploid wheat (Triticum aestivum, genomes AABBDD) is polymorphic for the 2 haplotypes. Bacterial artificial chromosome (BAC) clones hybridizing with PSR920 were isolated from Triticum urartu (genomes AA), Triticum monococcum (genomes AmAm), and T. turgidum ssp. durum (genomes AABB) and sequenced. PSR920 is a fragment of a putative ATP binding cassette (ABC) transporter gene (designated ABCT-1). The wheat ABCT-1 gene is more similar to the T. urartu gene than to the T. monococcum gene and diverged from the T. urartu gene about 0.7 MYA. The comparison of the sequence of the wheat A genome BAC clone with that of the T. urartu BAC clone provides the first insight into the microsynteny of the wheat A genome with that of T. urartu. Within 103 kb of orthologous intergenic space, 37 kb of new DNA has been inserted and 36 kb deleted leaving 49.7% of the region syntenic between the clones. The nucleotide substitution rate in the syntenic intergenic space has been 1.6 x 10(-8) nt(-1) year(-1), which is, respectively, 4 and 3 times as great as nucleotide substitution rates in the introns and the third codon positions of the juxtaposed gene. The RFLP is caused by a miniature inverted transposable element (MITE) insertion into intron 18 of the ABCT-A1 gene. Polymerase chain reaction primers were developed for the amplification of the MITE insertion site and its sequencing. The T. aestivum ABCT-A1a haplotype is identical to the haplotype of domesticated tetraploid wheat, and the ABCT-A1b haplotype is identical to that of wild tetraploid wheat. This finding shows for the first time that wild tetraploid wheat participated in the evolution of hexaploid wheat. A cline of the 2 haplotype frequencies exists across Euro-Asia in T. aestivum. It is suggested that T. aestivum in eastern Asia conserved the gene pool of the original T. aestivum more than wheat elsewhere.     

Puroindoline A-gene expression is involved in association of puroindolines to starch

Puroindolines largely influence cereal grain hardness. In order to understand how they exert this influence, we carried out a molecular analysis of the pina and pinb genes of many Italian wheat cultivars. On the basis of their pin genotypes they could be divided into three groups: Pina-D1a/Pinb-D1a; Pina-D1a/Pinb-D1b; and Pina-D1b/Pinb-D1a. Five cultivars from each group were chosen to be studied to examine the quantity of puroindolines associated with starch (friabilin) and the amount not associated with starch. In addition, the level of pina expression was measured using RT-PCR. Soft cultivars (Pina-D1a/Pinb-D1a) exhibited the highest level of expression of pina; among the hard cultivars, those with the Pina-D1a/Pinb-D1b genotype showed a lower level of expression, while those with the Pina-D1b/Pinb-D1a genotype did not express pina. Total puroindoline and friabilin content was then measured by flow cytometry. Soft Pina-D1a/Pinb-D1a cultivars displayed high puroindoline content that was primarily starch associated. Hard Pina-D1b/Pinb-D1a cultivars had very low puroindoline content with no puroindoline bound to starch. Hard Pina-D1a/Pinb-D1b cultivars were highly heterogeneous with respect to both the content of puroindolines and the level of association with starch. The accurate quantification of puroindolines in starch-bound and not starch-bound forms in association with molecular analysis, indicates that pina expression and presence controls the abundance of total puroindoline and its association with starch.Communicated by H.F. Linskens     

Intraspecific sequence comparisons reveal similar rates of non-collinear gene insertion in the B and D genomes of bread wheat

ABSTRACT: BACKGROUND: Polyploidization is considered one of the main mechanisms of plant genome evolution. The presence of multiple copies of the same gene reduces selection pressure and permits sub-functionalization and neo-functionalization leading to plant diversification, adaptation and speciation. In bread wheat, polyploidization and the prevalence of transposable elements resulted in massive gene duplication and movement. As a result, the number of genes which are non-collinear to genomes of related species seems markedly increased in wheat. RESULTS: We used new-generation sequencing (NGS) to generate sequence of a Mb-sized region from wheat chromosome arm 3DS. Sequence assembly of 24 BAC clones resulted in two scaffolds of 1,264,820 and 333,768 bases. The sequence was annotated and compared to the homoeologous region on wheat chromosome 3B and orthologous loci of Brachypodium distachyon and rice. Among 39 coding sequences in the 3DS scaffolds, 32 have a homoeolog on chromosome 3B. In contrast, only fifteen and fourteen orthologs were identified in the corresponding regions in rice and Brachypodium, respectively. Interestingly, five pseudogenes were identified among the non-collinear coding sequences at the 3B locus, while none was found at the 3DS locus. CONCLUSION: Direct comparison of two Mb-sized regions of the B and D genomes of bread wheat revealed similar rates of non-collinear gene insertion in both genomes with a majority of gene duplications occurring before their divergence. Relatively low proportion of pseudogenes was identified among non-collinear coding sequences. Our data suggest that the pseudogenes did not originate from insertion of non-functional copies, but were formed later during the evolution of hexaploid wheat. Some evidence was found for gene erosion along the B genome locus.     

Detailed comparison between the wheat chromosome group 7 short arms and the rice chromosome arms 6S and 8L with special reference to genes involved in starch biosynthesis

Rice bacterial artificial chromosome (BAC) clones have been identified that contain sequences orthologous to each EST localized to wheat chromosome 7AS deletion stocks by Southern blot hybridization. This information has been used to relate the DNA sequence included in each wheat deletion stock to a complement of rice BACs. A virtual contig was used that covered 90 cM (21 Mb) of DNA sequence (with a gap for the 6S/8L junction). Comparison of the positions of orthologous genes on the rice virtual contig and on wheat chromosome 7AS showed that there was an unexpectedly low level of synteny (31.4%) and a high level of chromosome rearrangements (68.6%). The non-syntenous loci were of two classes: wheat and rice genes found at different locations in the genome (32.6%), and ESTs in wheat not present in rice (36.0%). Four starch synthetic genes, GBSSI, SSI, SSIIa and DBEI, were located at similar positions on wheat chromosome 7AS and the virtual rice contig covering wheat chromosome 7AS. A preliminary comparison between the short arms of chromosome 7A and 7D in wheat showed that both chromosomes had a similar level of sequence synteny with rice. Therefore, there appears to be considerable variation in gene order between wheat chromosome 7S and rice chromosome 6S and 8L.     

Macro- and microcolinearity between the genomic region of wheat chromosome 5B containing the Tsn1 gene and the rice genome

The Tsn1 gene in wheat confers sensitivity to a proteinaceous host-selective toxin (Ptr ToxA) produced by the tan spot fungus (Pyrenophora tritici-repentis) and lies within a gene-rich region of chromosome 5B. To use the rice genome sequence information for the map-based cloning of Tsn1, colinearity between the wheat genomic region containing Tsn1 and the rice genome was determined at the macro- and microlevels. Macrocolinearity was determined by testing 28 expressed sequence markers (ESMs) spanning a 25.5-cM segment and encompassing Tsn1 for similarity to rice sequences. Twelve ESMs had no similarity to rice sequences, and 16 had similarity to sequences on seven different rice chromosomes. Segments of colinearity with rice chromosomes 3 and 9 were identified, but frequent rearrangements and disruptions occurred. Microcolinearity was determined by testing the sequences of 26 putative genes identified from BAC contigs of 205 and 548 kb in length and flanking Tsn1 for similarity to rice genomic sequences. Fourteen of the predicted genes detected orthologous sequences on six different rice chromosomes, whereas the remaining 12 had no similarity with rice sequences. Four genes were colinear on rice chromosome 9, but multiple disruptions, rearrangements, and duplications were observed in wheat relative to rice. The data reported provide a detailed analysis of a region of wheat chromosome 5B that is highly rearranged relative to rice.     

High-resolution mapping of the S and Z loci of Phalaris coerulescens.

Self incompatibility (SI) in Phalaris coerulescens is gametophytically determined by two unlinked multi allelic loci (S and Z). Neither the S nor Z genes have yet been cloned. As part of a map-based cloning strategy, high-resolution maps of the S and Z regions were generated from distorted segregating populations using RFLP probes from wheat, barley, oat, and Phalaris. The S locus was delimited to 0.26 cM with two boundary markers (Xwg811 and Xpsr168) and cosegregated with Xbm2 and Xbcd762. Xbcd266 was the closest marker linked to Z (0.9 cM). A high level of colinearity in the S and Z regions was found in both self-incompatible and -compatible species. The S locus was localized to the subcentromere region of chromosome 1 and the Z locus to the long arm end of chromosome 2. Several rice BAC clones orthologous to the S and Z locus regions were identified. This opens the possibility of using the rice genome sequence data to generate more closely linked markers and identify SI candidate genes. These results add further support to the conservation of gene order in the S and Z regions of the grass genomes.     

A sodium transporter (HKT7) is a candidate for Nax1, a gene for salt tolerance in durum wheat

Durum wheat (Triticum turgidum subsp. durum) is more salt sensitive than bread wheat (Triticum aestivum). A novel source of Na(+) exclusion conferring salt tolerance to durum wheat is present in the durum wheat Line 149 derived from Triticum monococcum C68-101, and a quantitative trait locus contributing to low Na(+) concentration in leaf blades, Nax1, mapped to chromosome 2AL. In this study, we used the rice (Oryza sativa) genome sequence and data from the wheat expressed sequence tag deletion bin mapping project to identify markers and construct a high-resolution map of the Nax1 region. Genes on wheat chromosome 2AL and rice chromosome 4L had good overall colinearity, but there was an inversion of a chromosomal segment that includes the Nax1 locus. Two putative sodium transporter genes (TmHKT7) related to OsHKT7 were mapped to chromosome 2AL. One TmHKT7 member (TmHKT7-A1) was polymorphic between the salt-tolerant and -sensitive lines, and cosegregated with Nax1 in the high-resolution mapping family. The other TmHKT7 member (TmHKT7-A2) was located within the same bacterial artificial chromosome contig of approximately 145 kb as TmHKT7-A1. TmHKT7-A1 and -A2 showed 83% amino acid identity. TmHKT7-A2, but not TmHKT7-A1, was expressed in roots and leaf sheaths of the salt-tolerant durum wheat Line 149. The expression pattern of TmHKT7-A2 was consistent with the physiological role of Nax1 in reducing Na(+) concentration in leaf blades by retaining Na(+) in the sheaths. TmHKT7-A2 could control Na(+) unloading from xylem in roots and sheaths.     

Structural characterization of Brachypodium genome and its syntenic relationship with rice and wheat

Brachypodium distachyon (Brachypodium) has been recently recognized as an emerging model system for both comparative and functional genomics in grass species. In this study, 55,221 repeat masked Brachypodium BAC end sequences (BES) were used for comparative analysis against the 12 rice pseudomolecules. The analysis revealed that ~26.4% of BES have significant matches with the rice genome and 82.4% of the matches were homologous to known genes. Further analysis of paired-end BES and ~1.0 Mb sequences from nine selected BACs proved to be useful in revealing conserved regions and regions that have undergone considerable genomic changes. Differential gene amplification, insertions/deletions and inversions appeared to be the common evolutionary events that caused variations of microcolinearity at different orthologous genomic regions. It was found that ~17% of genes in the two genomes are not colinear in the orthologous regions. Analysis of BAC sequences also revealed higher gene density (~9 kb/gene) and lower repeat DNA content (~13.1%) in Brachypodium when compared to the orthologous rice regions, consistent with the smaller size of the Brachypodium genome. The 119 annotated Brachypodium genes were BLASTN compared against the wheat EST database and deletion bin mapped wheat ESTs. About 77% of the genes retrieved significant matches in the EST database, while 9.2% matched to the bin mapped ESTs. In some cases, genes in single Brachypodium BACs matched to multiple ESTs that were mapped to the same deletion bins, suggesting that the Brachypodium genome will be useful for ordering wheat ESTs within the deletion bins and developing specific markers at targeted regions in the wheat genome.