Complex evolution of orthologous and paralogous decarboxylase genes

The decarboxylases are involved in neurotransmitter synthesis in animals, and in pathways of secondary metabolism in plants. Different decarboxylase proteins are characterized for their different substrate specificities, but are encoded by homologous genes. We study, within a maximum-likelihood framework, the evolutionary relationships among dopa decarboxylase (Ddc), histidine decarboxylase (Hdc) and alpha-methyldopa hypersensitive (amd) in animals, and tryptophan decarboxylase (Wdc) and tyrosine decarboxylase (Ydc) in plants. The evolutionary rates are heterogeneous. There are differences between paralogous genes in the same lineages: 4.13 x 10(-10) nucleotide substitutions per site per year in mammalian Ddc vs. 1.95 in Hdc; between orthologous genes in different lineages, 7.62 in dipteran Ddc vs. 4.13 in mammalian Ddc; and very large temporal variations in some lineages, from 3.7 up to 54.9 in the Drosophila Ddc lineage. Our results are inconsistent with the molecular clock hypothesis.     

Exploiting the yeast Saccharomyces cerevisiae for the study of the organization and evolution of complex genomes

Yeast artificial chromosome (YAC) cloning systems have advanced the analysis of complex genomes considerably. They permit the cloning of larger fragments than do bacterial artificial chromosome systems, and the cloned material is more easily modified. We recently developed a novel YAC cloning system called transformation-associated recombination (TAR) cloning. Using in vivo recombination in yeast, TAR cloning selectively isolates, as circular YACs, desired chromosome segments or entire genes from complex genomes. The ability to do that without constructing a representative genomic library of random clones greatly facilitates analysis of gene function and its role in disease. In this review, we summarize how recombinational cloning techniques have advanced the study of complex genome organization, gene expression, and comparative genomics.     

Intervening sequences in paralogous genes: a comparative genomic approach to study the evolution of X chromosome introns

The enlargement of the genome size and the decrease in genome compactness with increase in the number and size of introns is a general pattern during the evolution of eukaryotes. Among the possible mechanisms for modifying intron size, it has been suggested that the insertion of transposable elements might have an important role in driving intron evolution. The analysis of large portions of the human genome demonstrated that a relatively recent (50 to 100 MYA) accumulation of transposable elements appears to be biased, favoring a preferential insertion of LINE1 transposons into sex chromosomes rather than into autosomes. In the present work, the effect of chromosomal location on the increase in size of introns was evaluated with a comparative analysis performed on pairs of human paralogous genes, one located on the X chromosome and the second on an autosome. A phylogenetic analysis was also performed on the X-encoded proteins and their paralogs to confirm orthology-paralogy and to approximately estimate the time of gene duplication. Statistical analysis of total intron length for each pair of paralogous genes provided no evidence for a larger size of introns in the gene copies located on the X chromosome. On the opposite, introns of autosomal genes were found to be significantly longer than introns of their X-linked paralogs. Likewise, LINE1 elements were not significantly more frequent in X-chromosome introns, whereas the frequency of SINE elements showed a marginally significant bias toward autosomal introns.     

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.     

Retrotransposons: central players in the structure, evolution and function of plant genomes

One class of mobile DNA, the retrotransposons, constitutes a major portion of all eukaryotic nuclear genomes (often half of the total DNA), and has made a significant contribution towards the genome rearrangements that give rise to altered gene order and novel gene regulation. Research on plant retrotransposons is supported in only a handful of laboratories worldwide, but their progress continues to be striking, as shown at a recent workshop¹.     

Characterizing the composition and evolution of homoeologous genomes in hexaploid wheat through BAC-end sequencing on chromosome 3B

