Recurrent Evolution of Melanism in South American Felids

Morphological variation in natural populations is a genomic test bed for studying the interface between molecular evolution and population genetics, but some of the most interesting questions involve non-model organisms that lack well annotated reference genomes. Many felid species exhibit polymorphism for melanism but the relative roles played by genetic drift, natural selection, and interspecies hybridization remain uncertain. We identify mutations of Agouti signaling protein (ASIP) or the Melanocortin 1 receptor (MC1R) as independent causes of melanism in three closely related South American species: the pampas cat (Leopardus colocolo), the kodkod (Leopardus guigna), and Geoffroy’s cat (Leopardus geoffroyi). To assess population level variation in the regions surrounding the causative mutations we apply genomic resources from the domestic cat to carry out clone-based capture and targeted resequencing of 299 kb and 251 kb segments that contain ASIP and MC1R, respectively, from 54 individuals (13–21 per species), achieving enrichment of ~500–2500-fold and ~150x coverage. Our analysis points to unique evolutionary histories for each of the three species, with a strong selective sweep in the pampas cat, a distinctive but short melanism-specific haplotype in the Geoffroy’s cat, and reduced nucleotide diversity for both ancestral and melanism-bearing chromosomes in the kodkod. These results reveal an important role for natural selection in a trait of longstanding interest to ecologists, geneticists, and the lay community, and provide a platform for comparative studies of morphological variation in other natural populations.     

A Deep-Learning Approach for Inference of Selective Sweeps from the Ancestral Recombination Graph

Detecting signals of selection from genomic data is a central problem in population genetics. Coupling the rich information in the ancestral recombination graph (ARG) with a powerful and scalable deep-learning framework, we developed a novel method to detect and quantify positive selection: Selection Inference using the Ancestral recombination graph (SIA). Built on a Long Short-Term Memory (LSTM) architecture, a particular type of a Recurrent Neural Network (RNN), SIA can be trained to explicitly infer a full range of selection coefficients, as well as the allele frequency trajectory and time of selection onset. We benchmarked SIA extensively on simulations under a European human demographic model, and found that it performs as well or better as some of the best available methods, including state-of-the-art machine-learning and ARG-based methods. In addition, we used SIA to estimate selection coefficients at several loci associated with human phenotypes of interest. SIA detected novel signals of selection particular to the European (CEU) population at the MC1R and ABCC11 loci. In addition, it recapitulated signals of selection at the LCT locus and several pigmentation-related genes. Finally, we reanalyzed polymorphism data of a collection of recently radiated southern capuchino seedeater taxa in the genus Sporophila to quantify the strength of selection and improved the power of our previous methods to detect partial soft sweeps. Overall, SIA uses deep learning to leverage the ARG and thereby provides new insight into how selective sweeps shape genomic diversity.     

Mutations in ASIP and MC1R: dominant black and recessive black alleles segregate in native Swedish sheep populations

By studying genes associated with coat colour, we can understand the role of these genes in pigmentation but also gain insight into selection history. North European short‐tailed sheep, including Swedish breeds, have variation in their coat colour, making them good models to expand current knowledge of mutations associated with coat colour in sheep. We studied ASIP and MC1R, two genes with known roles in pigmentation, and their association with black coat colour. We did this by sequencing the coding regions of ASIP in 149 animals and MC1R in 129 animals from seven native Swedish sheep breeds in individuals with black, white or grey fleece. Previously known mutations in ASIP [recessive black allele: g.100_105del (D₅) and/or g.5172T>A] were associated with black coat colour in Klövsjö and Roslag sheep breeds and mutations in both ASIP and MC1R (dominant black allele: c.218T>A and/or c.361G>A) were associated with black coat colour in Swedish Finewool. In Gotland, Gute, Värmland and Helsinge sheep breeds, coat colour inheritance was more complex: only 11 of 16 individuals with black fleece had genotypes that could explain their black colour. These breeds have grey individuals in their populations, and grey is believed to be a result of mutations and allelic copy number variation within the ASIP duplication, which could be a possible explanation for the lack of a clear inheritance pattern in these breeds. Finally, we found a novel missense mutation in MC1R (c.452G>A) in Gotland, Gute and Värmland sheep and evidence of a duplication of MC1R in Gotland sheep.     

How the Leopard Hides Its Spots: ASIP Mutations and Melanism in Wild Cats

