Phylogenetic analysis of nitrite, nitric oxide, and nitrous oxide respiratory enzymes reveal a complex evolutionary history for denitrification

Denitrification is a facultative respiratory pathway in which nitrite (NO2(-)), nitric oxide (NO), and nitrous oxide (N2O) are successively reduced to nitrogen gas (N(2)), effectively closing the nitrogen cycle. The ability to denitrify is widely dispersed among prokaryotes, and this polyphyletic distribution has raised the possibility of horizontal gene transfer (HGT) having a substantial role in the evolution of denitrification. Comparisons of 16S rRNA and denitrification gene phylogenies in recent studies support this possibility; however, these results remain speculative as they are based on visual comparisons of phylogenies from partial sequences. We reanalyzed publicly available nirS, nirK, norB, and nosZ partial sequences using Bayesian and maximum likelihood phylogenetic inference. Concomitant analysis of denitrification genes with 16S rRNA sequences from the same organisms showed substantial differences between the trees, which were supported by examining the posterior probability of monophyletic constraints at different taxonomic levels. Although these differences suggest HGT of denitrification genes, the presence of structural variants for nirK, norB, and nosZ makes it difficult to determine HGT from other evolutionary events. Additional analysis using phylogenetic networks and likelihood ratio tests of phylogenies based on full-length sequences retrieved from genomes also revealed significant differences in tree topologies among denitrification and 16S rRNA gene phylogenies, with the exception of the nosZ gene phylogeny within the data set of the nirK-harboring genomes. However, inspection of codon usage and G + C content plots from complete genomes gave no evidence for recent HGT. Instead, the close proximity of denitrification gene copies in the genomes of several denitrifying bacteria suggests duplication. Although HGT cannot be ruled out as a factor in the evolution of denitrification genes, our analysis suggests that other phenomena, such gene duplication/divergence and lineage sorting, may have differently influenced the evolution of each denitrification gene.     

Sequence variation in ITS spacers and 5.8S rDNA and relationship of E, St, P, Ns, Xm, and H genomes in the genera of Agropyron, Elytrigia, Leymus, Pascopyrum, Psathyrostachys, and Hordeum

Understanding species evolution and improvement requires information of their genome origin and differentiation. Among the species in the family Gramineae, genome identities of Agropyron-Elytrigia-Leymus group are still ambiguous. In order to delineate the genome relationship, nucleotide sequence analysis in the rDNA ITS regions was carried out among the species in the genera Elytrigia, Agropyron, Psathyrostachys, Leymus, and Psacopyrum containing E, St, P, Ns, and Xm genomes. The ITS-1 and ITS-2 showed a narrow range of variation in length except for the presence of a pentanucleotide, TGGGG, in/del in some haplotypes, whereas higher numbers of nucleotide substitutions were observed in most genera. There were 187 variable sites in the ITS-1, 5.8S, and ITS-2 regions, in which a few genome specific mutations were observed. While the level of variation was similar between ITS-1 and ITS-2, the rate of transition mutation versus transversion mutations was different among the ITS-1, 5.8S, and ITS-2 segments. GC contents of the ITS regions ranged between 55–65% between genomes and the haplotypes of P and H genomes were slightly higher than others. In phylogenetic analysis, the ITS haplotypes were classified into two groups; one containing H, Ns, NsXm genomes, and another containing P, St, and E genomes, which are congruous to the genome affinities from other studies. Among the four genomes in Pascopyrum smithii (2n=8x=56, StStNsNsHHXmXm), the haplotypes of H and St genomes were identified with the reference diploid species, but the haplotypes having Ns and Xm genomes were not found in the present analysis.     

Arachnid relationships based on mitochondrial genomes: asymmetric nucleotide and amino acid bias affects phylogenetic analyses

