Genomic basis of melanin-associated phenotypes suggests colour-specific environmental adaptations in tawny owls

Feathers comprise a series of evolutionary innovations but also harbour colour, a key biological trait known to co-vary with life history or complex traits. Those relationships are particularly true in melanin-based pigmentation species due to known pleiotropic effects of the melanocortin pathway - originating from melanin-associated phenotypes. Here, we explore the molecular basis of melanin colouration and expected co-variation at the molecular level in the melanin-based, colour polymorphic system of the tawny owl (Strix aluco). An extensive body of literature has revealed that grey and brown tawny owl colour morphs differ in a series of life history and behavioural traits. Thus, it is plausible to expect co-variation also at molecular level between colour morphs. To investigate this possibility, we assembled the first draft genome of the species against which we mapped ddRADseq reads from 220 grey and 150 brown morphs - representing 10 years of pedigree data from a population in Southern Finland - and explored genome-wide associations with colour phenotype. Our results revealed putative molecular signatures of cold adaptation strongly associated with the grey phenotype, namely, a non-synonymous substitution in MCHR1, plus 2 substitutions in non-coding regions of FTCD and FAM135A whose genotype combinations obtained a predictive power of up to 100% (predicting grey colour). These suggest a molecular basis of cold environment adaptations predicted to be grey-morph specific. Our results potentially reveal part of the molecular machinery of melanin-associated phenotypes and provide novel insights towards understanding the functional genomics of colour polymorphism in melanin-based pigmented species.

Characterization of the diversity of barn owl's mitochondrial genome reveals high copy number variations in the control region

Mitochondria are known to play an essential role in the cell. These organelles contain their own DNA, which is divided in a coding and non-coding region (NCR). While much of the NCR's function is unknown, tandem repeats have been observed in several vertebrates, with extreme intra-individual, intraspecific and interspecific variation. Taking advantage of a new complete reference for the mitochondrial genome of the Afro-European Barn Owl (Tyto alba), as well as 172 whole genome-resequencing; we (i) describe the reference mitochondrial genome with a special focus on the repeats in the NCR, (ii) quantify the variation in number of copies between individuals, and (iii) explore the possible factors associated with the variation in the number of repetitions. The reference mitochondrial genome revealed a long (256bp) and a short (80bp) tandem repeat in the NCR region. The re-sequenced genomes showed a great variation in number of copies between individuals, with 4 to 38 copies of the Long and 6 to 135 copies of the short repeat. Among the factors associated with this variation between individuals, the tissue used for extraction was the most significant. The exact mechanisms of the formations of these repeats are still to be discovered and understanding them will help explain the maintenance of the polymorphism in the number of copies, as well as their interactions with the metabolism, the aging and health of the individuals.

The genomic architecture of continuous plumage colour variation in the European barn owl (Tyto alba)

The maintenance of colour variation in wild populations has long fascinated evolutionary biologists, although most studies have focused on discrete traits exhibiting rather simple inheritance patterns and genetic architectures. However, the study of continuous colour traits and their potentially oligo- or polygenic genetic bases remains rare in wild populations. We studied the genetics of the continuously varying white-to-rufous plumage coloration of the European barn owl (Tyto alba) using a genome-wide association approach on the whole-genome data of 75 individuals. We confirmed a mutation at the melanocortin-1-receptor gene (MC1R) is involved in the coloration and identified two new regions, located in super-scaffolds 9 and 42. The combination of the three regions explains most of the colour variation (80.37%, 95% credible interval 58.45-100%). One discovered region, located in the sex chromosome, differs between the most extreme colorations in owls sharing a specific MC1R genotype. This region may play a role in the colour sex dimorphism of this species, possibly in interaction with the autosomal MC1R. We thus provide insights into the genetic architecture of continuous colour variation, pointing to an oligogenic basis with potential epistatic effects among loci that should aid future studies understanding how continuous colour variation is maintained in nature.

Description of Amblyomma monteiroae n. sp. (Acari: Ixodidae), a parasite of the great horned owl (Strigiformes: Strigidae) in southern Brazil

In 2020, adult hard ticks (males and females) were collected from great horned owls [Bubo virginianus (Gmelin, 1788)] in the coastal region in southern Brazil. The engorged females were allowed to oviposit in the laboratory and hatched larvae could be obtained. Analyses of the external morphology of the adult ticks revealed that they represent a new species, which was named Amblyomma monteiroae n. sp. Partial sequences of the mitochondrial 16S rRNA gene and the nuclear second internal transcribed spacer (ITS2) were generated from a male and a female. Their 16S rRNA haplotypes were identical to each other and closest (96% identity) to corresponding sequences of Amblyomma parvitarsum Neumann, 1901, and 90% identical to Amblyomma neumanni Ribaga, 1902. Their ITS2 haplotypes were 95.8 to 96.0 identical to the single ITS-2 partial sequence of A. parvitarsum available in GenBank. In the phylogenetic trees inferred by both 16S rRNA and ITS2 partial sequences, A. monteiroae n. sp. formed a clade with A. parvitarsum, with A. neumanni branching sister to this clade. Amblyomma monteiroae n. sp. is genetically and morphologically related to A. parvitarsum. Both tick species are unique in combining the following morphological characters: scutum extensively ornate; eyes rounded and bulging; coxa I with two moderate pointed spurs, the external longer than the internal; a single triangular short spur on coxae II-III; presence of two spines on the tibia of legs II-IV; hypostomal dentition 3/3, trochanters without spurs. However, the males of the two species can be separated by specific features in palps and festoons, whereas the females differ in specific features of the coxal spurs. The larva of A. monteiroae n. sp. can be morphologically distinguished from A. parvitarsum only by morphometry, with the former species being slightly smaller. Currently, A. monteiroae n. sp. is restricted to southern Brazil, and the only known host is B. virginianus (Strigiformes: Strigidae). The present study increases the Amblyomma Brazilian fauna to 34 species.

Range-wide phylogenomics of the Great Horned Owl (Bubo virginianus) reveals deep north-south divergence in northern Peru