The occurrence of melanism (darkening of the background coloration) is documented in 13 felid species, in some cases reaching high frequencies at the population level. Recent analyses have indicated that it arose multiple times in the Felidae, with three different species exhibiting unique mutations associated with this trait. The causative mutations in the remaining species have so far not been identified, precluding a broader assessment of the evolutionary dynamics of melanism in the Felidae. Among these, the leopard (Panthera pardus) is a particularly important target for research, given the iconic status of the ‘black panther’ and the extremely high frequency of melanism observed in some Asian populations. Another felid species from the same region, the Asian golden cat (Pardofelis temminckii), also exhibits frequent records of melanism in some areas. We have sequenced the coding region of the Agouti Signaling Protein (ASIP) gene in multiple leopard and Asian golden cat individuals, and identified distinct mutations strongly associated with melanism in each of them. The single nucleotide polymorphism (SNP) detected among the P. pardus individuals was caused by a nonsense mutation predicted to completely ablate ASIP function. A different SNP was identified in P. temminckii, causing a predicted amino acid change that should also induce loss of function. Our results reveal two additional cases of species-specific mutations implicated in melanism in the Felidae, and indicate that ASIP mutations may play an important role in naturally-occurring coloration polymorphism.     

Melanism in guinea fowl (Numida meleagris) is associated with a deletion of Phenylalanine‐256 in the MC1R gene

We have characterized a deletion in the MC1R gene causing the loss of one amino acid (p.Phe256del), which is perfectly associated with melanism in guinea fowl (Numida meleagris). Co‐segregation of the p.Phe256del with melanism was confirmed in 25 offspring born from a cross of two heterozygote birds; therefore we suggest that this mutation is responsible for the black phenotype. Interestingly, this is the first case of recessive melanism linked to MC1R.     

The evolution of black plumage from blue in Australian fairy‐wrens (Maluridae): genetic and structural evidence

Genetic variation in the melanocortin‐1 receptor (MC1R) locus is responsible for color variation, particularly melanism, in many groups of vertebrates. Fairy‐wrens, Maluridae, are a family of Australian and New Guinean passerines with several instances of dramatic shifts in plumage coloration, both intra‐ and inter‐specifically. A number of these color changes are from bright blue to black plumage. In this study, we examined sequence variation at the MC1R locus in most genera and species of fairy‐wrens. Our primary focus was subspecies of the white‐winged fairy‐wren Malurus leucopterus in which two subspecies, each endemic to islands off the western Australian coast, are black while the mainland subspecies is blue. We found fourteen variable amino acid residues within M. leucopterus, but at only one position were alleles perfectly correlated with plumage color. Comparison with other fairy‐wren species showed that the blue mainland subspecies, not the black island subspecies, had a unique genotype. Examination of MC1R protein sequence variation across our sample of fairy‐wrens revealed no correlation between plumage color and sequence in this group. We thus conclude that amino acid changes in the MC1R locus are not directly responsible for the black plumage of the island subspecies of M. leucopterus. Our examination of the nanostructure of feathers from both black and blue subspecies of M. leucopterus and other black and blue fairy‐wren species clarifies the evolution of black plumage in this family. Our data indicate that the black white‐winged fairy‐wrens evolved from blue ancestors because vestiges of the nanostructure required for the production of blue coloration exist within their black feathers. Based on our phylogeographic analysis of M. leucopterus, in which the two black subspecies do not appear to be each other's closest relatives, we infer that there have been two independent evolutionary transitions from blue to black plumage. A third potential transition from blue to black appears to have occurred in a sister clade.     

Skin Transcriptome Profiles Associated with Skin Color in Chickens

