I have a question about assembling mixed samples, where the organisms sampled are very similar. We usually sequence pure isolates but in a particular study ended up Illumina sequencing (paired ends) a sample that included a number of E. coli strains. These are all commensal E. coli but we know from separate work they have different phenotype and PFGE results. If each bacteria in the sample was a different organisms (e.g. E. coli, Salmonella, Acinetobacter) I know we could attempt to assemble the data and pull out separate genomes. But in this case, much of the genomes are identical or very similar. Yet other regions will include gene additions / deletions or genome re-arrangments. I was thinking that this would be analogous to assembling a plant or animal polyploid genome, where parts of the chromosome could be assembled to a haploid consensus but others would have multiple solutions, reflective of difference in chromosome copies. Can anyone suggest an assembly approach for our data or perhaps a polyploid assembler that may applied? Thanks!
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Genome assembly of mixed E. coli samples
-----
Andrew G. McArthur, Ph.D.
Associate Professor & Cisco Research Chair in Bioinformatics
Department of Biochemistry & Biomedical Sciences
McMaster University, Hamilton, Ontario, Canada
w. http://mcarthurbioinformatics.ca | e. [email protected] -
I don't know if the polyploid assemblers will work with your sample. Ignoring repetitive elements (which are difficult enough in haploids), polyploid genomes have (approximately) integer multiples of the haploid genome for variant representation. The data corresponding to one haploid genome, as well as the ploidy, can be estimated from kmer frequencies. For example, a tetraploid genome will typically have the tallest kmer peak at 4X the frequency of the first observable peak (the former represents the non-variant sequences, while the latter represents variants unique to one haploid genome), with peaks at 2X and 3X as well. Your sample has an unknown number of strains with non-integer representations, and you expect the majority of their genomes to be identical. Your predicted distribution of kmer frequencies would be one large peak corresponding to your depth of coverage, with no obvious peaks at lower coverage.
Perhaps the best approach would be to align the data to an E. coli reference with stringent parameters (e.g., minimal mismatches and gaps) to remove the bulk of the data, then try to assemble the unaligned reads with a metagenome assembler. Note that I have not attempted such an approach, so it may not work.
-
Comment
Latest Articles
Collapse
-
by seqadmin
Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.
Somatic Genomics
“We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...-
Channel: Articles
05-24-2024, 01:16 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 06-03-2024, 06:55 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
06-03-2024, 06:55 AM
|
||
Started by seqadmin, 05-30-2024, 03:16 PM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
05-30-2024, 03:16 PM
|
||
Comprehensive Sequencing of Great Ape Sex Chromosomes Yields Insights into Evolution and Genetic Variability
by seqadmin
Started by seqadmin, 05-29-2024, 01:32 PM
|
0 responses
29 views
0 likes
|
Last Post
by seqadmin
05-29-2024, 01:32 PM
|
||
Started by seqadmin, 05-24-2024, 07:15 AM
|
0 responses
216 views
0 likes
|
Last Post
by seqadmin
05-24-2024, 07:15 AM
|
Comment