Dataset Information

Genome assembly has a major impact on gene content: a comparison of annotation in two Bos taurus assemblies.

ABSTRACT: Gene and SNP annotation are among the first and most important steps in analyzing a genome. As the number of sequenced genomes continues to grow, a key question is: how does the quality of the assembled sequence affect the annotations? We compared the gene and SNP annotations for two different Bos taurus genome assemblies built from the same data but with significant improvements in the later assembly. The same annotation software was used for annotating both sequences. While some annotation differences are expected even between high-quality assemblies such as these, we found that a staggering 40% of the genes (>9,500) varied significantly between assemblies, due in part to the availability of new gene evidence but primarily to genome mis-assembly events and local sequence variations. For instance, although the later assembly is generally superior, 660 protein coding genes in the earlier assembly are entirely missing from the later genome's annotation, and approximately 3,600 (15%) of the genes have complex structural differences between the two assemblies. In addition, 12-20% of the predicted proteins in both assemblies have relatively large sequence differences when compared to their RefSeq models, and 6-15% of bovine dbSNP records are unrecoverable in the two assemblies. Our findings highlight the consequences of genome assembly quality on gene and SNP annotation and argue for continued improvements in any draft genome sequence. We also found that tracking a gene between different assemblies of the same genome is surprisingly difficult, due to the numerous changes, both small and large, that occur in some genes. As a side benefit, our analyses helped us identify many specific loci for improvement in the Bos taurus genome assembly.

SUBMITTER: Florea L

PROVIDER: S-EPMC3120881 | biostudies-literature | 2011

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Genome assembly has a major impact on gene content: a comparison of annotation in two Bos taurus assemblies.

Florea Liliana L Souvorov Alexander A Kalbfleisch Theodore S TS Salzberg Steven L SL

PloS one 20110622 6

Gene and SNP annotation are among the first and most important steps in analyzing a genome. As the number of sequenced genomes continues to grow, a key question is: how does the quality of the assembled sequence affect the annotations? We compared the gene and SNP annotations for two different Bos taurus genome assemblies built from the same data but with significant improvements in the later assembly. The same annotation software was used for annotating both sequences. While some annotation dif ...[more]

PMID: 21731731

Similar Datasets

Project description:BackgroundReference genomes are essential in the analysis of genomic data. As the cost of sequencing decreases, multiple reference genomes are being produced within species to alleviate problems such as low mapping accuracy and reference allele bias in variant calling that can be associated with the alignment of divergent samples to a single reference individual. The latest reference sequence adopted by the scientific community for the analysis of cattle data is ARS_UCD1.2, built from the DNA of a Hereford cow (Bos taurus taurus-B. taurus). A complementary genome assembly, UOA_Brahman_1, was recently built to represent the other cattle subspecies (Bos taurus indicus-B. indicus) from a Brahman cow haplotype to further support analysis of B. indicus data. In this study, we aligned the sequence data of 15 B. taurus and B. indicus breeds to each of these references.ResultsThe alignment of B. taurus individuals against UOA_Brahman_1 detected up to five million more single-nucleotide variants (SNVs) compared to that against ARS_UCD1.2. Similarly, the alignment of B. indicus individuals against ARS_UCD1.2 resulted in one and a half million more SNVs than that against UOA_Brahman_1. The number of SNVs with nearly fixed alternative alleles also increased in the alignments with cross-subspecies. Interestingly, the alignment of B. taurus cattle against UOA_Brahman_1 revealed regions with a smaller than expected number of counts of SNVs with nearly fixed alternative alleles. Since B. taurus introgression represents on average 10% of the genome of Brahman cattle, we suggest that these regions comprise taurine DNA as opposed to indicine DNA in the UOA_Brahman_1 reference genome. Principal component and admixture analyses using genotypes inferred from this region support these taurine-introgressed loci. Overall, the flagged taurine segments represent 13.7% of the UOA_Brahman_1 assembly. The genes located within these segments were previously reported to be under positive selection in Brahman cattle, and include functional candidate genes implicated in feed efficiency, development and immunity.ConclusionsWe report a list of taurine segments that are in the UOA_Brahman_1 assembly, which will be useful for the interpretation of interesting genomic features (e.g., signatures of selection, runs of homozygosity, increased mutation rate, etc.) that could appear in future re-sequencing analysis of indicine cattle.

Project description:BACKGROUND:Methane emissions by methanogen from livestock ruminants have significantly contributed to the agricultural greenhouse gas effect. It is worthwhile to compare methanogen from "energy-saving" animal (yak) and normal animal (cattle) in order to investigate the link between methanogen structure and low methane production. RESULTS:Diversity of methanogens from the yak and cattle rumen was investigated by analysis of 16S rRNA gene sequences from rumen digesta samples from four yaks (209 clones) and four cattle (205 clones) from the Qinghai-Tibetan Plateau area (QTP). Overall, a total of 414 clones (i.e. sequences) were examined and assigned to 95 operational taxonomic units (OTUs) using MOTHUR, based upon a 98% species-level identity criterion. Forty-six OTUs were unique to the yak clone library and 34 OTUs were unique to the cattle clone library, while 15 OTUs were found in both libraries. Of the 95 OTUs, 93 putative new species were identified. Sequences belonging to the Thermoplasmatales-affiliated Linage C (TALC) were found to dominate in both libraries, accounting for 80.9% and 62.9% of the sequences from the yak and cattle clone libraries, respectively. Sequences belonging to the Methanobacteriales represented the second largest clade in both libraries. However, Methanobrevibacter wolinii (QTPC 110) was only found in the cattle library. The number of clones from the order Methanomicrobiales was greater in cattle than in the yak clone library. Although the Shannon index value indicated similar diversity between the two libraries, the Libshuff analysis indicated that the methanogen community structure of the yak was significantly different than those from cattle. CONCLUSION:This study revealed for the first time the molecular diversity of methanogen community in yaks and cattle in Qinghai-Tibetan Plateau area in China. From the analysis, we conclude that yaks have a unique rumen microbial ecosystem that is significantly different from that of cattle, this may also help to explain why yak produce less methane than cattle.

Dataset Information

Genome assembly has a major impact on gene content: a comparison of annotation in two Bos taurus assemblies.

Publications

Genome assembly has a major impact on gene content: a comparison of annotation in two Bos taurus assemblies.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets