Historical biogeography of tits (Aves: Paridae, Remizidae)

Positive selection and convergent evolution shape molecular phenotypic traits of innate immunity receptors in tits (Paridae)

Despite widespread variability and redundancy abounding animal immunity, little is currently known about the rate of evolutionary convergence (functionally analogous traits not inherited from a common ancestor) in host molecular adaptations to parasite selective pressures. Toll-like receptors (TLRs) provide the molecular interface allowing hosts to recognize pathogenic structures and trigger early danger signals initiating an immune response. Using a novel combination of bioinformatic approaches, here we explore genetic variation in ligand-binding regions of bacteria-sensing TLR4 and TLR5 in 29 species belonging to the tit family of passerine birds (Aves: Paridae). Three out of the four consensual positively selected sites in TLR4 and six out of 14 positively selected positions in TLR5 were located on the receptor surface near the functionally important sites, and based on the phylogenetic pattern evolved in a convergent (parallel) manner. This type of evolution was also seen at one N-glycosylation site and two positively selected phosphorylation sites, providing the first evidence of convergence in post-translational modifications in evolutionary immunology. Finally, the overall mismatch between phylogeny and the clustering of surface charge distribution demonstrates that convergence is common in overall TLR4 and TLR5 molecular phenotypes involved in ligand binding. Our analysis did not reveal any broad ecological traits explaining the convergence observed in electrostatic potentials, suggesting that information on microbial symbionts may be needed to explain TLR evolution. Adopting state-of-the-art predictive structural bionformatics, we have outlined a new broadly applicable methodological approach to estimate the functional significance of positively selected variation and test for the adaptive molecular convergence in protein-coding polymorphisms.

Elevational divergence in the great tit complex revealed by major hemoglobin genes

Gene flow and demographic history can play important roles in the adaptive genetic differentiation of species, which is rarely understood in the high-altitude adaptive evolution of birds. To elucidate genetic divergence of populations in the great tit complex (Parus major, P. minor and P. cinereus) at different elevations, we compared the genetic structure and gene flow in hemoglobin genes with neutral loci. Our results revealed the elevationally divergent structure of α^A-globin gene, distinctive from that of the β^A-globin gene and neutral loci. We further investigated gene flow patterns among the populations in the central-northern (> 1,000 m a.s.l.), south-eastern (< 1,000 m a.s.l.) regions and the Southwest Mountains (> 2,000 m a.s.l.) in China. The high-altitude (> 1,000 m a.s.l.) diverged α^A-globin genetic structure coincided with higher α^A-globin gene flow between highland populations, in contrast to restricted neutral gene flow concordant with the phylogeny. The higher α^A-globin gene flow suggests the possibility of adaptive evolution during population divergence, contrary to the lower α^A-globin gene flow homogenized by neutral loci during population expansion. In concordance with patterns of historical gene flow, genotypic and allelic profiles provide distinctive patterns of fixation in different high-altitude populations. The fixation of alleles at contrasting elevations may primarily due to highland standing variants α^A49Asn/72Asn/108Ala originating from the south-western population. Our findings demonstrate a pattern of genetic divergence with gene flow in major hemoglobin genes depending on population demographic history.

Divergent and parallel routes of biochemical adaptation in high-altitude passerine birds from the Qinghai-Tibet Plateau

When different species experience similar selection pressures, the probability of evolving similar adaptive solutions may be influenced by legacies of evolutionary history, such as lineage-specific changes in genetic background. Here we test for adaptive convergence in hemoglobin (Hb) function among high-altitude passerine birds that are native to the Qinghai-Tibet Plateau, and we examine whether convergent increases in Hb-O₂ affinity have a similar molecular basis in different species. We documented that high-altitude parid and aegithalid species from the Qinghai-Tibet Plateau have evolved derived increases in Hb-O₂ affinity in comparison with their closest lowland relatives in East Asia. However, convergent increases in Hb-O₂ affinity and convergence in underlying functional mechanisms were seldom attributable to the same amino acid substitutions in different species. Using ancestral protein resurrection and site-directed mutagenesis, we experimentally confirmed two cases in which parallel substitutions contributed to convergent increases in Hb-O₂ affinity in codistributed high-altitude species. In one case involving the ground tit (Parus humilis) and gray-crested tit (Lophophanes dichrous), parallel amino acid replacements with affinity-enhancing effects were attributable to nonsynonymous substitutions at a CpG dinucleotide, suggesting a possible role for mutation bias in promoting recurrent changes at the same site. Overall, most altitude-related changes in Hb function were caused by divergent amino acid substitutions, and a select few were caused by parallel substitutions that produced similar phenotypic effects on the divergent genetic backgrounds of different species.

The reanalysis of biogeography of the Asian tree frog, Rhacophorus (Anura: Rhacophoridae): geographic shifts and climatic change influenced the dispersal process and diversification