Phylogenetic analyses based on mitochondrial DNA have yielded widely differing relationships among members of the arthropod lineage Arachnida, depending on the nucleotide coding schemes and models of evolution used. We enhanced taxonomic coverage within the Arachnida greatly by sequencing seven new arachnid mitochondrial genomes from five orders. We then used all 13 mitochondrial protein-coding genes from these genomes to evaluate patterns of nucleotide and amino acid biases. Our data show that two of the six orders of arachnids (spiders and scorpions) have experienced shifts in both nucleotide and amino acid usage in all their protein-coding genes, and that these biases mislead phylogeny reconstruction. These biases are most striking for the hydrophobic amino acids isoleucine and valine, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. To improve phylogenetic accuracy based on amino acid differences, we tested two recoding methods: (1) removing all isoleucine and valine sites and (2) recoding amino acids based on their physiochemical properties. We find that these methods yield phylogenetic trees that are consistent in their support of ancient intraordinal divergences within the major arachnid lineages. Further refinement of amino acid recoding methods may help us better delineate interordinal relationships among these diverse organisms.     

REVIEW—CONTRASTING MITOCHONDRIAL GENOME ORGANIZATIONS AND SEQUENCE AFFILIATIONS AMONG GREEN ALGAE: POTENTIAL FACTORS, MECHANISMS, AND EVOLUTIONARY SCENARIOS

The three green algal mitochondrial genomes completely sequenced to date — those of Chlamydomonas reinhardtii Dangeard, Chlamydomonas eugametos Gerloff, and Prototheca wickerhamii Soneda & Tubaki — revealed very different mitochondrial genome organizations and sequence affiliations. The Chlamydomonas genomes resemble the ciliate / fungal / animal counterparts, and the Prototheca genome resembles land plant homologues. This review points out that all the green algal mitochondrial genomes examined to date resemble either the Chlamydomonas or the Prototheca mitochondrial genome; the Chlamydomonas- like mitochondrial genomes are small and have a reduced gene content (no ribosomal protein or 5S rRNA genes and only a few protein-coding and tRNA genes) and fragmented and scrambled rRNA coding regions, whereas the Prototheca- like mitochondrial genomes are larger and have a larger set of protein-coding genes (including ribosomal protein genes), more tRNA genes, and 5S rRNA and conventional continuous small-subunit (SSU) and large-subunit (LSU) rRNA coding regions. It appears, therefore, that the differences previously observed between the mitochondrial genomes of C. reinhardtii and P. wickerhamii extend to the two green algal mitochondrial lineages to which they belong and are significant enough to raise questions about the causes and mechanisms responsible for such contrasting evolutionary strategies among green algae. This review suggests an integrative approach in explaining the occurrence of distinct evolutionary strategies and apparent phylogenetic affiliations among the known green algal mitochondrial lineages. The observed differences could be the result of distinct genetic potentials differentiated during the previous evolutionary history of the flagellate ancestors and / or of subsequent changes in habitat and life history of the more advanced green algal lineages.     

Distribution of genes and recombination in wheat and other eukaryotes

The genome sizes of eukaryotes may differ as much as 10,400-fold. A part of these differences may be attributed to polyploidy, and increase in gene number and size. Most of the genome size disparity is due to non-transcribed repeated DNA including retrotransposons and pseudogenes. Only a small fraction of the larger genomes such as those of many crop plants, contain genes. Genes are distributed unevenly along the chromosomes, often organized in clusters of varying sizes and gene-densities (gene-rich regions). The regions corresponding to gene-clusters in smaller genome plants such as rice may be divided into many ‘mini’ gene-clusters in the related larger genomes. The range of gene-density within the ‘mini2019; gene-clusters is about the same among plants with varying genome sizes. Recombination per chromosome is similar among eukaryotes, and thus is considerably independent of DNA content and chromosome size. Relatively little recombination occurs outside the gene-rich regions. Recombination varies dramatically among various gene regions, and is highly uneven within gene regions as well. Consequently, a significant number of genes may be inaccessible to recombination-based manipulations such as map-based cloning.     

