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
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-25-2024, 11:49 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
04-25-2024, 11:49 AM
|
||
Started by seqadmin, 04-24-2024, 08:47 AM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
04-24-2024, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
62 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
Comment