Rapid uplifts of the Tibetan Plateau and climate change in Asia are thought to have profoundly modulated the diversification of most of the species distributed throughout Asia. The ranoid tree frog genus Rhacophorus, the largest genus in the Rhacophoridae, is widely distributed in Asia and especially speciose in the areas south and east of the Tibetan Plateau. Here, we infer phylogenetic relationships among species and estimate divergence times, asking whether the spatiotemporal characteristics of diversification within Rhacophorus were related to rapid uplifts of the Tibetan Plateau and concomitant climate change. Phylogenetic analysis recovered distinct lineage structures in Rhacophorus, which indicated a clear distribution pattern from Southeast Asia toward East Asia and India. Molecular dating suggests that the first split within the genus date back to the Middle Oligocene (approx. 30 Ma). The Rhacophorus lineage through time (LTT) showed that there were periods of increased speciation rate: 14–12 Ma and 10–4 Ma. In addition, ancestral area reconstructions supported Southeast Asia as the ancestral area of Rhacophorus. According to the results of molecular dating, ancestral area reconstructions and LTT we think the geographic shifts, the staged rapid rises of the Tibetan Plateau with parallel climatic changes and reinforcement of the Asian monsoons (15 Ma, 8 Ma and 4–3 Ma), possibly prompted a burst of diversification in Rhacophorus.

Improved sampling at the subspecies level solves a taxonomic dilemma - A case study of two enigmatic Chinese tit species (Aves, Passeriformes, Paridae, Poecile)

A recent full species-level phylogeny of tits, titmice and chickadees (Paridae) has placed the Chinese endemic black-bibbed tit (Poecile hypermelaenus) as the sister to the Palearctic willow tit (P. montanus). Because this sister-group relationship is in striking disagreement with the traditional affiliation of P. hypermelaenus close to the marsh tit (P. palustris) we tested this phylogenetic hypothesis in a multi-locus analysis with an extended taxon sampling including sixteen subspecies of willow tits and marsh tits. As a taxonomic reference we included type specimens in our analysis. The molecular genetic study was complemented with an analysis of biometric data obtained from museum specimens. Our phylogenetic reconstructions, including a comparison of all GenBank data available for our target species, clearly show that the genetic lineage previously identified as P. hypermelaenus actually refers to P. weigoldicus because sequences were identical to that of a syntype of this taxon. The close relationship of P. weigoldicus and P. montanus - despite large genetic distances between the two taxa - is in accordance with current taxonomy and systematics. In disagreement with the previous phylogenetic hypothesis but in accordance with most taxonomic authorities, all our P. hypermelaenus specimens fell in the sister clade of all western and eastern Palearctic P. palustris. Though shared haplotypes among the Chinese populations of the two latter species might indicate mitochondrial introgression in this part of the breeding range, further research is needed here due to the limitations of our own sampling.

Evolution of body morphology and beak shape revealed by a morphometric analysis of 14 Paridae species

BACKGROUND

Morphological characters of birds reflect their adaptive evolution and ecological requirements and are also relevant to phylogenetic relationships within a group of related species. The tits (Paridae) are known to be outwardly homogeneous in shape, with one aberrant member, the Ground Tit (Pseudopodoces humilis), which is quite different from its relatives in both body morphology and beak shape. We combined traditional measurements and geometric morphometrics to quantify the variation in body morphology and beak shape of 14 Paridae species distributed in China. Based on these results, we sought to assess the contribution of phylogeny, altitude and species interactions to the evolution of morphological traits.

RESULTS

The basic features for discriminating among the 14 species studied here were overall body size, the ratio of body and tail length to culmen and tarsus length, and beak shape (long/slender/pointy vs. short/robust/blunt). These dimensions clearly separate Ps. humilis and Melanochlora sultanea from the other species in shape space. Body length and PC3 of beak shape (round outline vs. straight outline) show significant phylogenetic signals. Across 14 species, altitude is related to tarsus, culmen length and PC1 of beak shape. Within Parus major, altitude is related to body weight, body length, culmen length and PC1 of body morphology. Morphological distances and geographic distances among species are positively correlated.

CONCLUSIONS

The body morphology of Paridae species shows extensive evolutionary changes, while their beak has mainly evolved along the long/slender/pointy vs. short/robust/blunt dimension. Only body length and beak curvature show a phylogenetic signal. Altitude correlates with multiple traits both across and within species, suggesting that altitude is an important factor in promoting morphological divergence. The deviant appearance of Ps. humilis corresponds to its foraging and feeding adaptations to high-altitude steppe habitats. Our results also show a higher level of morphological divergence with greater difference in distribution ranges among the Paridae species involved in this study.

Extreme variation in patterns of tandem repeats in mitochondrial control region of yellow-browed tits (Sylviparus modestus, Paridae)

To investigate the evolutionary pattern and origins of tandem repeats in the mitochondrial control region of the yellow-browed tit (Sylviparus modestus), the control region and another four mitochondrial loci from fifteen individuals were analyzed. A 117-bp tandem repeat unit that repeated once, twice or three times in different individuals was found, and a rarely reported arrangement for this tandem repeats region that a 5' imperfect copy at its downstream and a 3' imperfect copy at its upstream was observed. The haplotype network, phylogenetic trees, and ancestral state reconstruction of the combined dataset of five loci suggested multiple origins of the same repeat number. The turnover model via slipped-strand mispairing was introduced to interpret the results, because mispairing occurred so frequently that multiple origins of certain repeat number were observed. Insertion via recombination should be a better explanation for the origin of this tandem repeat unit, considering characteristics of the combined sequence of the 3' and 5' imperfect copy, including identification of its homolog in other passerines and its predicted secondary structure.