The Great Horned Owl (Bubo virginianus) inhabits myriad habitats throughout the Americas and shows complex patterns of individual and geographic morphological variation. The owl family Strigidae is known to follow ecogeographic rules, such as Gloger's rule. Although untested at the species level, these ecogeographic rules may affect B. virginianus plumage coloration and body size. Previous studies have indicated that, despite this species' morphological variability, little genetic differentiation exists across parts of their range. This study uses reduced representation genome-wide nuclear and complete mitochondrial DNA sequence data to assess range-wide relationships among B. virginianus populations and the disputed species status of B. v. magellanicus (Magellanic or Lesser Horned Owl) of the central and southern Andes. We found shallow phylogenetic relationships generally structured latitudinally to the north of the central Andes, and a deep divergence between a southern and northern clade close to the Marañón Valley in the central Andes, a common biogeographic barrier. We identify evidence of gene flow between B. v. magellanicus and other subspecies based on mitonuclear discordance and F-branch statistics. Overall differences in morphology, plumage coloration, voice, and genomic divergence support species status for B. v. magellanicus.

Genomic Variation and Recent Population Histories of Spotted (Strix occidentalis) and Barred (Strix varia) Owls

Spotted owls (SOs, Strix occidentalis) are a flagship species inhabiting old-growth forests in western North America. In recent decades, their populations have declined due to ongoing reductions in suitable habitat caused by logging, wildfires, and competition with the congeneric barred owl (BO, Strix varia). The northern spotted owl (S. o. caurina) has been listed as "threatened" under the Endangered Species Act since 1990. Here, we use an updated SO genome assembly along with 51 high-coverage whole-genome sequences to examine population structure, hybridization, and recent changes in population size in SO and BO. We found that potential hybrids identified from intermediate plumage morphology were a mixture of pure BO, F1 hybrids, and F1 × BO backcrosses. Also, although SO underwent a population bottleneck around the time of the Pleistocene-Holocene transition, their population sizes rebounded and show no evidence of any historical (i.e., 100-10,000 years ago) population decline. This suggests that the current decrease in SO abundance is due to events in the past century. Finally, we estimate that western and eastern BOs have been genetically separated for thousands of years, instead of the previously assumed recent (i.e., <150 years) divergence. Although this result is surprising, it is unclear where the ancestors of western BO lived after the separation. In particular, although BO may have colonized western North America much earlier than the first recorded observations, it is also possible that the estimated divergence time reflects unsampled BO population structure within central or eastern North America.

Cranial evolution in the extinct Rodrigues Island owl Otus murivorus (Strigidae), associated with unexpected ecological adaptations

Island birds that were victims of anthropic extinctions were often more specialist species, having evolved their most distinctive features in isolation, making the study of fossil insular birds most interesting. Here we studied a fossil cranium of the 'giant' extinct scops owl Otus murivorus from Rodrigues Island (Mascarene Islands, southwestern Indian Ocean), to determine any potential unique characters. The fossil and extant strigids were imaged through X-ray microtomography, providing 3D views of external and internal (endocast, inner ear) cranial structures. Geometric morphometrics and analyses of traditional measurements yielded new information about the Rodrigues owl's evolution and ecology. Otus murivorus exhibits a 2-tier "lag behind" phenomenon for cranium and brain evolution, both being proportionately small relative to increased body size. It also had a much more developed olfactory bulb than congeners, indicating an unexpectedly developed olfactory sense, suggesting a partial food scavenging habit. In addition, O. murivorus had the eyes placed more laterally than O. sunia, the species from which it was derived, probably a side effect of a small brain; rather terrestrial habits; probably relatively fearless behavior; and a less vertical posture (head less upright) than other owls (this in part an allometric effect of size increase). These evolutionary features, added to gigantism and wing reduction, make the extinct Rodrigues owl's evolution remarkable, and with multiple causes.

Identification of two novel adenoviruses in smooth-billed ani and tropical screech owl

Avian adenoviruses (AdVs) are a very diverse group of pathogens causing diseases in poultry and wild birds. Wild birds, endangered by habitat loss and habitat fragmentation in the tropical forests, are recognised to play a role in the transmission of various AdVs. In this study, two novel, hitherto unknown AdVs were described from faecal samples of smooth-billed ani and tropical screech owl. The former was classified into genus Aviadenovirus, the latter into genus Atadenovirus, and both viruses most probably represent new AdV species as well. These results show that there is very limited information about the biodiversity of AdVs in tropical wild birds, though viruses might have a major effect on the population of their hosts or endanger even domesticated animals. Surveys like this provide new insights into the diversity, evolution, host variety, and distribution of avian AdVs.

Mitogenome of the little owl Athene noctua and phylogenetic analysis of Strigidae

New advances in molecular approaches for DNA analysis have enhanced our understanding of the phylogenetic relationship of birds. The Little Owl (Athene noctua) is of great significance for the integrated management of forest diseases and control of regional pests. Here, we sequenced and annotated the 17,772 bp complete mitogenome of A. noctua. The mitogenome encoded 37 typical mitochondrial genes: 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA, and one non-coding control region (D-loop). The organization and location of genes in the A. noctua mitogenome were consistent with those reported for other Strigidae birds. Phylogenetic relationships based on Bayesian inference and Maximum likelihood methods showed that A. noctua has close relationships with Athene brama and Glaucidium cuculoides, confirming that A. noctua belongs to the Strigidae family. The phylogenetic relationships among seven genera of the Strigidae family used in this study were: Ninox and the other six genera were far apart, Otus and the clade ((Bubo + Strix) + Asio) were clustered into one branch, and Athene and Glaucidium were clustered into one branch. This phylogenetic classification is consistent with prior taxonomic studies on the Strigidae family. Our results provide new mitogenomic data to support further phylogenetic and taxonomic studies of Strigidae.

Expression of glucocorticoid and mineralocorticoid receptor genes co-varies with a stress-related colour signal in barn owls