Nutritional and medicinal benefits have been attributed to the consumption of tissues from the black-boned chickens in oriental countries. Lueyang black-boned chicken is one of the native chicken breeds. However, some birds may instead have white or lighter skin, which directly causes economic losses every year. Previous studies of pigmentation have focused on a number of genes that may play important roles in coat color regulation. Illumina2000 sequencing technology was used to catalog the global gene expression profiles in the skin of the Lueyang chicken with white versus black skin. A total of 18,608 unigenes were assembled from the reads obtained from the skin of the white and black chickens. A total of 649 known genes were differentially expressed in the black versus white chickens, with 314 genes that were up regulated and 335 genes that were down-regulated, and a total of 162 novel genes were differentially expressed in the black versus white chickens, consisting of 73 genes that were up-regulated (including 4 highly expressed genes that were expressed exclusively in the skin of the black chickens) and 89 genes that were down-regulated. There were also a total of 8 known coat-color genes expressed in previous studies (ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R). In this study, 4 of which showed greater expression in the black chickens, and several were up-regulated, such as KIT, ASIP, TYR and OCA2. To our surprise, KITLG, MITF and MC1R showed no significant difference in expression between the black- and white-skinned chickens, and the expression of TYRP1 was not detected in either skin color. The expression of ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R was validated by real-time quantitative polymerase chain reaction (qPCR), and the results of the qPCR were consistent with the RNA-seq. This study provides several candidate genes that may be associated with the development of black versus white skin. More importantly, the fact that the MC1R gene showed no significant difference in expression between the black and white chickens is of particular interest for future studies that aim to elucidate its functional role in the regulation of skin color.     

Evolution of dark colour in toucans (Ramphastidae): a case of molecular adaptation?

In the last decades, researchers have been able to determine the molecular basis of some phenotypes, to test for evidence of natural selection upon them, and to demonstrate that the same genes or genetic pathways can be associated with convergent traits. Colour traits are often subject to natural selection because even small changes in these traits can have a large effect on fitness via camouflage, sexual selection or other mechanisms. The melanocortin‐1 receptor locus (MC1R) is frequently associated with intraspecific coat colour variation in vertebrates, but it has been far harder to demonstrate that this locus is involved in adaptive interspecific colour differences. Here, we investigate the contribution of the MC1R gene to the colour diversity found in toucans (Ramphastidae). We found divergent selection on MC1R in the clade represented by the genus Ramphastos and that this coincided with the evolution of darker plumage in members of this genus. Using phylogenetically corrected correlations, we show significant and specific relationships between the rate of nonsynonymous change in MC1R (dN) and plumage darkness across Ramphastidae, and also between the rate of functionally significant amino acid changes in MC1R and plumage darkness. Furthermore, three of the seven amino acid changes in MC1R that occurred in the ancestral Ramphastos branch are associated with melanism in other birds. Taken together, our results suggest that the dark colour of Ramphastos toucans was related to nonsynonymous substitutions in MC1R that may have been subject to positive selection or to a relaxation of selective pressure. These results also demonstrate a quantitative relationship between gene and phenotype evolution, representing an example of how MC1R molecular evolution may affect macroevolution of plumage phenotypes.     

Genomic evidence for convergent evolution of a key trait underlying divergence in island birds

Reproductive isolation can be initiated by changes in one or a few key traits that prevent random mating among individuals in a population. During the early stages of speciation, when isolation is often incomplete, there will be a heterogeneous pattern of differentiation across regions of the genome between diverging populations, with loci controlling these key traits appearing the most distinct as a result of strong diversifying selection. In this study, we used Illumina‐sequenced ddRAD tags to identify genomewide patterns of differentiation in three recently diverged island populations of the Monarcha castaneiventris flycatcher of the Solomon Islands. Populations of this species have diverged in plumage colour, and these differences in plumage colour, in turn, are used in conspecific recognition and likely important in reproductive isolation. Previous candidate gene sequencing identified point mutations in MC1R and ASIP, both known pigmentation genes, to be associated with the difference in plumage colour between islands. Here, we show that background levels of genomic differentiation based on over 70,000 SNPs are extremely low between populations of distinct plumage colour, with no loci reaching the level of differentiation found in either candidate gene. Further, we found that a phylogenetic analysis based on these SNPs produced a taxonomy wherein the two melanic populations appear to have evolved convergently, rather than from a single common ancestor, in contrast to their original classification as a single subspecies. Finally, we found evidence that the pattern of low genomic differentiation is the result of both incomplete lineage sorting and gene flow between populations.     

Investigation of the role of the agouti signaling protein gene (ASIP) in coat color evolution in primates

