Synteny between Arabidopsis thaliana and rice at the genome level: a tool to identify conservation in the ongoing rice genome sequencing project

BLASTX alignment between 189.5 Mb of rice genomic sequence and translated Arabidopsis thaliana annotated coding sequences (CDS) identified 60 syntenic regions involving 4–22 rice orthologs covering ≤3.2 cM (centiMorgan). Most regions are <3 cM in length. A detailed and updated version of a table representing these regions is available on our web site. Thirty-five rice loci match two distinct A.thaliana loci, as expected from the duplicated nature of the A.thaliana genome. One A.thaliana locus matches two distinct rice regions, suggesting that rice chromosomal sequence duplications exist. A high level of rearrangement characterizing the 60 syntenic regions illustrates the ancient nature of the speciation between A.thaliana and rice. The apparent reduced level of microcollinearity implies the dispersion to new genomic locations, via transposon activity, of single or small clusters of genes in the rice genome, which represents a significant additional effector of plant genome evolution.     

Long-Range Genomic Enrichment,Sequencing, and Assembly to Determine Unknown Sequences Flanking a Known microRNA

Conserved plant microRNAs (miRNAs) modulate important biological processes but little is known about conserved cis-regulatory elements (CREs) surrounding MIRNA genes. We developed a solution-based targeted genomic enrichment methodology to capture, enrich, and sequence flanking genomic regions surrounding conserved MIRNA genes with a locked-nucleic acid (LNA)-modified, biotinylated probe complementary to the mature miRNA sequence. Genomic DNA bound by the probe is captured by streptavidin-coated magnetic beads, amplified, sequenced and assembled de novo to obtain genomic DNA sequences flanking MIRNA locus of interest. We demonstrate the sensitivity and specificity of this enrichment methodology in Arabidopsis thaliana to enrich targeted regions spanning 10–20 kb surrounding known MIR166 and MIR165 loci. Assembly of the sequencing reads successfully recovered all targeted loci. While further optimization for larger, more complex genomes is needed, this method may enable determination of flanking genomic DNA sequence surrounding a known core (like a conserved mature miRNA) from multiple species that currently don''t have a full genome assembly available.     

FLAGdb/FST: a database of mapped flanking insertion sites (FSTs) of Arabidopsis thaliana T-DNA transformants