Bread wheat (Triticum aestivum) is one of the most important crops worldwide. However, because of its large, hexaploid, highly repetitive genome it is a challenge to develop efficient means for molecular analysis and genetic improvement in wheat. To better understand the composition and molecular evolution of the hexaploid wheat homoeologous genomes and to evaluate the potential of BAC-end sequences (BES) for marker development, we have followed a chromosome-specific strategy and generated 11 Mb of random BES from chromosome 3B, the largest chromosome of bread wheat. The sequence consisted of about 86% of repetitive elements, 1.2% of coding regions, and 13% remained unknown. With 1.2% of the sequence length corresponding to coding sequences, 6000 genes were estimated for chromosome 3B. New repetitive sequences were identified, including a Triticineae-specific tandem repeat (Fat) that represents 0.6% of the B-genome and has been differentially amplified in the homoeologous genomes before polyploidization. About 10% of the BES contained junctions between nested transposable elements that were used to develop chromosome-specific markers for physical and genetic mapping. Finally, sequence comparison with 2.9 Mb of random sequences from the D-genome of Aegilops tauschii suggested that the larger size of the B-genome is due to a higher content in repetitive elements. It also indicated which families of transposable elements are mostly responsible for differential expansion of the homoeologous wheat genomes during evolution. Our data demonstrate that BAC-end sequencing from flow-sorted chromosomes is a powerful tool for analysing the structure and evolution of polyploid and highly repetitive genomes.     

AFLP-derived SCARs facilitate construction of a 1.1 Mb sequence-ready map of a region that spans the Vf locus in the apple genome

The availability of high-density anchored markers is a prerequisite for reliable construction of a deep coverage BAC contig, which leads to creation of a sequence-ready map in the target chromosomal region. Unfortunately, such markers are not available for most plant species, including woody perennial plants. Here, we report on an efficient approach to build a megabase-size sequence-ready map in the apple genome for the Vf region containing apple scab resistance gene(s) by targeting AFLP-derived SCAR markers to this specific genomic region. A total of 11 AFLP-derived SCAR markers, previously tagged to the Vf locus, along with three other Vf-linked SCAR markers have been used to screen two apple genome BAC libraries. A single BAC contig which spans the Vf region at a physical distance of approximately 1,100 kb has been constructed by assembling the recovered BAC clones, followed by closure of inter-contig gaps. The contig is 4 ×deep, and provides a minimal tiling path of 16 contiguous and overlapping BAC clones, thus generating a sequence-ready map. Within the Vf region, duplication events have occurred frequently, and the Vf locus is restricted to the ca. 290 kb region covered by a minimum of three overlapping BAC clones.     

Multiplex sequencing of bacterial artificial chromosomes for assembling complex plant genomes

Hierarchical shotgun sequencing remains the method of choice for assembling high‐quality reference sequences of complex plant genomes. The efficient exploitation of current high‐throughput technologies and powerful computational facilities for large‐insert clone sequencing necessitates the sequencing and assembly of a large number of clones in parallel. We developed a multiplexed pipeline for shotgun sequencing and assembling individual bacterial artificial chromosomes (BACs) using the Illumina sequencing platform. We illustrate our approach by sequencing 668 barley BACs (Hordeum vulgare L.) in a single Illumina HiSeq 2000 lane. Using a newly designed parallelized computational pipeline, we obtained sequence assemblies of individual BACs that consist, on average, of eight sequence scaffolds and represent >98% of the genomic inserts. Our BAC assemblies are clearly superior to a whole‐genome shotgun assembly regarding contiguity, completeness and the representation of the gene space. Our methods may be employed to rapidly obtain high‐quality assemblies of a large number of clones to assemble map‐based reference sequences of plant and animal species with complex genomes by sequencing along a minimum tiling path.     

两种HMW-GS基因启动子片段的克隆及分析

利用聚合酶链反应(PCR)技术从小偃6号中获得400bp左右的扩增产物,将其与pGEM-T Easy载体连接后转入大肠杆菌,经过筛选获得HMW-8-P和HMW-38-P两种类型克隆。序列分析表明：HMW-38-P包括了HMW-GS14基因上游启动子及信号肽对应编码区,而另一段(HMW-8-P)为一未知HMW-GS基因启动子区及信号肽对应的编码区。将两序列和GenBank中已知的35种HWM-GS基因启动子区序列进行多序列比对,最后获得HMW-GS启动子的系统发生树。通过系统发生树可以清晰地看出位于不同染色体上的不同亚基类型的HMW-GS基因的进化关系,并可确定HMW-8-P为Glu-D-1类型HMW-GS的启动子区,小偃6号中Glu-D-1类型的亚基为2亚基,所以HMW-6-P为2亚基启动子区序列。     

稻属不同染色体组着丝粒BAC克隆的分离和鉴定

