Characterization of the Complete Mitochondrial Genome of Sea Duck Mergus serrator and Comparison with other Anseriformes Species
Characterization of the Complete Mitochondrial Genome of Sea Duck Mergus serrator and Comparison with other Anseriformes Species
Peng Chen1, Jiaqi Li1, Hongbo Li2, Qin Lu3, Wei Liu1,4,5* and Jianliang Zhang1*
1Nanjing Institute of Environmental Sciences, Ministry of Ecology and Environment, Nanjing 210042, China
2National Marine Environmental Monitoring Center, Ministry of Ecology and Environment, Dalian 116023, China.
3Nantong Museum, Nantong 226001, China
4College of Marine Life Sciences and Frontiers Science Center for Deep Multispheres and Earth System, Ocean University of China, Qingdao 266003, China
5Institute of Evolution and Marine Biodiversity, Ocean University of China, Qingdao 266003, China.
ABSTRACT
As a key group in the Anatidae, the Mergini tribe (sea ducks) is strongly structured, with clear genetic assignments and phylogenetic relationships. The tribes also differ in terms of breeding strategy, longevity, and bionomic strategy. The mitochondrial genome (mitogenome) of Mergus serrator was 16,627 bp in length, and its structure was circular. The structure and gene arrangement of the genome were basically the same as those of Anseriformes species. The mean base composition of the mitogenome of Anseriformes was T (22.31 ± 0.51%), C (32.63 ± 0.64%), A (29.36 ± 0.64%), and G (15.71 ± 0.52%), indicating a slight specific bias towards A and C. AT content ranged of the mitogenome was from 50.27% to 55.31%, with an average value of 51.67 ± 1.10%, higher than the GC content and similar to that of birds in general (50.5% to 57.7%). In addition, the start and stop codons, the mitogenome consists of 3662 codons. The most commonly used amino acid was leucine (13.63%) in the use of M. serrator. The analyses indicated that the Anseriformes include the families Anseranatidae, Anhimidae, and Anatidae. Furthermore, Anatinae is composed of Aythyini, Anatini, Somaterini, and Mergini. M serrator was sister to M. merganser and M. squamatus, and this group belongs to Mergini. The ω value of the ND3 gene in the Mergini tribe is lower than those for other tribes. The phylogenetic relationships were analyzed and M. serrator was sister to M. merganser and M. squamatus, and formed a closely evolved Mergini clade. Different evolutionary rates between the Mergini tribe and other tribes were found.
Article Information
Received 14 March 2022
Revised 05 April 2022
Accepted 23 April 2022
Available online 09 August 2022
(early access)
Published 01 September 2023
Authors’ Contribution
CP and LW conceived and designed the study. LW performed the statistical analysis, under supervision of LJ, LH, and ZJ. CP wrote the first draft of the manuscript, and the final version included edits from all authors. QL, CP and WL revised this manuscript.
Key words
Anseriformes, Mitochondrial DNA, Phylogenetic construction, Seabird, Mergus serrator
DOI: https://dx.doi.org/10.17582/journal.pjz/20220314030351
* Corresponding author: lwecology@126.com
0030-9923/2023/0005-2365 $ 9.00/0
Copyright 2023 by the authors. Licensee Zoological Society of Pakistan.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
INTRODUCTION
Mitochondrial genome (mitogenome) harbors useful molecular information that can be used to infer phylogenetic relationships among Anseriformes different taxa (Livezey, 1997; Donne-Goussé et al., 2002; Jiang et al., 2009; Sun et al., 2017), one of the best-studied groups of birds (Livezey, 1986). Anseriformes includes 3 families, 8 tribes, and more than 150 species (Johnsgard, 2010; Ringelman et al., 2015). Due to gene rearrangement, genetic introgression, gene flow, and hybridization, the taxonomic statuses and phylogenetic placements within the order remain in dispute (Scribner et al., 2001; Gonzalez et al., 2009; Lavretsky et al., 2021). The Mergini tribe includes 20 species of which two are extinct (Johnsgard, 1966). As a key group in the Anatidae, the Mergini tribe (sea ducks) is strongly structured, with clear genetic assignments and phylogenetic relationships (Johnsgard, 2010). Whereas, the relationships among the Mergini tribe species are not certain based on multiple approaches and different genes so far. For example, the phylogenetic status of Mergus serrator is inconsistent based on different genes (Liu et al., 2012, 2014; Lavretsky et al., 2021). There is a remarkable difference in the life history and behavioral traits among different tribes in the Anseriformes (Peters et al., 2012). Most tribes tend to live in inland streams, rivers, or ponds, such as the Anatini, Aythyini, and Anserini tribes (Debela et al., 2021; Ponomarenko et al., 2021). In contrast, the Mergini tribe roosts and forages in sea habitats (Johnsgard, 2010). The tribes also differ in terms of breeding strategy, longevity, and bionomic strategy (Savard et al., 1998; Johnsgard, 2010). The Mergini tribe as an excellent tribe to investigate molecular adaptions to elevate metabolic burden associated with other life-history characteristics. Phylogenetic analyses can help investigate the relationships between mitochondrial DNA (mtDNA) molecular evolution and metabolic performance (Sruoga et al., 2008; Sonsthagen al., 2011).
MtDNA is typically double-stranded and maternally (Wolstenholme, 1992; Roques et al., 2004), it is valuable for studies of avian taxonomy and phylogeny (Shen et al., 2009; Jetz et al., 2012; Ren et al., 2014), genetic structure (Ruan et al., 2018), biological identification (Lecroy and Barker, 2006; Zhu et al., 2021), and conservation genetics (Allendorf et al., 2010), as it is a neutral molecular marker (Torroni et al., 2001; Nabholz et al., 2008). As the main energy suppliers of the cell, mitochondria produce approximately 90% of cellular ATP reserves through the highly efficient chemical osmotic coupling of electron and proton transfer to ATP synthesis (Saraste, 1999; Starkov, 2008). Evaluating the selective pressure of environmental temperature and oxygen availability on changes in mtDNA molecular could provide key insights for the adaptive evolution of mitogenomes (Luo et al., 2013). The key to successful migration to a new habitat is whether it can adapt to different energetic demands. Energetic demand is essentially related to the ability of mitochondrial to generate energy through the process of oxidative phosphorylation, such as in certain positions encoded by the NADH dehydrogenase genes (ND4 and ND5) (Janssen et al., 2004; Liu et al., 2018). However, single or combined genes have obvious shortcomings in phylogenetic analysis, such as insufficient information (Liu et al., 2012). These are overcome when using complete sequences of mitochondria, and it is becoming the first-choice marker for resolving controversial species relationships (García et al., 2014; Peng et al., 2017).
In this study, we sequenced the mitogenome of a typical sea duck, M. serrator, to analyze the genetic structure of and phylogenetic relationships among the Mergini tribe. In recent years, an increasing number of mitogenomes of the Mergini tribe have been sequenced. This study aims to shed increased light on the phylogenetic status of M. serrator in the Mergini tribe and to elucidate the evolutionary rates and molecular signatures of natural selection for all 13 mtDNA protein coding genes (PCGs) in sea ducks.
MATERIALS AND METHODS
Ethics approval and consent to participate
The sample collection was strictly conducted under national ethical guidelines (Regulations for Administration of Affairs Concerning Experimental Animals, China, 1988) for animal husbandry and humane treatment.
Sample collection and genome sequencing
We collected a dead specimen of M. serrator from Dalian (38°51′52.22″N 121°39′39.11″E; Liaoning, China). A muscle sample was transferred to a 5 mL centrifuge tube, 1 mL 75% ethyl alcohol was added, and the sample was stored in a refrigerator at –20°C. The sample (specimen voucher: NIES20210320HXQSYS01) was deposited in the Museum of Laboratory of Biodiversity Investigation of Nanjing Institute of Environment Science.
Total genomic DNA was extracted from muscle tissue using the Easy Pure genomic DNA kit (Trans Gen Biotech Co, Beijing, China). The original sequence data (15.65 G) was deposited in NCBI’s Sequence Read Archive (SRA; accession: SRR13516389), and complete mitogenome sequencing (15.62 G) was performed on an Illumina NovaSeq6000 platform (Novogene Bioinformatics Technology Co. Ltd., Tianjin, China) (Tang et al., 2021). The mitogenome was assembled using the program NOVO Plasty 3.7 (Dierckxsens et al., 2017) and then adjusted manually. The annotated using the MITOS Web Server (Bernt et al., 2013). The assembled mitogenome sequences have been deposited in GenBank under accession numbers MZ365040.
Data acquisition and analysis
All sequences were aligned on Clustal X using the default options (Thompson et al., 1997). PCGs were identified by comparison with the corresponding known complete mtDNA sequence of M. merganser and M. squamatus using Sequin 11.0. The 22 tRNA genes of cloverleaf secondary structures and anticodon sequences were identified from MITOS WebServer (http://mitos.bioinf.uni-leipzig.de/index.py). The control region (CR) was identified by sequence homology analyses (Cadahia et al., 2009; Barker et al., 2012). The start and stop codons of 13 PCGs were identified in the M. serrator mitogenome by ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/) and use the CGView Server software to generate a graphical map of the mitogenome and modified it manually (Table I). The codon number and relative synonymous codon usage (RSCU) of three species of Mergus mitochondrial PCGs were obtained using MEGA X and PhyloSuite.
Table I. Organization of the complete mtDNAs of Mergus serrator.
Gene |
Location |
Length (bp) |
Start codon |
End codon |
Spacer (+) or overlap (-) |
Strand |
tRNAPhe |
1-72 |
72 |
-1 |
H |
||
12S rRNA |
72-1066 |
995 |
0 |
H |
||
tRNAVal |
1067-1137 |
71 |
0 |
H |
||
16S rRNA |
1138-2737 |
1600 |
0 |
H |
||
tRNALeu(UUR) |
2738-2811 |
74 |
6 |
H |
||
ND1 |
2818-3795 |
978 |
ATA |
AGG |
-2 |
H |
tRNAIle |
3794-3866 |
73 |
7 |
H |
||
tRNAGln |
3874-3944 |
71 |
-1 |
L |
||
tRNAMet |
3944-4012 |
69 |
0 |
H |
||
ND2 |
4013-5051 |
1039 |
ATG |
T-- |
0 |
H |
tRNATrp |
5052-5127 |
76 |
3 |
H |
||
tRNAAla |
5131-5199 |
69 |
1 |
L |
||
tRNAAsn |
5201-5273 |
73 |
0 |
L |
||
tRNACys |
5274-5337 |
64 |
0 |
L |
||
tRNATyr |
5338-5408 |
71 |
1 |
L |
||
COI |
5410-6960 |
1551 |
GTG |
AGG |
-9 |
H |
tRNASer(UCN) |
6952-7024 |
73 |
2 |
L |
||
tRNAAsp |
7027-7095 |
69 |
1 |
H |
||
COⅡ |
7097-7783 |
687 |
GTG |
TAA |
1 |
H |
tRNALys |
7785-7853 |
69 |
1 |
H |
||
ATP8 |
7855-8022 |
168 |
ATG |
TAA |
-10 |
H |
ATP6 |
8013-8696 |
684 |
ATG |
TAA |
-1 |
H |
COIII |
8696-9479 |
784 |
ATG |
T-- |
0 |
H |
tRNAGly |
9480-9548 |
69 |
0 |
H |
||
ND3 |
9549-9722 |
174 |
ATG |
TAG |
1 |
H |
9724-9900 |
177 |
1 |
H |
|||
tRNAArg |
9902-9971 |
70 |
0 |
H |
||
ND4L |
9972-10268 |
297 |
ATG |
TAA |
-7 |
H |
ND4 |
10262-11639 |
1378 |
ATG |
T-- |
0 |
H |
tRNAHis |
11640-11708 |
69 |
0 |
H |
||
tRNASer(AGY) |
11709-11774 |
66 |
-1 |
H |
||
tRNALeu(CUN) |
11774-11844 |
71 |
0 |
H |
||
ND5 |
11845-13668 |
1824 |
GTG |
TAA |
-1 |
H |
CYTB |
13668-14810 |
1143 |
ATG |
TAA |
2 |
H |
tRNAThr |
14813-14880 |
68 |
10 |
H |
||
tRNAPro |
14891-14960 |
70 |
10 |
L |
||
ND6 |
14971-15492 |
522 |
ATG |
TAG |
0 |
L |
tRNAGlu |
15493-15560 |
68 |
0 |
L |
||
CR |
15561-16627 |
1067 |
0 |
H |
Negative values indicate overlapping; H, heavy strand; L, light strand.
Table II. Nucleotide composition of the mitochondrial genomes of 53 Anseriformes species.
Family/Genus |
Species |
T |
C |
A |
G |
AT skew |
GC skew |
A+T |
G+C |
Total |
Accession number |
||||||||||
Anatidae |
|||||||||||||||||||||
Aix |
Aix galericulata |
22.35 |
32.77 |
29.21 |
15.66 |
0.13 |
-0.35 |
51.57 |
48.43 |
16605 |
NC023969 |
||||||||||
Anas |
Anas acuta |
22.12 |
32.86 |
28.78 |
16.24 |
0.13 |
-0.34 |
50.90 |
49.10 |
16599 |
NC024631 |
||||||||||
A. crecca |
22.35 |
32.64 |
29.05 |
15.96 |
0.13 |
-0.34 |
51.40 |
48.60 |
16601 |
NC022452 |
|||||||||||
A. platyrhynchos |
22.21 |
32.81 |
29.20 |
15.78 |
0.14 |
-0.35 |
51.41 |
48.59 |
16604 |
NC009684 |
|||||||||||
A. poecilorhyncha |
22.21 |
32.82 |
29.17 |
15.81 |
0.14 |
-0.35 |
51.37 |
48.63 |
16608 |
NC022418 |
|||||||||||
Anser |
Anser albifrons |
22.63 |
32.05 |
30.15 |
15.18 |
0.14 |
-0.36 |
52.78 |
47.22 |
16737 |
NC004539 |
||||||||||
A. anser |
22.58 |
32.14 |
30.18 |
15.10 |
0.14 |
-0.36 |
52.77 |
47.23 |
16738 |
NC011196 |
|||||||||||
A. cygnoides |
22.49 |
32.24 |
30.21 |
15.06 |
0.15 |
-0.36 |
52.70 |
47.30 |
16739 |
NC023832 |
|||||||||||
A. fabalis |
22.73 |
31.84 |
30.07 |
15.36 |
0.14 |
-0.35 |
52.80 |
47.20 |
16688 |
NC016922 |
|||||||||||
A. indicus |
22.50 |
32.23 |
30.17 |
15.10 |
0.15 |
-0.36 |
52.67 |
47.33 |
16730 |
NC025654 |
|||||||||||
Asarcornis |
Asarcornis scutulata |
21.77 |
33.22 |
28.88 |
16.13 |
0.14 |
-0.35 |
50.64 |
49.36 |
16539 |
NC052827 |
||||||||||
Aythya |
Aythya americana |
22.24 |
32.75 |
29.39 |
15.62 |
0.14 |
-0.35 |
51.62 |
48.38 |
16616 |
NC000877 |
||||||||||
A. ferina |
22.19 |
32.83 |
29.42 |
15.56 |
0.14 |
-0.36 |
51.61 |
48.39 |
16616 |
NC024602 |
|||||||||||
A. fuligula |
22.18 |
32.87 |
29.42 |
15.52 |
0.14 |
-0.36 |
51.61 |
48.39 |
16616 |
NC024595 |
|||||||||||
A. nyroca |
22.28 |
32.71 |
29.66 |
15.35 |
0.14 |
-0.36 |
51.94 |
48.06 |
16623 |
MW287344 |
|||||||||||
Branta |
Branta bernicla |
22.68 |
31.99 |
29.91 |
15.42 |
0.14 |
-0.35 |
52.59 |
47.41 |
16747 |
KJ680301 |
||||||||||
B. canadensis |
22.60 |
32.07 |
30.18 |
15.14 |
0.14 |
-0.36 |
52.79 |
47.21 |
16760 |
NC007011 |
|||||||||||
Bucephala |
Bucephala albeola |
21.82 |
33.36 |
28.52 |
16.30 |
0.13 |
-0.34 |
50.34 |
49.66 |
16614 |
MW849286 |
||||||||||
B. clangula |
21.65 |
33.46 |
28.68 |
16.21 |
0.14 |
-0.35 |
50.33 |
49.67 |
16625 |
MW849283 |
|||||||||||
B. islandica |
21.65 |
33.43 |
28.68 |
16.24 |
0.14 |
-0.35 |
50.32 |
49.68 |
16627 |
MW849281 |
|||||||||||
Cairina |
Cairina moschata |
21.93 |
32.95 |
29.01 |
16.12 |
0.14 |
-0.34 |
50.93 |
49.07 |
16610 |
NC010965 |
||||||||||
Clangula |
Clangula hyemalis |
21.83 |
33.31 |
29.18 |
15.68 |
0.14 |
-0.36 |
51.01 |
48.99 |
16610 |
MW849278 |
||||||||||
Cygnus |
Cygnus atratus |
22.20 |
32.55 |
29.52 |
15.73 |
0.14 |
-0.35 |
51.72 |
48.28 |
16748 |
NC012843 |
||||||||||
C. columbianus |
22.79 |
31.89 |
30.09 |
15.23 |
0.14 |
-0.35 |
52.88 |
47.12 |
16728 |
NC007691 |
|||||||||||
C. cygnus |
22.79 |
31.97 |
29.97 |
15.27 |
0.14 |
-0.35 |
52.76 |
47.24 |
16724 |
NC027095 |
|||||||||||
C. olor |
22.59 |
32.12 |
29.64 |
15.65 |
0.13 |
-0.34 |
52.23 |
47.77 |
16739 |
NC027096 |
|||||||||||
Dendrocygna |
Dendrocygna javanica |
23.67 |
30.44 |
30.44 |
15.45 |
0.13 |
-0.33 |
54.11 |
45.89 |
16753 |
NC012844 |
||||||||||
Heteronetta |
Heteronetta atricapilla |
23.00 |
31.93 |
30.01 |
15.06 |
0.13 |
-0.36 |
53.00 |
47.00 |
16723 |
CM021836 |
||||||||||
Histrionicus |
Histrionicus histrionicus |
21.79 |
33.29 |
29.32 |
15.61 |
0.15 |
-0.36 |
51.11 |
48.89 |
16634 |
MW849288 |
||||||||||
Lophodytes |
Lophodytes cucullatus |
21.96 |
33.02 |
29.01 |
16.01 |
0.14 |
-0.35 |
50.97 |
49.03 |
16549 |
MW849287 |
||||||||||
Mareca |
Mareca falcata |
22.35 |
32.65 |
28.88 |
16.11 |
0.13 |
-0.34 |
51.24 |
48.76 |
16601 |
NC023352 |
||||||||||
M. penelope |
22.29 |
32.71 |
28.89 |
16.11 |
0.13 |
-0.34 |
51.18 |
48.82 |
16596 |
NC050973 |
|||||||||||
M. strepera |
22.19 |
32.81 |
28.84 |
16.16 |
0.13 |
-0.34 |
51.03 |
48.97 |
16600 |
NC045373 |
|||||||||||
Melanitta |
Melanitta deglandi |
21.58 |
33.45 |
28.74 |
16.22 |
0.14 |
-0.35 |
50.33 |
49.67 |
16452 |
MW849279 |
||||||||||
M. nigra |
21.77 |
33.32 |
28.83 |
16.08 |
0.14 |
-0.35 |
50.60 |
49.40 |
16552 |
MW849282 |
|||||||||||
M. perspicillata |
21.61 |
33.49 |
28.66 |
16.24 |
0.14 |
-0.35 |
50.27 |
49.73 |
16499 |
MW849280 |
|||||||||||
Mergus |
Mergus merganser |
21.88 |
33.24 |
28.73 |
16.16 |
0.14 |
-0.35 |
50.60 |
49.40 |
16630 |
NC040986 |
||||||||||
M. serrator |
21.72 |
33.38 |
28.77 |
16.13 |
0.14 |
-0.35 |
50.49 |
49.51 |
16627 |
MZ365040 |
|||||||||||
M. squamatus |
22.26 |
32.76 |
29.03 |
15.95 |
0.13 |
-0.35 |
51.29 |
48.71 |
16595 |
NC016723 |
|||||||||||
Table continues on next page.................... |
|||||||||||||||||||||
Family/Genus |
Species |
T |
C |
A |
G |
AT skew |
GC skew |
A+T |
G+C |
Total |
Accession number |
||||||||||
Netta |
Netta rufina |
22.44 |
32.59 |
29.61 |
15.37 |
0.14 |
-0.36 |
52.04 |
47.96 |
16625 |
NC024922 |
||||||||||
Nettapus |
Nettapus auritus |
22.83 |
32.07 |
29.79 |
15.31 |
0.13 |
-0.35 |
52.62 |
47.38 |
16643 |
CM021833 |
||||||||||
Oxyura |
Oxyura jamaicensis |
22.94 |
32.30 |
28.81 |
15.94 |
0.11 |
-0.34 |
51.75 |
48.25 |
15914 |
CM021834 |
||||||||||
Polysticta |
Polysticta stelleri |
21.89 |
33.04 |
28.79 |
16.28 |
0.14 |
-0.34 |
50.68 |
49.32 |
16612 |
MW849289 |
||||||||||
Sibirionetta |
Sibirionetta formosa |
22.51 |
32.44 |
29.52 |
15.53 |
0.13 |
-0.35 |
52.03 |
47.97 |
16592 |
NC015482 |
||||||||||
Somateria |
Somateria fischeri |
21.90 |
33.12 |
28.74 |
16.25 |
0.14 |
-0.34 |
50.64 |
49.36 |
16601 |
MW849290 |
||||||||||
S. mollissima |
21.97 |
33.02 |
28.62 |
16.39 |
0.13 |
-0.34 |
50.59 |
49.41 |
16626 |
MW849292 |
|||||||||||
S. spectabilis |
21.93 |
33.07 |
28.71 |
16.29 |
0.13 |
-0.34 |
50.65 |
49.35 |
16627 |
MW849291 |
|||||||||||
Spatula |
Spatula clypeata |
22.47 |
32.49 |
29.39 |
15.65 |
0.13 |
-0.35 |
51.86 |
48.14 |
16599 |
NC028346 |
||||||||||
Stictonetta |
Stictonetta naevosa |
22.39 |
32.80 |
28.45 |
16.35 |
0.12 |
-0.33 |
50.84 |
49.16 |
16778 |
CM021835 |
||||||||||
Tadorna |
Tadorna ferruginea |
21.99 |
32.96 |
29.57 |
15.48 |
0.15 |
-0.36 |
51.56 |
48.44 |
16639 |
NC024640 |
||||||||||
T. tadorna |
22.21 |
32.82 |
29.20 |
15.78 |
0.14 |
-0.35 |
51.40 |
48.60 |
16604 |
NC024750 |
|||||||||||
Anhimidae |
|||||||||||||||||||||
Chauna |
Chauna torquata |
24.01 |
30.82 |
31.31 |
13.87 |
0.13 |
-0.38 |
55.31 |
44.69 |
16766 |
NC052807 |
||||||||||
Anseranatidae |
|||||||||||||||||||||
Anseranas |
Anseranas semipalmata |
23.49 |
31.38 |
30.92 |
14.20 |
0.14 |
-0.38 |
54.41 |
45.59 |
16868 |
NC005933 |
||||||||||
Averages |
22.31 |
32.63 |
29.36 |
15.71 |
0.14 |
-0.35 |
51.67 |
48.33 |
16632 |
||||||||||||
Standard deviation |
0.51 |
0.64 |
0.64 |
0.52 |
0.01 |
0.01 |
1.10 |
1.10 |
127.12 |
Retrieved and downloaded the complete mitogenomes of 53 Anseriformes species from GenBank, and including 2 families and 46 genera (Table II). We used MEGA5.0 (Tamura et al., 2011) and DNASTAR (Burland, 2000) to compare the sequences of these birds and to analyze their base content. The calculation for the skewing of nucleotide composition was AT skew = (A-T) / (A + T), GC skew = (G-C) / (G + C) (Perna and Kocher, 1995).
Analysis of phylogeny
We collected the complete mitogenome of 53 Anseriformes species for phylogenetic analysis (Table II). Corresponding Gallus gallus, Hydrophasianus chirurgus, and Porzana fusca sequence were used as outgroups. Phylogenetic trees were estimated using ML and BI methods. The optimal parameter model selected by MrModel Test 3.06 (Nylander et al., 2004) and PAUP 4.0b10 software (Swofford, 2003). The GTR + I + G model that was selected as the best fit model for nucleotide phylogenetic analysis. The BI trees of Anseriformes birds were constructed using MrBayes 3.1.2 software (Ronquist et al., 2012). The ML tree was constructed using the MEGA5.0 software (Tamura et al., 2011), as follows: run the files in Phy format, select the best model using the Mr Model Test 3.06 software filter, set the ML in bootstrap mode, repeat bootstraps 50 times, and run the whole process 10,000 times to obtain the ML tree. Finally, Treeview32 (Saldanha, 2004) software was used to create and annotate the evolutionary tree diagram.
Analysis of adaptive molecular evolution
In order to test the possible influence of topology on the inference of selected locations, estimated ω (dN/dS) rate using the topology inferred from genomic data (Nadeau et al., 2007; Stoletzki and Eyre-Walker, 2011). dN is the non-synonymous substitution rate, or the rate at which changes in nucleotide sites lead to changes in new amino acid chains. Things are the opposite for dS (synonymous substitution). The dN/dS ratio provides evidence for selective restraint and pressures acting on PCGs (Shen et al., 2009; Botero-Castro et al., 2017). If natural selection is favoring amino acid changes, ω > 1 indicates positive selection. A ratio below 1 indicates purification selection, and a ratio at 1 indicates neutral evolution.
The phylogenetic tree has many branches, and test the role of selection on the Mergini branch. We employed CodeML branch models that estimate the selective pressure between Mergini tribes and other tribes, which allow the estimation of dN/dS of specific branches and clades of interest related to other clades and other parts of the tree (Zhang et al., 2005). In the dual scale and all branch site models described below, we label these species in the Mergini tribes as the foreground branches in a separate analysis. Use branch-site analysis to determine whether a portion of sites was subject to positive selection. To further evaluate these sites, we performed a positive selection analysis using the site model in PAML4. Selection Test was carried out using likelihood-ratio tests (LRT). The model was compared with a null model, in which the ω of the foreground partition class of the selected location was restricted to 1 (Zhang et al., 2005). LRT was used to identify the significance of selective pressure between different pedigrees.
RESULTS AND DISCUSSION
Genome organization and composition
M. serrator is a typical sea duck in the Mergini tribe and its complete mitogenome was found to be 16,627 bp (Table I). So far, the length of the mitogenome of birds that has been sequenced is from 16,300 bp to 23,500 bp (Gao et al., 2009; Lei, 2015). Among Anseriformes species, the shortest mitogenome belongs to Oxyura jamaicensis (15,914 bp) and the longest is that of Anseranas semipalmata (16,868 bp). In the Anseriformes taxa, the length and gene arrangement of the mitogenome are quite conserved. The genomes of M. serrator contained a typical set of 37 genes (Fig. 1), including 13 PCGs, 2 ribosomal RNAs, 22 transfer RNAs, and a non-coding putative CR (Table I, Fig. 2). The structure and gene arrangement were basically the same as most birds (Gao et al., 2009; Lei, 2015). An ND6 gene and eight tRNA genes (tRNAGln, tRNAAla, tRNAAsn, tRNACys, tRNATyr, tRNASer (UCN), tRNAPro, and tRNAGlu) were encoded by the L-strand, whereas other genes were encoded by the H-strand (Table I).
The overall information about the mitogenome of M. serrator was as follows: T: 21.72%, C: 33.38%, A: 28.77%, and G: 16.13% and the content of GC is 49.51% with A+T-rich feature, same as in other species of Anseriformes. The mean base composition of Anseriformes is T (22.31 ± 0.51%), C (32.63 ± 0.64%), A (29.36 ± 0.64%), and G (15.71 ± 0.52%). Additionally, in our study, the sequence of base richness was mainly C > A > T > G, and only Chauna torquata had base contents of A > C >T > G. The AT content range is 50.27% ~ 55.31%, with an average value is 51.67 ± 1.10%, which is higher than the GC content and similar to that of birds in general (50.5% to 57.7%) (Gao et al., 2009).
Vertebrate mtDNA usually shows significant chain bias in nucleotide composition, and this chain bias can be measured as AT and GC skew. We measured the AT skew (0.14) and GC skew (-0.35) of the H-strand in the M. serrator mitogenome. The mean AT skew was 0.14 ± 0.01, and the GC skew was -0.35 ± 0.01. It shows that the nucleotide composition of the complete mitogenome of Anseriformes has a slight specific bias against A and C.
Analysis of PCGs and codon usage
The geneme of M. serrator had 13 PCGs, and 12 PCGs were encoded on the H-strand except ND6. The overall length of 13 PCGs was 11,229 bp, which accounted for 67.53% of the whole mitogenome (Table I). The longest was the ND5 gene (1,824 bp) and the shortest was the ATP8 gene (168 bp). Among the 13 PCGs, the start codon of 9 PCGs was ATG, and the nonstandard start codon GTG was found in COI, COII, and ND5 genes. In addition, the start codon of ND1 gene was ATA (Table I). The most common stop codon was TAA, which has been used six times. TAG and AGG were used twice as stop codons, and the terminal T probably served as the stop signal in the ND2, ND4, and COIII genes. In the ND3 gene, the base A of the stop codon was not translated.
The RSCU values of 13 PCGs are shown in Figure 3. In addition, the start and stop codons, the mitogenome consists of 3662 codons. The most commonly used amino acid was leucine (13.63%) in the use of 13 PCGs of M. serrator, followed by proline (11.47%), serine (9.59%), and threonine (8.06%). In three Mergini species, eight amino acids (Pro, Thr, Leu1, Arg, Ala, Ser1, Val, and Gly) were encoded by four different codons, and the others were encoded by two different codons. The overall AT-skew and GC-skew of the 13 PCGs was 0.06 and -0.36, respectively. The base composition of PCGs has a higher A + C contents, and a preferred to use A or C nucleotides in the first codon position (Fig. 3). This phenomenon also occurs in other species of Mergini, especially M. squamatus and M. merganser.
Ribosomal RNAs, transfer RNAs, and CR
In the M. serrator mitogenome, ribosomal RNA genes include 12S rRNA and 16S rRNA, which are 995 bp and 1600 bp in length, respectively. They are located between the tRNAPhe and tRNALeu genes, and separated by the tRNAVal, a feature often found in the mitogenomes of birds (Cheng et al., 2011; Jin et al., 2012). Two ribosomal RNAs were identified on the H-strand. The A + T content of both rRNA genes were collectively 52.5% and they had a positive AT-skew (0.28).
A total of 22 tRNAs were found in the mitogenome, and they were scattered between rRNAs and PCGs. Their sizes is between 64 bp and 76 bp (Table I, Supplementary Fig. 1). Three tRNA clusters, namely, IQM (tRNAIle tRNAGln tRNAMet), WANCY (tRNATrp tRNAAla tRNAAsn tRNACys tRNATyr), and HSL (tRNAHis tRNASer (AGY) tRNALeu (CUN)), were determined. Sequences of the tRNA genes could be folded into a canonical cloverleaf secondary structure.
The CR of mtDNAs of M. serrator was 1,067 bp, which was located between the tRNAGlu and tRNAPhe genes. The CR could be divided into three domains according to the base composition (5′-peripheral domain, Domain I; central conserved domain, Domain II; and 3′-peripheral domains, Domain III). In Domain I, extended termination-associated sequences 1 (ETAS1) was observed at position nt 77 - 98 (Fig. 2). A sequence block was also found in domain I (nt 298–323), which is similar to the conserved sequence block (CSB1-like). Four conserved sequence boxes (C, D, E, and F) were identified in Domain II. CBS1 (nt 407–589) and the H- and L-strand transcription promoter sites (nt 806–822) occurred in Domain III.
Phylogenetic analysis
The complete sequences should be suitable for the phylogenetic analysis and research of Anseriformes. Reconstruction of BI and ML trees using complete mitogenome sequences of 56 species. Use the best-fit model of GTR+G+I (-lnL = 258460.5951, AIC= 516937.2376) and calculated with the Modeltest and Mrmodeltest programs for analysis. Since the generated BI and ML trees show the same topological structure, only the BI tree is displayed (the numbers on the branches represent the bootstrap support for BI and ML trees).
Traditional morphological analysis divides Anatinae into Tadornini, Tachyerini, Cairinini, Anatini, Aythyini, Mergini, Oxyurini, and others. In the literature, Anseriformes species are relatively clearly classified at the level of the family but not the subfamily and genus (Liu et al., 2012; García et al., 2014; Hu et al., 2017). The replacement rate and mutation rate of the mitochondrial sequence of Anseriformes species are high (Fain et al., 2007; Yang et al., 2010). The BI and ML trees showed the complete mitogenome has the ability to distinguish 53 Anseriformes species. The phylogenetic analysis based on the complete mitogenome strongly supports the monophyleticity of the order Anseriformes. The analyses indicated that the Anseriformes include the families Anseranatidae, Anhimidae, and Anatidae. Anseranas semipalmata was included in the branch of Anseranatidae, and Chauna torquata was included in the branch of Anhimidae, and other species belong to the branch of Anatidae (Fig. 4). The branch of Anatidae included Dendrocygninae, Oxyurinae, Stictonettinae, Anserinae, Tadorninae, and Anatinae. The Anatinae had the largest number (32) of species. The Anatinae included Aythyini, Anatini, Somaterini, and Mergini. The results suggest that M. serrator is sister to M. merganser and M. squamatus, which is slightly different from a previous study (Liu et al., 2014). However, both studies suggest that this group belongs to the Mergini.
Selection analysis
Most mitochondrial genes besides the ND6 gene are conserved and evolved under purifying selection (Fig. 5). The branch model detected significant signals of positive selection for the ND6 gene tested in all tribes of Anseriformes. It also fitted better for most PCGs when the Mergini tribe was labeled as a foreground branch (Fig. 5) (LRT p value < 0.05). In addition, the ND6 gene showed the highest ω values (ω > 1) both for the foreground (ω = 4.59) and background (ω = 3.28) branches, compared to other genes. The ω value of the ND3 gene in the Mergini tribe is lower than those for other tribes. The ω values of other genes in the Mergini tribe are higher than those for other tribes. The free-ratio model fitting the data of 13 PCGs is significantly better than the null hypothesis (one-ratio model), indicating that there are different evolution rates among the tribes included in our data set.
CONCLUSION
The mitogenome (16,627 bp) was larger than those of most Anseriformes species. A positive AT skew consistent with other Anseriformes species was found, which implies a slightly specific bias towards A and C in the Anseriformes. PCGs accounted for 67.53% of the whole mitogenome and most of these genes were encoded by the L-strand. ATG and TAA were the most frequent start codon and stop codon, respectively. We summarized the RSCU values for the 13 PCGs and 3788 codons and high A+C contents were found, which also occurred in other species of Mergini. Two ribosomal RNA genes, 12S rRNA and 16S rRNA, were encoded on the H-strand and had a positive A-T skew. Sequences of 22 tRNA genes ranged from 64 bp to 76 bp and could be folded into a canonical cloverleaf secondary structure. The phylogenetic relationships were analyzed and M. serrator was sister to M. merganser and M. squamatus, and formed a closely evolved Mergini clade. Compared to other species of Anseriformes, most PCGs evolved under purifying selection except for the ND6 gene. Different evolutionary rates between the Mergini tribe and other tribes were found.
ACKNOWLEDGEMENTS
We thank researcher Yu Dandan for her help in writing this article. This study was supported by the Basic Research Business Project of National Scientific Institute (GYZX210405) and the National Science Foundation of China under Grant (32001111).
There is supplementary material associated with this article. Access the material online at: https://dx.doi.org/10.17582/journal.pjz/20220314030351
Conflict of interest statement
The authors have declared no conflict of interest.
REFERENCES
Allendorf, F.W., Hohenlohe, P.A. and Luikart, G., 2010. Genomics and the future of conservation genetics. Nat. Rev. Genet., 11: 697-709. https://doi.org/10.1038/nrg2844
Barker, F.K., Be-Nesh, M.K., Vandergon, A.J. and Lanyon, S.M., 2012. Contrasting evolutionary dynamics and information content of the avian mitochondrial control region and ND2 gene. PLoS One, 7: e46403. https://doi.org/10.1371/journal.pone.0046403
Bernt, M., Donath, A., Juhling, F., Externbrink, F., Florentz, C., Fritzsch, G., Putz, J., Middendorf, M. and Stadler, P.F., 2013. MITOS: Improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol., 69: 313-319. https://doi.org/10.1016/j.ympev.2012.08.023
Botero-Castro, F., Figuet, E., Tilak, M.K., Nabholz, B. and Galtier, N., 2017. Avian genomes revisited: Hidden genes uncovered and the rates versus traits paradox in birds. Mol. Biol. Evol., 34: 3123-3131. https://doi.org/10.1016/j.ympev.2012.08.023
Burland, T.G., 2000. Dnastar’s Lasergene sequence analysis software. Method Mol. Boil., (Clifton, N.J.), 132: 71-91. https://doi.org/10.1385/1-59259-192-2:71
Cadahia, L., Pinsker, W., Negro, J.J., Pavlicev, M., Urios, V. and Haring, E., 2009. Repeated sequence homogenization between the control and pseudo-control regions in the mitochondrial genomes of the subfamily Aquilinae. J. exp. Zool. B Mol. Dev. Evol., 312B: 171-185. https://doi.org/10.1002/jez.b.21282
Cheng, Y.Z., Wang, R.X., amd Xu, T.J., 2011. The mitochondrial genome of the spinyhead croaker Collichthys lucida: Genome organization and phylogenetic consideration. Mar. Genom., 4: 17-23. https://doi.org/10.1016/j.margen.2010.12.001
Debela, M.T., Wu, Q., Li, Z., Sun, X. and Li, Y., 2021. Habitat suitability assessment of wintering herbivorous Anseriformes in Poyang Lake, China. Diversity (Basel), 13: 171. https://doi.org/10.3390/d13040171
Dierckxsens, N., Mardulyn, P. and Smits, G., 2017. NOVOPlasty: De novo assembly of organelle genomes from whole genome data. Nucl. Acids Res., 45: e18. https://doi.org/10.1093/nar/gkw955
Donne-Goussé, Carole, Laudet, Vincent and Hänni, Catherine, 2002. A molecular phylogeny of anseriformes based on mitochondrial DNA analysis. Mol. Phylogenet. Evol., 23: 339-356. https://doi.org/10.1016/S1055-7903(02)00019-2
Fain, M.G., Krajewski, C. and Houde, P., 2007. Phylogeny of core Gruiformes (Aves: Grues) and resolution of the Limpkin–Sungrebe problem. Mol. Phylogenet. Evol., 43: 515-529. https://doi.org/10.1016/j.ympev.2007.02.015
Gao, Y., Miao, Y., Su, X., Chi, Z., Yu, B., and Jiang, F., 2009. A comprehensive analysis on 74 avian mitochondrial genome base compositions. J. Yunnan Agric. Univ., 24: 51-58. https://doi.org/10.1360/972009-1142
García, R.J.C., Gibb, G.C., and Trewick, S.A., 2014. Eocene diversification of crown group rails (Aves, Gruiformes, Rallidae). PLoS One, 9: e109635-e109635. https://doi.org/10.1371/journal.pone.0109635
Gonzalez, J., Düttmann, H., and Wink, M., 2009. Phylogenetic relationships based on two mitochondrial genes and hybridization patterns in Anatidae. J. Zool., 279: 310-318. https://doi.org/10.1111/j.1469-7998.2009.00622.x
Hu, C.C., Zhang, C.L., Sun, L., Zhang, Y., Xie, W.L., Zhang, B.W. and Chang, Q., 2017. The mitochondrial genome of pin-tailed snipe Gallinago stenura, and its implications for the phylogeny of Charadriiformes. PLoS One, 12: e0175244-e0175244. https://doi.org/10.1371/journal.pone.0175244
Janssen, R.J., Van Den Heuvel, L.P., and Smeitink, J.A., 2004. Genetic defects in the oxidative phosphorylation (OXPHOS) system. Exp. Rev. mol. Diagn., 4: 143-156. https://doi.org/10.1586/14737159.4.2.143
Jetz, W., Thomas, G.H., Joy, J.B., Hartmann, K., and Mooers, A.O., 2012. The global diversity of birds in space and time. Nature, 491: 444-448. https://doi.org/10.1038/nature11631
Jiang, F., Miao, Y., Liang, W., Ye, H., Liu, H. and Liu, B., 2009. The complete mitochondrial genomes of the whistling duck (Dendrocygna javanica) and black swan (Cygnus atratus): dating evolutionary divergence in Galloanserae. Mol. Boil. Rep., 37: 3001-3015. https://doi.org/10.1360/972009-1142
Jin, X., Wang, R., Xu, T. and Shi, G., 2012. Complete mitochondrial genome of Oxuderces dentatus (Perciformes, Gobioidei). Mitochondrial DNA, 23: 142-144. https://doi.org/10.3109/19401736.2012.660930
Johnsgard, P.A., 1966. Handbook of waterfowl behavior. Handb. Waterf. Behavior, 37: 860. https://doi.org/10.2307/4511272
Johnsgard, P.A., 2010. Ducks, geese, and swans of the world: Tribe Mergini (Sea Ducks). Paul johnsgard.
Lavretsky, P., Wilson, Robert, E., Talbot, Sandra, L. and Sonsthagen, S.A., 2021. Phylogenomics reveals ancient and contemporary gene flow contributing to the evolutionary history of sea ducks (Tribe Mergini). Mol. Phylogenet. Evol., 161: 107164. https://doi.org/10.1016/j.ympev.2021.107164
Lecroy, M. and Barker, F.K., 2006. A new species of bush-warbler from Bougainville Island and a monophyletic origin for Southwest Pacific Cettia. Am. Mus. Novit., 3511: 1-20. https://doi.org/10.1206/0003-0082(2006)3511[1:ANSOBF]2.0.CO;2
Lei, Z., 2015. Advances on the mitochondrial genome and phylogeny of Aves. J. Wuhu Inst. Technol.,
Liu, G., Zhou, L.Z., and Gu, C.M., 2012. Complete sequence and gene organization of the mitochondrial genome of scaly-sided merganser (Mergus squamatus) and phylogeny of some Anatidae species. Mol. Biol. Rep., 39: 2139-2145. https://doi.org/10.1007/s11033-011-0961-5
Liu, G., Zhou, L., Li, B., and Zhang, L., 2014. The complete mitochondrial genome of Aix galericulata and Tadorna ferruginea: Bearings on their phylogenetic position in the Anseriformes. PLoS One, 9: e109701. https://doi.org/10.1371/journal.pone.0109701
Liu, W., Hu, C., Xie, W., Chen, P., Zhang, Y., Yao, R., Li, K., and Chang, Q., 2018. The mitochondrial genome of red-necked phalarope Phalaropus lobatus (Charadriiformes: Scolopacidae) and phylogeny analysis among Scolopacidae. Genes Genomics, 40: 455-463. https://doi.org/10.1007/s13258-017-0632-6
Livezey, B.C., 1986. A phylogenetic analysis of recent anseriform genera using morphological characters. Auk, 103: 737-754. https://doi.org/10.1371/journal.pone.0184529
Livezey, B.C., 1997. A phylogenetic classification of waterfowl (Aves, Anseriformes), including selected fossil species. Annls Carnegie Mus., 66: 457-496. Available: https://www.biodiversitylibrary.org/part/215141
Luo, Y., Yang, X., and Gao, Y., 2013. Mitochondrial DNA response to high altitude: a new perspective on high-altitude adaptation. Mitochondrial DNA, 24: 313-319. https://doi.org/10.3109/19401736.2012.760558
Nabholz, B., Glémin, S., and Galtier, N., 2008. Extreme variation of mtDNA neutral substitution rate across mammalian species the longevity hypothesis. Mol. Biol. Evol., 25: 795-795. https://doi.org/10.1186/1471-2148-8-20
Nadeau, N.J., Burke, T., and Mundy, N.I., 2007. Evolution of an avian pigmentation gene correlates with a measure of sexual selection. Proc. R. Soc. B Biol. Sci., 274: 1807-1813. Available: http://www.jstor.org/stable/25249253, https://doi.org/10.1098/rspb.2007.0174
Nylander, J.A., Ronquist, F., Huelsenbeck, J.P., and Nieves-Aldrey, J., 2004. Bayesian phylogenetic analysis of combined data. Syst. Boil., 53: 47-67. https://doi.org/10.1080/10635150490264699
Peng, C., Yuqing, H., Chaoying, Z., Bin, G., and Luzhang, R., 2017. Complete mitochondrial genome of Porzana fusca and Porzana pusilla and phylogenetic relationship of 16 Rallidae species. Genetica. https://doi.org/10.1371/journal.pone.0109635
Perna, N.T. and Kocher, T.D., 1995. Unequal base frequencies and the estimation of substitution rates. Mol. Biol. Evol., 2: 359-361. https://doi.org/10.1093/oxfordjournals.molbev.a040211
Peters, J.L., Bolender, K.A., and Pearce, J.M., 2012. Behavioural vs. molecular sources of conflict between nuclear and mitochondrial DNA: the role of male-biased dispersal in a Holarctic sea duck. Mol. Ecol., 21: 3562-3575. https://doi.org/10.1111/j.1365-294X.2012.05612.x
Ponomarenko, O., Banik, M., and Zhukov, O., 2021. Assessing habitat suitability for the common pochard, (Anseriformes, Anatidae) at different spatial scales in Orel’River Valley, Ukraine. Ekológia (Bratislava), 40: 154-162. https://doi.org/10.2478/eko-2021-0018
Ren, Q., Yuan, J., Ren, L., Zhang, L., Zhang, L., Jiang, L., Chen, D., Kan, X., and Zhang, B., 2014. The complete mitochondrial genome of the yellow-browed bunting, Emberiza chrysophrys (Passeriformes: Emberizidae), and phylogenetic relationships within the genus Emberiza. J. Genet., 93: 699-707. https://doi.org/10.1007/s12041-014-0428-2
Ringelman, Kevin, and M., 2015. Ducks, geese, and swans of North America. J. Field Ornithol. J. Ornithol. Invest., https://doi.org/10.1093/sysbio/sys029
Ronquist, F., Teslenko, M., Van Der Mark, P., Ayres, D.L., Darling, A., Höhna, S., Larget, B., Liu, L., Suchard, M.A., and Huelsenbeck, J.P., 2012. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol., 61: 539-542. https://doi.org/10.1093/sysbio/sys029
Roques, S., Godoy, J.A., Negro, J.J., and Hiraldo, F., 2004. Organization and variation of the mitochondrial control region in two vulture species, Gypaetus barbatus and Neophron percnopterus. J. Hered., 95: 332-337. https://doi.org/10.1093/jhered/esh047
Ruan, L., Xu, W., Han, Y., Zhu, C., Guan, B., Xu, C., Gao, B., and Zhao, D., 2018. Gene flow from multiple sources maintains high genetic diversity and stable population history of Common Moorhen Gallinula chloropus in China. Ibis, 160: 855-869. https://doi.org/10.1111/ibi.12579
Saldanha, A.J., 2004. Java Treeview-extensible visualization of microarray data. Bioinformatics, 20: 3246-3248. https://doi.org/10.1002/acp.1211
Saraste, M., 1999. Oxidative phosphorylation at the fin de siecle. Science, 283: 1488-1493. https://doi.org/10.1002/acp.1211
Savard, J.P.L., Reed, A., and Lesage, L., 1998. Brood amalgamation in surf scoters Melanitta perspicillata and other Mergini. Wildfowl, 49: 129-138. https://doi.org/10.2173/tbna.363.p
Scribner, K.T., Petersen, M.R., Fields, R.L., Talbot, S.L., Pearce, J.M., and Chesser, R.K., 2001. Sex-biased gene flow in spectacled eiders (Anatidae): Inferences from molecular markers with contrasting modes of inheritance. Evolution, 55: 2105-2115. https://doi.org/10.1111/j.0014-3820.2001.tb01325.x
Shen, Y.Y., Shi, P., Sun, Y.B., and Zhang, Y.P., 2009. Relaxation of selective constraints on avian mitochondrial DNA following the degeneration of flight ability. Genome Res., 19: 1760-1765. https://doi.org/10.1101/gr.093138.109
Sonsthagen, S.A., Talbot, S.L., Scribner, K.T., and McCracken, K.G., 2011. Multilocus phylogeography and population structure of common eiders breeding in North America and Scandinavia. J. Biogeogr., 38: 1368-1380. https://doi.org/10.1111/j.1365-2699.2011.02492.x
Sruoga, A., Slavenaite, S., Butkauskas, D., and Grazulevicius, G., 2008. Cross species applicability of microsatellite markers for investigation of sea ducks (Mergini) genetic differentiation. Proc. Latvian Acad. Sci. De Gruyter Poland. https://doi.org/10.1111/j.1365-2699.2011.02492.x
Starkov, A.A., 2008. The role of mitochondria in reactive oxygen species metabolism and signaling. Annls N. Y. Acad. Sci., 1147: 37-52. https://doi.org/10.1196/annals.1427.015
Stoletzki, N. and Eyre-Walker, A., 2011. The positive correlation between d N/d S and d S in mammals is due to runs of adjacent substitutions. Mol. Biol. Evol., 28: 1371-1380. https://doi.org/10.1093/molbev/msq320
Sun, Z., P. Tao, C. Hu, S. Lu, and B. Zhang, 2017. Rapid and recent diversification patterns in Anseriformes birds: Inferred from molecular phylogeny and diversification analyses. PLoS One, 12: e0184529. https://doi.org/10.1371/journal.pone.0184529
Swofford, D. L., 2003. PAUP*. Phylogenetic analysis using parsimony (*and other methods) version 4. Sinauer, Sunderland, Massachusetts, USA. Nat. Biotechnol., 18: 233-234.
Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S., 2011. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol., 28: 2731-2739. https://doi.org/10.1093/molbev/msr121
Tang, Y., Zhang, X.Y., Ma, Y.Y. and Zheng, X.D., 2021. Descriptive study of the mitogenome of the diamondback squid (Thysanoteuthis rhombus Troschel, 1857) and the evolution of mitogenome arrangement in oceanic squids. J. Zool. Syst. Evol. Res., 59: 981-991. https://doi.org/10.1111/jzs.12478
Thompson, J.D., Gibson, T.J., Plewniak, F., Jeanmougin, F., and Higgins, D.G., 1997. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl. Acids Res., 25: 4876-4882. https://doi.org/10.1093/nar/25.24.4876
Torroni, A., Rengo, C., Guida, V., Cruciani, F., Sellitto, D., Coppa, A., Calderon, F.L., Simionati, B., Valle, G., and Richards, M., 2001. Do the four clades of the mtDNA haplogroup L2 evolve at different rates? Am. J. Hum. Genet., 69: 1348-1356. https://doi.org/10.1086/324511
Wolstenholme, D.R., 1992. Genetic novelties in mitochondrial genomes of multicellular animals. Curr. Opin. Genet. Dev., 2: 918-925. https://doi.org/10.1016/S0959-437X(05)80116-9
Yang, R., Wu, X., Yan, P., Su, X., and Yang, B., 2010. Complete mitochondrial genome of Otis tarda (Gruiformes: Otididae) and phylogeny of Gruiformes inferred from mitochondrial DNA sequences. Mol. Biol. Rep., 37: 3057-3066. https://doi.org/10.1007/s11033-009-9878-7
Zhang, J., Nielsen, R., and Yang, Z., 2005. Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol. Biol. Evol., 22: 2472-2479. https://doi.org/10.1093/molbev/msi237
Zhu, B.R., Verkuil, Y.I., Conklin, J.R., Yang, A., Lei, W., Alves, J.A., Hassell, C.J., Dorofeev, D., Zhang, Z., and Piersma, T., 2021. Discovery of a morphologically and genetically distinct population of black-tailed Godwits in the East Asian-Australasian flyway. Ibis, 163: 448-462. https://doi.org/10.1111/ibi.12890
To share on other social networks, click on any share button. What are these?