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     

Dynamics of the evolution of orthologous and paralogous portions of a complex locus region in two genomes of allopolyploid wheat

Two overlapping bacterial artificial chromosome (BAC) clones from the B genome of the tetraploid wheat Triticum turgidum were identified, each of which contains one of the two high-molecular-weight (HMW) glutenin genes, comprising the complex Glu-B1 locus. The complete sequence (285 506 bp of DNA) of this chromosomal region was determined. The two paralogous x-type ( Glu-1-1 ) and y-type ( Glu-1-2 ) HMW-glutenin genes of the complex Glu-B1 locus were found to be separated by ca. 168 000 bp instead of the 51 000 bp separation previously reported for the orthologous Glu-D1 locus of Aegilops tauschii, the D-genome donor of hexaploid wheat. This difference in intergene spacing is due almost entirely to be the insertion of clusters of nested retrotransposons. Otherwise, the orientation and order of the HMW glutenins and adjacent genes were identical in the two genomes. A comparison of these orthologous regions indicates modes and patterns of sequence divergence, with implications for the overall Triticeae genome structure and evolution. A duplicate globulin gene, found 5' of each HMW-glutenin gene, assists to tentatively define the original duplication event leading to the paralogous x- and y-type HMW-glutenin genes. The intergenic regions of the two loci are composed of different patterns and classes of retrotransposons, indicating that insertion times of these retroelements were after the divergence of the two wheat genomes. In addition, a putative receptor kinase gene near the y-type HMW-glutenin gene at the Glu-B1 locus is likely active as it matches recently reported ESTs from germinating barley endosperm. The presence of four genes represented only in the Triticeae endosperm ESTs suggests an endosperm-specific chromosome domain.     

Rapid evolutionary dynamics in a 2.8‐Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi‐gene families, the DNA sequence of a 2.8‐Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance‐like gene families, was analyzed in the diploid grass Aegilops tauschii, the D‐genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non‐syntenic genes and only 13 ancestral genes are shared among these grass species. These non‐syntenic genes, including prolamin and resistance‐like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non‐syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin‐related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance‐like genes and subsequent differential expansion of each R gene family. The high frequency of non‐syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8‐Mb region with nine‐fold higher than the genome‐wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.     

Analysis of the barley and rice genomes by comparative RFLP linkage mapping

Comparative genetic mapping of rice and barley, both major crop species with extensive genetic resources, offers the possibility of uniting two well-established and characterized genetic systems. In the present study, we screened 229 molecular markers and utilized 110 polymorphic orthologous loci to construct comparative maps of the rice and barley genomes. While extensive chromosomal rearrangements, including inversions and intrachromosomal translocations, differentiate the rice and barley genomes, several syntenous chromosomes are evident. Indeed, several chromosomes and chromosome arms appear to share nearly identical gene content and gene order. Seventeen regions of conserved organization were detected, spanning 287 cM (24%) and 321 cM (31%) of the rice and barley genomes, respectively. The results also indicate that most (72%) of the single-copy sequences in barley are also single copy in rice, suggesting that the large barley genome arose by unequal crossing over and amplification of repetitive DNA sequences and not by the duplication of single-copy sequences. Combining these results with those previously reported for comparative analyses of rice and wheat identified nine putatively syntenous chromosomes among barley, wheat and rice. The high degree of gene-order conservation as detected by comparative mapping has astonishing implications for interpreting genetic information among species and for elucidating chromosome evolution and speciation.     

Localization of anchor loci representing five hundred annotated rice genes to wheat chromosomes using PLUG markers

PCR-based Landmark Unique Gene (PLUG) markers are EST-PCR markers developed based on the orthologous gene conservation between rice and wheat, and on the intron polymorphisms among the three orthologous genes derived from the A, B and D genomes of wheat. We designed a total of 960 primer sets from wheat ESTs that showed high similarity with 951 single-copy rice genes. When genomic DNA of Chinese Spring wheat was used as a template, 872 primer sets amplified one to five distinct products. Out of these 872 PLUG markers, 531 were assigned to one or more chromosomes by nullisomic-tetrasomic analysis. For each wheat chromosome, the number of loci detected ranged from 32 for chromosome 6A to 73 for chromosome 7D, with an average of 48 loci per chromosome. Several novel synteny perturbations were identified using deletion bin-mapping of markers. Furthermore, we demonstrated that PLUG markers can be used as probes to simultaneously identify BAC clones that contain homoeologous regions from all three genomes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Conserved globulin gene across eight grass genomes identify fundamental units of the loci encoding seed storage proteins

