Multilocus phylogeny of the avian family Alaudidae (larks) reveals complex morphological evolution, non-monophyletic genera and hidden species diversity

Selected Lark Mitochondrial Genomes Provide Insights into the Evolution of Second Control Region with Tandem Repeats in Alaudidae (Aves, Passeriformes)

The control region (CR) regulates the replication and transcription of the mitochondrial genome (mitogenome). Some avian mitogenomes possess two CRs, and the second control region (CR2) may enhance replication and transcription; however, the CR2 in lark mitogenome appears to be undergoing loss and is accompanied by tandem repeats. Here, we characterized six lark mitogenomes from Alaudala cheleensis, Eremophila alpestris, Alauda razae, and Calandrella cinerea and reconstructed the phylogeny of Passerida. Through further comparative analysis among larks, we traced the evolutionary process of CR2. The mitochondrial gene orders were conserved in all published lark mitogenomes, with Cytb-trnT-CR1-trnP-ND6-trnE-remnant CR2 with tandem repeat-trnF-rrnS. Phylogenetic analysis revealed Alaudidae and Panuridae are sister groups at the base of Sylvioidea, and sporadic losses of CR2 may occur in their common ancestor. CR sequence and phylogeny analysis indicated CR2 tandem repeats were generated within CR2, originating in the ancestor of all larks, rather than inherited from CR1. The secondary structure comparison of tandem repeat units within and between species suggested slipped-strand mispairing and DNA turnover as suitable models for explaining the origin and evolution of these repeats. This study reveals the evolutionary process of the CR2 containing tandem repeat in Alaudidae, providing reference for understanding the evolutionary characteristics and dynamics of tandem repeats.

Type specimens and type localities of birds (Aves) collected by Gustav Adolf Fischer (18481886) in East Africa

Gustav Adolf Fischer (18481886) was an important German naturalist and ornithological collector in Africa. The extensive ornithological collections from his several expeditions were sent to at least two different museums (Zoological Museum Berlin and Zoological Museum Hamburg), and to a smaller extent, some private sales are assumed. Together with Anton Reichenow, Fischer described many species new to science, so the many types in the collections make them especially important. We located 616 bird specimens in the collections of the Museum fr Naturkunde Berlin and the Museum der Natur Hamburg along with a few further specimens in the collections of the Museum Heineanum Halberstadt, the Naturalis Biodiversity Center Leiden, The Natural History Museum-UK (Tring), the American Museum of Natural History, New York and the National Museum of Natural History, Washington, D.C. Among these, type specimens for 133 species and subspecies of birds were traced. For all of them, an update of the collecting localities with particular emphasis on some previously unknown locations is provided. A taxonomic update of the described species is given and a designation of lectotypes for the following four taxa was necessary: Francolinus Altumi G.A. Fischer & Reichenow, 1884, Accipiter minullus tropicalis Reichenow, 1898, Hypochera ultramarina var. orientalisReichenow, 1894 and Spinus citrinelloides frontalis Reichenow, 1904b.

Cloacal microbiota are biogeographically structured in larks from desert, tropical and temperate areas

BACKGROUND

In contrast with macroorganisms, that show well-documented biogeographical patterns in distribution associated with local adaptation of physiology, behavior and life history, strong biogeographical patterns have not been found for microorganisms, raising questions about what determines their biogeography. Thus far, large-scale biogeographical studies have focused on free-living microbes, paying little attention to host-associated microbes, which play essential roles in physiology, behavior and life history of their hosts. Investigating cloacal gut microbiota of closely-related, ecologically similar free-living songbird species (Alaudidae, larks) inhabiting desert, temperate and tropical regions, we explored influences of geographical location and host species on α-diversity, co-occurrence of amplicon sequence variants (ASVs) and genera, differentially abundant and dominant bacterial taxa, and community composition. We also investigated how geographical distance explained differences in gut microbial community composition among larks.

RESULTS