Multidisciplinary approach to genome analysis in the diploid species,Thinopyrum bessarabicum and Th. elongatum (Lophopyrum elongatum), of the Triticeae

Summary The J and E genome species of the Triticeae are invaluable sources of salt tolerance. The evidence concerning the phyletic relatedness of the J genome of diploid Thinopyrum bessarabicum and the E genome of diploid Th. elongatum (=Lophopyrum elongatum) is discussed. Low level of chromosome pairing between J and E at different ploidy levels, suppression of J-E pairing by the Ph1 pairing regulator that inhibits homoeologous pairing, complete sterility of the diploid hybrids (JE), karyotypic divergence of the two genomes, differences in total content and distribution of heterochromatin along their chromosomes, and marked differences in gliadin proteins, isozymes, 5S DNA, and rDNA indicate that J and E are distinct genomes. Well-defined biochemical markers have been identified in the two genomes and may be useful in plant breeding. The level of distinction between J and E is comparable to that among the universally accepted homoeologous genomes A, B, and D of wheat. Therefore, the J and E genomes are homoeologous and not homologous, although some workers continue to call them homologous. The previous workers' data on chromosome pairing in diploid hybrids and/ or karyotypic differences in the conventionally stained chromosomes do not provide sufficient evidence for the proposed merger of J and E genomes (and, hence, of the genera Thinopyrum and Lophopyrum) specifically and for establishing genome relationships generally. Extra precautions should be exercised before changing the designation of an established genome and before merging two genera. A uniform, standardized system of genomic nomenclature for the entire Triticeae is proposed, which should benefit cytogeneticists, plant breeders, taxonomists, and evolutionists.Cooperative investigations of the USDA-Agricultural Research Service and the Utah Agricultural Experiment Station, Logan, UT 84322, USA. Approved as Journal Paper no. 3832     

Using single-nucleotide polymorphisms to discriminate disease-associated from carried genomes of Neisseria meningitidis

Neisseria meningitidis is one of the main agents of bacterial meningitis, causing substantial morbidity and mortality worldwide. However, most of the time N. meningitidis is carried as a commensal not associated with invasive disease. The genomic basis of the difference between disease-associated and carried isolates of N. meningitidis may provide critical insight into mechanisms of virulence, yet it has remained elusive. Here, we have taken a comparative genomics approach to interrogate the difference between disease-associated and carried isolates of N. meningitidis at the level of individual nucleotide variations (i.e., single nucleotide polymorphisms [SNPs]). We aligned complete genome sequences of 8 disease-associated and 4 carried isolates of N. meningitidis to search for SNPs that show mutually exclusive patterns of variation between the two groups. We found 63 SNPs that distinguish the 8 disease-associated genomes from the 4 carried genomes of N. meningitidis, which is far more than can be expected by chance alone given the level of nucleotide variation among the genomes. The putative list of SNPs that discriminate between disease-associated and carriage genomes may be expected to change with increased sampling or changes in the identities of the isolates being compared. Nevertheless, we show that these discriminating SNPs are more likely to reflect phenotypic differences than shared evolutionary history. Discriminating SNPs were mapped to genes, and the functions of the genes were evaluated for possible connections to virulence mechanisms. A number of overrepresented functional categories related to virulence were uncovered among SNP-associated genes, including genes related to the category "symbiosis, encompassing mutualism through parasitism."     

Cucumber: a model angiosperm for mitochondrial transformation?

