Adaptive evidence of mitochondrial genomes in Dolycoris baccarum (Hemiptera: Pentatomidae) to divergent altitude environments

Given mitochondrion is the ‘energy and oxygen usage factories’, adaptive signatures of mitochondrial genes have been extensively investigated in vertebrates from different altitudes, but few studies focus on insects. Here, we sequenced the complete mitochondrial genome (mitogenome) of Dolycoris. baccarum living in the Tibetan Plateau (DBHC, ∼3200 m above sea level (asl)) and conducted a detailed comparative analysis with another D. baccarum mitogenome (DBQY) from relatively low altitude (∼1300 m asl). All the 37 mitochondrial genes were highly conserved and under purifying selection, except for two mitochondrial protein-coding genes (MPCGs) (atp6 and nad5) that showed positively selected signatures. We therefore further examined non-synonymous substitutions in atp6 and nad5, by sequencing more individuals from three populations with different altitudes. We found that these non-synonymous substitutions were polymorphic in these populations, likely due to relaxed selection constraints in different altitudes. Purifying selection in all mitochondrial genes may be due to their functional importance for the precision of ATP production usually. Length difference in mitochondrial control regions between DBHC and DBQY was also conversed at the population level, indicating that sequence size adjustments in control region may be associated with adaptation to divergent altitudes. Quantitatively real-time PCR analysis for 12 MPCGs showed that gene expression patterns had a significant change between the two populations, suggesting that expression levels of MPCGs could be modulated by divergent environmental pressures (e.g. oxygen content and ambient temperature). These results provided an important guide for further uncovering genetic mechanisms of ecological adaptation in insects.