Geographic location did not explain variation in richness and Shannon diversity of cloacal microbiota in larks. Out of 3798 ASVs and 799 bacterial genera identified, 17 ASVs (< 0.5%) and 43 genera (5%) were shared by larks from all locations. Desert larks held fewer unique ASVs (25%) than temperate zone (31%) and tropical larks (34%). Five out of 33 detected bacterial phyla dominated lark cloacal gut microbiomes. In tropical larks three bacterial classes were overrepresented. Highlighting the distinctiveness of desert lark microbiota, the relative abundances of 52 ASVs differed among locations, which classified within three dominant and 11 low-abundance phyla. Clear and significant phylogenetic clustering in cloacal microbiota community composition (unweighted UniFrac) showed segregation with geography and host species, where microbiota of desert larks were distinct from those of tropical and temperate regions. Geographic distance was nonlinearly associated with pairwise unweighted UniFrac distances.

CONCLUSIONS

We conclude that host-associated microbiota are geographically structured in a group of widespread but closely-related bird species, following large-scale macro-ecological patterns and contrasting with previous findings for free-living microbes. Future work should further explore if and to what extent geographic variation in host-associated microbiota can be explained as result of co-evolution between gut microbes and host adaptive traits, and if and how acquisition from the environmental pool of bacteria contributes to explaining host-associated communities.

Climatic refugia and reduced extinction correlate with underdispersion in mammals and birds in Africa

Macroevolutionary patterns, often inferred from metrics of community relatedness, are often used to ascertain major evolutionary processes shaping communities. These patterns have been shown to be informative of biogeographic barriers, of habitat suitability and invasibility (especially with regard to environmental filtering), and of regions that function as evolutionary cradles (i.e., sources of diversification) or museums (i.e., regions of reduced extinction). Here, we analyzed continental datasets of mammal and bird distributions to identify primary drivers of community evolution on the African continent for mostly endothermic vertebrates. We find that underdispersion (i.e., relatively low phylogenetic diversity compared to species richness) closely correlates with specific ecoregions that have been identified as climatic refugia in the literature, regardless of whether these specific regions have been touted as cradles or museums. Using theoretical models of identical communities that differ only with respect to extinction rates, we find that even small suppressions of extinction rates can result in underdispersed communities, supporting the hypothesis that climatic stability can lead to underdispersion. We posit that large-scale patterns of under- and overdispersion between regions of similar species richness are more reflective of a particular region's extinction potential, and that the very nature of refugia can lead to underdispersion via the steady accumulation of species richness through diversification within the same ecoregion during climatic cycles. Thus, patterns of environmental filtering can be obfuscated by environments that coincide with biogeographic refugia, and considerations of regional biogeographic history are paramount for inferring macroevolutionary processes.

The complete mitochondrial genome and the phylogenetic position of Alauda gulgula (Aves: Passeriformes)

The Oriental Skylark (Alauda gulgula) is a small songbird in the Alaudidae. Here we assembled the complete mitochondrial genome of Alauda gulgula which is 17,055 bp in length and consisting of 13 protein-coding genes (PCGs), 2 ribosomal RNA, 22 transfer RNA, and 2 extensive heteroplasmic control regions. The overall A + T content of the mitogenome is 52.3%The maximum-likelihood (ML) tree based on the complete mitochondrial genome of A. gulgula revealed the close genetic relationship between A. gulgula and A. arvensis, but separate from A. cheleensis.

Molecular Species Delimitation of Larks (Aves: Alaudidae), and Integrative Taxonomy of the Genus Calandrella, with the Description of a Range-Restricted African Relic Taxon

Larks constitute an avian family of exceptional cryptic diversity and striking examples of convergent evolution. Therefore, traditional morphology-based taxonomy has recurrently failed to reflect evolutionary relationships. While taxonomy ideally should integrate morphology, vocalizations, behaviour, ecology, and genetics, this can be challenging for groups that span several continents including areas that are difficult to access. Here, we combine morphometrics and mitochondrial DNA to evaluate the taxonomy of Calandrella larks, with particular focus on the African C. cinerea and the Asian C. acutirostris complexes. We describe a new range-restricted West African taxon, Calandrella cinerea rufipecta ssp. nov. (type locality: Jos, Plateau State, Nigeria), with an isolated relic population 3000 km from its closest relative in the Rift Valley. We performed molecular species delimitation, employing coalescence-based multi-rate Poisson Tree Processes (mPTP) on cytochrome b sequences across 52 currently recognized lark species, including multiple taxa currently treated as subspecies. Three species-level splits were inferred within the genus Calandrella and another 13 across other genera, primarily among fragmented sub-Saharan taxa and taxa distributed from Northwest Africa to Arabia or East Africa. Previously unknown divergences date back as far as to the Miocene, indicating the presence of currently unrecognized species. However, we stress that taxonomic decisions should not be based on single datasets, such as mitochondrial DNA, although analyses of mitochondrial DNA can be a good indicator of taxa in need of further integrative taxonomic assessment.

