Reduction of genetic diversity in 'Alala (Hawaiian crow; Corvus hawaiiensis) between the late 1800s and the late 1900s.

Genetic and genomic data are increasingly used to aid conservation management of endangered species by providing insights into evolutionary histories, factors associated with extinction risks, and potential for future adaptation. For the 'Alala, or Hawaiian crow (Corvus hawaiiensis), genetic concerns include negative correlations between inbreeding and hatching success. However, it is unclear if low genetic diversity and inbreeding depression are consequences of a historical population bottleneck, or if 'Alala had historically low genetic diversity that predated human influence, perhaps as a result of earlier declines or founding events. In this study, we applied a hybridization-based sequence capture to generate a genome-wide single nucleotide polymorphism (SNP) dataset for comparing historical specimens collected in the 1890s, when 'Alala were more numerous, to samples taken between 1973 and 1998, when 'Alala population densities were near the lowest documented levels in the wild, prior to all individuals being collected for captive rearing. We found low genome-wide diversity in both sample groups, however, the modern sample group (1973 to 1998 cohort) exhibited relatively fewer polymorphic alleles, a lower proportion of polymorphic loci, and lower observed heterozygosity, consistent with a population decline and potential bottleneck effects. These results combined with a current low population size highlight the importance of continued efforts by conservation managers to mitigate inbreeding and maintain founder representation to preserve what genetic diversity remains.

Chromatin accessibility, not 5mC methylation covaries with partial dosage compensation in crows.