着丝粒在真核生物有丝分裂和减数分裂染色体正常的分离和传递中起着重要的作用。通过构建5个稻属二倍体野生种的基因组BAC文库, 采用菌落杂交和FISH技术, 筛选和鉴定了各染色体组着丝粒克隆, 并且分析了这些克隆在不同基因组间的共杂交情况, 结果表明: (1) C染色体组的野生种O. officinalis 和F染色体组的野生种O. brachyantha具有各自着丝粒特异的卫星DNA序列, 并且O. brachyantha着丝粒还具有特异的逆转座子序列; (2) A、B和E染色体组的野生稻O. glaberrima、O. punctata和O. australiensis着丝粒区域都含有与栽培稻着丝粒重复序列CentO和CRR同源的序列; (3) C染色体组野生稻O. officinalis的2条体细胞染色体着丝粒具有CentO的同源序列, 同时也发现其所有着丝粒区域都包含栽培稻CRR的同源序列。这些结果对克隆稻属不同染色体组的着丝粒序列、研究不同染色体组间着丝粒的进化关系和稻属不同着丝粒DNA序列与功能之间的关系均具有重要意义。     

Closing a gap in the physical map of the ovine major histocompatibility complex

A 184 kb gap in an ovine MHC physical map was successfully closed by identification of two overlapping clones (304C7 and 222G18) from a Chinese fine wool merino sheep BAC library. The location and tiling path of the two clones were confirmed by BAC‐end sequencing and PCR amplification of loci in overlapping regions. Full‐length sequencing of the clones identified 13 novel ovine genes in the gap between loci Notch4 and Btnl2, and eight of them belonging to the Butyrophilin‐like (Btn‐like or Btnl) gene family. The scattered distribution of the Btnl gene cluster at the gap provided a clue to explain the difficulties previously experienced in closing the gap. Completed BAC contigs of the ovine MHC will facilitate sequencing of the entire ovine leukocyte antigen (OLA) region, providing detailed information for comparative studies of MHC evolution.     

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

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

Synthetic hexaploid wheat enhances variation and adaptive evolution of bread wheat in breeding processes

Synthetic hexaploid wheat (SHW) that combines novel and elite genes from the tetraploid wheat Triticum turgidum L. and wild ancestor Aegilops tauschii Coss., has been used to genetically improve hexaploid common wheat. The abundant genetic diversity in SHW can effectively make breakthroughs in wheat genetic improvement through the inclusion of increased variation. In this paper, we reviewed the current advances in research and utilization of the primary SHW lines and SHW-derived wheat varieties that have enhanced evolution of modern wheat under conditions of natural and artificial selection in southwestern China. Using primary SHW lines, four high-yielding wheat varieties have been developed. In addition, using the SHW-derived varieties as breeding parents, 12 new wheat varieties were also developed. Results of genotype–phenotype and fingerprint analysis showed that the introgressed alleles from SHW lines have contributed a great number of elite characters to the new wheat varieties, and these elite characters include disease resistance, more spikes per plant, more grains per spike, larger grains, and higher grain-yield potential. We found that the primary SHW lines and SHW-derived varieties have identifiable effects to enhance genetic variation and adaptive evolution of modern hexaploid wheat, which significantly increased the grain yields of hexaploid wheat in recent years. These findings have significant implications in the breeding of high-yielding wheat varieties resistant to biotic and abiotic stresses using SHW as genetic resources.     

PCR analysis of x- and y-type genes present at the complex Glu-A1 locus in durum and bread wheat

Genes (x-type) corresponding to different high-molecular-weight glutenin subunits encoded at the Glu-A1 locus present in bread- and durum-wheat cultivars have been selectively amplified by the polymerase chain reaction (PCR). DNA fragments corresponding to an unexpressed x-type gene were also amplified. As unexpressed y-type genes may or may not contain an 8-kb transposon-like insertion, two different sets of primers were designed to obtain amplification of DNA fragments corresponding to these genes. Amplified DNA fragments were also digested with restriction enzymes. The digestion patterns of amplified fragments corresponding to unusual x-type subunits showed similarities with genes encoding the most common subunits 2^* and 1. The unexpressed amplified x-type gene showed a restriction pattern similar to the one obtained with the allelic gene encoding high-molecular-weight glutenin subunit 1; homologies were also found within the repetitive region of the linked y-type genes. On the basis of these observations it is postulated that an ancestral active x-type gene, most likely corresponding to subunit 1, was silenced following the insertion of the 8-kb transposon-like fragment into the linked y-type gene. Received: 8 April 1996 / Accepted: 30 August 1996     