Multiple species delimitation approaches applied to the avian lark genus Alaudala

Species delimitation has advanced from a purely phenotypic exercise to a branch of science that integrates multiple sources of data to identify independently evolving lineages that can be treated as species. We here test species limits in the avian Lesser Short-toed Lark Alaudala rufesens-Sand Lark A. raytal complex, which has an intricate taxonomic history, ranging from a single to three recognised species, with different inclusiveness in different treatments. Our integrative taxonomic approach is based on a combination of DNA sequences, plumage, biometrics, songs, song-flights, geographical distributions, habitat, and bioclimatic data, and using various methods including a species delimitation program (STACEY) based on the multispecies coalescent model. We propose that four species should be recognised: Lesser Short-toed Lark A. rufescens (sensu stricto), Heine's Short-toed Lark A. heinei, Asian Short-toed Lark A. cheleensis and Sand Lark A. raytal. There is also some evidence suggesting lineage separation within A. cheleensis and A. raytal, but additional data are required to evaluate this. The species delimitation based on STACEY agrees well with the non-genetic data. Although computer-based species delimitation programs can be useful in identifying independently evolving lineages, we stress that whenever possible, species hypotheses proposed by these programs should be tested by independent, non-genetic data. Our results highlight the difficulty and subjectivity of delimiting lineages and species, especially at early stages in the speciation process.

High MHC gene copy number maintains diversity despite homozygosity in a Critically Endangered single-island endemic bird, but no evidence of MHC-based mate choice

Small population sizes can, over time, put species at risk due to the loss of genetic variation and the deleterious effects of inbreeding. Losing diversity in the major histocompatibility complex (MHC) could be particularly harmful, given its key role in the immune system. Here, we assess MHC class I (MHC-I) diversity and its effects on mate choice and survival in the Critically Endangered Raso lark Alauda razae, a species restricted to the 7 km² islet of Raso, Cape Verde, since ~1460, whose population size has dropped as low as 20 pairs. Exhaustively genotyping 122 individuals, we find no effect of MHC-I genotype/diversity on mate choice or survival. However, we demonstrate that MHC-I diversity has been maintained through extreme bottlenecks by retention of a high number of gene copies (at least 14), aided by cosegregation of multiple haplotypes comprising 2-8 linked MHC-I loci. Within-locus homozygosity is high, contributing to low population-wide diversity. Conversely, each individual had comparably many alleles, 6-16 (average 11), and the large and divergent haplotypes occur at high frequency in the population, resulting in high within-individual MHC-I diversity. This functional immune gene diversity will be of critical importance for this highly threatened species' adaptive potential.

Biomolecular analyses reveal the age, sex and species identity of a near-intact Pleistocene bird carcass

Ancient remains found in permafrost represent a rare opportunity to study past ecosystems. Here, we present an exceptionally well-preserved ancient bird carcass found in the Siberian permafrost, along with a radiocarbon date and a reconstruction of its complete mitochondrial genome. The carcass was radiocarbon dated to approximately 44–49 ka BP, and was genetically identified as a female horned lark. This is a species that usually inhabits open habitat, such as the steppe environment that existed in Siberia at the time. This near-intact carcass highlights the potential of permafrost remains for evolutionary studies that combine both morphology and ancient nucleic acids.

Nicolas Dussex et al. identify a 44,000–49,000 year old bird found in Siberian permafrost as a female horned lark using ancient DNA. This exceptionally well-preserved specimen illustrates the potential contribution to science of permafrost deposits, such as the study of ecology and evolution of ancient ecosystems, calibration of molecular clocks, and furthering our understanding of processes such as biological regulation and gene expression in relation to climate change.

