Identifying a cohort of ovary-specific genes across drosophilid species.

A, Volcano plot for one example species, D. primaeva (Dpri), showing the results of a differential gene expression analysis comparing the ovary to the carcass. The x-axis shows the log2 fold change of expression across transcripts, and the y-axis shows the adjusted p-value, log10 transformed. Points that are significantly differentially expressed are shown in black. B, Jitter plots showing the results of the same analysis across the twelve species studied here. The x-axis shows the log2 fold change of expression across transcripts, and points are arranged with random jitter within species on the y-axis. C, The same jitter plots as in B, now colored according to whether or not transcripts belong to a cohort of core ovary genes. These are defined as genes, grouped by BLAST similarity to D. melanogaster transcripts, for which at least one transcript is upregulated in the ovary of ten or more of the twelve species. D, A comparison of mean expression change across Hawaiian species to reported ovary-enrichment values from D. melanogaster, as reported in FlyAtlas2 [29]. Core ovary genes are marked in magenta. E, The boxed region shown in D, magnified and now showing only core ovary genes, annotated with the gene symbol from D. melanogaster.