The wheat high molecular weight (HMW) glutenins are important seed storage proteins that determine bread-making quality in hexaploid wheat (Triticum aestivum). In this study, detailed comparative sequence analyses of large orthologous HMW glutenin genomic regions from eight grass species, representing a wide evolutionary history of grass genomes, reveal a number of lineage-specific sequence changes. These lineage-specific changes, which resulted in duplications, insertions, and deletions of genes, are the major forces disrupting gene colinearity among grass genomes. Our results indicate that the presence of the HMW glutenin gene in Triticeae genomes was caused by lineage-specific duplication of a globulin gene. This tandem duplication event is shared by Brachypodium and Triticeae genomes, but is absent in rice, maize, and sorghum, suggesting the duplication occurred after Brachypodium and Triticeae genomes diverged from the other grasses ~35 Ma ago. Aside from their physical location in tandem, the sequence similarity, expression pattern, and conserved cis-acting elements responsible for endosperm-specific expression further support the paralogous relationship between the HMW glutenin and globulin genes. While the duplicated copy in Brachypodium has apparently become nonfunctional, the duplicated copy in wheat has evolved to become the HMW glutenin gene by gaining a central prolamin repetitive domain.     

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.     

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.     

Analysis of a contiguous 211 kb sequence in diploid wheat (Triticum monococcum L.) reveals multiple mechanisms of genome evolution

In plant species with large genomes such as wheat or barley, genome organization at the level of DNA sequence is largely unknown. The largest sequences that are publicly accessible so far from Triticeae genomes are two 60 kb and 66 kb intervals from barley. Here, we report on the analysis of a 211 kb contiguous DNA sequence from diploid wheat (Triticum monococcum L.). Five putative genes were identified, two of which show similarity to disease resistance genes. Three of the five genes are clustered in a 31 kb gene-enriched island while the two others are separated from the cluster and from each other by large stretches of repetitive DNA. About 70% of the contig is comprised of several classes of transposable elements. Ten different types of retrotransposons were identified, most of them forming a pattern of nested insertions similar to those found in maize and barley. Evidence was found for major deletion, insertion and duplication events within the analysed region, suggesting multiple mechanisms of genome evolution in addition to retrotransposon amplification. Seven types of foldback transposons, an element class previously not described for wheat genomes, were characterized. One such element was found to be closely associated with genes in several Triticeae species and may therefore be of use for the identification of gene-rich regions in these species.     

Comparative sequence analysis of a single-gene conserved segment in mouse and human

Comparative mapping and sequencing of the mouse and human genomes have defined large, conserved chromosomal segments in which gene content and order are highly conserved. These regions span megabase-sized intervals and together comprise the vast majority of both genomes. However, the evolutionary relationships among the small remaining portions of these genomes are not as well characterized. Here we describe the sequencing and annotation of a 341-kb region of mouse Chr 2 containing nine genes, including biliverdin reductase A (Blvra), and its comparison with the orthologous regions of the human and rat genomes. These analyses reveal that the known conserved synteny between mouse Chromosome (Chr) 2 and human Chr 7 reflects an interval containing one gene (Blvra/BLVRA) that is, at most, just 34 kb in the mouse genome. In the mouse, this segment is flanked proximally by genes orthologous to human chromosome 15q21 and distally by genes orthologous to human Chr 2q11. The observed differences between the human and mouse genomes likely resulted from one or more rearrangements in the rodent lineage. In addition to the resulting changes in gene order and location, these rearrangements also appear to have included genomic deletions that led to the loss of at least one gene in the rodent lineage. Finally, we also have identified a recent mouse-specific segmental duplication. These finding illustrate that small genomic regions outside the large mouse–human conserved segments can contain a single gene as well as sequences that are apparently unique to one genome. The nucleotide sequence data reported in this paper have been submitted to GenBank and assigned the accession numbers AC074224 and AC074041.     

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.     

Novel rearrangements of herpes simplex virus DNA sequences resulting from duplication of a sequence within the unique region of the L component.

We constructed insertion mutants of herpes simplex virus type 1 that contained a duplication of DNA sequences from the BamHI-L fragment (map units 0.706 to 0.744), which is located in the unique region of the L component (UL) of the herpes simplex virus type 1 genome. The second copy of the BamHI-L sequence was inserted in inverted orientation into the viral thymidine kinase gene (map units 0.30 to 0.32), also located within UL. A significant fraction of the progeny produced by these insertion mutants had genomes with rearranged DNA sequences, presumably resulting from intramolecular or intermolecular recombination between the BamHI-L sequences at the two different genomic locations. The rearranged genomes either had an inversion of the DNA sequence flanked by the duplication or were recombinant molecules in which different regions of the genome had been duplicated and deleted. Genomic rearrangements similar to those described here have been reported previously but only for herpes simplex virus insertion mutants containing an extra copy of the repetitive a sequence. Such rearrangements have not been reported for insertion mutants that contain duplications of herpes simplex virus DNA sequences from largely unique regions of the genome. The implications of these results are discussed.     