Repeated sex chromosome evolution in vertebrates supported by expanded avian sex chromosomes

Sex chromosomes have evolved from the same autosomes multiple times across vertebrates, suggesting that selection for recombination suppression has acted repeatedly and independently on certain genetic backgrounds. Here, we perform comparative genomics of a bird clade (larks and their sister lineage; Alaudidae and Panuridae) where multiple autosome-sex chromosome fusions appear to have formed expanded sex chromosomes. We detected the largest known avian sex chromosome (195.3 Mbp) and show that it originates from fusions between parts of four avian chromosomes: Z, 3, 4A and 5. Within these four chromosomes, we found evidence of five evolutionary strata where recombination had been suppressed at different time points, and show that stratum age explained the divergence rate of Z-W gametologs. Next, we analysed chromosome content and found that chromosome 3 was significantly enriched for genes with predicted sex-related functions. Finally, we demonstrate extensive homology to sex chromosomes in other vertebrate lineages: chromosomes Z, 3, 4A and 5 have independently evolved into sex chromosomes in fish (Z), turtles (Z, 5), lizards (Z, 4A), mammals (Z, 4A) and frogs (Z, 3, 4A, 5). Our results provide insights into and support for repeated evolution of sex chromosomes in vertebrates.

DNA barcoding a unique avifauna: an important tool for evolution, systematics and conservation

BACKGROUND

DNA barcoding utilises a standardised region of the cytochrome c oxidase I (COI) gene to identify specimens to the species level. It has proven to be an effective tool for identification of avian samples. The unique island avifauna of New Zealand is taxonomically and evolutionarily distinct. We analysed COI sequence data in order to determine if DNA barcoding could accurately identify New Zealand birds.

RESULTS

We sequenced 928 specimens from 180 species. Additional Genbank sequences expanded the dataset to 1416 sequences from 211 of the estimated 236 New Zealand species. Furthermore, to improve the assessment of genetic variation in non-endemic species, and to assess the overall accuracy of our approach, sequences from 404 specimens collected outside of New Zealand were also included in our analyses. Of the 191 species represented by multiple sequences, 88.5% could be successfully identified by their DNA barcodes. This is likely a conservative estimate of the power of DNA barcoding in New Zealand, given our extensive geographic sampling. The majority of the 13 groups that could not be distinguished contain recently diverged taxa, indicating incomplete lineage sorting and in some cases hybridisation. In contrast, 16 species showed evidence of distinct intra-species lineages, some of these corresponding to recognised subspecies. For species identification purposes a character-based method was more successful than distance and phylogenetic tree-based methods.

CONCLUSIONS

DNA barcodes accurately identify most New Zealand bird species. However, low levels of COI sequence divergence in some recently diverged taxa limit the identification power of DNA barcoding. A small number of currently recognised species would benefit from further systematic investigations. The reference database and analysis presented will provide valuable insights into the evolution, systematics and conservation of New Zealand birds.

Thermal physiology of a range-restricted desert lark

Much recent work on avian physiological adaptation to desert environments has focused on larks (Passeriformes: Alaudidae). We tested the prediction that the threatened red lark (Calendulauda burra), a species restricted to very arid parts of South Africa and which is not known to drink, exhibits highly efficient evaporative cooling and makes pronounced use of facultative hyperthermia when exposed to high air temperatures (T_a). We also predicted that C. burra possesses similarly low basal metabolic rate (BMR) and total evaporative water loss (EWL) at moderate T_a as reported for species from the deserts of the Middle East. Rest-phase thermoregulation in C. burra was characterized by an unusually low lower critical limit of thermoneutrality at T_a = ~ 21 °C and a BMR of 0.317 ± 0.047 W, the lowest BMR relative to allometrically-expected values yet reported in any lark. During the diurnal active phase, red larks were able to tolerate T_a up to 50 °C, with the onset of panting occurring at T_a = 38 °C. Maximum EWL was 1.475 ± 0.107 g h^- 1 at T_a = 50 °C, equivalent to 620% of minimum EWL at thermoneutrality. The maximum ratio of evaporative heat dissipation to metabolic heat production was 1.58, a value towards the lower end of the range reported for passerines. Our data support the prediction that C. burra shows metabolic traits similar to those of other larks inhabiting extremely arid climates, but not the notion that evaporative cooling at high T_a in this species is more efficient than in most passerines.