Plants possess three major genomes, carried in the chloroplast, mitochondrion, and nucleus. The chloroplast genomes of higher plants tend to be of similar sizes and structure. In contrast both the nuclear and mitochondrial genomes show great size differences, even among closely related species. The largest plant mitochondrial genomes exist in the genus Cucumis at 1500 to 2300 kilobases, over 100 times the sizes of the yeast or human mitochondrial genomes. Biochemical and molecular analyses have established that the huge Cucumis mitochondrial genomes are due to extensive duplication of short repetitive DNA motifs. The organellar genomes of almost all organisms are maternally transmitted and few methods exist to manipulate these important genomes. Although chloroplast transformation has been achieved, no routine method exists to transform the mitochondrial genome of higher plants. A mitochondrial-transformation system for a higher plant would allow geneticists to use reverse genetics to study mitochondrial gene expression and to establish the efficacy of engineered mitochondrial genes for the genetic improvement of the mitochondrial genome. Cucumber possesses three unique attributes that make it a potential model system for mitochondrial transformation of a higher plant. Firstly, its mitochondria show paternal transmission. Secondly, microspores possess relatively few, huge mitochondria. Finally, there exists in cucumber unique mitochondrial mutations conditioning strongly mosaic (msc) phenotypes. The msc phenotypes appear after regeneration of plants from cell culture and sort with specific rearranged and deleted regions in the mitochondrial genome. These mitochondrial deletions may be a useful genetic tool to develop selectable markers for mitochondrial transformation of higher plants.     

Detecting outlier peptides in quantitative high-throughput mass spectrometry data

Quantitative high-throughput mass spectrometry has become an established tool to measure relative gene expression proteome-wide. The output of such an experiment usually consists of a list of expression ratios (fold changes) for several thousand proteins between two conditions. However, we observed that individual peptide fold changes may show a significantly different behavior than other peptides from the same protein and that these differences cannot be explained by imprecise measurements. Such outlier peptides can be the consequence of several technical (misidentifications, misquantifications) or biological (post-translational modifications, differential regulation of isoforms) reasons. We developed a method to detect outlier peptides in mass spectrometry data which is able to delineate imprecise measurements from real outlier peptides with high accuracy when the true difference is as small as 1.4 fold. We applied our method to experimental data and investigated the different technical and biological effects that result in outlier peptides. Our method will assist future research to reduce technical bias and can help to identify genes with differentially regulated protein isoforms in high throughput mass spectrometry data.     

Inheritance and recombination of mitochondrial genomes in plants, fungi and animals

It is generally assumed that mitochondrial genomes are uniparentally transmitted, homoplasmic and nonrecombining. However, these assumptions draw largely from early studies on animal mitochondrial DNA (mtDNA). In this review, we show that plants, animals and fungi are all characterized by episodes of biparental inheritance, recombination among genetically distinct partners, and selfish elements within the mitochondrial genome, but that the extent of these phenomena may vary substantially across taxa. We argue that occasional biparental mitochondrial transmission may allow organisms to achieve the best of both worlds by facilitating mutational clearance but continuing to restrict the spread of selfish genetic elements. We also show that methodological biases and disproportionately allocated study effort are likely to have influenced current estimates of the extent of biparental inheritance, heteroplasmy and recombination in mitochondrial genomes from different taxa. Despite these complications, there do seem to be discernible similarities and differences in transmission dynamics and likelihood of recombination of mtDNA in plant, animal and fungal taxa that should provide an excellent opportunity for comparative investigation of the evolution of mitochondrial genome dynamics.     

Genetic relatedness of the genomes of equine herpesvirus types 1, 2, and 3.

Genomic DNAs of equine herpesvirus type 1 (EHV-1), EHV-2 (equine cytomegalovirus), and EHV-3 were examined by reassociation kinetic and thermal denaturation analyses to determine the extent and degree of homology among the three viral DNAs. Results of reassociation analyses indicated a limited homology among the three EHV genomes. Homologous DNA sequences equivalent to 1.8 to 3.7 megadaltons between EHV-1 and equine cytomegalovirus, 7.6 to 8.2 megadaltons between EHV-1 and EHV-3, and 1.3 to 1.9 megadaltons between equine cytomegalovirus and EHV-3 were detected. Examination by thermal denaturation of the DNA homoduplexes and heteroduplexes formed during reassociation revealed a high degree of base pairing within the duplexes, suggesting that closely related sequences may be conserved among the genomes of EHV.     

