Aphids live in symbiosis with a variety of bacteria, including the obligate symbiont Buchnera aphidicola and diverse facultative symbionts. The symbiotic associations for one aphid species, especially for polyphagous species, often differ across populations. In the present study, by using high-throughput 16S rRNA sequencing, we surveyed in detail the microbiota in natural populations of the cotton aphid Aphis gossypii in China and assessed differences in bacterial diversity with respect to host plant and geography. The microbial community of A. gossypii was dominated by a few heritable symbionts. Arsenophonus was the most dominant secondary symbiont, and Spiroplasma was detected for the first time. Statistical tests and ordination analyses showed that host plants rather than geography seemed to have shaped the associated symbiont composition. Special symbiont communities inhabited the Cucurbitaceae-feeding populations, which supported the ecological specialization of A. gossypii on cucurbits from the viewpoint of symbiotic bacteria. Correlation analysis suggested antagonistic interactions between Buchnera and coexisting secondary symbionts and more complicated interactions between different secondary symbionts. Our findings lend further support to an important role of the host plant in structuring symbiont communities of polyphagous aphids and will improve our understanding of the interactions among phytophagous insects, symbionts, and environments.

Symbiotic association is universal in nature, and an array of symbionts play a crucial part in host life history. Aphids and their diverse symbionts have become a good model system to study insect-symbiont interactions. Previous symbiotic diversity surveys have mainly focused on a few aphid clades, and the relative importance of different factors regulating microbial community structure is not well understood. In this study, we collected 65 colonies representing eight species of the aphid genus Mollitrichosiphum from different regions and plants in southern China and Nepal and characterized their microbial compositions using Illumina sequencing of the V3 − V4 hypervariable region of the 16S rRNA gene. We evaluated how microbiota varied across aphid species, geography and host plants and the correlation between microbial community structure and host aphid phylogeny. Heritable symbionts dominated the microbiota associated with Mollitrichosiphum, and multiple infections of secondary symbionts were prevalent. Ordination analyses and statistical tests highlighted the contribution of aphid species in shaping the structures of bacterial, symbiont and secondary symbiont communities. Moreover, we observed a significant correlation between Mollitrichosiphum aphid phylogeny and microbial community composition, providing evidence for a pattern of phylosymbiosis between natural aphid populations and their microbial associates.