We investigated variation in the gene encoding the agouti signaling protein (ASIP) in relation to coat color evolution in primates. We found little evidence that mutations in the coding region of ASIP have been involved in color changes among closely related primate species. Among many closely related species with differing coat color, the coding region of ASIP was identical. In two cases (Sulawesi macaque and black lion tamarin) where species with almost completely black coat color had derived point mutations in exon 4 of the ASIP coding sequence, the same mutations did not alter coloration in other mammals and so probably do not affect ASIP function. Evolutionary reconstructions of two key phenotypes that are typically related to ASIP function—transverse phaeomelanin bands on hairs and pale ventral coloration—showed that these usually evolved concurrently, suggesting that loci acting downstream of ASIP may be involved. Analysis of dN/dS ratios revealed a likely change in functional constraint on ASIP following loss of agouti-banded hairs + pale ventral coloration, particularly in catarrhine primates (humans, apes, and Old World monkeys). Together with previous results on a lack of association of coat color with MC1R variation, these results suggest that other loci probably have an important role in primate coat color evolution.     

Contrasting Mode of Evolution at a Coat Color Locus in Wild and Domestic Pigs

Despite having only begun ~10,000 years ago, the process of domestication has resulted in a degree of phenotypic variation within individual species normally associated with much deeper evolutionary time scales. Though many variable traits found in domestic animals are the result of relatively recent human-mediated selection, uncertainty remains as to whether the modern ubiquity of long-standing variable traits such as coat color results from selection or drift, and whether the underlying alleles were present in the wild ancestor or appeared after domestication began. Here, through an investigation of sequence diversity at the porcine melanocortin receptor 1 (MC1R) locus, we provide evidence that wild and domestic pig (Sus scrofa) haplotypes from China and Europe are the result of strikingly different selection pressures, and that coat color variation is the result of intentional selection for alleles that appeared after the advent of domestication. Asian and European wild boar (evolutionarily distinct subspecies) differed only by synonymous substitutions, demonstrating that camouflage coat color is maintained by purifying selection. In domestic pigs, however, each of nine unique mutations altered the amino acid sequence thus generating coat color diversity. Most domestic MC1R alleles differed by more than one mutation from the wild-type, implying a long history of strong positive selection for coat color variants, during which time humans have cherry-picked rare mutations that would be quickly eliminated in wild contexts. This pattern demonstrates that coat color phenotypes result from direct human selection and not via a simple relaxation of natural selective pressures.     

Dominant Red Coat Color in Holstein Cattle Is Associated with a Missense Mutation in the Coatomer Protein Complex,Subunit Alpha (COPA) Gene

Coat color in Holstein dairy cattle is primarily controlled by the melanocortin 1 receptor (MC1R) gene, a central determinant of black (eumelanin) vs. red/brown pheomelanin synthesis across animal species. The major MC1R alleles in Holsteins are Dominant Black (MC1R^D) and Recessive Red (MC1R^e). A novel form of dominant red coat color was first observed in an animal born in 1980. The mutation underlying this phenotype was named Dominant Red and is epistatic to the constitutively activated MC1R^D. Here we show that a missense mutation in the coatomer protein complex, subunit alpha (COPA), a gene with previously no known role in pigmentation synthesis, is completely associated with Dominant Red in Holstein dairy cattle. The mutation results in an arginine to cysteine substitution at an amino acid residue completely conserved across eukaryotes. Despite this high level of conservation we show that both heterozygotes and homozygotes are healthy and viable. Analysis of hair pigment composition shows that the Dominant Red phenotype is similar to the MC1R Recessive Red phenotype, although less effective at reducing eumelanin synthesis. RNA-seq data similarly show that Dominant Red animals achieve predominantly pheomelanin synthesis by downregulating genes normally required for eumelanin synthesis. COPA is a component of the coat protein I seven subunit complex that is involved with retrograde and cis-Golgi intracellular coated vesicle transport of both protein and RNA cargo. This suggests that Dominant Red may be caused by aberrant MC1R protein or mRNA trafficking within the highly compartmentalized melanocyte, mimicking the effect of the Recessive Red loss of function MC1R allele.     

Description of a New Galapagos Giant Tortoise Species (Chelonoidis; Testudines: Testudinidae) from Cerro Fatal on Santa Cruz Island

The taxonomy of giant Galapagos tortoises (Chelonoidis spp.) is currently based primarily on morphological characters and island of origin. Over the last decade, compelling genetic evidence has accumulated for multiple independent evolutionary lineages, spurring the need for taxonomic revision. On the island of Santa Cruz there is currently a single named species, C. porteri. Recent genetic and morphological studies have shown that, within this taxon, there are two evolutionarily and spatially distinct lineages on the western and eastern sectors of the island, known as the Reserva and Cerro Fatal populations, respectively. Analyses of DNA from natural populations and museum specimens, including the type specimen for C. porteri, confirm the genetic distinctiveness of these two lineages and support elevation of the Cerro Fatal tortoises to the rank of species. In this paper, we identify DNA characters that define this new species, and infer evolutionary relationships relative to other species of Galapagos tortoises.     