Glucocorticoid hormones are important intermediates between an organism and its environment. They enable an organism to adjust its behavioural and physiological processes in response to environmental changes by binding to mineralocorticoid receptors (MR) and glucocorticoid receptors (GR) expressed in many tissues, including the integument. The regulation of glucocorticoids co-varies with melanin-based colouration in numerous species, an association that might result from pleiotropic effects of genes in the melanocortin system and evolve within a signalling context. Most studies have focused on the circulating levels of glucocorticoids disregarding the receptors that mediate their action, and that might partly account for the covariation between the regulation of stress and melanin-based colouration. We investigated the association of the expression levels of GR and MR genes with melanin-based colouration in the growing feathers of nestling barn owls (Tyto alba). We also explored the association between GR and MR expression levels and the expression of genes related to the melanocortin system and melanogenesis to better understand the origin of the link between the expression of receptors to which corticosterone binds and melanin-based colouration. Nestling barn owls displaying larger eumelanic black feather spots expressed GR and MR at lower levels than smaller-spotted individuals. However, we found that the expression of the GR and MR genes was positively rather than negatively correlated with the expression of genes involved in the deposition of melanin pigments at the time we sampled the nestlings. This provides mixed evidence of the association between melanin-based traits and MR and GR gene expression. The finding that the expression of GR and MR was positively associated with the expression of the PCSK2 gene (encoding one of the protein convertase responsible for the production of hormone peptide ACTH and α-MSH) suggests that the melanocortin system may be implicated in the establishment of the covariation between melanin-based colour and the expression of receptors to which glucocorticoids bind. However, further studies investigating the expression of melanin-based traits with stress-related endpoints at different time points of feather development will be necessary to understand better the proximate mechanism linking melanin-based traits with stress.

Baseline and stress-induced corticosterone levels are heritable and genetically correlated in a barn owl population

The hypothalamic-pituitary-adrenal (HPA) axis is responsible for the regulation of corticosterone, a hormone that is essential in the mediation of energy allocation and physiological stress. As a continuous source of challenge and stress for organisms, the environment has promoted the evolution of physiological adaptations and led to a great variation in corticosterone profiles within or among individuals, populations and species. In order to evolve via natural selection, corticosterone levels do not only depend on the strength of selection exerted on them, but also on the extent to which the regulation of corticosterone is heritable. Nevertheless, the heritability of corticosterone profiles in wild populations is still poorly understood. In this study, we estimated the heritability of baseline and stress-induced corticosterone levels in barn owl (Tyto alba) nestlings from 8 years of data, using a multivariate animal model based on a behavioural pedigree. We found that baseline and stress-induced corticosterone levels are strongly genetically correlated (r = 0.68-0.80) and that the heritability of stress-induced corticosterone levels (h2 = 0.24-0.33) was moderate and similar to the heritability of baseline corticosterone levels (h2 = 0.19-0.30). These findings suggest that the regulation of stress-induced corticosterone and baseline levels evolves at a similar pace when selection acts with the same intensity on both traits and that contrary to previous studies, the evolution of baseline and stress-induced level is interdependent in barn owls, as they may be strongly genetically correlated.

Preliminary Findings of Structure and Expression of Opioid Receptor Genes in a Peregrine Falcon ( Falco peregrinus), a Snowy Owl ( Bubo scandiacus), and a Blue-fronted Amazon Parrot ( Amazona aestiva)

To further knowledge of the physiology of opioid receptors in birds, the structure and expression of the μ-, δ-, and κ-opioid receptor genes were studied in a peregrine falcon ( Falco peregrinus), a snowy owl ( Bubo scandiacus), and a blue-fronted Amazon parrot ( Amazona aestiva). Tissue samples were obtained from birds that had been euthanatized for poor release prognosis or medical reasons. Samples were taken from the brain (telencephalon, thalamus, pituitary gland, cerebellum, pons, medulla oblongata, mesencephalon), the spinal cord and dorsal root ganglions, and plantar foot skin. Messenger RNA was recovered, and reverse transcription polymerase chain reaction (RT-PCR) was performed to generate complementary DNA (cDNA) sequences. Gene structures were documented by directly comparing cDNA sequences with recently published genomic sequences for the peregrine falcon and the blue-fronted Amazon parrot or by comparisons with genomic sequences of related species for the snowy owl. Structurally, the avian μ-opioid receptor messenger RNA (mRNA) species were complex, displaying differential splicing, alternative stop codons, and multiple polyadenylation signals. In comparison, the structure of the avian κ-receptor mRNA was relatively simple. In contrast to what is seen in humans, the avian δ-receptor mRNA structure was found to be complex, demonstrating novel 3-prime coding and noncoding exons not identified in mammals. The role of the δ-opioid receptor merits further investigation in avian species.

Isolation by distance versus landscape resistance: Understanding dominant patterns of genetic structure in Northern Spotted Owls (Strix occidentalis caurina)

Landscape genetics investigations examine how the availability and configuration of habitat influence genetic structure of plants and animals. We used landscape genetics to evaluate the role that forest connectivity plays in determining genetic structure of the federally-threatened Northern Spotted Owl (Strix occidentalis caurina) using genotypes of 339 Northern Spotted Owls obtained for 10 microsatellite loci. Spatial clustering analyses identified a distinct genetic cluster at the southern extent of the region examined. This cluster could not be linked to landscape connectivity patterns and suggested that post-Pleistocene processes were involved with its development rather than contemporary landscape configuration. We also compared matrices of pairwise inter-individual genetic distances with resistance distances derived from a circuit-theory based framework. Resistance distances were obtained for an idealized raster map that reflected continuous unimpeded dispersal habitat across the landscape along with five empirically-derived raster maps reflecting the 1870’s, 1940’s, 1986, 1994, and 2012. Resistance distances from the idealized map served as surrogates for linear geographic distances. Relative to idealized conditions, resistance distances were ~250% higher in the 1940’s and ~200% higher from 1986 onward. Resistance distances from the 1870’s were ~40% higher than idealized conditions. Inter-individual genetic distances were most highly correlated with resistance distances from the idealized map rather than any of the empirical maps. Two hypotheses explain our results. First, our results may reflect temporal lags between the onset of large-scale habitat alterations and their novel effects on genetic structure in long-lived species such as Northern Spotted Owls. Second, because Northern Spotted Owls disperse over long distances, our results may indicate that forest habitat has never been sufficiently fragmented to the point where connectivity was disrupted. The second hypothesis could indicate that forest management practices mandated by the Northwest Forest Plan succeeded with one of its primary goals. However, our results do not represent a complete portrayal of the status of Northern Spotted Owls given detection of significant population declines and bottlenecks in other studies. Future investigations based on computer simulations may help distinguish between hypotheses.

