BACKGROUND
next-generation sequencing technologies, such as the  <dig>  and solexa  platforms, provide researchers working on emerging- and non-model species an affordable means of addressing a wide range of questions  <cit> .
while completely assembled genomes of non-model organisms may not be easily obtained, very low-coverage shotgun sequencing can be used for various applications, such as to identify microsatellites for population genetic analyses  <cit> .
low coverage genome survey sequences  can also provide information about gene content, polymorphisms, functional elements and repetitive elements  <cit> .
in support of the utility of low-coverage sequencing, simulations have shown that most of the coding sequence in a genome can be surveyed with less than 2× genomic coverage  <cit> .
as a case in point, the original  <dig> × assembly of the dog genome provided partial sequence of thousands of dog orthologs of human genes  <cit> .
at even lower coverage, wernersson et al.
<cit>  were able to recover 38% of the coding fraction of the mouse-human alignment with only  <dig> × coverage of the pig genome.
furthermore, low-coverage survey sequencing appears to be an efficient way of identifying common repetitive dna sequences.
for example, a large fraction of the repetitive dna sequences in the complex, highly repetitive barley genome of were computationally identified with only ~10% of the genome sequenced  <cit> .
these studies suggest that as next generation sequencing becomes more widespread, low-coverage genome surveys will play a prominent role in studies of non-model species.
to gauge the suitability of genome surveys at depths of coverage lower than those used in previous studies, we partially sequenced the genome of the scuttle fly megaselia scalaris  by  <dig> pyrosequencing.
while no genomic resources had been developed for this species prior to this point, m. scalaris has a rich history in biological research  and can serve as a representative of a potentially emerging model species for ecology and evolutionary biology.
to this end, one-fourth of a standard  <dig> run was conducted with randomly sheared genomic dna.
based on flow cytometry estimates of genome size, the depth of coverage across the genome was approximately  <dig> - <dig> x. however, our low-coverage likely reflects the position of many researchers using next-generation technologies for sequencing in non-model eukaryotes with large, complex genomes.
we therefore view our low-coverage as a bioinformatics challenge and focus on analyses that can be conducted with limited gss data.
even with only ~ <dig> × coverage, we were able to generate a considerable amount of biologically useful information and genomic resources for m. scalaris.
first, because of the substantial impact repetitive elements  can have on the structure and evolution of genomes  <cit> , we searched for res in m. scalaris based on homology to well-characterized res.
novel or lineage-specific tandem res were then identified by using a custom program we developed.
the effectiveness of these computational procedures was evaluated by performing the same analyses on low-coverage sequence simultaneously generated from drosophila pseudoobscura, for which a well-annotated reference genome sequence is available  <cit> .
secondly, we assembled a nearly complete sequence for the ~ <dig>  kb mitochondrial  genome of m. scalaris.
thirdly, we looked for microsatellite loci to develop molecular markers for m. scalaris.
finally, we identified coding regions and other functional elements in the m. scalaris genome by comparisons to the completed genomes of drosophila melanogaster and other dipterans.
a discussion of the utility and concerns raised by extremely low-coverage survey sequencing is also presented.
RESULTS
a total of  <dig>  sequence reads with a mean read length of  <dig> bp were generated from  <dig> sequencing of randomly sheared megaselia scalaris genomic dna.
if we assume the genome size to be approximately  <dig> megabases, as estimated from our flow cytometry, this coverage amounts to  <dig> ×.
if instead the genome size is  <dig> megabases  <cit> , we have  <dig> × coverage of the genome.
to expedite analysis, the full set of sequences was filtered to remove redundant reads that were shorter, internal fragments of longer reads.
the non-redundant set contained  <dig>  reads and is hereinafter referred to as the genome survey sequence  set.
table  <dig> summarizes the proportions of redundant, mitochondrial and identified coding and repetitive  sequence reads contained in the gss set.
the nucleotide content of the gss reads showed that the megaselia genome is g-c poor  compared to a-t content .
removing repetitive and mitochondrial sequences from the gss further increased the nucleotide usage bias to a g-c content of 30%.
*based on significant similarity to known elements in other dipteran genomes
survey of repetitive elements
two approaches were used to identify repetitive elements in the m. scalaris genome.
the first approach relied on homology of repeats to well-described res in publicly available databases.
a total of  <dig> res were identified on the basis of similarity to known elements in the drosophila and anopheles repbase libraries .
based on the number of reads with significant similarity to various classes of res, retroelements appear more abundant in copy number than dna transposons, as is the tendency in most eukaryotes  <cit> .
ltr retroelements were the most abundant class , with many copies of gypsy and copia elements identified.
the second most abundant class was non-ltr elements , representing 38% of identified res, including many jockey elements.
only 4% of identified res were dna transposons.
the second approach used our refinder.plx http://www.biology.duke.edu/noorlab/noorsoft.html program to identify tandem res de novo.
this method has the advantage of detecting previously uncharacterized or highly divergent res, but was not designed to detect interspersed elements.
before running the program on the m. scalaris gss reads, we tested the program on genome survey sequences from d. pseudoobscura bogotana, which were generated from the same  <dig> run and at the same coverage as the m. scalaris data.
this allowed us to optimize parameters in the refinder.plx program and to test the program's ability to detect previously described res in the d. pseudoobscura genome.
the contigs generated by the program matched single units of elements known in d. pseudoobscura such as leviathan  <cit> , a bilbo non-ltr, and its sequences.
furthermore, to confirm that refinder.plx correctly assembles repetitive sequences, contigs that the program generated were queried against the assembled d. pseudoobscura reference genome.
in most cases, our contigs aligned to stretches of the d. pseudoosbcura genome that were as long or just slightly shorter than the contigs generated by the program.
that the alignments spanned the length of multiple reads suggests that these contigs were correctly assembled rather than artificial hybrid sequences.
starting from a set of  <dig>  candidate repeat sequences that were present in at least two copies in the m. scalaris gss reads, refinder.plx identified  <dig> tandem repeats,  <dig> of which were unique repeats that did not show significant similarity to other computationally identified res .
the average coverage across all unique repeat motifs, calculated from coverage across  <dig> bp windows, was  <dig>   with a few repeat motifs reaching as great as 200× coverage.
this high level of redundancy in repeat coverage suggests that most of the identified res are present in high copy number in the m. scalaris genome.
furthermore, for  <dig> of the identified res, no significant hits were found by blast searches of genbank, suggesting that most of the identified res were unique to m. scalaris or previously uncharacterized in other species.
the m. scalaris mitochondrial genome
we initially assembled the m. scalaris mitochondrial  genome by aligning contiguous but overlapping sequence reads while using the d. melanogaster mt genome as a scaffold.
from this initial assembly,  <dig> mitochondrial gss sequences were identified and aligned into several large contigs.
gaps between contigs were amplified and sequenced by traditional sanger sequencing when possible.
there were two gaps for which we could not obtain high-quality sequence.
the first is located between cytochrome c oxidase subunit ii and atpase subunit  <dig> and is flanked by a-t homopolymer tracts .
based on the size of a pcr fragment that spans this gap, it is approximately  <dig> bp long.
the second gap is located between the 16s ribosomal rnas and nadh dehyrdrogenase subunit  <dig>  but we do not have an estimate for its length since it falls within a highly a-t rich region that varies in size among other dipterans and we could not successfully amplify a fragment spanning the gap.
although we did not obtain complete coverage, we did generate a  <dig>  kb consensus sequence for the m. scalaris mt genome .
coverage across the mt genome was high except in extremely a-t rich regions , where no  <dig> sequence traces were recovered.
the mean depth of coverage was greater than 10×, with smaller regions attaining depths greater than 20× .
overall, the m. scalaris mt genome is similar to that of other dipterans.
the a-t content is high, at ~80%.
the arrangement of protein-coding genes and the small and large ribosomal subunit rnas is identical to that of most other arthropods , although atpase subunit  <dig> was not present in our assembly .
the most parsimonious explanation for the absence of atpase subunit  <dig> is that the gene is located within gap  <dig> of our assembly between cytochrome c oxidase ii and atpase subunit  <dig>  where it is normally located in other arthropods  <cit> .
since mtdna sequences are commonly used in phylogenetic and population genetic analyses, we asked if our low-depth gss strategy would typically recover a complete mt genome sequence in other species.
therefore, using the same methods as for m. scalaris, we assembled the  <dig>  kb d. pseudoobscura bogotana mt genome from our gss reads .
for d. pseudoobscura,  <dig> % of the gss reads were mitochondrial as compared to  <dig> % for m. scalaris.
the higher relative abundance of mt reads in d. pseudoobscura allowed for the assembly of the complete mt genome with a mean of 20× coverage.
these results demonstrate that a relatively small number of sequences can recover complete mt sequences in gss studies of some species.
scan for microsatellites
gss reads were scanned for the presence of all possible perfect di- or tri-nucleotide repeat motifs with the program microscan.plx http://www.biology.duke.edu/noorlab/noorsoft.html to find potentially variable microsatellite loci.
only  <dig> microsatellites with  <dig> or more repeat units were identified through this approach in m. scalaris.
this was apparently not due to a failure of the program to detect microsatellites, as it recovered  <dig> such microsatellites when run with the d. pseudoobscura bogotana gss sequences, despite there being 25% fewer sequence reads than in the m. scalaris gss.
the m. scalaris genome therefore appears to harbor a much smaller number of microsatellites.
from these sequence reads, we were able to design primers and identify variability among strains at a subset of these microsatellites .
identification of gene homologs
homology searches were conducted to find m. scalaris reads with significant similarity to annotated genes and functional elements in other dipteran genomes.
evaluating the number of "true" matches between drosophila melanogaster gene sequences and m. scalaris gss reads is challenging because of multigene families within both species, including histones, trnas and rrnas.
hence, a single d. melanogaster gene may resemble many m. scalaris genes, and vice versa.
considering this issue, we examined the number of sequences with similarity between d. melanogaster and m. scalaris in two ways.
as a first approach, we used blast to identify which genes in d. melanogaster bore significant similarity to the nonrepetitive m. scalaris gss set.
with an e-value cutoff of 10- <dig>  we found  <dig> d. melanogaster coding genes with at least one significant match to a m. scalaris read .
however, this approach can still identify multiple similar d. melanogaster genes matching a single m. scalaris sequence trace.
to address this, we then eliminated all redundant hits to single m. scalaris traces, and identified  <dig> unique m. scalaris reads with significant sequence homology to an annotated d. melanogaster gene.
for comparison, we also used the annotated genomes of drosophila virilis and anopheles gambiae to conduct gene homology searches.
with d. virilis,  <dig> potentially homologous genes were identified, although most of the genes are uncharacterized in the d. virilis genome annotation .
as with d. melanogaster, we eliminated all redundant hits to a single m. scalaris sequence trace, which identified  <dig> unique m. scalaris reads with significant homology to a d. virilis gene.
with a. gambiae, only  <dig> potential homologs from  <dig> unique m. scalaris reads were identified , likely due to the combination of the smaller number of annotated genes in the a. gambiae genome and the longer divergence time between m. scalaris and anopheles relative to m. scalaris and drosophila.
in all,  <dig> m. scalaris reads were identified as potential coding sequences by homology searches performed with other dipteran genomes.
discussion
although our preliminary sequencing of the megaselia scalaris genome resulted in extremely low-coverage , we were able to perform a number of bioinformatic analyses that provided useful information for characterizing this genome as well as generating various genomic resources.
we were able to characterize numerous repetitive sequences in the genome, including some with homology to known elements and some that have not been characterized previously.
useful resources such as a nearly complete mitochondrial genome sequence and microsatellite markers were also easily developed from the gss data.
moreover, partial sequences for hundreds of orthologs of drosophila and anopheles genes were generated.
an assumption laden in some of our analyses is that the genome survey sequences studied are "random" segments from across the genome.
we cannot exclude the possibility that certain regions of the genome were more or less likely to be surveyed due to features such as gc-content.
indeed, we observed coverage of the mitochondrial genome was lower or missing in the most extremely a-t rich regions.
this bias may have resulted from the sequencing process itself or issues with sample preparation and/or library generation.
nonetheless, the approaches here provide a first, albeit imperfect, approximation of various features of a previously unexplored genome, and several of our conclusions do not depend upon a truly random sampling of the genome.
m. scalaris was chosen for partial genome sequencing because of its interesting natural history and potential to become a model species in ecology and evolutionary biology.
previous work on m. scalaris has already revealed much about its ecology, development, sex-determination system, and life cycle .
the species is widely distributed and many aspects of its ecology are peculiar.
for example, m. scalaris larvae are notable for the wide range of organic matter on which they can feed; reportedly the widest range of any insect  <cit> .
because larvae are also facultative parasites, they can enter open wounds and therefore pose some threat to human health, especially in the developing world  <cit> .
a complete m. scalaris genome sequence would also strengthen comparative and evolutionary genomic studies of the dipterans.
while there are completed genomes for  <dig> drosophila species  <cit>  and the mosquitoes anopheles gambiae  <cit>  and aedes aegypti  <cit> , no genome sequence is currently available for any dipteran species outside of the drosophilids and mosquitoes.
a phorid fly such as m. scalaris would also serve as a good outgroup in comparative genomic studies of the drosophilids.
for example, the genome of m. scalaris could facilitate the identification of regulatory elements and assessing patterns of evolution, as has been recently suggested also for tephritids  <cit> .
applications of low-coverage genome sequence
we anticipate that researchers studying a wide range of non-model taxa will be drawn to newer, less-expensive genome sequencing technologies, often for generating microsatellites  <cit>  or other markers  <cit>  to survey population variability and connectivity, phylogenetic position, and other questions.
based on our study of m. scalaris, using  <dig> pyrosequencing to sequence genomic dna appears to be an effective strategy for generating low-coverage sequence data, with read-lengths amenable for assembly or blast  <cit>  analyses.
sequence reads also appear to be distributed throughout the genome, allowing for partial coverage of many functional elements and hundreds of orthologs of known genes.
thus, low-depth sequencing provides mostly new sequence and avoids the high redundancy seen in large-scale genome projects.
the ability to find repetitive sequences is another important test of the applicability of survey sequencing since identifying and masking repetitive sequences can be crucial for accurately estimating genome coverage, identifying low-copy "gene space", and assembling large contigs.
we identified over  <dig> m. scalaris transposable element copies by homology searches, most of which were ltr retroelements and non-ltr retrotransposons.
these res could be masked in future genomic work in m, scalaris, facilitating assembly of the short sequence reads obtained through  <dig> or other short-read sequencing.
low-coverage genome surveys therefore appear to be an effective way to identify repetitive sequences, as several previous studies have successfully identified repetitive sequences with low genome coverage in other systems  <cit> .
while available programs like repeatmasker  and others can identify previously known res, identifying novel res in unassembled genomes remains problematic.
our refinder.plx program was designed to quickly identify as many novel res in unassembled genomes as possible.
we further validated this program by applying it to comparable gss from a species with a fully sequenced and assembled genome, drosophila pseudoobscura, and identifying known elements.
however, it was not designed to detect all classes of transposable elements and, because the program works by assembling and identifying potentially repetitive sequences in contigs, it can only identify res in tandem arrays.
it should also be noted that our program was not designed to identify higher-order repeats or identify the exact boundaries of res.
other programs for de novo detection of res, such as reas  <cit>  or rerep  <cit> , may provide better detection of other classes of repeats, such as interspersed elements, in low-coverage genome surveys.
it is also possible that some res we detected are hybrids of different elements or that some non-repetitive flanking ends of res were incorporated.
nonetheless, it provides a useful starting point for characterizing a novel genome of its repetitive element content.
since no attempt to remove mtdna from nuclear dna was made prior to sequencing, mtdna sequences were present in high copy number, which allowed us to assemble most of the m. scalaris mt genome.
even more encouraging was that we were able to assemble a complete mt genome at 20× coverage from the d. pseudoobscura bogotana gss reads.
this suggests that low-coverage genome surveys can also be an easy way of obtaining mtdna sequences for phylogenetic studies and markers for population genetic studies.
the proportion of mitochondrial traces was  <dig> %  for the m. scalaris gss and  <dig> %  for d. p. bogotana, consistent with the estimated greater nuclear genome size of the former .
while it would be helpful to know exactly how much sequence data is needed to completely cover a mt genome, this cannot be easily quantified.
based on a binomial distribution, the expected coverage of a target sequence given a certain depth of coverage or level of redundancy, r, can be approximated by the equation: e =  <dig> - e-r
based on this relationship, for a  <dig> kb mt genome and a mean sequence read length of  <dig> bp, approximately  <dig> reads of mitochondrial sequence are needed to obtain full coverage.
however, this approximation will not hold if sequence reads are nonrandomly distributed over the target sequence.
for instance, a bias towards sequence reads being in g-c rich regions across the m. scalaris mt genome likely explains why we did not obtain the sequence of the a-t rich mitochondrial control region even though we recovered  <dig> mt sequence reads, far more than theory suggests are necessary.
the amount of sequence required for full coverage of a mt genome therefore depends on biases in sequencing and dna preparation, as well as biological differences among organisms  in mt copy number.
the point raised above for mt genome sequencing brings up a more general caveat for researchers using low-coverage gss strategies.
with low depths of coverage, the probability of obtaining complete coverage of any target sequence becomes exceedingly low.
this holds true for coding sequences in the nuclear genome as well as organellar genomes.
if specific sequences are the ultimate goal of genome sequencing, then more directed approaches would be more appropriate than our random gss approach.
CONCLUSIONS
we were able to generate genomic resources for megaselia scalaris with very limited sequence data obtained through 454-pyrosequencing.
although this was a preliminary study, the data we have generated is both immediately useful and will be used to guide future larger-scale sequencing of the m. scalaris genome.
we have also developed scripts for facilitating bioinformatics analysis of gss data and made them available to the public.
our encouraging results suggest that low-coverage gss approaches will become more popular among researchers working on non-model organisms, especially as the cost of next-generation sequencing continues to decline.

 next-generation sequencing technologies, such as the and solexa platforms, provide researchers working on emerging- and non-model species an affordable means of addressing a wide range of questions. while completely assembled genomes of non-model organisms may not be easily obtained, very low-coverage shotgun sequencing can be used for various applications, such as to identify microsatellites for population genetic analyses. low coverage genome survey sequences can also provide information about gene content, polymorphisms, functional elements and repetitive elements. in support of the utility of low-coverage sequencing, simulations have shown that most of the coding sequence in a genome can be surveyed with less than 2× genomic coverage. as a case in point, the original × assembly of the dog genome provided partial sequence of thousands of dog orthologs of human genes. at even lower coverage, wernersson et al. were able to recover 38% of the coding fraction of the mouse-human alignment with only × coverage of the pig genome. furthermore, low-coverage survey sequencing appears to be an efficient way of identifying common repetitive dna sequences. for example, a large fraction of the repetitive dna sequences in the complex, highly repetitive barley genome of were computationally identified with only ~10% of the genome sequenced. these studies suggest that as next generation sequencing becomes more widespread, low-coverage genome surveys will play a prominent role in studies of non-model species. to gauge the suitability of genome surveys at depths of coverage lower than those used in previous studies, we partially sequenced the genome of the scuttle fly megaselia scalaris by pyrosequencing. while no genomic resources had been developed for this species prior to this point, m. scalaris has a rich history in biological research and can serve as a representative of a potentially emerging model species for ecology and evolutionary biology. to this end, one-fourth of a standard run was conducted with randomly sheared genomic dna. based on flow cytometry estimates of genome size, the depth of coverage across the genome was approximately - x. however, our low-coverage likely reflects the position of many researchers using next-generation technologies for sequencing in non-model eukaryotes with large, complex genomes. we therefore view our low-coverage as a bioinformatics challenge and focus on analyses that can be conducted with limited gss data. even with only ~ × coverage, we were able to generate a considerable amount of biologically useful information and genomic resources for m. scalaris. first, because of the substantial impact repetitive elements can have on the structure and evolution of genomes, we searched for res in m. scalaris based on homology to well-characterized res. novel or lineage-specific tandem res were then identified by using a custom program we developed. the effectiveness of these computational procedures was evaluated by performing the same analyses on low-coverage sequence simultaneously generated from drosophila pseudoobscura, for which a well-annotated reference genome sequence is available. secondly, we assembled a nearly complete sequence for the ~ kb mitochondrial genome of m. scalaris. thirdly, we looked for microsatellite loci to develop molecular markers for m. scalaris. finally, we identified coding regions and other functional elements in the m. scalaris genome by comparisons to the completed genomes of drosophila melanogaster and other dipterans. a discussion of the utility and concerns raised by extremely low-coverage survey sequencing is also presented.  a total of sequence reads with a mean read length of bp were generated from sequencing of randomly sheared megaselia scalaris genomic dna. if we assume the genome size to be approximately megabases, as estimated from our flow cytometry, this coverage amounts to ×. if instead the genome size is megabases, we have × coverage of the genome. to expedite analysis, the full set of sequences was filtered to remove redundant reads that were shorter, internal fragments of longer reads. the non-redundant set contained reads and is hereinafter referred to as the genome survey sequence set. table summarizes the proportions of redundant, mitochondrial and identified coding and repetitive sequence reads contained in the gss set. the nucleotide content of the gss reads showed that the megaselia genome is g-c poor compared to a-t content. removing repetitive and mitochondrial sequences from the gss further increased the nucleotide usage bias to a g-c content of 30%. *based on significant similarity to known elements in other dipteran genomes survey of repetitive elements two approaches were used to identify repetitive elements in the m. scalaris genome. the first approach relied on homology of repeats to well-described res in publicly available databases. a total of res were identified on the basis of similarity to known elements in the drosophila and anopheles repbase libraries. based on the number of reads with significant similarity to various classes of res, retroelements appear more abundant in copy number than dna transposons, as is the tendency in most eukaryotes. ltr retroelements were the most abundant class, with many copies of gypsy and copia elements identified. the second most abundant class was non-ltr elements, representing 38% of identified res, including many jockey elements. only 4% of identified res were dna transposons. the second approach used our refinder.plx http://www.biology.duke.edu/noorlab/noorsoft.html program to identify tandem res de novo. this method has the advantage of detecting previously uncharacterized or highly divergent res, but was not designed to detect interspersed elements. before running the program on the m. scalaris gss reads, we tested the program on genome survey sequences from d. pseudoobscura bogotana, which were generated from the same run and at the same coverage as the m. scalaris data. this allowed us to optimize parameters in the refinder.plx program and to test the program's ability to detect previously described res in the d. pseudoobscura genome. the contigs generated by the program matched single units of elements known in d. pseudoobscura such as leviathan, a bilbo non-ltr, and its sequences. furthermore, to confirm that refinder.plx correctly assembles repetitive sequences, contigs that the program generated were queried against the assembled d. pseudoobscura reference genome. in most cases, our contigs aligned to stretches of the d. pseudoosbcura genome that were as long or just slightly shorter than the contigs generated by the program. that the alignments spanned the length of multiple reads suggests that these contigs were correctly assembled rather than artificial hybrid sequences. starting from a set of candidate repeat sequences that were present in at least two copies in the m. scalaris gss reads, refinder.plx identified tandem repeats, of which were unique repeats that did not show significant similarity to other computationally identified res. the average coverage across all unique repeat motifs, calculated from coverage across bp windows, was with a few repeat motifs reaching as great as 200× coverage. this high level of redundancy in repeat coverage suggests that most of the identified res are present in high copy number in the m. scalaris genome. furthermore, for of the identified res, no significant hits were found by blast searches of genbank, suggesting that most of the identified res were unique to m. scalaris or previously uncharacterized in other species. the m. scalaris mitochondrial genome we initially assembled the m. scalaris mitochondrial genome by aligning contiguous but overlapping sequence reads while using the d. melanogaster mt genome as a scaffold. from this initial assembly, mitochondrial gss sequences were identified and aligned into several large contigs. gaps between contigs were amplified and sequenced by traditional sanger sequencing when possible. there were two gaps for which we could not obtain high-quality sequence. the first is located between cytochrome c oxidase subunit ii and atpase subunit and is flanked by a-t homopolymer tracts. based on the size of a pcr fragment that spans this gap, it is approximately bp long. the second gap is located between the 16s ribosomal rnas and nadh dehyrdrogenase subunit but we do not have an estimate for its length since it falls within a highly a-t rich region that varies in size among other dipterans and we could not successfully amplify a fragment spanning the gap. although we did not obtain complete coverage, we did generate a kb consensus sequence for the m. scalaris mt genome. coverage across the mt genome was high except in extremely a-t rich regions, where no sequence traces were recovered. the mean depth of coverage was greater than 10×, with smaller regions attaining depths greater than 20×. overall, the m. scalaris mt genome is similar to that of other dipterans. the a-t content is high, at ~80%. the arrangement of protein-coding genes and the small and large ribosomal subunit rnas is identical to that of most other arthropods, although atpase subunit was not present in our assembly. the most parsimonious explanation for the absence of atpase subunit is that the gene is located within gap of our assembly between cytochrome c oxidase ii and atpase subunit where it is normally located in other arthropods. since mtdna sequences are commonly used in phylogenetic and population genetic analyses, we asked if our low-depth gss strategy would typically recover a complete mt genome sequence in other species. therefore, using the same methods as for m. scalaris, we assembled the kb d. pseudoobscura bogotana mt genome from our gss reads. for d. pseudoobscura, % of the gss reads were mitochondrial as compared to % for m. scalaris. the higher relative abundance of mt reads in d. pseudoobscura allowed for the assembly of the complete mt genome with a mean of 20× coverage. these results demonstrate that a relatively small number of sequences can recover complete mt sequences in gss studies of some species. scan for microsatellites gss reads were scanned for the presence of all possible perfect di- or tri-nucleotide repeat motifs with the program microscan.plx http://www.biology.duke.edu/noorlab/noorsoft.html to find potentially variable microsatellite loci. only microsatellites with or more repeat units were identified through this approach in m. scalaris. this was apparently not due to a failure of the program to detect microsatellites, as it recovered such microsatellites when run with the d. pseudoobscura bogotana gss sequences, despite there being 25% fewer sequence reads than in the m. scalaris gss. the m. scalaris genome therefore appears to harbor a much smaller number of microsatellites. from these sequence reads, we were able to design primers and identify variability among strains at a subset of these microsatellites. identification of gene homologs homology searches were conducted to find m. scalaris reads with significant similarity to annotated genes and functional elements in other dipteran genomes. evaluating the number of "true" matches between drosophila melanogaster gene sequences and m. scalaris gss reads is challenging because of multigene families within both species, including histones, trnas and rrnas. hence, a single d. melanogaster gene may resemble many m. scalaris genes, and vice versa. considering this issue, we examined the number of sequences with similarity between d. melanogaster and m. scalaris in two ways. as a first approach, we used blast to identify which genes in d. melanogaster bore significant similarity to the nonrepetitive m. scalaris gss set. with an e-value cutoff of 10- we found d. melanogaster coding genes with at least one significant match to a m. scalaris read. however, this approach can still identify multiple similar d. melanogaster genes matching a single m. scalaris sequence trace. to address this, we then eliminated all redundant hits to single m. scalaris traces, and identified unique m. scalaris reads with significant sequence homology to an annotated d. melanogaster gene. for comparison, we also used the annotated genomes of drosophila virilis and anopheles gambiae to conduct gene homology searches. with d. virilis, potentially homologous genes were identified, although most of the genes are uncharacterized in the d. virilis genome annotation. as with d. melanogaster, we eliminated all redundant hits to a single m. scalaris sequence trace, which identified unique m. scalaris reads with significant homology to a d. virilis gene. with a. gambiae, only potential homologs from unique m. scalaris reads were identified, likely due to the combination of the smaller number of annotated genes in the a. gambiae genome and the longer divergence time between m. scalaris and anopheles relative to m. scalaris and drosophila. in all, m. scalaris reads were identified as potential coding sequences by homology searches performed with other dipteran genomes. discussion although our preliminary sequencing of the megaselia scalaris genome resulted in extremely low-coverage, we were able to perform a number of bioinformatic analyses that provided useful information for characterizing this genome as well as generating various genomic resources. we were able to characterize numerous repetitive sequences in the genome, including some with homology to known elements and some that have not been characterized previously. useful resources such as a nearly complete mitochondrial genome sequence and microsatellite markers were also easily developed from the gss data. moreover, partial sequences for hundreds of orthologs of drosophila and anopheles genes were generated. an assumption laden in some of our analyses is that the genome survey sequences studied are "random" segments from across the genome. we cannot exclude the possibility that certain regions of the genome were more or less likely to be surveyed due to features such as gc-content. indeed, we observed coverage of the mitochondrial genome was lower or missing in the most extremely a-t rich regions. this bias may have resulted from the sequencing process itself or issues with sample preparation and/or library generation. nonetheless, the approaches here provide a first, albeit imperfect, approximation of various features of a previously unexplored genome, and several of our conclusions do not depend upon a truly random sampling of the genome. m. scalaris was chosen for partial genome sequencing because of its interesting natural history and potential to become a model species in ecology and evolutionary biology. previous work on m. scalaris has already revealed much about its ecology, development, sex-determination system, and life cycle. the species is widely distributed and many aspects of its ecology are peculiar. for example, m. scalaris larvae are notable for the wide range of organic matter on which they can feed; reportedly the widest range of any insect. because larvae are also facultative parasites, they can enter open wounds and therefore pose some threat to human health, especially in the developing world. a complete m. scalaris genome sequence would also strengthen comparative and evolutionary genomic studies of the dipterans. while there are completed genomes for drosophila species and the mosquitoes anopheles gambiae and aedes aegypti, no genome sequence is currently available for any dipteran species outside of the drosophilids and mosquitoes. a phorid fly such as m. scalaris would also serve as a good outgroup in comparative genomic studies of the drosophilids. for example, the genome of m. scalaris could facilitate the identification of regulatory elements and assessing patterns of evolution, as has been recently suggested also for tephritids. applications of low-coverage genome sequence we anticipate that researchers studying a wide range of non-model taxa will be drawn to newer, less-expensive genome sequencing technologies, often for generating microsatellites or other markers to survey population variability and connectivity, phylogenetic position, and other questions. based on our study of m. scalaris, using pyrosequencing to sequence genomic dna appears to be an effective strategy for generating low-coverage sequence data, with read-lengths amenable for assembly or blast analyses. sequence reads also appear to be distributed throughout the genome, allowing for partial coverage of many functional elements and hundreds of orthologs of known genes. thus, low-depth sequencing provides mostly new sequence and avoids the high redundancy seen in large-scale genome projects. the ability to find repetitive sequences is another important test of the applicability of survey sequencing since identifying and masking repetitive sequences can be crucial for accurately estimating genome coverage, identifying low-copy "gene space", and assembling large contigs. we identified over m. scalaris transposable element copies by homology searches, most of which were ltr retroelements and non-ltr retrotransposons. these res could be masked in future genomic work in m, scalaris, facilitating assembly of the short sequence reads obtained through or other short-read sequencing. low-coverage genome surveys therefore appear to be an effective way to identify repetitive sequences, as several previous studies have successfully identified repetitive sequences with low genome coverage in other systems. while available programs like repeatmasker and others can identify previously known res, identifying novel res in unassembled genomes remains problematic. our refinder.plx program was designed to quickly identify as many novel res in unassembled genomes as possible. we further validated this program by applying it to comparable gss from a species with a fully sequenced and assembled genome, drosophila pseudoobscura, and identifying known elements. however, it was not designed to detect all classes of transposable elements and, because the program works by assembling and identifying potentially repetitive sequences in contigs, it can only identify res in tandem arrays. it should also be noted that our program was not designed to identify higher-order repeats or identify the exact boundaries of res. other programs for de novo detection of res, such as reas or rerep, may provide better detection of other classes of repeats, such as interspersed elements, in low-coverage genome surveys. it is also possible that some res we detected are hybrids of different elements or that some non-repetitive flanking ends of res were incorporated. nonetheless, it provides a useful starting point for characterizing a novel genome of its repetitive element content. since no attempt to remove mtdna from nuclear dna was made prior to sequencing, mtdna sequences were present in high copy number, which allowed us to assemble most of the m. scalaris mt genome. even more encouraging was that we were able to assemble a complete mt genome at 20× coverage from the d. pseudoobscura bogotana gss reads. this suggests that low-coverage genome surveys can also be an easy way of obtaining mtdna sequences for phylogenetic studies and markers for population genetic studies. the proportion of mitochondrial traces was % for the m. scalaris gss and % for d. p. bogotana, consistent with the estimated greater nuclear genome size of the former. while it would be helpful to know exactly how much sequence data is needed to completely cover a mt genome, this cannot be easily quantified. based on a binomial distribution, the expected coverage of a target sequence given a certain depth of coverage or level of redundancy, r, can be approximated by the equation: e = - e-r based on this relationship, for a kb mt genome and a mean sequence read length of bp, approximately reads of mitochondrial sequence are needed to obtain full coverage. however, this approximation will not hold if sequence reads are nonrandomly distributed over the target sequence. for instance, a bias towards sequence reads being in g-c rich regions across the m. scalaris mt genome likely explains why we did not obtain the sequence of the a-t rich mitochondrial control region even though we recovered mt sequence reads, far more than theory suggests are necessary. the amount of sequence required for full coverage of a mt genome therefore depends on biases in sequencing and dna preparation, as well as biological differences among organisms in mt copy number. the point raised above for mt genome sequencing brings up a more general caveat for researchers using low-coverage gss strategies. with low depths of coverage, the probability of obtaining complete coverage of any target sequence becomes exceedingly low. this holds true for coding sequences in the nuclear genome as well as organellar genomes. if specific sequences are the ultimate goal of genome sequencing, then more directed approaches would be more appropriate than our random gss approach.  we were able to generate genomic resources for megaselia scalaris with very limited sequence data obtained through 454-pyrosequencing. although this was a preliminary study, the data we have generated is both immediately useful and will be used to guide future larger-scale sequencing of the m. scalaris genome. we have also developed scripts for facilitating bioinformatics analysis of gss data and made them available to the public. our encouraging results suggest that low-coverage gss approaches will become more popular among researchers working on non-model organisms, especially as the cost of next-generation sequencing continues to decline. 
