The Korea brassica genome project: a glimpse of the brassica genome based on comparative genome analysis with Arabidopsis

A complete genome sequence provides unlimited information in the sequenced organism as well as in related taxa. According to the guidance of the Multinational Brassica Genome Project (MBGP), the Korea Brassica Genome Project (KBGP) is sequencing chromosome 1 (cytogenetically oriented chromosome #1) of Brassica rapa. We have selected 48 seed BACs on chromosome 1 using EST genetic markers and FISH analyses. Among them, 30 BAC clones have been sequenced and 18 are on the way. Comparative genome analyses of the EST sequences and sequenced BAC clones from Brassica chromosome 1 revealed their homeologous partner regions on the Arabidopsis genome and a syntenic comparative map between Brassica chromosome 1 and Arabidopsis chromosomes. In silico chromosome walking and clone validation have been successfully applied to extending sequence contigs based on the comparative map and BAC end sequences. In addition, we have defined the (peri)centromeric heterochromatin blocks with centromeric tandem repeats, rDNA and centromeric retrotransposons. In-depth sequence analyses of five homeologous BAC clones and an Arabidopsis chromosomal region reveal overall co-linearity, with 82% sequence similarity. The data indicate that the Brassica genome has undergone triplication and subsequent gene losses after the divergence of Arabidopsis and Brassica. Based on in-depth comparative genome analyses, we propose a comparative genomics approach for conquering the Brassica genome. In 2005 we intend to construct an integrated physical map, including sequence information from 500 BAC clones and integration of fingerprinting data and end sequence data of more than 100 000 BAC clones. The sequences have been submitted to GenBank with accession numbers: 10 204 BAC ends of the KBrH library (CW978640-CW988843); KBrH138P04, AC155338; KBrH117N09, AC155337; KBrH097M21, AC155348; KBrH093K03, AC155347; KBrH081N08, AC155346; KBrH080L24, AC155345; KBrH077A05, AC155343; KBrH020D15, AC155340; KBrH015H17, AC155339; KBrH001H24, AC155335; KBrH080A08, AC155344; KBrH004D11, AC155341; KBrH117M18, AC146875; KBrH052O08, AC155342.     

Sequence-level analysis of the diploidization process in the triplicated FLOWERING LOCUS C region of Brassica rapa

Strong evidence exists for polyploidy having occurred during the evolution of the tribe Brassiceae. We show evidence for the dynamic and ongoing diploidization process by comparative analysis of the sequences of four paralogous Brassica rapa BAC clones and the homologous 124-kb segment of Arabidopsis thaliana chromosome 5. We estimated the times since divergence of the paralogous and homologous lineages. The three paralogous subgenomes of B. rapa triplicated 13 to 17 million years ago (MYA), very soon after the Arabidopsis and Brassica divergence occurred at 17 to 18 MYA. In addition, a pair of BACs represents a more recent segmental duplication, which occurred approximately 0.8 MYA, and provides an exception to the general expectation of three paralogous segments within the B. rapa genome. The Brassica genome segments show extensive interspersed gene loss relative to the inferred structure of the ancestral genome, whereas the Arabidopsis genome segment appears little changed. Representatives of all 32 genes in the Arabidopsis genome segment are represented in Brassica, but the hexaploid complement of 96 has been reduced to 54 in the three subgenomes, with compression of the genomic region lengths they occupy to between 52 and 110 kb. The gene content of the recently duplicated B. rapa genome segments is identical, but intergenic sequences differ.     

Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication

Background

Brassica rapa is one of the most economically important vegetable crops worldwide. Owing to its agronomic importance and phylogenetic position, B. rapa provides a crucial reference to understand polyploidy-related crop genome evolution. The high degree of sequence identity and remarkably conserved genome structure between Arabidopsis and Brassica genomes enables comparative tiling sequencing using Arabidopsis sequences as references to select the counterpart regions in B. rapa, which is a strong challenge of structural and comparative crop genomics.

Results

We assembled 65.8 megabase-pairs of non-redundant euchromatic sequence of B. rapa and compared this sequence to the Arabidopsis genome to investigate chromosomal relationships, macrosynteny blocks, and microsynteny within blocks. The triplicated B. rapa genome contains only approximately twice the number of genes as in Arabidopsis because of genome shrinkage. Genome comparisons suggest that B. rapa has a distinct organization of ancestral genome blocks as a result of recent whole genome triplication followed by a unique diploidization process. A lack of the most recent whole genome duplication (3R) event in the B. rapa genome, atypical of other Brassica genomes, may account for the emergence of B. rapa from the Brassica progenitor around 8 million years ago.

Conclusions

This work demonstrates the potential of using comparative tiling sequencing for genome analysis of crop species. Based on a comparative analysis of the B. rapa sequences and the Arabidopsis genome, it appears that polyploidy and chromosomal diploidization are ongoing processes that collectively stabilize the B. rapa genome and facilitate its evolution.     

New computational tools for brassica genome research

With the increasing quantities of Brassica genomic data being entered into the public domain and in preparation for the complete Brassica genome sequencing effort, there is a growing requirement for the structuring and detailed bioinformatic analysis of Brassica genomic information within a user-friendly database. At the Plant Biotechnology Centre, Melbourne, Australia, we have developed a series of tools and computational pipelines to assist in the processing and structuring of genomic data, to aid its application to agricultural biotechnology research. These tools include a sequence database, ASTRA, a sequence processing pipeline incorporating annotation against GenBank, SwissProt and Arabidopsis Gene Ontology (GO) data and tools for molecular marker discovery and comparative genome analysis. All sequences are mined for simple sequence repeat (SSR) molecular markers using 'SSR primer' and mapped onto the complete Arabidopsis thaliana genome by sequence comparison. The database may be queried using a text-based search of sequence annotation or GO terms, BLAST comparison against resident sequences, or by the position of candidate orthologues within the Arabidopsis genome. Tools have also been developed and applied to the discovery of single nucleotide polymorphism (SNP) molecular markers and the in silico mapping of Brassica BAC end sequences onto the Arabidopsis genome. Planned extensions to this resource include the integration of gene expression data and the development of an EnsEMBL-based genome viewer.     

Comparative physical mapping of segments of the genome of Brassica oleracea var. alboglabra that are homoeologous to sequenced regions of chromosomes 4 and 5 of Arabidopsis thaliana

Due to their relatedness to Arabidopsis thaliana (Arabidopsis), the cultivated Brassica species represent the first group of crops with which to evaluate comparative genomics approaches to understanding biological processes and manipulating traits. We have constructed a high-quality binary BAC library (JBo) from genomic DNA of Brassica oleracea var. alboglabra, in order to underpin such investigations. Using the Arabidopsis genome sequence and clones from the JBo library, we have analysed aspects of gene conservation and microsynteny between a 222 kb region of the genome of Arabidopsis and homoeologous segments of the genome of B. oleracea. All 19 predicted genes tested were found to hybridize to clones in the JBo library, indicating a high level of gene conservation. Further analyses and physical mapping with the BAC clones identified allowed us to construct clone contig maps and analyse in detail the gene content and organization in the set of paralogous segments identified in the genome of B. oleracea. Extensive divergence of gene content was observed, both between the B. oleracea paralogous segments and between them and their homoeologous segment within the genome of Arabidopsis. However, the genes present show highly conserved collinearity with their orthologues in the genome of Arabidopsis. We have identified one example of a Brassica gene in a non-collinear position and one rearrangement. Some of the genes not present in the discernible homoeologous regions appear to be located elsewhere in the B. oleracea genome. The implications of our findings for comparative map-based cloning of genes from crop species are discussed.     

Conservation of the microstructure of genome segments in Brassica napus and its diploid relatives

The cultivated Brassica species are the group of crops most closely related to Arabidopsis thaliana (Arabidopsis). They represent models for the application in crops of genomic information gained in Arabidopsis and provide an opportunity for the investigation of polyploid genome formation and evolution. The scientific literature contains contradictory evidence for the dynamics of the evolution of polyploid genomes. We aimed at overcoming the inherent complexity of Brassica genomes and clarify the effects of polyploidy on the evolution of genome microstructure in specific segments of the genome. To do this, we have constructed bacterial artificial chromosome (BAC) libraries from genomic DNA of B. rapa subspecies trilocularis (JBr) and B. napus var Tapidor (JBnB) to supplement an existing BAC library from B. oleracea. These allowed us to analyse both recent polyploidization (under 10,000 years in B. napus) and more ancient polyploidization events (ca. 20 Myr for B. rapa and B. oleracea relative to Arabidopsis), with an analysis of the events occurring on an intermediate time scale (over the ca. 4 Myr since the divergence of the B. rapa and B. oleracea lineages). Using the Arabidopsis genome sequence and clones from the JBr library, we have analysed aspects of gene conservation and microsynteny between six regions of the genome of B. rapa with the homoeologous regions of the genomes of B. oleracea and Arabidopsis. Extensive divergence of gene content was observed between the B. rapa paralogous segments and their homoeologous segments within the genome of Arabidopsis. A pattern of interspersed gene loss was identified that is similar, but not identical, to that observed in B. oleracea. The conserved genes show highly conserved collinearity with their orthologues across genomes, but a small number of species-specific rearrangements were identified. Thus the evolution of genome microstructure is an ongoing process. Brassica napus is a recently formed polyploid resulting from the hybridization of B. rapa (containing the Brassica A genome) and B. oleracea (containing the Brassica C genome). Using clones from the JBnB library, we have analysed the microstructure of the corresponding segments of the B. napus genome. The results show that there has been little or no change to the microstructure of the analysed segments of the Brassica A and C genomes as a consequence of the hybridization event forming natural B. napus. The observations indicate that, upon polyploid formation, these segments of the genome did not undergo a burst of evolution discernible at the scale of microstructure.     

Development of a high density integrated reference genetic linkage map for the multinational Brassica rapa Genome Sequencing Project

We constructed a high-density Brassica rapa integrated linkage map by combining a reference genetic map of 78 doubled haploid lines derived from Chiifu-401-42?× Kenshin (CKDH) and a new map of 190 F2 lines derived from Chiifu-401-42?× rapid cycling B. rapa (CRF2). The integrated map contains 1017 markers and covers 1262.0 cM of the B. rapa genome, with an average interlocus distance of 1.24 cM. High similarity of marker order and position was observed among the linkage groups of the maps with few short-distance inversions. In total, 155 simple sequence repeat (SSR) markers, anchored to 102 new bacterial artificial chromosomes (BACs) and 146 intron polymorphic (IP) markers were mapped in the integrated map, which would be helpful to align the sequenced BACs in the ongoing multinational Brassica rapa Genome Sequencing Project (BrGSP). Further, comparison of the B. rapa consensus map with the 10 B. juncea A-genome linkage groups by using 98 common IP markers showed high-degree colinearity between the A-genome linkage groups, except for few markers showing inversion or translocation. Suggesting that chromosomes are highly conserved between these Brassica species, although they evolved independently after divergence. The sequence information coming out of BrGSP would be useful for B. juncea breeding. and the identified Arabidopsis chromosomal blocks and known quantitative trait loci (QTL) information of B. juncea could be applied to improve other Brassica crops including B. rapa.     

Comparative fluorescence in situ hybridization mapping of a 431-kb Arabidopsis thaliana bacterial artificial chromosome contig reveals the role of chromosomal duplications in the expansion of the Brassica rapa genome

Comparative genome studies are important contributors to our understanding of genome evolution. Most comparative genome studies in plants have been based on genetic mapping of homologous DNA loci in different genomes. Large-scale comparative physical mapping has been hindered by the lack of efficient and affordable techniques. We report here the adaptation of fluorescence in situ hybridization (FISH) techniques for comparative physical mapping between Arabidopsis thaliana and Brassica rapa. A set of six bacterial artificial chromosomes (BACs) representing a 431-kb contiguous region of chromosome 2 of A. thaliana was mapped on both chromosomes and DNA fibers of B. rapa. This DNA fragment has a single location in the A. thaliana genome, but hybridized to four to six B. rapa chromosomes, indicating multiple duplications in the B. rapa genome. The sizes of the fiber-FISH signals from the same BACs were not longer in B. rapa than those in A. thaliana, suggesting that this genomic region is duplicated but not expanded in the B. rapa genome. The comparative fiber-FISH mapping results support that chromosomal duplications, rather than regional expansion due to accumulation of repetitive sequences in the intergenic regions, played the major role in the evolution of the B. rapa genome.     

Toward unraveling the structure of Brassica rapa genome

Genomic research in any organism encompasses understanding structure of the target genome and genes, their function, and evolution. Brassica rapa , which is phylogenetically related to Arabidopsis thaliana , is an important species with respect to its uses as vegetable, oil, and fodder. The availability of suitable genetic and genomic resources is a prerequisite to undertake genomic research in B. rapa . We have developed reference mapping populations of Chinese cabbage ( B. rapa ssp. pekinensis ) comprising 78 doubled haploid lines and over 250 recombinant inbred lines. Two Bacterial Artificial Chromosome (BAC) libraries, generated by restriction enzymes Hin dIII (KBrH) and Bam HI (KBrB), comprise 56 592 and 50 688 clones, respectively. We have also constructed 22 cDNA libraries from different plant tissues consisting of 104 914 clones with an average length of 575 bp. Initial BAC-end sequence analysis of 1473 clones of the KBrH library led us to understand the structure of B. rapa genome with respect to extent of genic sequences and their annotation, and relative abundance of different types of repetitive DNAs. Full-length sequence analysis of BAC clones revealed extensive triplication of B. rapa DNA segments coupled with variable gene losses within the segments. The formulation of the 'Multinational Brassica Genome Project' has laid the foundation to sequence the complete genome of B. rapa ssp. pekinensis by the international Brassica research community. It has been proposed to undertake BAC-to-BAC sequencing of genetically mapped seed BACs. In recent years, development of bioinformatics tools in Brassica has given a boost to structural genomics research in Brassica species. The research undertaken with the availability of various genomic resources in the public domain has added to our understanding of the structure of B. rapa .     

Structural and functional comparative mapping between the Brassica A genomes in allotetraploid Brassica napus and diploid Brassica rapa

Brassica napus (AACC genome) is an important oilseed crop that was formed by the fusion of the diploids B. rapa (AA) and B. oleracea (CC). The complete genomic sequence of the Brassica A genome will be available soon from the B. rapa genome sequencing project, but it is not clear how informative the A genome sequence in B. rapa (A(r)) will be for predicting the structure and function of the A subgenome in the allotetraploid Brassica species B. napus (A(n)). In this paper, we report the results of structural and functional comparative mapping between the A subgenomes of B. napus and B. rapa based on genetic maps that were anchored with bacterial artificial chromosomes (BACs)-sequence of B. rapa. We identified segmental conservation that represented by syntenic blocks in over one third of the A genome; meanwhile, comparative mapping of quantitative trait loci for seed quality traits identified a dozen homologous regions with conserved function in the A genome of the two species. However, several genomic rearrangement events, such as inversions, intra- and inter-chromosomal translocations, were also observed, covering totally at least 5% of the A genome, between allotetraploid B. napus and diploid B. rapa. Based on these results, the A genomes of B. rapa and B. napus are mostly functionally conserved, but caution will be necessary in applying the full sequence data from B. rapa to the B. napus as a result of genomic rearrangements in the A genome between the two species.     

The CACTA transposon Bot1 played a major role in Brassica genome divergence and gene proliferation

We isolated and characterized a Brassica C genome-specific CACTA element, which was designated Bot1 (Brassica oleracea transposon 1). After analysing phylogenetic relationships, copy numbers and sequence similarity of Bot1 and Bot1 analogues in B. oleracea (C genome) versus Brassica rapa (A genome), we concluded that Bot1 has encountered several rounds of amplification in the oleracea genome only, and has played a major role in the recent rapa and oleracea genome divergence. We performed in silico analyses of the genomic organization and internal structure of Bot1, and established which segment of Bot1 is C-genome specific. Our work reports a fully characterized Brassica repetitive sequence that can distinguish the Brassica A and C chromosomes in the allotetraploid Brassica napus, by fluorescent in situ hybridization. We demonstrated that Bot1 carries a host S locus-associated SLL3 gene copy. We speculate that Bot1 was involved in the proliferation of SLL3 around the Brassica genome. The present study reinforces the assumption that transposons are a major driver of genome and gene evolution in higher plants.     

Isolation of circadian-associated genes in Brassica rapa by comparative genomics with Arabidopsis thaliana

Elucidation of the roles of circadian associated factors requires a better understanding of the molecular mechanisms of circadian rhythms, control of flowering time through photoperiodic pathways, and photosensory signal transduction. In Arabidopsis, the APRR1 quintet, APRRs 1, 3, 5, 7, and 9, are known as central oscillator genes. Other plants may share the molecular mechanism underlying the circadian rhythm. To identify and characterize these circadian response genes in Brassica crops whose genome was triplicated after divergence from Arabidopsis, we identified B. rapa BAC clones containing these genes by BLAST analysis of B. rapa BAC end sequences against the five corresponding Arabidopsis regions. Subsequent fingerprinting, Southern hybridization, and PCR allowed identification of five BAC clones, one for each of the five circadian-related genes. By draft shotgun sequencing of the BAC clones, we identified the complete gene sequences and cloned the five expressed B. rapa circadian-associated gene members, BrPRRs 1, 3, 5, 7, and 9. Phylogenetic analysis revealed that each BrPRR was orthologous to the corresponding APRR at the sequence level. Northern hybridization revealed that the five genes were transcribed at distinct points in the 24 hour period, and Southern hybridization revealed that they are present in 2, 1, 2, 2, and 1 copies, respectively in the B. rapa genome, which was triplicated and then diploidized during the last 15 million years.     

Exploiting comparative mapping among Brassica species to accelerate the physical delimitation of a genic male-sterile locus (BnRf) in Brassica napus

The recessive genic male sterility (RGMS) line 9012AB has been used as an important pollination control system for rapeseed hybrid production in China. Here, we report our study on physical mapping of one male-sterile locus (BnRf) in 9012AB by exploiting the comparative genomics among Brassica species. The genetic maps around BnRf from previous reports were integrated and enriched with markers from the Brassica A7 chromosome. Subsequent collinearity analysis of these markers contributed to the identification of a novel ancestral karyotype block F that possibly encompasses BnRf. Fourteen insertion/deletion markers were further developed from this conserved block and genotyped in three large backcross populations, leading to the construction of high-resolution local genetic maps where the BnRf locus was restricted to a less than 0.1-cM region. Moreover, it was observed that the target region in Brassica napus shares a high collinearity relationship with a region from the Brassica rapa A7 chromosome. A BnRf-cosegregated marker (AT3G23870) was then used to screen a B. napus bacterial artificial chromosome (BAC) library. From the resulting 16 positive BAC clones, one (JBnB089D05) was identified to most possibly contain the BnRf (c) allele. With the assistance of the genome sequence from the Brassica rapa homolog, the 13.8-kb DNA fragment covering both closest flanking markers from the BAC clone was isolated. Gene annotation based on the comparison of microcollinear regions among Brassica napus, B. rapa and Arabidopsis showed that five potential open reading frames reside in this fragment. These results provide a foundation for the characterization of the BnRf locus and allow a better understanding of the chromosome evolution around BnRf.     

An integrated physical and genetic map of the rice genome

Rice was chosen as a model organism for genome sequencing because of its economic importance, small genome size, and syntenic relationship with other cereal species. We have constructed a bacterial artificial chromosome fingerprint–based physical map of the rice genome to facilitate the whole-genome sequencing of rice. Most of the rice genome (~90.6%) was anchored genetically by overgo hybridization, DNA gel blot hybridization, and in silico anchoring. Genome sequencing data also were integrated into the rice physical map. Comparison of the genetic and physical maps reveals that recombination is suppressed severely in centromeric regions as well as on the short arms of chromosomes 4 and 10. This integrated high-resolution physical map of the rice genome will greatly facilitate whole-genome sequencing by helping to identify a minimum tiling path of clones to sequence. Furthermore, the physical map will aid map-based cloning of agronomically important genes and will provide an important tool for the comparative analysis of grass genomes.     

The ABC's of comparative genomics in the Brassicaceae: building blocks of crucifer genomes

In this review we summarize recent advances in our understanding of phylogenetics, polyploidization and comparative genomics in the family Brassicaceae. These findings pave the way for a unified comparative genomic framework. We integrate several of these findings into a simple system of 24 conserved chromosomal blocks (labeled A-X). The naming, order, orientation and color-coding of these blocks are based on their positions in a proposed ancestral karyotype (n=8), rather than by their position in the reduced genome of Arabidopsis thaliana (n=5). We show how these crucifer building blocks can be rearranged to model the genome structures of A. thaliana, Arabidopsis lyrata, Capsella rubella and Brassica rapa. A framework for comparison between species is timely because several crucifer genome-sequencing projects are underway.     

Analysis of expressed sequence tags from Brassica rapa L. ssp. pekinensis

Non-redundant expressed sequence tags (ESTs) were generated from six different organs at various developmental stages of Chinese cabbage, Brassica rapa L. ssp. pekinensis. Of the 1,295 ESTs, 915 (71%) showed significantly high homology in nucleotide or deduced amino acid sequences with other sequences deposited in databases, while 380 did not show similarity to any sequences. Briefly, 598 ESTs matched with proteins of identified biological function, 177 with hypothetical proteins or non-annotated Arabidopsis genome sequences, and 140 with other ESTs. About 82% of the top-scored matching sequences were from Arabidopsis or Brassica, but overall 558 (43%) ESTs matched with Arabidopsis ESTs at the nucleotide sequence level. This observation strongly supports the idea that gene-expression profiles of Chinese cabbage differ from that of Arabidopsis, despite their genome structures being similar to each other. Moreover, sequence analyses of 21 Brassica ESTs revealed that their primary structure is different from those of corresponding annotated sequences of Arabidopsis genes. Our data suggest that direct prediction of Brassica gene expression pattern based on the information from Arabidopsis genome research has some limitations. Thus, information obtained from the Brassica EST study is useful not only for understanding of unique developmental processes of the plant, but also for the study of Arabidopsis genome structure.     

Birth, death and subfunctionalization in the Arabidopsis genome

Arabidopsis thaliana is now a model system, not just for plant biology but also for comparative genomics. The completion of the sequences of two closely related species, Arabidopsis lyrata and Brassica rapa, is complemented by genomic comparisons among A. thaliana accessions and mutation accumulation lines. Together these genomic data document the birth of new genes via gene duplication, transposon exaptation and de novo formation of new genes from noncoding sequence. Most novel loci exhibit low expression, and are undergoing pseudogenization or subfunctionalization. Comparatively, A. thaliana has lost large amounts of sequence through deletion, particularly of transposable elements. Intraspecific genomic variation indicates high rates of deletion mutations and deletion polymorphisms across accessions, shedding light on the history of Arabidopsis genome architecture.     

SynBrowse: a synteny browser for comparative sequence analysis

MOTIVATION: The recent efforts of various sequence projects to sequence deeply into various phylogenies provide great resources for comparative sequence analysis. A generic and portable tool is essential for scientists to visualize and analyze sequence comparisons. RESULTS: We have developed SynBrowse, a synteny browser for visualizing and analyzing genome alignments both within and between species. It is intended to help scientists study macrosynteny, microsynteny and homologous genes between sequences. It can also aid with the identification of uncharacterized genes, putative regulatory elements and novel structural features of a species. SynBrowse is a GBrowse (the Generic Genome Browser) family software tool that runs on top of the open source BioPerl modules. It consists of two components: a web-based front end and a set of relational database back ends. Each database stores pre-computed alignments from a focus sequence to reference sequences in addition to the genome annotations of the focus sequence. The user interface lets end users select a key comparative alignment type and search for syntenic blocks between two sequences and zoom in to view the relationships among the corresponding genome annotations in detail. SynBrowse is portable with simple installation, flexible configuration, convenient data input and easy integration with other components of a model organism system. AVAILABILITY: The software is available at http://www.gmod.org CONTACT: vbrendel@iastate.edu     

Prediction and analysis of paralogous proteins in Trichomonas vaginalis genome

Trichomonas vaginalis causes trichomoniasis, second most sexually transmitted disease. The genome sequence draft of T. vaginalis was published by The Institute of Genomic Research reveals an abnormally large genome size of 160 Mb. It was speculated that a significant portion of the proteome contains paralogous proteins. The present study was aimed at identification and analysis of the paralogous proteins. The all against all search approach is used to identify the paralogous proteins. The dataset of proteins was retrieved from TIGR and TrichDB FTP server. The BLAST-P program performed all against all database searches against the protein database of Trichomonas vaginalis available at NCBI genome database. In the present study about 50,000 proteins were searched where 2,700 proteins were found to be paralogous under the rigid selection criteria. The Pfam database search has identified significant number of paralogous proteins which were further categorized among different 1496 paralogous protein in pfam families, 1027 paralogous protein contains domain, 60 proteins were having different repeats and 1092 paralogous protein sequences of clans. Such identification and functional annotation of paralogous proteins will also help in removing paralogous proteins from possible drug targets in future. Presence of huge number of paralogous proteins across wide range of gene families and domains may be one of the possible mechanisms involved in the T. vaginalis genome expansion and evolution.