DNA replication and strand asymmetry in prokaryotic and mitochondrial genomes

Different patterns of strand asymmetry have been documented in a variety of prokaryotic genomes as well as mitochondrial genomes. Because different replication mechanisms often lead to different patterns of strand asymmetry, much can be learned of replication mechanisms by examining strand asymmetry. Here I summarize the diverse patterns of strand asymmetry among different taxonomic groups to suggest that (1) the single-origin replication may not be universal among bacterial species as the endosymbionts Wigglesworthia glossinidia, Wolbachia species, cyanobacterium Synechocystis 6803 and Mycoplasma pulmonis genomes all exhibit strand asymmetry patterns consistent with the multiple origins of replication, (2) different replication origins in some archaeal genomes leave quite different patterns of strand asymmetry, suggesting that different replication origins in the same genome may be differentially used, (3) mitochondrial genomes from representative vertebrate species share one strand asymmetry pattern consistent with the strand-displacement replication documented in mammalian mtDNA, suggesting that the mtDNA replication mechanism in mammals may be shared among all vertebrate species, and (4) mitochondrial genomes from primitive forms of metazoans such as the sponge and hydra (representing Porifera and Cnidaria, respectively), as well as those from plants, have strand asymmetry patterns similar to single-origin or multi-origin replications observed in prokaryotes and are drastically different from mitochondrial genomes from other metazoans. This may explain why sponge and hydra mitochondrial genomes, as well as plant mitochondrial genomes, evolves much slower than those from other metazoans.     

Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates

Genome size and complexity vary tremendously among eukaryotic species and their organelles. Comparisons across deeply divergent eukaryotic lineages have suggested that variation in mutation rates may explain this diversity, with increased mutational burdens favoring reduced genome size and complexity. The discovery that mitochondrial mutation rates can differ by orders of magnitude among closely related angiosperm species presents a unique opportunity to test this hypothesis. We sequenced the mitochondrial genomes from two species in the angiosperm genus Silene with recent and dramatic accelerations in their mitochondrial mutation rates. Contrary to theoretical predictions, these genomes have experienced a massive proliferation of noncoding content. At 6.7 and 11.3 Mb, they are by far the largest known mitochondrial genomes, larger than most bacterial genomes and even some nuclear genomes. In contrast, two slowly evolving Silene mitochondrial genomes are smaller than average for angiosperms. Consequently, this genus captures approximately 98% of known variation in organelle genome size. The expanded genomes reveal several architectural changes, including the evolution of complex multichromosomal structures (with 59 and 128 circular-mapping chromosomes, ranging in size from 44 to 192 kb). They also exhibit a substantial reduction in recombination and gene conversion activity as measured by the relative frequency of alternative genome conformations and the level of sequence divergence between repeat copies. The evolution of mutation rate, genome size, and chromosome structure can therefore be extremely rapid and interrelated in ways not predicted by current evolutionary theories. Our results raise the hypothesis that changes in recombinational processes, including gene conversion, may be a central force driving the evolution of both mutation rate and genome structure.     

Reconciling nuclear microsatellite and mitochondrial marker estimates of population structure: breeding population structure of Chesapeake Bay striped bass (Morone saxatilis)