random pyrosequencing of the m. scalaris genome provided a depth of coverage much lower than typical gss studies. we also use our gss data to develop genomic resources such as a nearly complete mitochondrial genome sequence and microsatellite markers for m. scalaris. we conclude that low-coverage genome surveys are effective at generating useful information about organisms currently lacking genomic sequence data. we demonstrate that, even with extremely low-coverage sequencing, bioinformatics approaches can yield extensive information about functional and repetitive elements. while obtaining high-depth, completely assembled genome sequences for most non-model organisms remains challenging, low-coverage genome survey sequences can provide a wealth of biologically useful information at low cost. here, using a random pyrosequencing approach, we sequence the genome of the scuttle fly megaselia scalaris and evaluate the utility of our low-coverage gss approach.
 the declining cost of dna sequencing is making genome sequencing a feasible option for more organisms, including many of interest to ecologists and evolutionary biologists. while obtaining high-depth, completely assembled genome sequences for most non-model organisms remains challenging, low-coverage genome survey sequences can provide a wealth of biologically useful information at low cost. here, using a random pyrosequencing approach, we sequence the genome of the scuttle fly megaselia scalaris and evaluate the utility of our low-coverage gss approach.  random pyrosequencing of the m. scalaris genome provided a depth of coverage much lower than typical gss studies. we demonstrate that, even with extremely low-coverage sequencing, bioinformatics approaches can yield extensive information about functional and repetitive elements. we also use our gss data to develop genomic resources such as a nearly complete mitochondrial genome sequence and microsatellite markers for m. scalaris.  we conclude that low-coverage genome surveys are effective at generating useful information about organisms currently lacking genomic sequence data. 