The evolution of genetic sex determination is often accompanied by degradation of the sex-limited chromosome. Male heterogametic systems have evolved convergent, epigenetic mechanisms restoring the resulting imbalance in gene dosage between diploid autosomes (AA) and the hemizygous sex chromosome (X). Female heterogametic systems (AAf Zf, AAm ZZm) tend to only show partial dosage compensation (0.5 < Zf:AAf < 1) and dosage balance (0.5

Increased genetic diversity and immigration after West Nile virus emergence in American crows: No evidence for a genetic bottleneck.

Infectious diseases can cause steep declines in wildlife populations, leading to changes in genetic diversity that may affect the susceptibility of individuals to infection and the overall resilience of populations to pathogen outbreaks. Here, we examine evidence for a genetic bottleneck in a population of American crows (Corvus brachyrhynchos) before and after the emergence of West Nile virus (WNV). More than 50% of marked birds in this population were lost over the 2-year period of the epizootic, representing a 10-fold increase in adult mortality. Using analyses of single-nucleotide polymorphisms (SNPs) and microsatellite markers, we tested for evidence of a genetic bottleneck and compared levels of inbreeding and immigration in the pre- and post-WNV populations. Counter to expectations, genetic diversity (allelic diversity and the number of new alleles) increased after WNV emergence. This was likely due to increases in immigration, as the estimated membership coefficients were lower in the post-WNV population. Simultaneously, however, the frequency of inbreeding appeared to increase: Mean inbreeding coefficients were higher among SNP markers, and heterozygosity-heterozygosity correlations were stronger among microsatellite markers, in the post-WNV population. These results indicate that loss of genetic diversity at the population level is not an inevitable consequence of a population decline, particularly in the presence of gene flow. The changes observed in post-WNV crows could have very different implications for their response to future pathogen risks, potentially making the population as a whole more resilient to a changing pathogen community, while increasing the frequency of inbred individuals with elevated susceptibility to disease.

Satellite DNA evolution in Corvoidea inferred from short and long reads.

Satellite DNA (satDNA) is a fast-evolving portion of eukaryotic genomes. The homogeneous and repetitive nature of such satDNA causes problems during the assembly of genomes, and therefore it is still difficult to study it in detail in nonmodel organisms as well as across broad evolutionary timescales. Here, we combined the use of short- and long-read data to explore the diversity and evolution of satDNA between individuals of the same species and between genera of birds spanning ~40 millions of years of bird evolution using birds-of-paradise (Paradisaeidae) and crow (Corvus) species. These avian species highlighted the presence of a GC-rich Corvoidea satellitome composed of 61 satellite families and provided a set of candidate satDNA monomers for being centromeric on the basis of length, abundance, homogeneity and transcription. Surprisingly, we found that the satDNA of crow species rapidly diverged between closely related species while the satDNA appeared more similar between birds-of-paradise species belonging to different genera.

High heritability of telomere length and low heritability of telomere shortening in wild birds.

Telomere length and telomere shortening predict survival in many organisms. This raises the question of the contribution of genetic and environmental effects to variation in these traits, which is still poorly known, particularly for telomere shortening. We used experimental (cross-fostering) and statistical (quantitative genetic "animal models") means to disentangle and estimate genetic and environmental contributions to telomere length variation in pedigreed free-living jackdaws (Corvus monedula). Telomere length was measured twice in nestlings, at ages 4 (n = 715) and 29 days (n = 474), using telomere restriction fragment (TRF) analysis, adapted to exclude interstitial telomeric sequences. Telomere length shortened significantly over the nestling period (10.4 ± 0.3 bp day-1 ) and was highly phenotypically (rP  = 0.95 ± 0.01) and genetically (rG  > 0.99 ± 0.01) correlated within individuals. Additive genetic effects explained a major part of telomere length variation among individuals, with its heritability estimated at h2  = 0.74 on average. We note that TRF-based studies reported higher heritabilities than qPCR-based studies, and we discuss possible explanations. Parent-offspring regressions yielded similar heritability estimates for mothers and fathers when accounting for changes in paternal telomere length over life. Year effects explained a small but significant part of telomere length variation. Heritable variation for telomere shortening was low (h2  = 0.09 ± 0.11). The difference in heritability between telomere length (high) and telomere shortening (low) agrees with evolutionary theory, in that telomere shortening has stronger fitness consequences in this population. Despite the high heritability of telomere length, its evolvability, which scales the additive genetic variance by mean telomere length, was on average 0.48%. Hence, evolutionary change of telomere length due to selection is likely to be slow.

A genome-wide investigation of adaptive signatures in protein-coding genes related to tool behaviour in New Caledonian and Hawaiian crows.

Very few animals habitually manufacture and use tools. It has been suggested that advanced tool behaviour co-evolves with a suite of behavioural, morphological and life history traits. In fact, there are indications for such an adaptive complex in tool-using crows (genus Corvus species). Here, we sequenced the genomes of two habitually tool-using and ten non-tool-using crow species to search for genomic signatures associated with a tool-using lifestyle. Using comparative genomic and population genetic approaches, we screened for signals of selection in protein-coding genes in the tool-using New Caledonian and Hawaiian crows. While we detected signals of recent selection in New Caledonian crows near genes associated with bill morphology, our data indicate that genetic changes in these two lineages are surprisingly subtle, with little evidence at present for convergence. We explore the biological explanations for these findings, such as the relative roles of gene regulation and protein-coding changes, as well as the possibility that statistical power to detect selection in recently diverged lineages may have been insufficient. Our study contributes to a growing body of literature aiming to decipher the genetic basis of recently evolved complex behaviour.

Extra-pair paternity as a strategy to reduce the costs of heterospecific reproduction? Insights from the crow hybrid zone.

Within hybrid zones of socially monogamous species, the number of mating opportunities with a conspecific can be limited. As a consequence, individuals may mate with a heterospecific (social) partner despite possible fitness costs to their hybrid offspring. Extra-pair copulations with a conspecific may thus arise as a possible post hoc strategy to reduce the costs of hybridization. We here assessed the rate of extra-pair paternity in the hybrid zone between all-black carrion crows (Corvus (corone) corone) and grey hooded crows (C. (c.) cornix) and tested whether extra-pair paternity (EPP) was more likely in broods where parents differed in plumage colour. The proportion of broods with at least one extra-pair offspring and the proportion of extra-pair offspring were low overall (6.98% and 2.90%, respectively) with no evidence of hybrid broods having higher EPP rates than purebred nests.

Cryptic and extensive hybridization between ancient lineages of American crows.

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow (Corvus caurinus) and American Crow (Corvus brachyrhynchos) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%-divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late-generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human-modified habitats.

Polymorphism Data Assist Estimation of the Nonsynonymous over Synonymous Fixation Rate Ratio ω for Closely Related Species.

The ratio of nonsynonymous over synonymous sequence divergence, dN/dS, is a widely used estimate of the nonsynonymous over synonymous fixation rate ratio ω, which measures the extent to which natural selection modulates protein sequence evolution. Its computation is based on a phylogenetic approach and computes sequence divergence of protein-coding DNA between species, traditionally using a single representative DNA sequence per species. This approach ignores the presence of polymorphisms and relies on the indirect assumption that new mutations fix instantaneously, an assumption which is generally violated and reasonable only for distantly related species. The violation of the underlying assumption leads to a time-dependence of sequence divergence, and biased estimates of ω in particular for closely related species, where the contribution of ancestral and lineage-specific polymorphisms to sequence divergence is substantial. We here use a time-dependent Poisson random field model to derive an analytical expression of dN/dS as a function of divergence time and sample size. We then extend our framework to the estimation of the proportion of adaptive protein evolution α. This mathematical treatment enables us to show that the joint usage of polymorphism and divergence data can assist the inference of selection for closely related species. Moreover, our analytical results provide the basis for a protocol for the estimation of ω and α for closely related species. We illustrate the performance of this protocol by studying a population data set of four corvid species, which involves the estimation of ω and α at different time-scales and for several choices of sample sizes.

Population Genomics and Structure of the Critically Endangered Mariana Crow (Corvus kubaryi).

The Mariana Crow, or Åga (Corvus kubaryi), is a critically endangered species (IUCN -International Union for Conservation of Nature), endemic to the islands of Guam and Rota in the Mariana Archipelago. It is locally extinct on Guam, and numbers have declined dramatically on Rota to a historical low of less than 55 breeding pairs throughout the island in 2013. Because of its extirpation on Guam and population decline on Rota, it is of critical importance to assess the genetic variation among individuals to assist ongoing recovery efforts. We conducted a population genomics analysis comparing the Guam and Rota populations and studied the genetic structure of the Rota population. We used blood samples from five birds from Guam and 78 birds from Rota. We identified 145,552 candidate single nucleotide variants (SNVs) from a genome sequence of an individual from Rota and selected a subset of these to develop an oligonucleotide in-solution capture assay. The Guam and Rota populations were genetically differentiated from each other. Crow populations sampled broadly across their range on Rota showed significant genetic structuring ⁻ a surprising result given the small size of this island and the good flight capabilities of the species. Knowledge of its genetic structure will help improve management strategies to help with its recovery.

Epistatic mutations under divergent selection govern phenotypic variation in the crow hybrid zone.

The evolution of genetic barriers opposing interspecific gene flow is key to the origin of new species. Drawing from information on over 400 admixed genomes sourced from replicate transects across the European hybrid zone between all-black carrion crows and grey-coated hooded crows, we decipher the interplay between phenotypic divergence and selection at the molecular level. Over 68% of plumage variation was explained by epistasis between the gene NDP and a ~2.8-megabase region on chromosome 18 with suppressed recombination. Both pigmentation loci showed evidence for divergent selection resisting introgression. This study reveals how few, large-effect loci can govern prezygotic isolation and shield phenotypic divergence from gene flow.

In situ quantification of individual mRNA transcripts in melanocytes discloses gene regulation of relevance to speciation.

Functional validation of candidate genes involved in adaptation and speciation remains challenging. Here, we exemplify the utility of a method quantifying individual mRNA transcripts in revealing the molecular basis of divergence in feather pigment synthesis during early-stage speciation in crows. Using a padlock probe assay combined with rolling circle amplification, we quantified cell-type-specific gene expression in the histological context of growing feather follicles. Expression of Tyrosinase Related Protein 1 (TYRP1), Solute Carrier Family 45 member 2 (SLC45A2) and Hematopoietic Prostaglandin D Synthase (HPGDS) was melanocyte-limited and significantly reduced in follicles from hooded crow, explaining the substantially lower eumelanin content in grey versus black feathers. The central upstream Melanocyte Inducing Transcription Factor (MITF) only showed differential expression specific to melanocytes - a feature not captured by bulk RNA-seq. Overall, this study provides insight into the molecular basis of an evolutionary young transition in pigment synthesis, and demonstrates the power of histologically explicit, statistically substantiated single-cell gene expression quantification for functional genetic inference in natural populations.

Apparent inbreeding preference despite inbreeding depression in the American crow.

Although matings between relatives can have negative effects on offspring fitness, apparent inbreeding preference has been reported in a growing number of systems, including those with documented inbreeding depression. Here, we examined evidence for inbreeding depression and inbreeding preference in two populations (Clinton, New York, and Davis, California, USA) of the cooperatively breeding American crow (Corvus brachyrhynchos). We then compared observed inbreeding strategies with theoretical expectations for optimal, adaptive levels of inbreeding, given the inclusive fitness benefits and population-specific magnitude of inbreeding depression. We found that low heterozygosity at a panel of 33 microsatellite markers was associated with low survival probability (fledging success) and low white blood cell counts among offspring in both populations. Despite these costs, our data were more consistent with inbreeding preference than avoidance: The observed heterozygosity among 396 sampled crow offspring was significantly lower than expected if local adults were mating by random chance. This pattern was consistent across a range of spatial scales in both populations. Adaptive levels of inbreeding, given the magnitude of inbreeding depression, were predicted to be very low in the California population, whereas complete disassortative mating was predicted in the New York population. Sexual conflict might have contributed to the apparent absence of inbreeding avoidance in crows. These data add to an increasing number of examples of an "inbreeding paradox," where inbreeding appears to be preferred despite inbreeding depression.

Genomic evidence of speciation reversal in ravens.

Many species, including humans, have emerged via complex reticulate processes involving hybridisation. Under certain circumstances, hybridisation can cause distinct lineages to collapse into a single lineage with an admixed mosaic genome. Most known cases of such 'speciation reversal' or 'lineage fusion' involve recently diverged lineages and anthropogenic perturbation. Here, we show that in western North America, Common Ravens (Corvus corax) have admixed mosaic genomes formed by the fusion of non-sister lineages ('California' and 'Holarctic') that diverged ~1.5 million years ago. Phylogenomic analyses and concordant patterns of geographic structuring in mtDNA, genome-wide SNPs and nuclear introns demonstrate long-term admixture and random interbreeding between the non-sister lineages. In contrast, our genomic data support reproductive isolation between Common Ravens and Chihuahuan Ravens (C. cryptoleucus) despite extensive geographic overlap and a sister relationship between Chihuahuan Ravens and the California lineage. These data suggest that the Common Raven genome was formed by secondary lineage fusion and most likely represents a case of ancient speciation reversal that occurred without anthropogenic causes.

Population mitogenomics provides insights into evolutionary history, source of invasions and diversifying selection in the House Crow (Corvus splendens).

The House Crow (Corvus splendens) is a useful study system for investigating the genetic basis of adaptations underpinning successful range expansion. The species originates from the Indian subcontinent, but has successfully spread through a variety of thermal environments across Asia, Africa and Europe. Here, population mitogenomics was used to investigate the colonisation history and to test for signals of molecular selection on the mitochondrial genome. We sequenced the mitogenomes of 89 House Crows spanning four native and five invasive populations. A Bayesian dated phylogeny, based on the 13 mitochondrial protein-coding genes, supports a mid-Pleistocene (~630,000 years ago) divergence between the most distant genetic lineages. Phylogeographic patterns suggest that northern South Asia is the likely centre of origin for the species. Codon-based analyses of selection and assessments of changes in amino acid properties provide evidence of positive selection on the ND2 and ND5 genes against a background of purifying selection across the mitogenome. Protein homology modelling suggests that four amino acid substitutions inferred to be under positive selection may modulate coupling efficiency and proton translocation mediated by OXPHOS complex I. The identified substitutions are found within native House Crow lineages and ecological niche modelling predicts suitable climatic areas for the establishment of crow populations within the invasive range. Mitogenomic patterns in the invasive range of the species are more strongly associated with introduction history than climate. We speculate that invasions of the House Crow have been facilitated by standing genetic variation that accumulated due to diversifying selection within the native range.

Sequencing of the complete mitochondrial genome of the common raven Corvus corax (Aves: Corvidae) confirms mitogenome-wide deep lineages and a paraphyletic relationship with the Chihuahuan raven C. cryptoleucus.

Previous studies based on single mitochondrial markers have shown that the common raven (Corvus corax) consists of two highly diverged lineages that are hypothesised to have undergone speciation reversal upon secondary contact. Furthermore, common ravens are paraphyletic with respect to the Chihuahuan raven (C. cryptoleucus) based on mitochondrial DNA (mtDNA). Here we explore the causes of mtDNA paraphyly by sequencing whole mitochondrial genomes of 12 common ravens from across the Northern Hemisphere, in addition to three Chihuahuan ravens and one closely related brown-necked raven (C. ruficollis) using a long-range PCR protocol. Our raven mitogenomes ranged between 16925-16928 bp in length. GC content varied from 43.3% to 43.8% and the 13 protein coding genes, two rRNAs and 22 tRNAs followed a standard avian mitochondrial arrangement. The overall divergence between the two common raven clades was 3% (range 0.3-5.8% in 16 regions including the protein coding genes, rRNAs and the control region). Phylogenies constructed from whole mitogenomes recovered the previously found mitochondrial sister relationship between the common raven California clade and the Chihuahuan raven (overall divergence 1.1%), which strengthens the hypothesis that mtDNA paraphyly in the common raven results from speciation reversal of previously distinct Holarctic and California lineages.

Combination of short-read, long-read, and optical mapping assemblies reveals large-scale tandem repeat arrays with population genetic implications.

Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (ρ) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently.

Covariation in levels of nucleotide diversity in homologous regions of the avian genome long after completion of lineage sorting.

Closely related species may show similar levels of genetic diversity in homologous regions of the genome owing to shared ancestral variation still segregating in the extant species. However, after completion of lineage sorting, such covariation is not necessarily expected. On the other hand, if the processes that govern genetic diversity are conserved, diversity may potentially covary even among distantly related species. We mapped regions of conserved synteny between the genomes of two divergent bird species-collared flycatcher and hooded crow-and identified more than 600 Mb of homologous regions (66% of the genome). From analyses of whole-genome resequencing data in large population samples of both species we found nucleotide diversity in 200 kb windows to be well correlated (Spearman's ρ = 0.407). The correlation remained highly similar after excluding coding sequences. To explain this covariation, we suggest that a stable avian karyotype and a conserved landscape of recombination rate variation render the diversity-reducing effects of linked selection similar in divergent bird lineages. Principal component regression analysis of several potential explanatory variables driving heterogeneity in flycatcher diversity levels revealed the strongest effects from recombination rate variation and density of coding sequence targets for selection, consistent with linked selection. It is also possible that a stable karyotype is associated with a conserved genomic mutation environment contributing to covariation in diversity levels between lineages. Our observations imply that genetic diversity is to some extent predictable.

The origin and phylogenetic relationships of the New Zealand ravens.

The relationships of the extinct New Zealand ravens (Corvus spp.) are poorly understood. We sequenced the mitogenomes of the two currently recognised species and found they were sister-taxa to a clade comprising the Australian raven, little raven, and forest raven (C.coronoides, C. mellori and C. tasmanicus respectively). The divergence between the New Zealand ravens and Australian raven clade occurred in the latest Pliocene, which coincides with the onset of glacial deforestation. We also found that the divergence between the two putative New Zealand species C. antipodum and C. moriorum probably occurred in the late Pleistocene making their separation as species untenable. Consequently, we consider Corax antipodum (Forbes, 1893) to be a subspecies of Corvus moriorum Forbes, 1892. We re-examine the osteological evidence that led 19th century researchers to assign the New Zealand taxa to a separate genus, and re-assess these features in light of our new phylogenetic hypotheses. Like previous researchers, we conclude that the morphology of the palate of C. moriorum is unique among the genus Corvus, and suggest this may be an adaptation for a specialist diet.

Evolution of heterogeneous genome differentiation across multiple contact zones in a crow species complex.

Uncovering the genetic basis of species diversification is a central goal in evolutionary biology. Yet, the link between the accumulation of genomic changes during population divergence and the evolutionary forces promoting reproductive isolation is poorly understood. Here, we analysed 124 genomes of crow populations with various degrees of genome-wide differentiation, with parallelism of a sexually selected plumage phenotype, and ongoing hybridization. Overall, heterogeneity in genetic differentiation along the genome was best explained by linked selection exposed on a shared genome architecture. Superimposed on this common background, we identified genomic regions with signatures of selection specific to independent phenotypic contact zones. Candidate pigmentation genes with evidence for divergent selection were only partly shared, suggesting context-dependent selection on a multigenic trait architecture and parallelism by pathway rather than by repeated single-gene effects. This study provides insight into how various forms of selection shape genome-wide patterns of genomic differentiation as populations diverge.