Northern Spotted Owl (Strix occidentalis caurina) Genome: Divergence with the Barred Owl (Strix varia) and Characterization of Light-Associated Genes

We report here the assembly of a northern spotted owl (Strix occidentalis caurina) genome. We generated Illumina paired-end sequence data at 90× coverage using nine libraries with insert lengths ranging from ∼250 to 9,600 nt and read lengths from 100 to 375 nt. The genome assembly is comprised of 8,108 scaffolds totaling 1.26 × 109 nt in length with an N50 length of 3.98 × 106 nt. We calculated the genome-wide fixation index (FST) of S. o. caurina with the closely related barred owl (Strix varia) as 0.819. We examined 19 genes that encode proteins with light-dependent functions in our genome assembly as well as in that of the barn owl (Tyto alba). We present genomic evidence for loss of three of these in S. o. caurina and four in T. alba. We suggest that most light-associated gene functions have been maintained in owls and their loss has not proceeded to the same extent as in other dim-light-adapted vertebrates.

Circulating testosterone and feather-gene expression of receptors and metabolic enzymes in relation to melanin-based colouration in the barn owl

Knowledge of how and why secondary sexual characters are associated with sex hormones is important to understand their signalling function. Such a link can occur if i) testosterone participates in the elaboration of sex-traits, ii) the display of an ornament triggers behavioural response in conspecifics that induce a rise in testosterone, or iii) genes implicated in the elaboration of a sex-trait pleiotropically regulate testosterone physiology. To evaluate the origin of the co-variation between melanism and testosterone, we measured this hormone and the expression of enzymes involved in its metabolism in feathers of barn owl (Tyto alba) nestlings at the time of melanogenesis and in adults outside the period of melanogenesis. Male nestlings displaying smaller black feather spots had higher levels of circulating testosterone, potentially suggesting that testosterone could block the production of eumelanin pigments, or that genes involved in the production of small spots pleiotropically regulate testosterone production. In contrast, the enzyme 5α-reductase, that metabolizes testosterone to DHT, was more expressed in feathers of reddish-brown than light-reddish nestlings. This is consistent with the hypothesis that testosterone might be involved in the expression of reddish-brown pheomelanic pigments. In breeding adults, male barn owls displaying smaller black spots had higher levels of circulating testosterone, whereas in females the opposite result was detected during the rearing period, but not during incubation. The observed sex- and age-specific co-variations between black spottiness and testosterone in nestling and adult barn owls may not result from testosterone-dependent melanogenesis, but from melanogenic genes pleiotropically regulating testosterone, or from colour-specific life history strategies that influence testosterone levels.

Heritability of fear of humans in urban and rural populations of a bird species

Flight initiation distance (FID), a measure of an animal's tolerance to human disturbance and a descriptor of its fear of humans, is increasingly employed for conservation purposes and to predict the response of species to urbanization. However, most work devoted to understanding variability in FID has been conducted at the population level and little is still known about inter-individual variability in this behaviour. We estimated the heritability of FID, a factor fundamental to understanding the strength and evolutionary consequences of selection of particular phenotypes associated with human disturbances. We used a population of burrowing owls (Athene cunicularia) monitored long-term and for which FID was previously shown to be highly consistent across an individual's lifespan. Heritability estimates varied between 0.37 and 0.80, depending on the habitat considered (urban-rural) and method used (parent-offspring regressions or animal models). These values are unusually high compared with those previously reported for other behavioural traits. Although more research is needed to fully understand the underlying causes of this resemblance between relatives, selection pressures acting on this behaviour should be seriously considered as an important evolutionary force in animal populations increasingly exposed to human disturbance worldwide.

Genetic Signatures of Demographic Changes in an Avian Top Predator during the Last Century: Bottlenecks and Expansions of the Eurasian Eagle Owl in the Iberian Peninsula

The study of the demographic history of species can help to understand the negative impact of recent population declines in organisms of conservation concern. Here, we use neutral molecular markers to explore the genetic consequences of the recent population decline and posterior recovery of the Eurasian eagle owl (Bubo bubo) in the Iberian Peninsula. During the last century, the species was the object of extermination programs, suffering direct persecution by hunters until the 70's. Moreover, during the last decades the eagle owl was severely impacted by increased mortality due to electrocution and the decline of its main prey species, the European rabbit (Oryctolagus cuniculus). In recent times, the decrease of direct persecution and the implementation of some conservation schemes have allowed the species' demographic recovery. Yet, it remains unknown to which extent the past population decline and the later expansion have influenced the current species' pattern of genetic diversity. We used eight microsatellite markers to genotype 235 eagle owls from ten Spanish subpopulations and analyse the presence of genetic signatures attributable to the recent population fluctuations experienced by the species. We found moderate levels of differentiation among the studied subpopulations and Bayesian analyses revealed the existence of three genetic clusters that grouped subpopulations from central, south-western and south-eastern Spain. The observed genetic structure could have resulted from recent human-induced population fragmentation, a patchy distribution of prey populations and/or the philopatric behaviour and habitat selection of the species. We detected an old population bottleneck, which occurred approximately 10,000 years ago, and significant signatures of recent demographic expansions. However, we did not find genetic signatures for a recent bottleneck, which may indicate that population declines were not severe enough to leave detectable signals on the species genetic makeup or that such signals have been eroded by the rapid demographic recovery experienced by the species in recent years.

