Collembola are a basal group of Hexapoda renowned for both unique morphological characters and significant ecological roles. However, a robust and plausible phylogenetic relationship between its deeply divergent lineages has yet to be achieved. We carried out a mitophylogenomic study based on a so far the most comprehensive mitochondrial genome dataset. Our data matrix contained mitogenomes of 31 species from almost all major families of all four orders, with 16 mitogenomes newly sequenced and annotated. We compared the linear arrangements of genes along mitochondria across species. Then we conducted 13 analyses each under a different combination of character coding, partitioning scheme and heterotachy models, and assessed their performance in phylogenetic inference. Several hypothetical tree topologies were also tested. Mitogenomic structure comparison revealed that most species share the same gene order of putative ancestral pancrustacean pattern, while seven species from Onychiuridae, Poduridae and Symphypleona bear different levels of gene rearrangements, indicating phylogenetic signals. Tomoceroidea was robustly recovered for the first time in the presence of all its families and subfamilies. Monophyly of Onychiuroidea was supported using unpartitioned models alleviating LBA. Paronellidae was revealed polyphyletic with two subfamilies inserted independently into Entomobryidae. Although Entomobryomorpha has not been well supported, more than half of the analyses obtained convincing topologies by placing Tomoceroidea within or near remaining Entomobryomorpha. The relationship between elongate-shaped and sphericalshaped collembolans still remained ambiguous, but Neelipleona tend to occupy the basal position in most trees. This study showed that mitochondrial genomes could provide important information for reconstructing the relationships among Collembola when suitable analytical approaches are implemented. Of all the data refining and model selecting schemes used in this study, the combination of nucleotide sequences, partitioning model and

Colibacillosis, caused by Escherichia coli , is one of the most common bacterial diseases of chickens. The high incidence and considerable economic losses associated with colibacillosis make it a significant concern worldwide. In recent years, the efficacy of colistin has been severely impacted by the emergence of plasmid-mediated colistin resistance genes, especially mcr-1 . Therefore, monitoring of antibiotic resistance, particularly colistin resistance, amongst E. coli strains is vitally important to the future growth and sustainability of the poultry industry. In this study, a total of 130 E. coli strains were isolated from the livers of chickens displaying symptoms of colibacillosis in Tai'an, China. Isolates were screened for their susceptibility to various antibiotics and for the presence of mobile colistin resistance genes and other antibiotic resistance genes. Overall, 75 (57.7%) isolates showed resistance to colistin and were positive for mcr-1 . The mobile colistin resistance genes, mcr-2 , -3 , and -4 , were not detected in this study. Of the 75 mcr-1 -positive isolates, all (100%) also carried tetracycline resistance genes, 71 (94.7%) also contained genes associated with β-lactam resistance, 59 (78.7%) contained aminoglycoside resistance genes, and 57 (76%) contained sulfonamide resistance genes. This high prevalence of multidrug resistance among mcr-1 -positive E. coli isolates, including the production of extended-spectrum β-lactamases, is highly concerning. The surveillance findings presented here will be conducive to our understanding of the prevalence and characteristics of multidrug-resistance in E. coli in the Tai'an area and will provide a better scientific basis for the clinical treatment of colibacillosis in chickens.

A novel MDR plasmid, pXZ, harbouring four resistance genes (rmtB, fosA3, bla(TEM-1) and bla(CTX-M)) was identified. To our knowledge, this is the first report of an aminoglycoside resistance plasmid harbouring the fosA3 gene.