A large collection of T-DNA insertion transformants of Arabidopsis thaliana has been generated at the Institute of Agronomic Research, Versailles, France. The molecular characterisation of the insertion sites is currently performed by sequencing genomic regions flanking the inserted T-DNA (FST). The almost complete sequence of the nuclear genome of A.thaliana provides the framework for organising FSTs in a genome oriented database, FLAGdb/FST (http://genoplante-info.infobiogen.fr). The main scope of FLAGdb/FST is to help biologists to find the FSTs that interrupt the genes in which they are interested. FSTs are anchored to the genome sequences of A.thaliana and positions of both predicted genes and FSTs are shown graphically on sequences. Requests to locate the genomic position of a query sequence are made using BLAST programs. The response delivered by FLAGdb/FST is a graphical representation of the putative FSTs and of predicted genes in a 20 kb region.     

Genetic and molecular characterization of embryonic mutants identified following seed transformation in Arabidopsis

Over 5000 transgenic families of Arabidopsis thaliana produced following seed transformation with Agrobacterium tumefaciens were screened for embryonic lethals, defectives, and pattern mutants. One hundred and seventy-eight mutants with a wide range of developmental abnormalities were identified. Forty-one mutants appear from genetic studies to be tagged (36% of the 115 mutants examined in detail). Mapping with visible markers demonstrated that mutant genes were randomly distributed throughout the genome. Seven mutant families appeared to contain chromosomal translocations because the mutant genes exhibited linkage to visible markers on two different chromosomes. Chromosomal rearrangements may therefore be widespread following seed transformation. DNA gel blot hybridizations with 34 tagged mutants and three T-DNA probes revealed a wide range of insertion patterns. Models of T-DNA structure at each mutant locus were constructed to facilitate gene isolation. The value of such models was demonstrated by using plasmid rescue to clone flanking plant DNA from four tagged mutants. Further analysis of genes isolated from these insertional mutants should help to elucidate the relationship between gene function and plant embryogenesis.     

Comparative genomics and reverse genetics analysis reveal indispensable functions of the serine acetyltransferase gene family in Arabidopsis

Ser acetyltransferase (SERAT), which catalyzes O-acetyl-Ser (OAS) formation, plays a key role in sulfur assimilation and Cys synthesis. Despite several studies on SERATs from various plant species, the in vivo function of multiple SERAT genes in plant cells remains unaddressed. Comparative genomics studies with the five genes of the SERAT gene family in Arabidopsis thaliana indicated that all three Arabidopsis SERAT subfamilies are conserved across five plant species with available genome sequences. Single and multiple knockout mutants of all Arabidopsis SERAT gene family members were analyzed. All five quadruple mutants with a single gene survived, with three mutants showing dwarfism. However, the quintuple mutant lacking all SERAT genes was embryo-lethal. Thus, all five isoforms show functional redundancy in vivo. The developmental and compartment-specific roles of each SERAT isoform were also demonstrated. Mitochondrial SERAT2;2 plays a predominant role in cellular OAS formation, while plastidic SERAT2;1 contributes less to OAS formation and subsequent Cys synthesis. Three cytosolic isoforms, SERAT1;1, SERAT3;1, and SERAT3;2, may play a major role during seed development. Thus, the evolutionally conserved SERAT gene family is essential in cellular processes, and the substrates and products of SERAT must be exchangeable between the cytosol and organelles.     

Sequence and structure of Brassica rapa chromosome A3

Background

The species Brassica rapa includes important vegetable and oil crops. It also serves as an excellent model system to study polyploidy-related genome evolution because of its paleohexaploid ancestry and its close evolutionary relationships with Arabidopsis thaliana and other Brassica species with larger genomes. Therefore, its genome sequence will be used to accelerate both basic research on genome evolution and applied research across the cultivated Brassica species.

Results

We have determined and analyzed the sequence of B. rapa chromosome A3. We obtained 31.9 Mb of sequences, organized into nine contigs, which incorporated 348 overlapping BAC clones. Annotation revealed 7,058 protein-coding genes, with an average gene density of 4.6 kb per gene. Analysis of chromosome collinearity with the A. thaliana genome identified conserved synteny blocks encompassing the whole of the B. rapa chromosome A3 and sections of four A. thaliana chromosomes. The frequency of tandem duplication of genes differed between the conserved genome segments in B. rapa and A. thaliana, indicating differential rates of occurrence/retention of such duplicate copies of genes. Analysis of 'ancestral karyotype' genome building blocks enabled the development of a hypothetical model for the derivation of the B. rapa chromosome A3.

Conclusions

We report the near-complete chromosome sequence from a dicotyledonous crop species. This provides an example of the complexity of genome evolution following polyploidy. The high degree of contiguity afforded by the clone-by-clone approach provides a benchmark for the performance of whole genome shotgun approaches presently being applied in B. rapa and other species with complex genomes.     

The complete nucleotide sequence of the mitochondrial genome of sugar beet (Beta vulgaris L.) reveals a novel gene for tRNA(Cys)(GCA)

We determined the complete nucleotide sequence of the mitochondrial genome of an angiosperm, sugar beet (Beta vulgaris cv TK81-O). The 368 799 bp genome contains 29 protein, five rRNA and 25 tRNA genes, most of which are also shared by the mitochondrial genome of Arabidopsis thaliana, the only other completely sequenced angiosperm mitochondrial genome. However, four genes identified here (namely rps13, trnF-GAA, ccb577 and trnC2-GCA) are missing in Arabidopsis mitochondria. In addition, four genes found in Arabidopsis (ccb228, rpl2, rpl16 and trnY2-GUA) are entirely absent in sugar beet or present only in severely truncated form. Introns, duplicated sequences, additional reading frames and inserted foreign sequences (chloroplast, nuclear and plasmid DNA sequences) contribute significantly to the overall size of the sugar beet mitochondrial genome. Nevertheless, 55.6% of the genome has no obvious features of information. We identified a novel tRNA^Cys gene (trnC2-GCA) which shows no sequence homology with any tRNA^Cys genes reported so far in higher plants. Intriguingly, this tRNA gene is actually transcribed into a mature tRNA, whereas the native tRNA^Cys gene (trnC1-GCA) is most likely a pseudogene.     

MIPS Arabidopsis thaliana Database (MAtDB): an integrated biological knowledge resource based on the first complete plant genome

Arabidopsis thaliana is the first plant for which the complete genome has been sequenced and published. Annotation of complex eukaryotic genomes requires more than the assignment of genetic elements to the sequence. Besides completing the list of genes, we need to discover their cellular roles, their regulation and their interactions in order to understand the workings of the whole plant. The MIPS Arabidopsis thaliana Database (MAtDB; http://mips.gsf.de/proj/thal/db) started out as a repository for genome sequence data in the European Scientists Sequencing Arabidopsis (ESSA) project and the Arabidopsis Genome Initiative. Our aim is to transform MAtDB into an integrated biological knowledge resource by integrating diverse data, tools, query and visualization capabilities and by creating a comprehensive resource for Arabidopsis as a reference model for other species, including crop plants.     

Comparative genomics reveals expansion of the FLC region in the genus Arabidopsis

Mechanisms of genome evolution are poorly understood although recent genome sequencing is providing the tools to begin to illuminate such mechanisms. Using high-resolution molecular cytogenetic tools, we examined the structural evolution of 790 kb surrounding the evolutionarily important FLC locus of Arabidopsis thaliana in three of its relatives, Arabidopsis halleri, Arabidopsis neglecta and Arabidopsis arenosa. Sequenced BACs from A. thaliana were used as heterologous probes across these species and genome expansion was found in all three species relative to A. thaliana, ranging from 16 to 27%. Expansion was seen along the length of the entire region but molecular analyses revealed no characteristic pattern of either intra- or intergenic expansion among these species. Mapping of BACs on DNA fibers from A. thaliana revealed one possible error, ~14 kb missing from the reported sequence, indicating that for comparative studies it is important to confirm the reference sequence to which comparison will be made.     

LEA Gene Introns: is the Intron of Dehydrin Genes a Characteristic of the Serine-Segment?

Dehydrins, which belong to group 2 LEA proteins, are a family of intrinsically unstructured plant proteins that accumulate during the late stages of embryogenesis and in response to abiotic stresses. We have previously reported that the OpsDHN1 gene, encoding an SK₃-type acidic dehydrin protein from Opuntia streptacantha, contains an intron inserted within the sequence encoding the S-motif. Herein, we present an in silico analysis of intron sequences in dehydrin genes from mono- and dicotyledonous plants that reveals a preference for insertion within the nucleotide sequence encoding the S-motif. Sequence comparison of ten Dhn genes from Arabidopsis thaliana and the orthologous genes in Arabidopsis lyrata revealed that introns maintain considerable sequence identity and conserve the insertion pattern. Furthermore, syntenic regions were identified among eight orthologous genes of A. thaliana and A. lyrata, showing that correlated gene arrangements are conserved between these Arabidopsis species. Our study shows that most SKn-type dehydrins contain one intron that is conserved in phase and location; this intron is linked to the nucleotide sequence that encodes the S-motif.     

A first survey of the rye (Secale cereale) genome composition through BAC end sequencing of the short arm of chromosome 1R

Background

Thellungiella halophila (also known as Thellungiella salsuginea) is a model halophyte with a small plant size, short life cycle, and small genome. It easily undergoes genetic transformation by the floral dipping method used with its close relative, Arabidopsis thaliana. Thellungiella genes exhibit high sequence identity (approximately 90% at the cDNA level) with Arabidopsis genes. Furthermore, Thellungiella not only shows tolerance to extreme salinity stress, but also to chilling, freezing, and ozone stress, supporting the use of Thellungiella as a good genomic resource in studies of abiotic stress tolerance.

Results

We constructed a full-length enriched Thellungiella (Shan Dong ecotype) cDNA library from various tissues and whole plants subjected to environmental stresses, including high salinity, chilling, freezing, and abscisic acid treatment. We randomly selected about 20 000 clones and sequenced them from both ends to obtain a total of 35 171 sequences. CAP3 software was used to assemble the sequences and cluster them into 9569 nonredundant cDNA groups. We named these cDNAs "RTFL" (RIKEN Thellungiella Full-Length) cDNAs. Information on functional domains and Gene Ontology (GO) terms for the RTFL cDNAs were obtained using InterPro. The 8289 genes assigned to InterPro IDs were classified according to the GO terms using Plant GO Slim. Categorical comparison between the whole Arabidopsis genome and Thellungiella genes showing low identity to Arabidopsis genes revealed that the population of Thellungiella transport genes is approximately 1.5 times the size of the corresponding Arabidopsis genes. This suggests that these genes regulate a unique ion transportation system in Thellungiella.

Conclusion

As the number of Thellungiella halophila (Thellungiella salsuginea) expressed sequence tags (ESTs) was 9388 in July 2008, the number of ESTs has increased to approximately four times the original value as a result of this effort. Our sequences will thus contribute to correct future annotation of the Thellungiella genome sequence. The full-length enriched cDNA clones will enable the construction of overexpressing mutant plants by introduction of the cDNAs driven by a constitutive promoter, the complementation of Thellungiella mutants, and the determination of promoter regions in the Thellungiella genome.     

Loss/retention and evolution of NBS-encoding genes upon whole genome triplication of Brassica rapa

A genome triplication took place in the ancestor of Brassiceae species after the split of the Arabidopsis lineage. The postfragmentation and shuffling of the genome turned the ancestral hexaploid back to diploids and caused the radiation of Brassiceae species. The course of speciation was accompanied by the loss of duplicate genes and also influenced the evolution of retained genes. Of all the genes, those encoding NBS domains are typical R genes that confer resistance to invading pathogens. In this study, using the genome of Arabidopsis thaliana as a reference, we examined the loss/retention of orthologous NBS-encoding loci in the tripled Brassica rapa genome and discovered differential loss/retention frequencies. Further analysis indicated that loci of different retention ratios showed different evolutionary patterns. The loci of classesII and III (maintaining two and three syntenic loci, respectively, multi-loci) show sharper expansions by tandem duplications, have faster evolutionary rates and have more potential to be associated with novel gene functions. On the other hand, the loci that are retained at the minimal rate (keeping only one locus, class I, single locus) showed opposite patterns. Phylogenetic analysis indicated that recombination and translocation events were common among multi-loci in B. rapa, and differential evolutionary patterns between multi- and single-loci are likely the consequence of recombination. Investigations towards other gene families demonstrated different evolutionary characteristics between different gene families. The evolution of genes is more likely determined by the property of each gene family, and the whole genome triplication provided only a specific condition.