Large-scale Phylogenomic Analysis Provides New Insights Into the Phylogeny of the Order Gadiformes and Evolution of Freshwater Gadiform Species Burbot (Lota lota)

Our understanding of phylogenetic relationships among Gadiformes fish is obtained through the analysis of a small number of genes, but uncertainty remains around critical nodes. A series of phylogenetic controversial exists at the suborder, family, subfamily, and species levels. A total of 1105 orthologous exon sequences and translated amino acid sequences from 36 genomes and 12 transcriptomes covering 33 species were applied to investigate the phylogenetic relationships within Gadiformes and address these problems. Phylogenetic trees reconstructed with the amino acid data set using different tree-building methods (RAxML and MrBayes) showed consistent topology. The monophyly of Gadifromes was confirmed in our study. However, the three suborders Muraenolepidoidei, Macrouroidei, and Gadoidei were not well recovered by our phy-logenomic study, rejecting the validity of suborder Muraenolepidoidei. Four major lineages were revealed in this study. The family Bregmacerotidae forming clade I was the basal lineage of Gadiformes. The family Merluciidae formed clade II. Clade III contained families Melanonidae, Muraenolepididae, Macrouridae (with subfamilies Trachyrincinae, Macrourinae, and Bathygadinae), and Moridae. Clade IV contained at least three families of suborder Gadoidei, i.e., Gadidae, Phycidae, and Ranicipitidae. The subspecies of Lota lota from Amur River were confirmed, indicating that exon markers were a valid high-resolution method for delimiting subspecies or distinct lineages within species level. The PSMC analysis of different populations of L. lota suggests a continuous decline since 2 Myr.