Comparative analyses of nuclear and organelle genetic markers may help delineate evolutionarily significant units or management units, although population differentiation estimates from multiple genomes can also conflict. Striped bass (Morone saxatilis) are long-lived, highly migratory anadromous fish recently recovered from a severe decline in population size. Previous studies with protein, nuclear DNA and mitochondrial DNA (mtDNA) markers produced discordant results, and it remains uncertain if the multiple tributaries within Chesapeake Bay constitute distinct management units. Here, 196 young-of-the-year (YOY) striped bass were sampled from Maryland's Choptank, Potomac and Nanticoke Rivers and the north end of Chesapeake Bay in 1999 and from Virginia's Mataponi and Rappahannock Rivers in 2001. A total of 10 microsatellite loci exhibited between two and 27 alleles per locus with observed heterozygosities between 0.255 and 0.893. The 10-locus estimate of R(ST) among the six tributaries was -0.0065 (95% confidence interval -0.0198 to 0.0018). All R(ST) and all but one theta estimates for pairs of populations were not significantly different from zero. Reanalysis of Chesapeake Bay striped bass mtDNA data from two previous studies estimated population differentiation between theta=-0.002 and 0.160, values generally similar to mtDNA population differentiation predicted from microsatellite R(ST) after adjusting for reduced effective population size and uniparental inheritance in organelle genomes. Based on mtDNA differentiation, breeding sex ratios or gene flow may have been slightly male biased in some years. The results reconcile conflicting past studies based on different types of genetic markers, supporting a single Chesapeake Bay management unit encompassing a panmictic striped bass breeding population.     

Assignment of orthologous genes via genome rearrangement

The assignment of orthologous genes between a pair of genomes is a fundamental and challenging problem in comparative genomics. Existing methods that assign orthologs based on the similarity between DNA or protein sequences may make erroneous assignments when sequence similarity does not clearly delineate the evolutionary relationship among genes of the same families. In this paper, we present a new approach to ortholog assignment that takes into account both sequence similarity and evolutionary events at a genome level, where orthologous genes are assumed to correspond to each other in the most parsimonious evolving scenario under genome rearrangement. First, the problem is formulated as that of computing the signed reversal distance with duplicates between the two genomes of interest. Then, the problem is decomposed into two new optimization problems, called minimum common partition and maximum cycle decomposition, for which efficient heuristic algorithms are given. Following this approach, we have implemented a high-throughput system for assigning orthologs on a genome scale, called SOAR, and tested it on both simulated data and real genome sequence data. Compared to a recent ortholog assignment method based entirely on homology search (called INPARANOID), SOAR shows a marginally better performance in terms of sensitivity on the real data set because it is able to identify several correct orthologous pairs that are missed by INPARANOID. The simulation results demonstrate that SOAR, in general, performs better than the iterated exemplar algorithm in terms of computing the reversal distance and assigning correct orthologs.     

Plant proteome analysis by mass spectrometry: principles, problems, pitfalls and recent developments

The genome of several species has now been elucidated; these genomes indicate the proteomic potential of the cell. While identification of genomes has been, and continues to be, a technically and intellectually demanding process, the identification of the proteome contains inherently greater difficulties. The proteome of each living cell is dynamic, altering in response to the individual cell's metabolic state and reception of intracellular and extracellular signal molecules, and many of the proteins which are expressed will be post-translationally altered. Thus if the purpose of the proteome analysis is to aid the understanding of protein function and interaction, then it is identification of the proteins in their final state which is required: for this mass spectrometric identification of individual proteins, indicating site and nature of modifications, is essential. Here we review the principles of the methodologies involved in such analyses, give some indication of current achievements in plant proteomics, and indicate imminent and prospective technical developments.     

Recruitment of single-stranded recombinant adeno-associated virus vector genomes and intermolecular recombination are responsible for stable transduction of liver in vivo