Evolutionary conservation of Kv3.1 in the barn owl Tyto alba

For prey capture in the dark, the barn owl Tyto alba has evolved into an auditory specialist with an exquisite capability of sound localization. Adaptations include asymmetrical ears, enlarged auditory processing centers, the utilization of minute interaural time differences, and phase locking along the entire hearing range up to 10 kHz. Adaptations on the molecular level have not yet been investigated. Here, we tested the hypothesis that divergence in the amino acid sequence of the voltage-gated K(+) channel Kv3.1 contributes to the accuracy and high firing rates of auditory neurons in the barn owl. We therefore cloned both splice variants of Kcnc1, the gene encoding Kv3.1. Both splice variants, Kcnc1a and Kcnc1b, encode amino acids identical to those of the chicken, an auditory generalist. Expression analyses confirmed neural-restricted expression of the channel. In summary, our data reveal strong evolutionary conservation of Kcnc1 in the barn owl and point to other genes involved in auditory specializations of this animal. The data also demonstrate the feasibility to address neuroethological questions in organisms with no reference genome by molecular approaches. This will open new avenues for neuroethologists working in these organisms.

A new owl species of the genus Otus (aves: strigidae) from Lombok, Indonesia

The avifauna of Indonesia is one of the richest in the world but the taxonomic status of many species remains poorly documented. The sole species of scops owl known from Lombok has long been assigned to the widespread Moluccan Scops Owl Otus magicus on the basis of superficial similarities in morphology. Field work in 2003 has shown that the territorial song of the scops owls inhabiting the foothills of Gunung Rinjani differs dramatically from that of O. magicus and is more similar to those of Rufescent Scops Owl O. rufescens and Singapore Scops Owl O. cnephaeus. Detailed comparisons of sound recordings and museum specimens with those of other scops owls in Wallacea and the Indo-Malayan region have confirmed the distinctiveness of the Lombok population. We describe Otus jolandae as a new species, the Rinjani Scops Owl. It is locally common at elevations from 25–1350 m. and occurs within Gunung Rinjani National Park. The new species is known from seven specimens collected by Alfred Everett in 1896. Otus jolandae represents the first endemic bird species from Lombok.

Identification of avian wax synthases

BACKGROUND

Bird species show a high degree of variation in the composition of their preen gland waxes. For instance, galliform birds like chicken contain fatty acid esters of 2,3-alkanediols, while Anseriformes like goose or Strigiformes like barn owl contain wax monoesters in their preen gland secretions. The final biosynthetic step is catalyzed by wax synthases (WS) which have been identified in pro- and eukaryotic organisms.

RESULTS

Sequence similarities enabled us to identify six cDNAs encoding putative wax synthesizing proteins in chicken and two from barn owl and goose. Expression studies in yeast under in vivo and in vitro conditions showed that three proteins from chicken performed WS activity while a sequence from chicken, goose and barn owl encoded a bifunctional enzyme catalyzing both wax ester and triacylglycerol synthesis. Mono- and bifunctional WS were found to differ in their substrate specificities especially with regard to branched-chain alcohols and acyl-CoA thioesters. According to the expression patterns of their transcripts and the properties of the enzymes, avian WS proteins might not be confined to preen glands.

CONCLUSIONS

We provide direct evidence that avian preen glands possess both monofunctional and bifunctional WS proteins which have different expression patterns and WS activities with different substrate specificities.

Comparative analysis of complete mitochondrial DNA control region of four species of Strigiformes

The sequence of the whole mitochondrial (mt) DNA control region (CR) of four species of Strigiformes was obtained. Length of the CR was 3,290 bp, 2,848 bp, 2,444 bp, and 1,771 bp for Asio flammeus, Asio otus, Athene noctua, and Strix aluco, respectively. Interestingly, the length of the control region was maximum in Asio flammeus among all the avian mtDNA control regions sequenced thus far. In addition, the base composition and organization of mtDNA CR of Asio flammeus were identical to those reported for other birds. On the basis of the differential frequencies of base substitutions, the CR may be divided two variable domains, I and III, and a central conserved domain, II. The 3' end of the CR contained many tandem repeats of varying lengths and repeat numbers. In Asio flammeus, the repeated sequences consisted of a 126 bp sequence that was repeated seven times and a 78 bp sequence that was repeated 14 times. In Asio otus, there were also two repeated sequences, namely a 127 bp sequence that was repeated eight times and a 78 bp sequence that was repeated six times. The control region of Athene noctua contained three sets of repeats: a 89 bp sequence that was repeated three times, a 77 bp sequence that was repeated four times, and a 71 bp sequence that was repeated six times. Strix aluco, however, had only one repeated sequence, a 78 bp sequence that was repeated five times. The results of this study seem to indicate that these tandem repeats may have resulted from slipped-strand mispairing during mtDNA replication. Moreover, there are many conserved motifs within the repeated units. These sequences could form stable stem-loop secondary structures, which suggests that these repeated sequences play an important role in regulating transcription and replication of the mitochondrial genome.

Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco)

A recent study of mitochondrial phylogeography of tawny owls (Strix aluco) in western Europe suggested that this species survived the Pleistocene glaciations in three allopatric refugia located in Iberia, Italy, and the Balkans, and the latter was likely the predominant source of postglacial colonization of northern Europe. New data from seven microsatellite loci from 184 individual owls distributed among 14 populations were used to assess the genetic congruence between nuclear and mitochondrial DNA (mtDNA) markers. Microsatellites corroborated the major phylogeographical conclusions reached on the basis of the mtDNA sequences, but also showed important differences leading to novel inferences. Microsatellites corroborated the three major refugia and supported the Balkan origin of northern populations. When corrected for differences in effective population size, microsatellites and mtDNA yielded generally congruent overall estimates of population structure (N*ST=0.12 vs. RST=0.16); however, there was substantial heterogeneity in the RST among the seven nuclear loci that was not correlated with heterozygosity. Populations representing the Balkans postglacial expansion interact with populations from the other two refugia forming two clines near the Alps and the Pyrenees. In both cases, the apparent position of the contact zones differed substantially between markers due to the genetic composition of populations sampled in northern Italy and Madrid. Microsatellite data did not corroborate the lower genetic diversity of northern, recently populated regions as was found with mtDNA; this discrepancy was taken as evidence for a recent bottleneck recovery. Finally, this study suggests that congruence among genetic markers should be more likely in cases of range expansion into new areas than when populations interact across contact zones.