Family of melanocortin receptor (MCR) genes in mammals—mutations, polymorphisms and phenotypic effects

The melanocortin receptor gene family consists of five single-exon members, which are located on autosomes. Three genes (MC2R, MC4R and MC5R) are syntenic in the human, mouse, cattle and dog genomes, while in the pig, the syntenic group comprises MC1R, MC2R and MC5R. Two genes (MC1R and MC4R) have been extensively studied due to their function in melanogenesis (MC1R) and energy control (MC4R). Conservative organisation of these genes in five mammalian species (human, mouse, cattle, pig and dog), in terms of the encoded amino acid sequence, is higher in the case of MC4R compared to MC1R. Polymorphisms of these two genes are responsible or associated with variation of pigmentation (MC1R) and adipose tissue deposition (MC4R). Polymorphic variants in MC1R, causing coat colour variation, were described in humans and domestic mammals (cattle, horse, pig, sheep, dog), as well as farm red and arctic foxes. The MC4R gene is very polymorphic in humans and it is well known that some variants cause monogenic obesity or significantly contribute to the development of polygenic obesity. Such relationships are not so evident in domestic mammals; however, at least one missense substitution (298Asp?>?Asn) in the porcine MC4R significantly contributes, at least in some breeds, to fat tissue accumulation, feed conversion ratio and daily weight gain. Knowledge on the phenotypic effects of polymorphisms of MC2R, MC3R and MC5R in domestic mammals is scarce, probably due to the small number of reports addressing these genes. Thus, further studies focused on these genes should be undertaken.     

The 8818G allele of the agouti signaling protein (ASIP) gene is ancestral and is associated with darker skin color in African Americans

Skin color, a predictor of social interactions and risk factor for several types of cancer, is due to two contrasting forms of melanin, the darker eumelanin and lighter phaeomelanin. The lighter pigment phaeomelanin is the product of the antagonistic function of the agouti signaling protein (ASIP) on the -melanocyte stimulating hormone receptor (MC1R). Studies have shown that a single-nucleotide polymorphism (SNP) in the 3UTR of the ASIP gene is associated with dark hair and eyes; however, little is known about its role in inter-individual variation in skin color. Here we examine the relationship between the ASIP g.8818A>G SNP and skin color (M index) as assessed by reflectometry in 234 African Americans. Analyses of variance (ANOVA) were performed to evaluate the effects of ASIP genotypes, age, individual ancestry, and sex on skin color variation. Significant effects on M index variation were observed for ASIP genotypes (F(2,236)=4.37, P=0.01), ancestry (F(1,243)=37.2, P<0.001), and sex (F(1,244)=4.08, P=0.05). Subsequent analyses revealed a strong effect on M index from ASIP genotypes in African American females (P<0.001). Our study suggests that the ASIP G>A polymorphism exhibits a dominant effect leading to lighter skin color and that variation in the ASIP gene may have been one of several factors contributing to reductions in pigmentation in some populations. Further study is needed to reveal how interactions between ASIP and several other genes, such as MC1R and P, predict human pigmentation.     

Molecular population genetics of the melanic plumage polymorphism in Arctic skuas (Stercorarius parasiticus): evidence for divergent selection on plumage colour

The Arctic skua (Stercorarius parasiticus) is a classic example of an avian plumage polymorphism, with variation in melanin‐based ventral plumage coloration defining pale, intermediate and dark morphs in adults of both sexes. However, despite several decades of field research, there is an incomplete understanding of how the polymorphism in ventral plumage colour is maintained and the selective forces involved. Here, we investigate selection on a locus (MC1R) that is strongly associated with plumage colour variation in Arctic skuas using patterns of nucleotide variation and comparison to neutral loci (nuclear introns and mtDNA). We find that three linked nonsynonymous mutations in MC1R, including the single mutation described previously, are associated with plumage colour in the Arctic skua. The position of nonsynonymous mutations on a MC1R haplotype network implies that divergent selection drove the initial evolution of the colour morphs. Comparisons of F_STs of MC1R vs. nuclear introns among five skua populations differing in proportion of dark morphs along an approximate north–south cline reveal a signature of divergent selection on MC1R. In contrast, we find limited evidence for balancing selection on MC1R within populations, although the power is low. Our results provide strong evidence for both past and ongoing selection on MC1R, and, by implication, plumage colour in Arctic skuas. The results suggest that a fruitful avenue for future ecological studies will be analysis of selection on morphs in colonies at the extremes along the morph ratio cline.     

