Genetic diversity, demographic history and neo-sex chromosomes in the Critically Endangered Raso lark

Small effective population sizes could expose island species to inbreeding and loss of genetic variation. Here, we investigate factors shaping genetic diversity in the Raso lark, which has been restricted to a single islet for approximately 500 years, with a population size of a few hundred. We assembled a reference genome for the related Eurasian skylark and then assessed diversity and demographic history using RAD-seq data (75 samples from Raso larks and two related mainland species). We first identify broad tracts of suppressed recombination in females, indicating enlarged neo-sex chromosomes. We then show that genetic diversity across autosomes in the Raso lark is lower than in its mainland relatives, but inconsistent with long-term persistence at its current population size. Finally, we find that genetic signatures of the recent population contraction are overshadowed by an ancient expansion and persistence of a very large population until the human settlement of Cape Verde. Our findings show how genome-wide approaches to study endangered species can help avoid confounding effects of genome architecture on diversity estimates, and how present-day diversity can be shaped by ancient demographic events

De Novo Assembly of the Northern Cardinal (Cardinalis cardinalis) Genome Reveals Candidate Regulatory Regions for Sexually Dichromatic Red Plumage Coloration

Northern cardinals (Cardinalis cardinalis) are common, mid-sized passerines widely distributed in North America. As an iconic species with strong sexual dichromatism, it has been the focus of extensive ecological and evolutionary research, yet genomic studies investigating the evolution of genotype–phenotype association of plumage coloration and dichromatism are lacking. Here we present a new, highly-contiguous assembly for C. cardinalis. We generated a 1.1 Gb assembly comprised of 4,762 scaffolds, with a scaffold N50 of 3.6 Mb, a contig N50 of 114.4 kb and a longest scaffold of 19.7 Mb. We identified 93.5% complete and single-copy orthologs from an Aves dataset using BUSCO, demonstrating high completeness of the genome assembly. We annotated the genomic region comprising the CYP2J19 gene, which plays a pivotal role in the red coloration in birds. Comparative analyses demonstrated non-exonic regions unique to the CYP2J19 gene in passerines and a long insertion upstream of the gene in C. cardinalis. Transcription factor binding motifs discovered in the unique insertion region in C. cardinalis suggest potential androgen-regulated mechanisms underlying sexual dichromatism. Pairwise Sequential Markovian Coalescent (PSMC) analysis of the genome reveals fluctuations in historic effective population size between 100,000–250,000 in the last 2 millions years, with declines concordant with the beginning of the Pleistocene epoch and Last Glacial Period. This draft genome of C. cardinalis provides an important resource for future studies of ecological, evolutionary, and functional genomics in cardinals and other birds

Olfactory receptor subgenome and expression in a highly olfactory procellariiform seabird.

Procellariiform seabirds rely on their sense of smell for foraging and homing. Both genomes and transcriptomes yield important clues about how olfactory receptor (OR) subgenomes are shaped by natural and sexual selection, yet no transcriptomes have been made of any olfactory epithelium of any bird species thus far. Here, we assembled a high-quality genome and nasal epithelium transcriptome of the Leach's storm-petrel (Oceanodroma leucorhoa) to extensively characterize their OR repertoire. Using a depth-of-coverage-assisted counting method, we estimated over 160 intact OR genes (∼500 including OR fragments). This method reveals the highest number of intact OR genes and the lowest proportion of pseudogenes compared to other waterbirds studied, and suggests that rates of OR gene duplication vary between major clades of birds, with particularly high rates in passerines. OR expression patterns reveal two OR genes (OR6-6 and OR5-11) highly expressed in adults, and four OR genes (OR14-14, OR14-12, OR10-2, and OR14-9) differentially expressed between age classes of storm-petrels. All four genes differentially expressed between age classes were more highly expressed in chicks compared to adults, suggesting that OR genes may exhibit ontogenetic specializations. Three highly differentially expressed OR genes also had high copy number ratios, suggesting that expression variation may be linked to copy number in the genome. We provide better estimates of OR gene number by using a copy number-assisted counting method, and document ontogenetic changes in OR gene expression that may be linked to olfactory specialization. These results provide valuable insight into the expression, development, and macroevolution of olfaction in seabirds

Feather chemicals contain information about the major histocompatibility complex in a highly scented seabird.

