SCiMS: Sex Calling in Metagenomic Sequences

BackgroundHost sex is a critical determinant of microbial community structure across many host species, influenced by hormonal profiles, physiology, and sex-stratified behaviors. Despite its importance, sex metadata is frequently missing in microbiome studies, including for animal-associated samples. Host chromosomal sex can be inferred from the host-derived reads present in metagenomic data, but existing genomic sex prediction tools rely on fixed coverage thresholds calibrated for human XY chromosomes and require relatively high host reads, limiting their use on low host-biomass samples such as stool and on organisms with other sex-determination systems. ResultsHere, we present SCiMS (Sex Calling in Metagenomic Sequences), a bioinformatic tool that leverages host-derived DNA within shotgun metagenomic data to predict host chromosomal sex, even at low host coverage. SCiMS uses a multinomial likelihood computed from observed read counts under each sex and reports chromosomal sex calls. Because the expected read distribution is derived directly from chromosome lengths and ploidy under each candidate karyotype, SCiMS applies to any organism with a heterogametic sex-determination system. We benchmarked SCiMS against existing tools on simulated metagenomic data, human metagenomic samples spanning multiple body sites, and metagenomic samples from seven animal species. SCiMS matched or outperformed existing tools, with its noticeable advantage at low host read conditions. ConclusionsSCiMS provides an accurate, scalable, and cross-species generalizable solution for host chromosomal sex classification, even when host DNA is minimal. By enabling recovery of missing sex metadata, it serves as a quality-control tool analyses in microbiome research. SCiMS is freely available at http://github.com/davenport-lab/SCiMS.