Structural organization of the barley D-hordein locus in comparison with its orthologous regions of wheat genomes.

D hordein, a prolamin storage protein of barley endosperms, is highly homologous to the high molecular weight (HWM) glutenin subunits, which are the major determinants of bread-making quality in wheat flour. In hexaploid wheat (AABBDD), each genome contains two paralogous copies of HMW-glutenin genes that encode the x- and y-type HMW-glutenin subunits. Previously, we reported the sequence analysis of a 102-kb genomic region that contains the HMW-glutenin locus of the D genome from Aegilops tauschii, the donor of the D genome of hexaploid wheat. Here, we present the sequence analysis of a 120-kb D-hordein region of the barley genome, a more distantly related member of the Triticeae grass tribe. Comparative sequence analysis revealed that gene content and order are generally conserved. Genes included in both of these orthologous regions are arranged in the following order: a Xa21-like receptor kinase, an endosperm globulin, an HMW prolamin, and a serine (threonine) protein kinase. However, in the wheat D genome, a region containing both the globulin and HMW-glutenin gene was duplicated, indicating that this duplication event occurred after the separation of the wheat and barley genomes. The intergenic regions are divergent with regard to the sequence and structural organization. It was found that different types of retroelements are responsible for the intergenic structure divergence in the wheat and barley genomes. In the barley region, we identified 16 long terminal repeat (LTR) retrotransposons in three distinct nested clusters. These retroelements account for 63% of the contig sequence. In addition, barley D hordein was compared with wheat HMW glutenins in terms of cysteine residue conservation and repeat domain organization.     

Comparative Sequence Analysis of the <Emphasis Type="Italic">Phytochrome C</Emphasis> Gene and its Upstream Region in Allohexaploid Wheat Reveals New Data on the Evolution of its Three Constituent Genomes

Bread wheat is an allohexaploid with genome composition AABBDD. Phytochrome C is a gene involved in photomorphogenesis that has been used extensively for phylogenetic analyses. In wheat, the PhyC genes are single copy in each of the three homoeologous genomes and map to orthologous positions on the long arms of the group 5 chromosomes. Comparative sequence analysis of the three homoeologous copies of the wheat PhyC gene and of some 5 kb of upstream region has demonstrated a high level of conservation of PhyC, but frequent interruption of the upstream regions by the insertion of retroelements and other repeats. One of the repeats in the region under investigation appeared to have inserted before the divergence of the diploid wheat genomes, but was degraded to the extent that similarity between the A and D copies could only be observed at the amino acid level. Evidence was found for the differential presence of a foldback element and a miniature inverted-repeat transposable element (MITE) 5′ to PhyC in different wheat cultivars. The latter may represent the first example of an active MITE family in the wheat genome. Several conserved non-coding sequences were also identified that may represent functional regulatory elements. The level of sequence divergence (K_s) between the three wheat PhyC homoeologs suggests that the divergence of the diploid wheat ancestors occurred some 6.9 Mya, which is considerably earlier than the previously estimated 2.5–4.5 Mya. K_a/K_s ratios were <0.15 indicating that all three homoeologs are under purifying selection and presumably represent functional PhyC genes. RT-PCR confirmed expression of the A, B and D copies. The discrepancy in evolutionary age of the wheat genomes estimated using sequences from different parts of the genome may reflect a mosaic origin of some of the Triticeae genomes.     

MicroRNA genes derived from repetitive elements and expanded by segmental duplication events in mammalian genomes

MicroRNAs (miRNAs) are a class of small noncoding RNAs that regulate gene expression by targeting mRNAs for translation repression or mRNA degradation. Many miRNAs are being discovered and studied, but in most cases their origin, evolution and function remain unclear. Here, we characterized miRNAs derived from repetitive elements and miRNA families expanded by segmental duplication events in the human, rhesus and mouse genomes. We applied a comparative genomics approach combined with identifying miRNA paralogs in segmental duplication pair data in a genome-wide study to identify new homologs of human miRNAs in the rhesus and mouse genomes. Interestingly, using segmental duplication pair data, we provided credible computational evidence that two miRNA genes are located in the pseudoautosomal region of the human Y chromosome. We characterized all the miRNAs whether they were derived from repetitive elements or not and identified significant differences between the repeat-related miRNAs (RrmiRs) and non-repeat-derived miRNAs in (1) their location in protein-coding and intergenic regions in genomes, (2) the minimum free energy of their hairpin structures, and (3) their conservation in vertebrate genomes. We found some lineage-specific RrmiR families and three lineage-specific expansion families, and provided evidence indicating that some RrmiR families formed and expanded during evolutionary segmental duplication events. We also provided computational and experimental evidence for the functions of the conservative RrmiR families in the three species. Together, our results indicate that repetitive elements contribute to the origin of miRNAs, and large segmental duplication events could prompt the expansion of some miRNA families, including RrmiR families. Our study is a valuable contribution to the knowledge of evolution and function of non-coding region in genome.     