Genetic structure, introgression, and a narrow hybrid zone between northern and California spotted owls (Strix occidentalis)

The northern spotted owl (Strix occidentalis caurina) is a threatened subspecies and the California spotted owl (Strix occidentalis occidentalis) is a subspecies of special concern in the western United States. Concern for their continued viability has arisen because of habitat loss caused by timber harvesting. The taxonomic status of the northern subspecies has been the subject of continuing controversy. We investigated the phylogeographical and population genetic structure of northern and California spotted owls with special reference to their region of contact. Mitochondrial DNA (mtDNA) control region sequences confirmed the existence of two well-differentiated lineages connected by a narrow hybrid zone in a region of low population density in north central California. Maximum-likelihood estimates indicated bidirectional gene flow between the lineages but limited introgression outside the region of contact. The lengths of both the mtDNA hybrid zone and the reduced density patch were similar and slightly exceeded estimates of natal dispersal distances. This suggests that the two subspecies were in secondary contact in a hybrid zone trapped by a population density trough. Consequently, the zone of interaction is expected to be geographically stable. We discovered a third, rare clade of haplotypes, which we interpreted to be a result of incomplete lineage sorting; those haplotypes result in a paraphyletic northern spotted owl with respect to the California spotted owl. A congeneric species, the barred owl (Strix varia), occasionally hybridizes with spotted owls; our results indicated an upper bound for the frequency of barred owl mtDNA haplotypes in northern spotted owl populations of 3%.

The influence of Pleistocene glacial refugia on tawny owl genetic diversity and phylogeography in western Europe

The glacial refugia hypothesis indicates that during the height of the Pleistocene glaciations the temperate species that are today widespread in western Europe must have survived in small and climatically favourable areas located in the southern peninsulas of Iberia, Italy and Balkans. One such species is the tawny owl, a relatively sedentary, nonmigratory bird presently distributed throughout Europe. It is a tree-nesting species closely associated with deciduous and mixed coniferous woodlands. In this study I used control region mtDNA sequences from 187 individuals distributed among 14 populations to determine whether current genetic patterns in tawny owl populations were consistent with postglacial expansion from peninsular refugia. European, North African and Asian tawny owls were found to represent three distinct lineages, where North Africa is the sister clade to all European owls. Within Europe, I found three well-supported clades that correspond to each of the three allopatric refugia. Expansion patterns indicate that owls from the Balkan refugium repopulated most of northern Europe, while expansion out of Iberia and Italy had only regional effects leading to admixture in France. Estimates of population divergence times between refugia populations are roughly similar, but one order of magnitude smaller between Greece and northern Europe. Based on a wide range of mutation rates and generation times, divergence between refugia appears to date to the Pleistocene.

Convergent evolution of strigiform and caprimulgiform dark-activity is supported by phylogenetic analysis using the arylalkylamine N-acetyltransferase (Aanat) gene

Alternative hypotheses propose the sister order of owls (Strigiformes) to be either day-active raptors (Falconiformes) or dark-active nightjars and allies (Caprimulgiformes). In an effort to identify molecular characters distinguishing between these hypotheses we examined a gene, arylalkylamine N-acetyltransferase (Aanat), potentially associated with the evolution of avian dark-activity. Partial Aanat coding sequences, and two introns, were obtained from the genomic DNA of 16 species: Strigiformes (four species), Falconiformes (four species), Caprimulgiformes (five species), with outgroups: Ciconiiformes (one species), Passeriformes (one species), and Apterygiformes (one species). Phylogenetic trees derived from aligned, evolutionarily conserved Aanat regions did not consistently recover clades corresponding to orders Strigiformes and Falconiformes but did place a caprimulgiform clade more distant from the strigiform and falconiform species than the latter two groups are to each other. This finding was supported by spectral analysis. The taxonomic distribution of seven intronic indels is consistent with the Aanat derived phylogenetic trees and supports conventional family-level groupings within both Strigiformes and Caprimulgiformes. The phylogenetic analyses also indicate that Caprimulgiformes is a polyphyletic grouping. In conclusion the data support, but do not conclusively prove, the proposal that Falconiformes is the sister order to Strigiformes and therefore, that the dark-activity characteristic of Strigiformes and Caprimulgiformes arose by convergent evolution.

Single-generation estimates of individual fitness as proxies for long-term genetic contribution

Individual fitness is a central evolutionary concept, but the question of how it should be defined in empirical studies of natural selection remains contentious. Using founding cohorts from long-term population studies of two species of individually marked birds (collared flycatcher Ficedula albicollis and Ural owl Strix uralensis), we compared a rate-sensitive (lambdaind) and a rate-insensitive (lifetime reproductive success [LRS]) estimate of individual fitness with an estimate of long-term genetic fitness. The latter was calculated as the number of gene copies present in the population after more than two generations, as estimated by tracing genetic lineages and accounting for the fact that populations were not completely closed. When counting fledglings, rate-insensitive estimates of individual fitness correlated better than rate-sensitive estimates with estimated long-term genetic contribution. When counting recruits, both classes of estimates performed equally well. The results support the contention that simple, rate-insensitive measures of fitness, such as LRS, provide a valid and good estimate of fitness in evolutionary studies of natural populations.

A new family of satellite DNA sequences as a major component of centromeric heterochromatin in owls (Strigiformes)