Complete species-level phylogeny of the leaf warbler (Aves: Phylloscopidae) radiation

The leaf warbler radiation (Aves: Phylloscopidae) has undergone a c. 50% increase in the number of recognised species over the last three decades, mainly as a result of analyses of vocalisations and DNA. Using a multilocus dataset for all of the species in this family, and multispecies coalescent-based as well as concatenation methods, we provide the first complete species-level phylogeny for this important group, as well as an estimate of the timing of diversification. The most recent common ancestor for the family was dated at 11.7 million years ago (mya) (95% highest posterior density 9.8-13.7 mya), and divergence times between sister species ranged from 0.5 mya (0.3-0.8 mya) to 6.1 mya (4.8-7.5 mya). Based on our results, we support synonymising Seicercus with Phylloscopus, which results in a monogeneric Phylloscopidae. We discuss the pros and cons of this treatment,and we argue againstproliferation of taxonomic names,and conclude that a large monogeneric Phylloscopidae leads to the fewest taxonomic changes compared to traditional classifications. We briefly discuss morphological evolution in the light of the phylogeny. The time calibrated phylogeny is a major improvement compared to previous studies based on a smaller number of species and loci and can provide a basis for future studies of other aspects of phylloscopid evolution.

Rapid and recent diversification patterns in Anseriformes birds: Inferred from molecular phylogeny and diversification analyses

The Anseriformes is a well-known and widely distributed bird order, with more than 150 species in the world. This paper aims to revise the classification, determine the phylogenetic relationships and diversification patterns in Anseriformes by exploring the Cyt b, ND2, COI genes and the complete mitochondrial genomes (mito-genomes). Molecular phylogeny and genetic distance analyses suggest that the Dendrocygna species should be considered as an independent family, Dendrocygnidae, rather than a member of Anatidae. Molecular timescale analyses suggests that the ancestral diversification occurred during the Early Eocene Climatic Optimum (58 ~ 50 Ma). Furthermore, diversification analyses showed that, after a long period of constant diversification, the median initial speciation rate was accelerated three times, and finally increased to approximately 0.3 sp/My. In the present study, both molecular phylogeny and diversification analyses results support that Anseriformes birds underwent rapid and recent diversification in their evolutionary history, especially in modern ducks, which show extreme diversification during the Plio-Pleistocene (~ 5.3 Ma). Therefore, our study support that the Plio-Pleistocene climate fluctuations are likely to have played a significant role in promoting the recent diversification for Anseriformes.

Equilibrium Bird Species Diversity in Atlantic Islands

Half a century ago, MacArthur and Wilson proposed that the number of species on islands tends toward a dynamic equilibrium diversity around which species richness fluctuates [1]. The current prevailing view in island biogeography accepts the fundamentals of MacArthur and Wilson's theory [2] but questions whether their prediction of equilibrium can be fulfilled over evolutionary timescales, given the unpredictable and ever-changing nature of island geological and biotic features [3-7]. Here we conduct a complete molecular phylogenetic survey of the terrestrial bird species from four oceanic archipelagos that make up the diverse Macaronesian bioregion-the Azores, the Canary Islands, Cape Verde, and Madeira [8, 9]. We estimate the times at which birds colonized and speciated in the four archipelagos, including many previously unsampled endemic and non-endemic taxa and their closest continental relatives. We develop and fit a new multi-archipelago dynamic stochastic model to these data, explicitly incorporating information from 91 taxa, both extant and extinct. Remarkably, we find that all four archipelagos have independently achieved and maintained a dynamic equilibrium over millions of years. Biogeographical rates are homogeneous across archipelagos, except for the Canary Islands, which exhibit higher speciation and colonization. Our finding that the avian communities of the four Macaronesian archipelagos display an equilibrium diversity pattern indicates that a diversity plateau may be rapidly achieved on islands where rates of in situ radiation are low and extinction is high. This study reveals that equilibrium processes may be more prevalent than recently proposed, supporting MacArthur and Wilson's 50-year-old theory.