Identification of a new class of recombinant prolamin genes in wheat.

A novel storage protein gene with obvious [corrected] chimeric structure was isolated from an immature kernel-specific cDNA library prepared from the old Hungarian wheat [corrected] variety, Bánkúti 1201. This clone contains gamma-gliadin sequences in the 5' region and LMW-glutenin sequences on the 3' end. A frameshift mutation was also introduced by the putative recombination event. Hence, the amino acid sequence of the C-terminal region was transformed to a completely new polypeptide. Based on this finding, 7 additional recombinant prolamin genes of similar structure were isolated with specific PCR primers. The 8 chimeric clones seem to be derived from 4 individual gamma-gliadin and 3 LMW-glutenin sequences. These genes show remarkable diversity in size, gliadin:glutenin ratio, frameshift mutations, and sulphur content. The putative functional characteristics of the chimeric polypeptides and problems related to the origin of the encoding genes are discussed.     

Rapid genome evolution revealed by comparative sequence analysis of orthologous regions from four triticeae genomes

Bread wheat (Triticum aestivum) is an allohexaploid species, consisting of three subgenomes (A, B, and D). To study the molecular evolution of these closely related genomes, we compared the sequence of a 307-kb physical contig covering the high molecular weight (HMW)-glutenin locus from the A genome of durum wheat (Triticum turgidum, AABB) with the orthologous regions from the B genome of the same wheat and the D genome of the diploid wheat Aegilops tauschii (Anderson et al., 2003; Kong et al., 2004). Although gene colinearity appears to be retained, four out of six genes including the two paralogous HMW-glutenin genes are disrupted in the orthologous region of the A genome. Mechanisms involved in gene disruption in the A genome include retroelement insertions, sequence deletions, and mutations causing in-frame stop codons in the coding sequences. Comparative sequence analysis also revealed that sequences in the colinear intergenic regions of these different genomes were generally not conserved. The rapid genome evolution in these regions is attributable mainly to the large number of retrotransposon insertions that occurred after the divergence of the three wheat genomes. Our comparative studies indicate that the B genome diverged prior to the separation of the A and D genomes. Furthermore, sequence comparison of two distinct types of allelic variations at the HMW-glutenin loci in the A genomes of different hexaploid wheat cultivars with the A genome locus of durum wheat indicates that hexaploid wheat may have more than one tetraploid ancestor.     

Polyphenol oxidase (PPO) in wheat and wild relatives: molecular evidence for a multigene family

Wheat polyphenol oxidase (PPO) is the major cause of browning reactions that discolor Asian noodles and other wheat products. It has been hypothesized that genes encoding wheat PPOs may have evolved by gene duplication into a multigene family. Here we characterized PPO genomic sequences from diploid (Triticum monococcum, T. urartu, Aegilops tauschii, and Ae. speltoides), tetraploid (T. turgidum, subspecies dicoccoides and durum) and hexaploid (T. aestivum cultivars Klasic and ID377s) wheat species to gain a better understanding of the structure and organization of PPO genes. DNA fragments were amplified from a highly polymorphic and phylogenetic informative region of the gene. As a result, we obtained highly discriminative sequences. Three distinct PPOs, obtained from the A genome of T. monococcum, provided evidence for gene duplication events (paralogous loci). Furthermore, the number of sequences obtained for bread and durum wheat was higher than the expected number of orthologous loci. Sequence comparison revealed nucleotide and structural diversity, and detected five sequence intron types, all with a common insertion position. This was hypothesized to be homologous to that of intron 2 of previously reported wheat PPOs. A MITE of the Stowaway family accounted for the major difference between the five intervening sequences, and was unique to T. aestivum cv. Klasic. Nucleotide and structural diversity, together with well-resolved phylogenetic trees, provided molecular evidence to support the hypothesis of a PPO multigene family structure and organization. Mention of trademark or proprietary products does not constitute a guarantee or warranty of a product by the US Department of Agriculture and does not imply its approval to the exclusion of other products that may also be suitable. This article is in the public domain and not copyrightable.