Comparative analysis of complete orthologous centromeres from two subspecies of rice reveals rapid variation of centromere organization and structure

Centromeres are sites for assembly of the chromosomal structures that mediate faithful segregation at mitosis and meiosis. This function is conserved across species, but the DNA components that are involved in kinetochore formation differ greatly, even between closely related species. To shed light on the nature, evolutionary timing and evolutionary dynamics of rice centromeres, we decoded a 2.25‐Mb DNA sequence covering the centromeric region of chromosome 8 of an indica rice variety, ‘Kasalath’ (Kas‐Cen8). Analysis of repetitive sequences in Kas‐Cen8 led to the identification of 222 long terminal repeat (LTR)‐retrotransposon elements and 584 CentO satellite monomers, which account for 59.2% of the region. A comparison of the Kas‐Cen8 sequence with that of japonica rice ‘Nipponbare’ (Nip‐Cen8) revealed that about 66.8% of the Kas‐Cen8 sequence was collinear with that of Nip‐Cen8. Although the 27 putative genes are conserved between the two subspecies, only 55.4% of the total LTR‐retrotransposon elements in ‘Kasalath’ had orthologs in ‘Nipponbare’, thus reflecting recent proliferation of a considerable number of LTR‐retrotransposons since the divergence of two rice subspecies of indica and japonica within Oryza sativa. Comparative analysis of the subfamilies, time of insertion, and organization patterns of inserted LTR‐retrotransposons between the two Cen8 regions revealed variations between ‘Kasalath’ and ‘Nipponbare’ in the preferential accumulation of CRR elements, and the expansion of CentO satellite repeats within the core domain of Cen8. Together, the results provide insights into the recent proliferation of LTR‐retrotransposons, and the rapid expansion of CentO satellite repeats, underlying the dynamic variation and plasticity of plant centromeres.     

Potato flour viscosity improvement is associated with the expression of a wheat LMW-glutenin gene

It has been previously shown that expression of a high-molecular-weight glutenin (HMW-GS) in transgenic wheat seeds resulted in the improvement of flour functional properties. In this study, potato flour viscosity was improved through a specific expression of a low-molecular-weight glutenin (LMW-GS-MB1) gene in tuber. The resulting construct was introduced into potato leaf explants (Solanum tuberosum cv Kennebec) through Agrobacterium tumefaciens-mediated gene transfer. Southern and Northern analysis of transgenic potato confirmed that the integration of LMW-GS-MB1 in genomic DNA was stable and its mRNA was abundant in transgenic line 16 tubers. Western blot analysis of line 16 extract shows a LMW-GS subunit accumulation in tuber. To demonstrate the capacity of transgenic lines to produce tubers with improved flour functional properties, transgenic lines 9 and 16 exhibiting, respectively, moderate and high expression of LMW-GS-MB1 mRNA and nontransgenic plants were transferred to field plots. The mean viscosity value of flour obtained from the field-grown tubers of transgenic line 16 exhibited a 3-fold increase in viscosity at 23 degrees C when compared to flour from nontransgenic tubers.     

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.     

Structure of the catfish IGH locus: analysis of the region including the single functional IGHM gene

The catfish IGH locus is large (∼1 Mb) and complex, having undergone multiple internal duplications and transpositions. To define the structure of the locus that contains the single expressed IGHM gene, two overlapping bacterial-artificial-chromosome (BAC) clones spanning the most 3′ end of the channel catfish immunoglobulin heavy (IGH) chain locus have been completely sequenced. The analyses created a contig of 257,153 bp containing 55 VH, 6 D, 12 JH genes and the IGH constant region genes encoding the functional secreted and membrane forms of IgM and the membrane form of IgD. This analysis revealed three major features. First, no C-region genes were found aside from the previously described IGHM1 and IGHD1, with the latter gene being the most 3′ C-region gene of the catfish IGH locus. There was no evidence in the region sequenced for genes that could encode an Ig class similar to the IgZ/IgT described in zebrafish, trout and pufferfish. Second, there are a high number of VH pseudogenes, 28 out of 55 (51%). In contrast, the entire zebrafish IGH locus has 40 functional VH genes and eight pseudogenes (17%). Third, an internal duplication of a 52.4-kb block of VH genes has occurred. These observations suggest that the IGH locus of teleost fish varies significantly from species to species in the diversity of C-region genes as well as the numbers of genes encoding V regions.