Recombinant adeno-associated virus (rAAV) vectors stably transduce hepatocytes in experimental animals. Following portal-vein administration of rAAV vectors in vivo, single-stranded (ss) rAAV genomes become double stranded (ds), circularized, and/or concatemerized concomitant with a slow rise and, eventually, steady-state levels of transgene expression. Over time, at least some of the stabilized genomes become integrated into mouse chromosomal DNA. The mechanism(s) of formation of stable ds rAAV genomes from input ss DNA molecules has not been delineated, although second-strand synthesis and genome amplification by a rolling-circle model has been proposed. To begin to delineate a mechanism, we produced rAAV vectors in the presence of bacterial PaeR7 or Dam methyltransferase or constructed rAAV vectors labeled with different restriction enzyme recognition sites and introduced them into mouse hepatocytes in vivo. A series of molecular analyses demonstrated that second-strand synthesis and rolling-circle replication did not appear to be the major processes involved in the formation of stable ds rAAV genomes. Rather, recruitment of complementary plus and minus ss genomes and subsequent random head-to-head, head-to-tail, and tail-to-tail intermolecular joining were primarily responsible for the formation of ds vector genomes. These findings contrast with the previously described mechanism(s) of transduction based on in vitro studies. Understanding the mechanistic process responsible for vector transduction may allow the development of new strategies for improving rAAV-mediated gene transfer in vivo.     

Comparative genomics of small RNAs in bacterial genomes

In recent years, various families of small non-coding RNAs (sRNAs) have been discovered by experimental and computational approaches, both in bacterial and eukaryotic genomes. Although most of them await elucidation of their function, it has been reported that some play important roles in gene regulation. Here we carried out comparative genomics analysis of possible sRNAs that are computationally identified in 30 bacterial genomes from gamma- and alpha-proteobacteria and Deinococcus radiodurans. Identified sRNAs are clustered by a complete-linkage clustering method to see conservation among the organisms. On average, sRNAs are found in approximately 30% of intergenic regions of each genome sequence. Of these, 25.7% are conserved among three or more organisms. Approximately 60% of the conserved sRNAs do not locate in orthologous intergenic regions, implying that sRNAs may be shuffled their positions in genomes. The current study implies that sRNAs may be involved in a more extensive range of functions in bacteria.     

Shotgun metagenomics of soil invertebrate communities reflects taxonomy,biomass, and reference genome properties

• Metagenomics – shotgun sequencing of all DNA fragments from a community DNA extract – is routinely used to describe the composition, structure, and function of microorganism communities. Advances in DNA sequencing and the availability of genome databases increasingly allow the use of shotgun metagenomics on eukaryotic communities. Metagenomics offers major advances in the recovery of biomass relationships in a sample, in comparison to taxonomic marker gene‐based approaches (metabarcoding). However, little is known about the factors which influence metagenomics data from eukaryotic communities, such as differences among organism groups, the properties of reference genomes, and genome assemblies.
• We evaluated how shotgun metagenomics records composition and biomass in artificial soil invertebrate communities at different sequencing efforts. We generated mock communities of controlled biomass ratios from 28 species from all major soil mesofauna groups: mites, springtails, nematodes, tardigrades, and potworms. We shotgun sequenced these communities and taxonomically assigned them with a database of over 270 soil invertebrate genomes.
• We recovered over 95% of the species, and observed relatively high false‐positive detection rates. We found strong differences in reads assigned to different taxa, with some groups (e.g., springtails) consistently attracting more hits than others (e.g., enchytraeids). Original biomass could be predicted from read counts after considering these taxon‐specific differences. Species with larger genomes, and with more complete assemblies, consistently attracted more reads than species with smaller genomes. The GC content of the genome assemblies had no effect on the biomass–read relationships. Results were similar among different sequencing efforts.
• The results show considerable differences in taxon recovery and taxon specificity of biomass recovery from metagenomic sequence data. The properties of reference genomes and genome assemblies also influence biomass recovery, and they should be considered in metagenomic studies of eukaryotes. We show that low‐ and high‐sequencing efforts yield similar results, suggesting high cost‐efficiency of metagenomics for eukaryotic communities. We provide a brief roadmap for investigating factors which influence metagenomics‐based eukaryotic community reconstructions. Understanding these factors is timely as accessibility of DNA sequencing and momentum for reference genomes projects show a future where the taxonomic assignment of DNA from any community sample becomes a reality.