We isolated a new family of satellite DNA sequences from HaeIII- and EcoRI-digested genomic DNA of the Blakiston's fish owl ( Ketupa blakistoni). The repetitive sequences were organized in tandem arrays of the 174 bp element, and localized to the centromeric regions of all macrochromosomes, including the Z and W chromosomes, and microchromosomes. This hybridization pattern was consistent with the distribution of C-band-positive centromeric heterochromatin, and the satellite DNA sequences occupied 10% of the total genome as a major component of centromeric heterochromatin. The sequences were homogenized between macro- and microchromosomes in this species, and therefore intraspecific divergence of the nucleotide sequences was low. The 174 bp element cross-hybridized to the genomic DNA of six other Strigidae species, but not to that of the Tytonidae, suggesting that the satellite DNA sequences are conserved in the same family but fairly divergent between the different families in the Strigiformes. Secondly, the centromeric satellite DNAs were cloned from eight Strigidae species, and the nucleotide sequences of 41 monomer fragments were compared within and between species. Molecular phylogenetic relationships of the nucleotide sequences were highly correlated with both the taxonomy based on morphological traits and the phylogenetic tree constructed by DNA-DNA hybridization. These results suggest that the satellite DNA sequence has evolved by concerted evolution in the Strigidae and that it is a good taxonomic and phylogenetic marker to examine genetic diversity between Strigiformes species.

The complete mitochondrial genomes sequences of Asio flammeus and Asio otus and comparative analysis

The complete mitochondrial genomes of Asio flammeus and Asio otus were sequenced and found to span 18858 bp and 18493 bp, respectively. It is surprising to find the former to be the largest among all avian mitochondrial genomes sequenced so far. The two genomes have very similar gene order with that of Gallus gallus, neither contains the pseudo control region, but both have a single extra base, namely Cytidine, at position 174 in ND3 gene. The control regions of Asio flammeus and Asio otus' mitochondrial genomes span 3288 bp and 2926 bp respectively, which are the longest among vertebrates except for Myxine glutinosa and contribute to the large size of two genomes. The 3' end of the control region of Asio flammeus and Asio otus contains many tandemly repeated sequences, which are highly similar to a putative control element, i.e. Mt5, and may form stable stem-loop secondary structures. Such repeated sequences probably play an important role in regulating transcription and replication of mitochondrial genome. Our results may provide important clues for uncovering the origin and evolution mechanisms of mitochondrion genome.

Comparative analysis of Avian BMAL1 and CLOCK protein sequences: a search for features associated with owl nocturnal behaviour

Animals differ widely in the phasing of their daily rhythms with respect to daily environmental rhythms. While birds are predominantly day-active, nocturnal activity is a characteristic feature of the order Strigiformes (owls). To study the evolution of owl night-activity cDNA sequences encoding the circadian core oscillator (CCO) proteins BMAL1 and CLOCK were obtained from barn owl (Tyto alba). The predicted proteins showed high sequence identity with their Galliform homologues (BMAL1: 99%; CLOCK: 95.6%). A computer-predicted chicken BMAL1 casein kinase-1 phosphorylation site is absent from T. alba BMAL1, but also absent from homologues of other six bird species (5 orders) (night-active (n=2), day-active (n=4)) indicating no evolutionary association with night activity. Sequence differences between T. alba and Galliform CLOCK frequently involved serine and threonine residues suggesting potential differences in their phosphorylation. The length of a poly-glutamine string in the CLOCK C-terminus varied between and within 25 species (6 orders) examined, however, no discernible feature distinguishing day and night active species was found. No differences were found between day (n=5) and night (n=7)-active species (12 species, 6 orders) in a region of the PER2 protein implicated in altered rhythm phasing in humans. In conclusion the avian CCO components examined showed strong evolutionary conservation. Molecular evolution associated with owl night-activity may have involved alterations in the CCO relationship with 'output' genes rather than in the molecular structure of the CCO itself.

Genetic and environmental components of variation in eumelanin and phaeomelanin sex-traits in the barn owl

Knowledge of the mechanism underlying the expression of melanin-based sex-traits may help us to understand their signalling function. Potential sources of inter-individual variation are the total amount of melanins produced but also how biochemical precursors are allocated into the eumelanin and phaeomelanin pigments responsible for black and reddish-brown colours, respectively. In the barn owl (Tyto alba), a eumelanin trait (referred to as 'plumage spottiness') signals immunocompetence towards an artificially administrated antigen and parasite resistance in females, whereas a phaeomelanin trait ('plumage coloration') signals investment in reproduction in males. This raises the question whether plumage coloration and spottiness are expressed independent of each other. To investigate this question, we have studied the genetics of these two plumage traits. Crossfostering experiments showed that, for each trait, phenotypic variation has a strong genetic component, whereas no environmental component could be detected. Plumage coloration is autosomally inherited, as suggested by the similar paternal-to-maternal contribution to offspring coloration. In contrast, plumage spottiness may be sex-linked inherited (in birds, females are heterogametic). That proposition arises from the observation that sons resembled their mother more than their father and that daughters resembled only their father. Despite plumage coloration and spottiness signalling different qualities, these two traits are not inherited independent of each other, darker birds being spottier. This suggests that the extent to which coloration and spottiness are expressed depends on the total amount of melanin produced (with more melanin leading to a both darker and spottier plumage) rather than on differential allocation of melanin into plumage coloration and spottiness (in such a case, darker birds should have been less spotted). A gene controlling the production of melanin pigments may be located on sex-chromosomes, since the phenotypic correlation between coloration and spottiness was stronger in males than in females.

A female melanin ornament signals offspring fluctuating asymmetry in the barn owl

