不同食性雁形目鸟类线粒体基因组的适应性进化分析.pdf

生物技术进展生物技术进展 2023 年 第 13 卷 第 5 期 748 754Current Biotechnology ISSN 20952341研究论文研究论文Articles不同食性雁形目鸟类线粒体基因组的适应性进化分析张卫平1，2，邱冰滢1，2，张东升1，2*1.上海海洋大学水产与生命学院，水产科学国家级实验教学示范中心，上海 201306；2.上海海洋大学水产与生命学院，环境DNA技术与水生态健康评估工程中心，上海 201306摘要：以往的研究表明，动物的食性与能量代谢及生存适应密切相关。能量代谢主要发生在线粒体中，线粒体编码的13个蛋白质亚基是氧化磷酸化复合体的重要组成部分。雁形目鸟类的食性主要包括肉食性、杂食性和植食性3种类型。为了分析食性对鸟类基因组进化的影响，研究选取了20种雁形目鸟类，并根据食性分成3组，下载其线粒体基因组，通过适应性进化分析、放松性选择分析、多态性氨基酸位点检测以及3D结构预测分析，研究这3组鸟类的线粒体蛋白编码基因在进化上的表现。结果发现，食肉组鸟类线粒体基因组的进化速率高于食草组和杂食组鸟类，并且只有食肉组鸟类线粒体蛋白编码基因受到了放松性的选择压力作用。此外，多态性氨基酸位点检测表明，食肉组线粒体基因组编码的蛋白质中，多态性位点和有害位点数量远高于食草组和杂食组。而杂食组鸟类线粒体蛋白编码基因进化速率较低，大部分基因受到强化性选择作用，蛋白质序列中的多态性位点也较少。研究结果表明，不同食性的雁形目鸟类，其线粒体基因组编码的蛋白质受到了不同的选择压力，这为食性差异影响鸟类线粒体基因组适应性进化提供了分子依据。关键词：雁形目；食性；线粒体；适应性进化DOI：10.19586/j.20952341.2023.0054 中图分类号：Q951 文献标志码：AAdaptive Evolution Analysis of Mitochondrial Genomes in Anseriform Birds with Various Feeding HabitsZHANG Weiping1，2，QIU Bingying1，2，ZHANG Dongsheng1，2*1.National Demonstration Center for Experimental Fisheries Science Education，College of Fisheries and Life Science，Shanghai Ocean University，Shanghai 201306，China；2.Environmental DNA Technology and Water Ecological Health Assessment Engineering Center，College of Fisheries and Life Science，Shanghai Ocean University，Shanghai 201306，ChinaAbstract：Previous studies showed that，animals feeding habits are related to their energy metabolism and adaptation.Energy metabolism takes place in mitochondria，and 13 subunits of the oxidative phosphorylation complexes were encoded by mitochondrial genome.According to feeding habits，birds of Anteriformes could be divided into three groups，including carnivorous，herbivorous and omnivorous.To study the evolutionary effects of feeding habit on mitogenomes，we selected 20 anteriformes species and divided them into three groups.We downloaded their mitogenomes，and did adaptive analysis，relax analysis，polymorphism analysis and 3D structure prediction.Our results showed that，the carnivorous group had higher evolutionary rates than the other two groups，and only the carnivorous group had multiple genes under relaxed selection.Meanwhile，the carnivorous group had more polymorphism sites and deterious sites among their proteins sequences than the other two groups.The genes from the omnivorous group showed decreased evolutionary rates，and most of the genes were under intensification selection，furthermore，less polymorphic sites were detected in this group.The results indicated that the proteins encoded by mitochondrial genomes of anteriformes birds with different feeding habits are subject to different selective pressures.This study provides a molecular basis for the influence of dietary differences on the adaptive evolution of avian mitochondrial genomes.Key words：anteriformes；diets；mitochondria；adaptive evolution收稿日期：20230417；接受日期：20230517基金项目：国家重点研发计划项目（2022YFC2601301）。联系方式：张卫平 E-mail:；*通信作者 张东升 E-mail:张卫平，等：不同食性雁形目鸟类线粒体基因组的适应性进化分析世界上有10 000多种鸟类1，它们生活于不同的环境中，并演化出相应的解剖学和行为学特征，其中鸟类的食性是备受关注的特征之一。大量研究表明，鸟类通过适应不同的食物资源来推动其形态、生理和行为特征的进化2-3。例如，加拉帕戈斯群岛达尔文雀是摄食适应的典型代表，食物的选择影响喙大小和形状的进化4-5。鸟类食性多样，几乎包括所有类型的食物，如无脊椎动物、鱼虾、内脏腐肉、水果、花蜜、果实、植物草茎、杂草、苔藓等，许多鸟类吃一种或几种特定类型的食物，而有些鸟类食谱较广6-7。食物基质的不同导致这些物质被摄取后所需的酶、机制以及消化吸收的时间存在差异8。不少研究表明，鸟类和哺乳动物的基本代谢速率与自身的食性有关9-12。此外，鸟种的食物构成对物种的形成和灭绝产生深刻影响6，鸟类是研究饮食对物种形成和灭绝作用的一个很好的模型系统13-14。动物通过降解食物中的营养物质来获取能量，而线粒体是细胞产生能量的主要细胞器，为真核细胞提供95%的能量15-17。同时，线粒体自身的基因组编码13个氧化磷酸化途径中的重要亚基，因此线粒体及其基因组在动物的能量代谢中发挥重要作用18-19。线粒体基因组编码基因参与能量代谢，并在动物的长期演化过程中受到适应性进化压力20-22。研究发现，基础代谢率较高的哺乳动物比基础代谢率较低的哺乳动物对线粒体 DNA 的净化选择更强23。此外，丧失飞行能力的鸟类比飞行能力强的鸟类在mtDNA中积累了更多的非同义取代，这可能是因为强运动物种维持高效的能量代谢，mtDNA经历了强纯化选择24。雁形目鸟类是备受关注的水鸟种群之一，全世界共有 3科 178种 1。雁形目鸟类食谱较广，包括杂食性、草食性和肉食性物种。草食性雁形目主要以植物草茎、杂草、苔藓、真菌、豆类为食，肉食性雁形目主食鱼和无脊椎动物，专食化程度较高7。杂食性雁形目则以植物草茎、杂草、苔藓、真菌、豆类以及无脊椎动物为食7。本研究以 20 种雁形目鸟类为研究对象，通过分析不同食性组雁形目鸟类的线粒体编码蛋白基因，以探讨雁形目鸟类食性的差异与这些基因进化表现之间的关系。1材料与方法1.1数据收集与序列比对本研究选取了20种雁形目鸟类（Anseriformes birds）作为研究对象。将食性组成中食肉比例达80%及以上的鸟类划分为食肉组；食草比例达80%及以上的鸟类划分为食草组；两者皆食且食草或食肉比例均不超过80%的鸟类视为杂食组。依据食性差异，这些鸟类被划分为食肉组、食草组和杂食组 3 组。其中，食草组有 9 种：Anser albifrons（NC_004539.1），Anser anser（NC_011196.1），Anser cygnoides（NC_023832.1），Anser fabalis（HQ890328.1），Anser indicus（NC_025654.1），Branta canadensis（NC_007011.1），Cygnus columbianus（NC_007691.1），Cygnus cygnus（NC_027095.1），Cygnus olor（NC_027096.1）；食肉组有7 种：Bucephala clangula（MW849283.1），Clangula hyemalis（MW849278.1），Histrionicus histrionicus（MW849288.1），Mergus merganser（NC_040986.1），Mergus squamatus（NC_016723.1），Polta stelleri（MW849289.1），Mergus serrator（MZ365040.1）；杂食组有 4 种：Aythya baeri（NC_057048.1），Aythya fer-ina（NC_024602.1），Aythya fuligula（NC_024595.1），Aythya nyroca（NC_059801.1）。这些物种的序列信息均来源于NCBI数据库，分类信息来源于EltronTraits 1.0数据库7。1.2分子进化分析本研究以鸡形目雉鸡（Phasianus colchicus）作为外群，雉鸡与雁形目同属鸡雁小纲。将13个线粒体蛋白编码基因删除终止子，再拼接成全长编码序列，即拼接序列（spliced data），使用Mega1125程序的Muscle26比对方法，对拼接序列进行比对，生成.fas比对文件，并以比对好的拼接序列为基础，通过最大似然法构建系统发育树。同时，按照相同方法生成不含外群的.fas比对文件和.tree树文件，用于Paml和Relax的输入文件。1.3放松性选择分析借助 Relax27在线服务器（https：/www.datamonkey.org/RELAX）对各基因进行放松性选择分析，选择肉食性、杂食性和草食性的鸟类中其中1组作为前景支，以另外2组鸟类作为背景支，来检测这些分支经历的选择压力。当K值大于1，则表示749生物技术进展生物技术进展 Current Biotechnology该前景支受到了强化性选择；当K值小于1，则表示该前景支受到了放松性选择。1.4选择压力分析本研究使用了Paml 4.828中的Codeml程序评估3组鸟线粒体13个编码基因受到的选择压力，且以自由比率模型（free-ratio model）对拼接序列以及每一个编码基因进行适应性进化分析。当dN/dS1时，预示着该基因很可能受到了正选择作用，正选择也是物种适应性进化的重要基础。Codeml 程序中的支-位点模型（branch-site model）存在模型A和null假设，根据2种模型之间的似然比检验推测某些位点是否经历了正选择作用，当P0.05时表示有充足证据可以判断前景支中存在正选择位点。用位点模型（site model）探索每个位点间选择压力的差异，位点模型假定同一分支上不同位点的值（表示进化速率）都不同，但在进化树不同分支上却无差异。1.5氨基酸多态性及有害突变检测通过 Python 编程获取氨基酸多态性位点信息，并使用Provean程序29评估这些多态性位点的有害性。Score（预测的分值）=-2.5是判断有害位点或中性突变位点的临界值，当预测分值小于-2.5，可视为有害突变位点，反之视为中性突变。考虑到3组鸟类的物种数不一致（食肉组7；食草组 9；杂食组 4），本研究利用氨基酸多态性比例多态性位点数/（该基因氨基酸数该组物种数）和有害性比例 有害突变位点数/（该基因氨基酸数该组物种数）表示3组鸟类间的差异。1.6蛋白质三级结构预测利用Swiss-Model 30 对本研究涉及的所有线粒体蛋白质进行三维结构构建，并且使用PyMOL进行修饰，观察3组食性鸟类中这些位点的分布情况。2结果与分析2.1分子进化分析从NCBI数据库中下载已选定的20种雁形目鸟类的线粒体蛋白编码序列，并根据食性分为食肉组、杂食组和食草组。利用 Mega1125程序的Muscle26对编码线粒体蛋白质的 13 个基因拼接后的序列进行比对，构建的系统发育树（图1）显示实验选取的雁鸭类鸟类中，相同食性的鸟种在进化树中聚类为一个单系群。2.2放松性选择分析通过Relax在线分析线粒体蛋白编码基因在不同食性组中受到的放松性选择情况，结果发现（表 1），食 肉 组 有 COX2、COX3、CYTB、ND1、图1基于拼接序列构建的雁形目鸟类系统发育树Fig.1Phylogenetic tree based on the spliced data of Anseriformes birds750张卫平，等：不同食性雁形目鸟类线粒体基因组的适应性进化分析ND3、ND4 共 6 个基因受到了放松性选择，且没有任何基因受到强化选择；食草组有ND2和ND6共 2个基因受到了放松性选择，有 ATP6、COX2、COX3、ND3共4个基因受到了强化性选择；而杂食组有 COX1、COX2、CYTB、ND2、ND3 共 5 个基因受到了强化性选择，且没有任何基因受到放松性选择。2.3选择压力分析为了估算谱系特异性值，通过Codeml程序的自由比率模型（free-ratio model）来评估不同分支的进化速率，得到的编码数据包括每个谱系基因的N、S、dN、dS和值。对3组数据（dN，dS，）中的3类鸟（拼接序列）做两两之间的t检验，结果显示（图2），食肉组的dS值显著高于食草组和杂食组（P0.05）。对13个基因分别进行检测并作t检验，食肉组与杂食组进行t检验结果表明，食肉组有10个基 因 的 dS（ATP8、COX1、COX2、COX3、CYTB、ND1、ND2、ND4、ND5、ND6）、7 个 基 因 的 dN（ATP8、COX2、COX3、CYTB、ND4、ND5、ND6）显著高于杂食组；而食肉组与食草组进行比较表明，食肉 组 在 dS（ATP8、COX1、COX2、COX3、CYTB、ND1、ND2、ND4、ND5、ND6）、dN（ATP8、COX3、CYTB、ND4、ND5）显著高于食草组；但当食草组与杂食组相比较时，dS、dN均未表现出显著性，仅CYTB基因在杂食组中的值显著高于食草组。此外，本研究采用 Codeml 程序的支-位点（branch-site）模型和 site 模型探究哪些基因受到了正选择，结果均未发现正选择位点。2.4氨基酸多态性分析基于氨基酸序列的比对结果，通过Python编程得到氨基酸多态性位点信息，并使用Provean程序29评估这些多态性位点的有害性31。结果发现，在食肉组的这些蛋白质序列中，多态性位点和有害性位点出现的频率均高于另外两组鸟。对13个线粒体蛋白分别检测的结果（图3）表明，在食肉组的13个基因中，其中有11个基因的氨基酸多态性频率高于另外两组（图3B），有7个基因的有害突变频率高于另外两组（图3C）。2.5线粒体蛋白质三级结构分析本研究使用 Swiss-Model 在线网址预测了13种蛋白质的三级结构，利用Pymol对中性突变位点和有害位点进行标记：黄色标记中性突变，红表1结果中有显著放松性选择或显著强化性选择的基因Table 1The significant results of the relaxation and intensification选择类型放松性选择强化性选择测试支食肉组食草组食草组杂食组基因COX2COX3CYTBND1ND3ND4ND2ND6ATP6COX2COX3ND3COX1COX2CYTBND2ND3P值0.0270.0010.0270.0040.0450.00100.00200.01200.0420.0170.0020.0060.0340.028K值0.740.050.790.800.620.6400.076.5414.135.7022.0224.932.2341.174.0644.75注：dN非同义替换率；dS同义替换率；进化速率。*表示两组数据间差异显著(P0.01)图2食肉组、食草组、杂食组的dN、dS和比较Fig.2Comparisons of dN,dS and values among the groups of carnivore birds,herbivorous birds and omnivorous birds751生物技术进展生物技术进展 Current Biotechnology色标记有害突变。结果（图4）显示，3组食性鸟类的所有编码基因的蛋白质三维结构位点均呈现多态性，有害位点并未集中分布且主要存在于食肉组中。3讨论由于食物资源的时空分布不一，鸟类演化出了不同的饮食偏好和取食方式及形态特征，因此，普遍认为鸟类的食性选择会影响物种的进化。以往的研究主要集中在食性与生理结构8，32、行为特A：拼接序列（spliced data）的氨基酸多态性和有害性位点频率；B：13个蛋白质序列的多态性位点频率；C：13个蛋白质序列的有害性位点频率。图3食肉组、食草组、杂食组的氨基酸多态性和有害性的频率Fig.3Proportion of amino acid polymorphism sites and amino acid harmful sites in carnivore,herbivorous and omnivorous birdsA：三组食性鸟类中ND2的3D结构；B：三组食性鸟类ND2的氨基酸多态性位点频率；C：三组食性鸟类ND2的有害性位点比例。图4食肉组、食草组、杂食组中ND2的3D结构预测和氨基酸多态性、有害性位点的比例Fig.43D structure prediction and proportion of polymorphism sites and harmful sites of ND2 in carnivore,herbivorous and omnivorous birds752张卫平，等：不同食性雁形目鸟类线粒体基因组的适应性进化分析征2-3及特定基因33-34之间的关系等，而极少在基因组层面研究两者的关系。为了在基因组层面探讨两者的关系，本研究选取20种雁形目鸟类，利用Mega、Paml、Relax、Python、Swiss-Model、Pymol等生物信息学软件或方法对线粒体基因组进行进化速率分析、选择性分析、氨基酸多态性和有害性分析、蛋白质3D结构分析等。已有对鸟类等动物的研究中通常认为线粒体编码的基因在能量代谢相关的特征进化过程中受到强烈的选择压力18-19，35，纯化选择的放松可能与进化速度的加快有关36-38。本研究对鸟类食性与线粒体基因组的这些编码基因做了进化分析，结果也证实了这种现象。与杂食组和食草组鸟类相比，食肉组的雁形目鸟类大部分编码蛋白的线粒体基因处于放松性选择，而杂食组鸟类的强化性选择信号显著。利用Paml-Codeml对3种食性鸟类进行两两间的进化分析比较，虽然值没有表现出显著性，但是食肉组鸟类的dS和dN值显著高于另外两组鸟类，杂食组鸟类的dS和dN值显著偏低。此外，不同食性组鸟类线粒体基因组在氨基酸多态性及氨基酸有害性方面也表现出明显的差异，食肉组中有11个基因表现出较高的氨基酸多态性频率，有7个基因表现出较高的氨基酸有害性频率。研究结果表明，摄取肉食比例较高的雁形目鸟类，其线粒体基因组受到放松性选择作用。专门的蜜食动物和果食动物多在“四季常青”的地方生存，食虫动物则在全年都有昆虫的热带环境中生存39，因此资源的可获得性相对稳定。杂食动物在不稳定的资源和环境中更有优势40。但是，在生态竞争的竞技场中，杂食动物作为可食多种食物资源的“多面手”，在利用特定食物资源方面表现不佳41，因此在物种宏观进化方面，表现出低物种形成率或高灭绝率。Burin等7通过分析比较现存的近万中鸟类的进化关系和食性，发现杂食性鸟类是鸟类食性的“进化陷阱”，同其他食性的鸟类相比，杂食性鸟类的灭绝率较高，多样化率较低。但是，该研究并未涉及这些现象的分子机制。本研究从基因组层面揭示杂食性的雁形目鸟类线粒体基因组受到强化性的选择压力，氨基酸多态性下降，有害突变位点被淘汰，为杂食性鸟类物种的低多样化率提供了一定的分子依据。参 考 文 献 1 GILL F D D.IOC world bird list(version 13.1)EB/OL.2022-12-20 International Ornithologists Union,2022.https:/www.worldbirdnames.org/new/.2 ODONNELL S,LOGAN C J,CLAYTON N S.Specializations of birds that attend army ant raids:an ecological approach to cognitive and behavioral studiesJ.Behav.Processes,2012,91(3):267-274.3 ABRAHAMCZYK S,KESSLER M.Morphological and behavioural adaptations to feed on nectar:how feeding ecology determines the diversity and composition of hummingbird assemblagesJ.J.Ornithol.,2015,156(2):333-347.4 PR G.Ecology and evolution of Darwins finchesM.Princeton University Press,1999:565-566.5 OLSEN A M.Feeding ecology is the primary driver of beak shape diversification in waterfowlJ.Funct.Ecol.,2017,31(10):1985-1995.6 WILMAN H,BELMAKER J,SIMPSON J,et al.EltonTraits 1.0:species-level foraging attributes of the worlds birds and mammalsJ/OL.Ecology,2014,95(7):20272014-07-26.https:/doi.org/10.1890/13-1917.1.7 BURIN G,KISSLING W D,GUIMARES P R,et al.Omnivory in birds is a macroevolutionary sinkJ/OL.Nat.Commun.,2016,7:112502022-12-30.https:/org.doi/10.1038/ncomms11250.8 KARASOV W H,CARLOS M,CAVIEDES-VIDAL E.Ecological physiology of diet and digestive systemsJ/OL.Annu.Rev.Physiol.,2011,73(1):692022-12-30.https:/org.doi/10.1146/annurev-physiol-012110-142152.9 MCNAB B K.The economics of temperature regulation in neotropical batsJ.Comp.Biochem.Physiol.,1969,31(2):227-268.10 MCNAB B K.The influence of food habits on the energetics of eutherian mammalsJ.Ecol.Monogr.,1986,56(1):1-19.11 MCNAB B K.Food habits and the basal rate of metabolism in birdsJ.Oecologia,1988,77(3):343-349.12 MCNAB B K.Energetics,body size,and the limits to endothermyJ.J.Zool.,2010,199(1):1-29.13 BROWN W L,WILSON E O.Character displacementJ/OL.Syst.Zool.,1956,5(2):492023-07-26.https:/doi.org/10.2307/411924.14 GRANT P R,GRANT B R.Evolution of character displacement in Darwins finchesJ.Science,2006,313(5784):224-226.15 ROGGE M M.The role of impaired mitochondrial lipid oxidation in obesityJ.Biol.Res.Nurs.,2009,10(4):356-373.16 OSELLAME L D,BLACKER T S,DUCHEN M R.Cellular and molecular mechanisms of mitochondrial functionJ.Best Pract.Res.Clin.Endocrinol.Metab.,2012,26(6):711-723.17 KUSMINSKI C M,SCHERER P E.Mitochondrial dysfunction in white adipose tissueJ.Trends Endocrinol.Metab.,2012,23(9):435-443.18 BALLARD J W O,WHITLOCK M C.The incomplete natural history of mitochondriaJ.Mol.Ecol.,2004,13(4):729-744.19 RAND D M.The units of selection on mitochondrial DNAJ.Annu.Rev.Ecol.Syst.,2001,32:415-448.20 SHEN Y Y,LIANG L,ZHU Z H,et al.Adaptive evolution of 753生物技术进展生物技术进展 Current Biotechnologyenergy metabolism genes and the origin of flight in batsJ.Proc.Natl.Acad.Sci.USA,2010,107(19):8666-8671.21 YANG Y,XU S,XU J,et al.Adaptive evolution of mitochondrial energy metabolism genes associated with increased energy demand in flying insectsJ/OL.PLoS ONE,2014,9(6):e991202023-07-26.http:/org.doi/10.1371/journal.pone.0099120.22 RAMIREZ J,FOLKOW L P,BLIX A S.Hypoxia tolerance in mammals and birds:from the wilderness to the clinicJ.Annu.Rev.Physiol.,2007,69:113-143.23 FENG P,ZHAO H,LU X.Evolution of mitochondrial DNA and its relation to basal metabolic rateJ.Mitochon.DNA,2015,26(4):566-571.24 SHEN Y Y,SHI P,SUN Y B,et al.Relaxation of selective constraints on avian mitochondrial DNA following the degeneration of flight abilityJ.Genome Res.,2009,19(10):1760-1765.25 KOICHIRO T,GLEN S,SUDHIR K.MEGA11:molecular evolutionary genetics analysis version 11J/OL.Mol.Biol.Evol.,2021,(7):72022-12-30.https:/org.doi/10.1093/molbev/msab120.26 EDGAR R C.MUSCLE:a multiple sequence alignment method with reduced time and space complexityJ/OL.BMC Bioinform.,2004,5:1132022-12-30.https:/org.doi/10.1186/1471-2105-5-113.27 WERTHEIM J O,MURRELL B,SMITH M D,et al.LAX:detecting relaxed selection in a phylogenetic frameworkJ.Mol.Biol.Evol.,2015,32(3):820-832 28 YANG Z.PAML 4:phylogenetic analysis by maximum likelihoodJ.Mol.Biol.Evol.,2007,24(8):1586-1591.29 CHOI Y,SIMS G E,MURPHY S,et al.Predicting the functional effect of amino acid substitutions and indelsJ/OL.PLoS ONE,2012,7(10):e466882022-12-30.https:/org.doi/10.1371/journal.pone.0046688.30 BIASINI M,BIENERT S,WATERHOUSE A,et al.SWISS-MODEL:mdelling protein tertiary and quaternary structure using evolutionary informationJ.Nucleic Acids Res.,2014,42:252-258.31 胡田伟,王世航,张东升.线粒体基因组在红血南极鱼和白血南极鱼中的适应性进化J.基因组学与应用生物学,2021,40(Z4):3448-3457.32 李福来.鸟类食性与消化道的特点J.生物学通报,1984,19(2):24-26+10.33 ZHANG G,LI C,LI Q,et al.Comparative genomics reveals insights into avian genome evolution and adaptationJ.Science,2014,346(6215):1311-1320.34 周明.鸟类TLRs基因的适应性进化研究D.雅安:四川农业大学,2019.35 SUN Y B,SHEN Y Y,IRWIN D M,et al.Evaluating the roles of energetic functional constraints on teleost mitochondrial-encoded protein evolutionJ.Mol.Biol.Evol.,2011,28(1):39-44.36 OMETTO L,LI M,BRESADOLA L,et al.Rates of evolution in stress-related genes are associated with habitat preference in two Cardamine lineagesJ/OL.BMC Evol.Biol.,2012,12:72022-12-30.https:/org.doi/10.1093/molbev/msq256.37 STROHM J H T,GWIAZDOWSKI R A,HANNER R.Fast fish face fewer mitochondrial mutations:patterns of dN/dS across fish mitogenomesJ.Gene,2015,572(1):27-34.38 LU Y,XING H,ZHANG D.Evidence for relaxed selection of mitogenome in rapid-flow cyprinidsJ.Genes Genom.,2019,41(7):863-869.39 KISSLING W D,SEKERCIOGLU C H,JETZ W.Bird dietary guild richness across latitudes,environments and biogeographic regionsJ.Glob.Ecol.Biogeogr.,2012,21(3):328-340.40 MACARTHUR R,LEVINS R.The limiting similarity,convergence,and divergence of coexisting speciesJ.Am.Nat.,1967,101(921):377-385.41 MACARTHUR R H,PIANKA E R.On optimal use of a patchy environmentJ.Am.Nat.,1966,100(916):603-609.754