Variation in insulative feather structure in songbirds replacing each other along a tropical elevation gradient

High‐elevation organisms are expected to evolve physiological adaptations to cope with harsh environmental conditions. Yet, evidence for such adaptive differences, especially compared to closely related lowland taxa occurring along the same elevational gradient, is rare. Revisiting an anecdotal natural history observation by O. Bangs from 1899 and based on new measurements of museum specimens, we confirmed that the high‐elevation hermit wood wren (Henicorhina anachoreta) from the Sierra Nevada de Santa Marta, Colombia, has longer, more insulative feathers on the chest and back, than its lower‐elevation counterpart the grey‐breasted wood wren (H. leucophrys). However, we did not find evidence for the same specializations in subspecies of H. leucophrys that live at high elevations on other elevational gradients in the Colombian Andes, although similar adaptive solutions have arisen in separate mountain systems like the Himalayas. Adaptations in plumage may be associated with the recurrence of elevational species replacements throughout the tropics.     

Micronized Copper Wood Preservatives: Efficacy of Ion,Nano, and Bulk Copper against the Brown Rot Fungus Rhodonia placenta

Recently introduced micronized copper (MC) formulations, consisting of a nanosized fraction of basic copper (Cu) carbonate (CuCO₃·Cu(OH)₂) nanoparticles (NPs), were introduced to the market for wood protection. Cu NPs may presumably be more effective against wood-destroying fungi than bulk or ionic Cu compounds. In particular, Cu- tolerant wood-destroying fungi may not recognize NPs, which may penetrate into fungal cell walls and membranes and exert their impact. The objective of this study was to assess if MC wood preservative formulations have a superior efficacy against Cu-tolerant wood-destroying fungi due to nano effects than conventional Cu biocides. After screening a range of wood-destroying fungi for their resistance to Cu, we investigated fungal growth of the Cu-tolerant fungus Rhodonia placenta in solid and liquid media and on wood treated with MC azole (MCA). In liquid cultures we evaluated the fungal response to ion, nano and bulk Cu distinguishing the ionic and particle effects by means of the Cu²⁺ chelator ammonium tetrathiomolybdate (TTM) and measuring fungal biomass, oxalic acid production and laccase activity of R. placenta. Our results do not support the presence of particular nano effects of MCA against R. placenta that would account for an increased antifungal efficacy, but provide evidence that attribute the main effectiveness of MCA to azoles.     

Genetic analysis of the endemic island loggerhead shrike, <Emphasis Type="Italic">Lanius ludovicianus anthonyi</Emphasis>

Based on limited research, the island loggerhead shrike, Lanius ludovicianus anthonyi has been considered a distinct subspecies endemic to the northern California Channel Islands. We used mtDNA control region sequences and microsatellite genotyping to compare loggerhead shrikes from the southern California mainland (L. l. gambeli), San Clemente Island (L. l. mearnsi), and the northern islands (L. l. anthonyi). Habitats on the islands are recovering due to the removal of non-native ungulates on the islands, but may be transitioning to habitats less supportive of loggerhead shrikes, so this evaluation comes at a critical time. We utilized 96 museum specimens that were collected over a century to evaluate both spatial and temporal genetic patterns. Analysis of multilocus microsatellite genotypes indicated that historical specimens of loggerhead shrikes (collected between 1897 and 1986) from the two northern islands of Santa Rosa and Santa Cruz are genetically distinct from adjacent mainland and island shrikes. Birds from Santa Catalina Island showed mixed ancestry and did not cluster with the northern island birds. Historical specimens of L. l. mearnsi from San Clemente Island also showed mixed ancestry. Our study provides evidence that a genetically distinct form of loggerhead shrikes, L. l. anthonyi, occurred on the islands of Santa Rosa and Santa Cruz.