Sexual selection theory predicts that males advertise quality by displaying extravagant ornaments. By contrast, whether phenotypic variation in females has a signalling function remains an open question. Here, to our knowledge, we provide the first evidence that a female plumage trait can signal fluctuating asymmetry in the offspring. We experimentally demonstrate in wild barn owls (Tyto alba) that the extent to which females display black spots on their plumage does not only signal offspring parasite resistance as shown in a previous study but also developmental homeostasis in the offspring. A greater number of spotted females produced offspring that had more symmetrical feathers during the period of growth. Males, that pair non-randomly with respect to female plumage spottiness therefore appear to gain substantial benefits by mating with heavily spotted females. Genetic variation in plumage spottiness is nevertheless maintained as the covariation between offspring body mass and mother plumage spottiness varies annually depending on environmental conditions.

Unravelling a biogeographical knot: origin of the 'leapfrog' distribution pattern of Australo-Papuan sooty owls (Strigiformes) and logrunners (Passeriformes)

Molecular analysis of two Australo-Papuan rainforest birds exhibiting correlated 'leapfrog' patterns were used to elucidate the evolutionary origin of this unusual pattern of geographical differentiation. In both sooty owls (Tyto) and logrunners (Orthonyx), phenotypically similar populations occupy widely disjunct areas (central-eastern Australia and upland New Guinea) with a third, highly distinctive population, occurring between them in northeastern Queensland. Two mechanisms have been proposed to explain the origin of leapfrog patterns in avian distributions: recent shared ancestry of terminal populations and unequal rates or phenotypic change among populations. As the former should generate correlated patterns of phenotypic and genetic differentiation, we tested for a sister relationship between populations from New Guinea and central-eastern Australia using nuclear and mitochondrial DNA sequences. The resulting phylogenies not only refute recent ancestry as an explanation for the leapfrog pattern, but provide evidence of vastly different spatio-temporal histories for sooty owls and logrunners within the Australo-Papuan rainforests. This incongruence indicates that the evolutionary processes responsible for generating leapfrog patterns in these co-distributed taxa are complex, possibly involving a combination of selection and drift in sooty owls and convergence or retention of ancestral characteristics in logrunners.

A modified protocol for sex identification of in ovo avian embryos and its application as a management tool for endangered species conservation programs

A simple, reliable, and safe protocol was developed for the collection of small amounts of blood from avian eggs of variable size and at early stages of development. Fifty eggs were used in the study; 40 were common chicken (Gallus domesticus) eggs, six were homing pigeon (Columba livia domestica) eggs, and four were burrowing owl (Athene cunicularia) eggs. Collection was attempted approximately halfway through incubation. The success rate for collection of blood or blood-tinged fluid from eggs was high, averaging 68% in the chicken eggs, 100% in the homing pigeon eggs, and 75% in the burrowing owl eggs. Collection did not affect subsequent hatchability. This blood could then be used to determine the sex of the embryo by utilizing a DNA probe or restriction fragment length polymorphism technique. Sex identification in ovo allowed the demographic management of small populations of birds within our institution.

Molecular evolution of the avian CHD1 genes on the Z and W sex chromosomes

Genes shared between the nonrecombining parts of the two types of sex chromosomes offer a potential means to study the molecular evolution of the same gene exposed to different genomic environments. We have analyzed the molecular evolution of the coding sequence of the first pair of genes found to be shared by the avian Z (present in both sexes) and W (female-specific) sex chromosomes, CHD1Z and CHD1W. We show here that these two genes evolve independently but are highly conserved at nucleotide as well as amino acid levels, thus not indicating a female-specific role of the CHD1W gene. From comparisons of sequence data from three avian lineages, the frequency of nonsynonymous substitutions (K(a)) was found to be higher for CHD1W (1.55 per 100 sites) than for CHD1Z (0.81), while the opposite was found for synonymous substitutions (K(s), 13.5 vs. 22.7). We argue that the lower effective population size and the absence of recombination on the W chromosome will generally imply that nonsynonymous substitutions accumulate faster on this chromosome than on the Z chromosome. The same should be true for the Y chromosome relative to the X chromosome in XY systems. Our data are compatible with a male-biased mutation rate, manifested by the faster rate of neutral evolution (synonymous substitutions) on the Z chromosome than on the female-specific W chromosome.

Female barn owls (Tyto alba) advertise good genes

The good genes hypothesis of sexual selection postulates that ornamentation signals superior genetic quality to potential mates. Support for this hypothesis comes from studies on male ornamentation only, while it remains to be shown that female ornamentation may signal genetic quality as well. Female barn owls (Tyto alba) display more black spots on their plumage than males. The expression of this plumage trait has a genetic basis and it has been suggested that males prefer to mate with females displaying more black spots. Given the role of parasites in the evolution of sexually selected traits and of the immune system in parasite resistance, we hypothesize that the extent of female plumage 'spottiness' reflects immunological defence. We assessed the genetic variation in specific antibody production against a non-pathogenic antigen among cross-fostered nestlings and studied its covariation with the plumage spottiness of genetic parents. The magnitude of the antibody response was positively correlated with the plumage spottiness of the genetic mother but not of the genetic father. Our study thereby provides the first experimental support, to our knowledge, for the hypothesis that female ornamentation signals genetic quality.

Sex ratio and fledging success of supplementary-fed Tengmalm's owl broods

A nest box population of Tengmalm's owls (Aegolius funereus) in northern Sweden was studied to investigate the effects of extra food on the sex ratio between hatching and fledging in this sexually size-dimorphic species. The brood size and brood sex ratio of supplementary-fed and control broods were compared. Newly hatched nestlings were blood sampled and sexed by polymerase chain reaction (PCR) amplification of the sex-linked CHD1Z and CHD1W genes. The brood sex ratio at hatching was strongly male biased (65%); this was also the case in broods where all eggs hatched (72%). There was no relationship between hatch order and sex ratio, and hatching sex ratio did not vary significantly with laying date. Brood size decreased between hatching and fledging, but did not differ between fed and control broods at either stage. Brood sex ratio did not differ between hatching and fledging, and fledging sex ratio did not differ between fed and control broods. It was concluded that, at least during the year in which the study was carried out, feeding had no effect on brood reduction, and that male and female nestlings did not show any differential mortality. The mechanisms behind the male-biased sex ratio at hatching, and any possible adaptive reasons for it, are not known.