Mate choice informed by the immune genes of the major histocompatibility complex (MHC) may provide fitness benefits including offspring with increased immunocompetence. Olfactory cues are considered the primary mechanism organisms use to evaluate the MHC of potential mates, yet this idea has received limited attention in birds. Motivated by a finding of MHC-dependent mate choice in the Leach's storm-petrel (Oceanodroma leucorhoa), we examined whether the chemical profiles of this highly scented seabird contain information about MHC genes. Whereas previous studies in birds examined non-volatile compounds, we used gas chromatography-mass spectrometry to measure the volatile compounds emitted from feathers that potentially serve as olfactory infochemicals about MHC and coupled this with locus-specific genotyping of MHC IIB genes. We found that feather chemicals reflected individual MHC diversity through interactions with sex and breeding status. Furthermore, similarity in MHC genotype was correlated with similarity in chemical profiles within female-female and male-female dyads. We provide the first evidence that volatile chemicals from bird feathers can encode information about the MHC. Our findings suggest that olfaction likely aids MHC-based mate choice in this species and highlight a role for chemicals in mediating genetic mate choice in birds where this mode of communication has been largely overlooked

Data from: Ecology can inform genetics: disassortative mating contributes to MHC polymorphism in Leach’s storm-petrels (Oceanodroma leucorhoa)

Studies of MHC-based mate choice in wild populations generally test hypotheses by assuming female choice and male-male competition, whether or not mate choice dynamics have been previously determined for the species under study. Here we examined mate choice patterns in a small burrow-nesting seabird, the Leach’s storm-petrel (Oceanodroma leucorhoa), using the Major Histocompatibility Complex (MHC). The life history and ecology of this species is extreme: both partners work together to fledge a single chick during the breeding season, a task that requires regularly traveling hundreds of kilometers to and from foraging grounds over a six to eight-week provisioning period. Using a five-year dataset unprecedented for this species (n=1027 adults and 925 chicks), we found a positive relationship between variation in female reproductive quality and heterozygosity at Ocle-DAB2, a MHC class IIB locus. Contrary to previous reports rejecting disassortative mating as a mechanism for maintaining genetic polymorphism in this species, here we show that males make significant disassortative mate choice decisions. Variability in female reproductive success suggests that the most common homozygous females (Ocle-DAB201 / Ocle-DAB201) may be physiologically disadvantaged and, therefore, less preferred as lifelong partners for choosy males. The results from this study support the role of mate choice in maintaining high levels of MHC variability in a wild seabird species, and highlights the need to incorporate a broader ecological framework and sufficient sample sizes into studies of MHC-based mating patterns in wild populations in general

Data from: Neo-sex chromosomes and demography shape genetic diversity in the critically endangered Raso lark

Generally small effective population sizes expose island species to inbreeding and loss of genetic variation. The Raso lark has been restricted to a single islet for ~500 years, with a population size of a few hundred. To investigate the factors shaping genetic diversity in the species, we assembled a reference genome for the related Eurasian skylark and then assessed genomic diversity and demographic history using RAD-seq data (26 Raso lark samples and 52 samples from its two most closely related mainland species). Genetic diversity in the Raso lark is lower than in its mainland relatives, but is nonetheless considerably higher than anticipated given its recent population size. This is partly explained by an unusual and dramatic effect of enlarged neo-sex chromosomes, which preserve high heterozygosity across 13% of the genome in females, and account for half of the overall genetic diversity in the population. In addition, by reconstructing past demography we find that genetic signatures of the recent population contraction are overshadowed by an ancient expansion and persistence of a very large population until the human settlement of Cape Verde. Nevertheless, relatedness analyses suggest that the population is at risk of inbreeding depression. Our findings are particularly important in that they reveal the hidden effects of genome architecture in shaping diversity estimates, and hence demonstrate the value of a reference genome and population genomic analyses over conventional metrics to study diversity in non-model and endangered species

DAB2_Amino_Acid_Matrix

Number of amino acid substitutions between Ocle-DAB2 allele pairings

Adult_morphometric_data

Year, mass, tarsus length, wing chord, for adults used in morphometrics*genotype analyse

Stairwayplot best model

Output from Stairwayplot for the best fitting model of population size over time for Raso lark