Complete mitochondrial genome of the Mongolian lark, Melanocorypha mongolica (Aves: Passeriformes)

We sequenced the entire mitochondrial genome of Melanocorypha mongolica for the first time. The mitogenome is 17,358 bp in length, which includes 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA (rRNA) genes, the control region (CR1), and the control region 2 (CR2). Gene order follows a pattern similar to those of Eurasian Skylark. Using mitochondrial genomes of Melanocorypha mongolica and other seven reference birds in Sylvioidea, we preformed Bayesian analysis based on concatenated protein-coding genes. The results reveal that Alaudidae and Acroccphalidae are clustered together, which is sister to the branches included Sylviidae and Leiothrichidae. Further sequencing of mitochondrial genomes in Alaudidae is useful to advance phylogenetic relationship of species in the family.

Are antimicrobial defences in bird eggs related to climatic conditions associated with risk of trans-shell microbial infection?

Introduction

All bird eggs are exposed to microbes in the environment, which if transmitted to the developing embryo, could cause hatching failure. However, the risk of trans-shell infection varies with environmental conditions and is higher for eggs laid in wetter environments. This might relate to generally higher microbial abundances and diversity in more humid environments, including on the surface of eggshells, as well as the need for moisture to facilitate microbial penetration of the eggshell. To protect against microbial infection, the albumen of avian eggs contains antimicrobial proteins, including lysozyme and ovotransferrin. We tested whether lysozyme and ovotransferrin activities varied in eggs of larks (Alaudidae) living along an arid-mesic gradient of environmental aridity, which we used as a proxy for risk of trans-shell infection.

Results

Contrary to expectations, lysozyme activity was highest in eggs from hotter, more arid locations, where we predicted the risk of trans-shell infection would be lower. Ovotransferrin concentrations did not vary with climatic factors. Temperature was a much better predictor of antimicrobial protein activity than precipitation, a result inconsistent with studies stressing the importance of moisture for trans-shell infection.

Conclusions

Our study raises interesting questions about the links between temperature and lysozyme activity in eggs, but we find no support for the hypothesis that antimicrobial protein deposition is higher in eggs laid in wetter environments.

Limited phylogeographic signal in sex-linked and autosomal loci despite geographically, ecologically, and phenotypically concordant structure of mtDNA variation in the Holarctic avian genus Eremophila

Phylogeographic studies of Holarctic birds are challenging because they involve vast geographic scale, complex glacial history, extensive phenotypic variation, and heterogeneous taxonomic treatment across countries, all of which require large sample sizes. Knowledge about the quality of phylogeographic information provided by different loci is crucial for study design. We use sequences of one mtDNA gene, one sex-linked intron, and one autosomal intron to elucidate large scale phylogeographic patterns in the Holarctic lark genus Eremophila. The mtDNA ND2 gene identified six geographically, ecologically, and phenotypically concordant clades in the Palearctic that diverged in the Early - Middle Pleistocene and suggested paraphyly of the horned lark (E. alpestris) with respect to the Temminck's lark (E. bilopha). In the Nearctic, ND2 identified five subclades which diverged in the Late Pleistocene. They overlapped geographically and were not concordant phenotypically or ecologically. Nuclear alleles provided little information on geographic structuring of genetic variation in horned larks beyond supporting the monophyly of Eremophila and paraphyly of the horned lark. Multilocus species trees based on two nuclear or all three loci provided poor support for haplogroups identified by mtDNA. The node ages calculated using mtDNA were consistent with the available paleontological data, whereas individual nuclear loci and multilocus species trees appeared to underestimate node ages. We argue that mtDNA is capable of discovering independent evolutionary units within avian taxa and can provide a reasonable phylogeographic hypothesis when geographic scale, geologic history, and phenotypic variation in the study system are too complex for proposing reasonable a priori hypotheses required for multilocus methods. Finally, we suggest splitting the currently recognized horned lark into five Palearctic and one Nearctic species.