Assembly and annotation of 2 high-quality columbid reference genomes from sequencing of a Columba livia × Columba guinea F1 hybrid.

Pigeons and doves (family Columbidae) are one of the most diverse extant avian lineages, and many species have served as key models for evolutionary genomics, developmental biology, physiology, and behavioral studies. Building genomic resources for columbids is essential to further many of these studies. Here, we present high-quality genome assemblies and annotations for 2 columbid species, Columba livia and Columba guinea. We simultaneously assembled C. livia and C. guinea genomes from long-read sequencing of a single F1 hybrid individual. The new C. livia genome assembly (Cliv_3) shows improved completeness and contiguity relative to Cliv_2.1, with an annotation incorporating long-read IsoSeq data for more accurate gene models. Intensive selective breeding of C. livia has given rise to hundreds of breeds with diverse morphological and behavioral characteristics, and Cliv_3 offers improved tools for mapping the genomic architecture of interesting traits. The C. guinea genome assembly is the first for this species and is a new resource for avian comparative genomics. Together, these assemblies and annotations provide improved resources for functional studies of columbids and avian comparative genomics in general.

Assembly and annotation of two high-quality columbid reference genomes from sequencing of a Columba livia x Columba guinea F1 hybrid.

Pigeons and doves (family Columbidae) are one of the most diverse extant avian lineages, and many species have served as key models for evolutionary genomics, developmental biology, physiology, and behavioral studies. Building genomic resources for colubids is essential to further many of these studies. Here, we present high-quality genome assemblies and annotations for two columbid species, Columba livia and C. guinea. We simultaneously assembled C. livia and C. guinea genomes from long-read sequencing of a single F1 hybrid individual. The new C. livia genome assembly (Cliv_3) shows improved completeness and contiguity relative to Cliv_2.1, with an annotation incorporating long-read IsoSeq data for more accurate gene models. Intensive selective breeding of C. livia has given rise to hundreds of breeds with diverse morphological and behavioral characteristics, and Cliv_3 offers improved tools for mapping the genomic architecture of interesting traits. The C. guinea genome assembly is the first for this species and is a new resource for avian comparative genomics. Together, these assemblies and annotations provide improved resources for functional studies of columbids and avian comparative genomics in general.

An allelic series at the EDNRB2 locus controls diverse piebalding patterns in the domestic pigeon.

Variation in pigment patterns within and among vertebrate species reflects underlying changes in cell migration and function that can impact health, reproductive success, and survival. The domestic pigeon (Columba livia) is an exceptional model for understanding the genetic changes that give rise to diverse pigment patterns, as selective breeding has given rise to hundreds of breeds with extensive variation in plumage color and pattern. Here, we map the genetic architecture of a suite of pigmentation phenotypes known as piebalding. Piebalding is characterized by patches of pigmented and non-pigmented feathers, and these plumage patterns are often breed-specific and stable across generations. Using a combination of quantitative trait locus mapping in F2 laboratory crosses and genome-wide association analysis, we identify a locus associated with piebalding across many pigeon breeds. This shared locus harbors a candidate gene, EDNRB2, that is a known regulator of pigment cell migration, proliferation, and survival. We discover multiple distinct haplotypes at the EDNRB2 locus in piebald pigeons, which include a mix of protein-coding, noncoding, and structural variants that are associated with depigmentation in specific plumage regions. These results identify a role for EDNRB2 in pigment patterning in the domestic pigeon, and highlight how repeated selection at a single locus can generate a diverse array of stable and heritable pigment patterns.

An allelic series at the EDNRB2 locus controls diverse piebalding patterns in the domestic pigeon.

Variation in pigment patterns within and among vertebrate species reflects underlying changes in cell migration and function that can impact health, reproductive success, and survival. The domestic pigeon (Columba livia) is an exceptional model for understanding the genetic changes that give rise to diverse pigment patterns, as selective breeding has given rise to hundreds of breeds with extensive variation in plumage color and pattern. Here, we map the genetic architecture of a suite of pigmentation phenotypes known as piebalding. Piebalding is characterized by patches of pigmented and non-pigmented feathers, and these plumage patterns are often breed-specific and stable across generations. Using a combination of quantitative trait locus mapping in F2 laboratory crosses and genome-wide association analysis, we identify a locus associated with piebalding across many pigeon breeds. This shared locus harbors a candidate gene, EDNRB2, that is a known regulator of pigment cell migration, proliferation, and survival. We discover multiple distinct haplotypes at the EDNRB2 locus in piebald pigeons, which include a mix of protein-coding, noncoding, and structural variants that are associated with depigmentation in specific plumage regions. These results identify a role for EDNRB2 in pigment patterning in the domestic pigeon, and highlight how repeated selection at a single locus can generate a diverse array of stable and heritable pigment patterns.

Activation of Xist by an evolutionarily conserved function of KDM5C demethylase.

XX female and XY male therian mammals equalize X-linked gene expression through the mitotically-stable transcriptional inactivation of one of the two X chromosomes in female somatic cells. Here, we describe an essential function of the X-linked homolog of an ancestral X-Y gene pair, Kdm5c-Kdm5d, in the expression of Xist lncRNA, which is required for stable X-inactivation. Ablation of Kdm5c function in females results in a significant reduction in Xist RNA expression. Kdm5c encodes a demethylase that enhances Xist expression by converting histone H3K4me2/3 modifications into H3K4me1. Ectopic expression of mouse and human KDM5C, but not the Y-linked homolog KDM5D, induces Xist in male mouse embryonic stem cells (mESCs). Similarly, marsupial (opossum) Kdm5c but not Kdm5d also upregulates Xist in male mESCs, despite marsupials lacking Xist, suggesting that the KDM5C function that activates Xist in eutherians is strongly conserved and predates the divergence of eutherian and metatherian mammals. In support, prototherian (platypus) Kdm5c also induces Xist in male mESCs. Together, our data suggest that eutherian mammals co-opted the ancestral demethylase KDM5C during sex chromosome evolution to upregulate Xist for the female-specific induction of X-inactivation.

Complex genetic architecture of three-dimensional craniofacial shape variation in domestic pigeons.

Deciphering the genetic basis of vertebrate craniofacial variation is a longstanding biological problem with broad implications in evolution, development, and human pathology. One of the most stunning examples of craniofacial diversification is the adaptive radiation of birds, in which the beak serves essential roles in virtually every aspect of their life histories. The domestic pigeon (Columba livia) provides an exceptional opportunity to study the genetic underpinnings of craniofacial variation because of its unique balance of experimental accessibility and extraordinary phenotypic diversity within a single species. We used traditional and geometric morphometrics to quantify craniofacial variation in an F2 laboratory cross derived from the straight-beaked Pomeranian Pouter and curved-beak Scandaroon pigeon breeds. Using a combination of genome-wide quantitative trait locus scans and multi-locus modeling, we identified a set of genetic loci associated with complex shape variation in the craniofacial skeleton, including beak shape, braincase shape, and mandible shape. Some of these loci control coordinated changes between different structures, while others explain variation in the size and shape of specific skull and jaw regions. We find that in domestic pigeons, a complex blend of both independent and coupled genetic effects underlie three-dimensional craniofacial morphology.

A ROR2 coding variant is associated with craniofacial variation in domestic pigeons.

Vertebrate craniofacial morphogenesis is a highly orchestrated process that is directed by evolutionarily conserved developmental pathways.1,2 Within species, canalized development typically produces modest morphological variation. However, as a result of millennia of artificial selection, the domestic pigeon displays radical craniofacial variation within a single species. One of the most striking cases of pigeon craniofacial variation is the short-beak phenotype, which has been selected in numerous breeds. Classical genetic experiments suggest that pigeon beak length is regulated by a small number of genetic factors, one of which is sex linked (Ku2 locus).3-5 However, the genetic underpinnings of pigeon craniofacial variation remain unknown. Using geometric morphometrics and quantitative trait locus (QTL) mapping on an F2 intercross between a short-beaked Old German Owl (OGO) and a medium-beaked Racing Homer (RH), we identified a single Z chromosome locus that explains a majority of the variation in beak morphology in the F2 population. Complementary comparative genomic analyses revealed that the same locus is strongly differentiated between breeds with short and medium beaks. Within the Ku2 locus, we identified an amino acid substitution in the non-canonical Wnt receptor ROR2 as a putative regulator of pigeon beak length. The non-canonical Wnt pathway serves critical roles in vertebrate neural crest cell migration and craniofacial morphogenesis.6,7 In humans, ROR2 mutations cause Robinow syndrome, a congenital disorder characterized by skeletal abnormalities, including a widened and shortened facial skeleton.8,9 Our results illustrate how the extraordinary craniofacial variation among pigeons can reveal genetic regulators of vertebrate craniofacial diversity.

Two Genomic Loci Control Three Eye Colors in the Domestic Pigeon (Columba livia).

The iris of the eye shows striking color variation across vertebrate species, and may play important roles in crypsis and communication. The domestic pigeon (Columba livia) has three common iris colors, orange, pearl (white), and bull (dark brown), segregating in a single species, thereby providing a unique opportunity to identify the genetic basis of iris coloration. We used comparative genomics and genetic mapping in laboratory crosses to identify two candidate genes that control variation in iris color in domestic pigeons. We identified a nonsense mutation in the solute carrier SLC2A11B that is shared among all pigeons with pearl eye color, and a locus associated with bull eye color that includes EDNRB2, a gene involved in neural crest migration and pigment development. However, bull eye is likely controlled by a heterogeneous collection of alleles across pigeon breeds. We also found that the EDNRB2 region is associated with regionalized plumage depigmentation (piebalding). Our study identifies two candidate genes for eye colors variation, and establishes a genetic link between iris and plumage color, two traits that vary widely in the evolution of birds and other vertebrates.

A copy number variant is associated with a spectrum of pigmentation patterns in the rock pigeon (Columba livia).

Rock pigeons (Columba livia) display an extraordinary array of pigment pattern variation. One such pattern, Almond, is characterized by a variegated patchwork of plumage colors that are distributed in an apparently random manner. Almond is a sex-linked, semi-dominant trait controlled by the classical Stipper (St) locus. Heterozygous males (ZStZ+ sex chromosomes) and hemizygous Almond females (ZStW) are favored by breeders for their attractive plumage. In contrast, homozygous Almond males (ZStZSt) develop severe eye defects and often lack plumage pigmentation, suggesting that higher dosage of the mutant allele is deleterious. To determine the molecular basis of Almond, we compared the genomes of Almond pigeons to non-Almond pigeons and identified a candidate St locus on the Z chromosome. We found a copy number variant (CNV) within the differentiated region that captures complete or partial coding sequences of four genes, including the melanosome maturation gene Mlana. We did not find fixed coding changes in genes within the CNV, but all genes are misexpressed in regenerating feather bud collar cells of Almond birds. Notably, six other alleles at the St locus are associated with depigmentation phenotypes, and all exhibit expansion of the same CNV. Structural variation at St is linked to diversity in plumage pigmentation and gene expression, and thus provides a potential mode of rapid phenotypic evolution in pigeons.

Experimental Analysis of Imprinted Mouse X-Chromosome Inactivation.

X-chromosome inactivation is a dosage compensation mechanism that equalizes X-linked gene expression between male and female mammals through the transcriptional silencing of most genes on one of the two X-chromosomes in females. With a few key exceptions, once the X-chromosome is inactivated replicated copies of that X-chromosome are maintained as inactive in all descendant cells. X-inactivation is therefore a paradigm of epigenetic inheritance. Imprinted X-inactivation is a specialized form of X-inactivation that results in the silencing of the paternally derived X-chromosome. Due to its parent-of-origin-specific pattern of inactivation, imprinted X-inactivation is a model of mitotic as well as meiotic, i.e., transgenerational, epigenetic inheritance. All cells of the early mouse embryo undergo imprinted X-inactivation, a pattern that is subsequently maintained in extraembryonic cell types in vivo and in vitro. Here, we describe both high- and low-throughput approaches to interrogate imprinted X-inactivation in the mouse embryo as well in cultured extraembryonic stem cells.

Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon.

Birds and other vertebrates display stunning variation in pigmentation patterning, yet the genes controlling this diversity remain largely unknown. Rock pigeons (Columba livia) are fundamentally one of four color pattern phenotypes, in decreasing order of melanism: T-check, checker, bar (ancestral), or barless. Using whole-genome scans, we identified NDP as a candidate gene for this variation. Allele-specific expression differences in NDP indicate cis-regulatory divergence between ancestral and melanistic alleles. Sequence comparisons suggest that derived alleles originated in the speckled pigeon (Columba guinea), providing a striking example of introgression. In contrast, barless rock pigeons have an increased incidence of vision defects and, like human families with hereditary blindness, carry start-codon mutations in NDP. In summary, we find that both coding and regulatory variation in the same gene drives wing pattern diversity, and post-domestication introgression supplied potentially advantageous melanistic alleles to feral populations of this ubiquitous urban bird.

Improved Genome Assembly and Annotation for the Rock Pigeon (Columba livia).

The domestic rock pigeon (Columba livia) is among the most widely distributed and phenotypically diverse avian species. C. livia is broadly studied in ecology, genetics, physiology, behavior, and evolutionary biology, and has recently emerged as a model for understanding the molecular basis of anatomical diversity, the magnetic sense, and other key aspects of avian biology. Here we report an update to the C. livia genome reference assembly and gene annotation dataset. Greatly increased scaffold lengths in the updated reference assembly, along with an updated annotation set, provide improved tools for evolutionary and functional genetic studies of the pigeon, and for comparative avian genomics in general.

PRC2 represses transcribed genes on the imprinted inactive X chromosome in mice.

BACKGROUND: Polycomb repressive complex 2 (PRC2) catalyzes histone H3K27me3, which marks many transcriptionally silent genes throughout the mammalian genome. Although H3K27me3 is associated with silenced gene expression broadly, it remains unclear why some but not other PRC2 target genes require PRC2 and H3K27me3 for silencing. RESULTS: Here we define the transcriptional and chromatin features that predict which PRC2 target genes require PRC2/H3K27me3 for silencing by interrogating imprinted mouse X-chromosome inactivation. H3K27me3 is enriched at promoters of silenced genes across the inactive X chromosome. To abrogate PRC2 function, we delete the core PRC2 protein EED in F1 hybrid trophoblast stem cells (TSCs), which undergo imprinted inactivation of the paternally inherited X chromosome. Eed -/- TSCs lack H3K27me3 and Xist lncRNA enrichment on the inactive X chromosome. Despite the absence of H3K27me3 and Xist RNA, only a subset of the inactivated X-linked genes is derepressed in Eed -/- TSCs. Unexpectedly, in wild-type (WT) TSCs these genes are transcribed and are enriched for active chromatin hallmarks on the inactive-X, including RNA PolII, H3K27ac, and H3K36me3, but not the bivalent mark H3K4me2. By contrast, PRC2 targets that remain repressed in Eed -/- TSCs are depleted for active chromatin characteristics in WT TSCs. CONCLUSIONS: A comparative analysis of transcriptional and chromatin features of inactive X-linked genes in WT and Eed -/- TSCs suggests that PRC2 acts as a brake to prevent induction of transcribed genes on the inactive X chromosome, a mode of PRC2 function that may apply broadly.

Sex-specific silencing of X-linked genes by Xist RNA.

X-inactive specific transcript (Xist) long noncoding RNA (lncRNA) is thought to catalyze silencing of X-linked genes in cis during X-chromosome inactivation, which equalizes X-linked gene dosage between male and female mammals. To test the impact of Xist RNA on X-linked gene silencing, we ectopically induced endogenous Xist by ablating the antisense repressor Tsix in mice. We find that ectopic Xist RNA induction and subsequent X-linked gene silencing is sex specific in embryos and in differentiating embryonic stem cells (ESCs) and epiblast stem cells (EpiSCs). A higher frequency of X(ΔTsix)Y male cells displayed ectopic Xist RNA coating compared with X(ΔTsix)X female cells. This increase reflected the inability of X(ΔTsix)Y cells to efficiently silence X-linked genes compared with X(ΔTsix)X cells, despite equivalent Xist RNA induction and coating. Silencing of genes on both Xs resulted in significantly reduced proliferation and increased cell death in X(ΔTsix)X female cells relative to X(ΔTsix)Y male cells. Thus, whereas Xist RNA can inactivate the X chromosome in females it may not do so in males. We further found comparable silencing in differentiating X(ΔTsix)Y and 39,X(ΔTsix) (X(ΔTsix)O) ESCs, excluding the Y chromosome and instead implicating the X-chromosome dose as the source of the sex-specific differences. Because X(ΔTsix)X female embryonic epiblast cells and EpiSCs harbor an inactivated X chromosome prior to ectopic inactivation of the active X(ΔTsix) X chromosome, we propose that the increased expression of one or more X-inactivation escapees activates Xist and, separately, helps trigger X-linked gene silencing.

Visualizing Long Noncoding RNAs on Chromatin.

Fluorescence in situ hybridization (FISH) enables the detection of specific nucleic acid sequences within single cells. For example, RNA FISH provides information on both the expression level and localization of RNA transcripts and, when combined with detection of associated proteins and chromatin modifications, can lend essential insights into long noncoding RNA (lncRNA) function. Epigenetic effects have been postulated for many lncRNAs, but shown for only a few. Advances in in situ techniques and microscopy, however, now allow for visualization of lncRNAs that are expressed at very low levels or are not very stable. FISH-based detections of RNA and DNA coupled with immunological staining of proteins/histone modifications offer the possibility to connect lncRNAs to epigenetic effects. Here, we describe an integrated set of protocols to detect, individually or in combination, specific RNAs, DNAs, proteins, and histone modifications in single cells at a high level of sensitivity using conventional fluorescence microscopy.

An Xist-activating antisense RNA required for X-chromosome inactivation.

The transcriptional imbalance due to the difference in the number of X chromosomes between male and female mammals is remedied through X-chromosome inactivation, the epigenetic transcriptional silencing of one of the two X chromosomes in females. The X-linked Xist long non-coding RNA functions as an X inactivation master regulator; Xist is selectively upregulated from the prospective inactive X chromosome and is required in cis for X inactivation. Here we discover an Xist antisense long non-coding RNA, XistAR (Xist Activating RNA), which is encoded within exon 1 of the mouse Xist gene and is transcribed only from the inactive X chromosome. Selective truncation of XistAR, while sparing the overlapping Xist RNA, leads to a deficiency in Xist RNA expression in cis during the initiation of X inactivation. Thus, the Xist gene carries within its coding sequence an antisense RNA that drives Xist expression.

A Primary Role for the Tsix lncRNA in Maintaining Random X-Chromosome Inactivation.

Differentiating pluripotent epiblast cells in eutherians undergo random X-inactivation, which equalizes X-linked gene expression between the sexes by silencing one of the two X-chromosomes in females. Tsix RNA is believed to orchestrate the initiation of X-inactivation, influencing the choice of which X remains active by preventing expression of the antisense Xist RNA, which is required to silence the inactive-X. Here we profile X-chromosome activity in Tsix-mutant (X(ΔTsix)) mouse embryonic epiblasts, epiblast stem cells, and embryonic stem cells. Unexpectedly, we find that Xist is stably repressed on the X(ΔTsix) in both sexes in undifferentiated epiblast cells in vivo and in vitro, resulting in stochastic X-inactivation in females despite Tsix-heterozygosity. Tsix is instead required to silence Xist on the active-X as epiblast cells differentiate in both males and females. Thus, Tsix is not required at the onset of random X-inactivation; instead, it protects the active-X from ectopic silencing once X-inactivation has commenced.

Differentiation-dependent requirement of Tsix long non-coding RNA in imprinted X-chromosome inactivation.

Imprinted X-inactivation is a paradigm of mammalian transgenerational epigenetic regulation resulting in silencing of genes on the paternally inherited X-chromosome. The preprogrammed fate of the X-chromosomes is thought to be controlled in cis by the parent-of-origin-specific expression of two opposing long non-coding RNAs, Tsix and Xist, in mice. Exclusive expression of Tsix from the maternal-X has implicated it as the instrument through which the maternal germline prevents inactivation of the maternal-X in the offspring. Here, we show that Tsix is dispensable for inhibiting Xist and X-inactivation in the early embryo and in cultured stem cells of extra-embryonic lineages. Tsix is instead required to prevent Xist expression as trophectodermal progenitor cells differentiate. Despite induction of wild-type Xist RNA and accumulation of histone H3-K27me3, many Tsix-mutant X-chromosomes fail to undergo ectopic X-inactivation. We propose a novel model of lncRNA function in imprinted X-inactivation that may also apply to other genomically imprinted loci.

Long nonoding RNAs in the X-inactivation center.

The X-inactivation center is a hotbed of functional long noncoding RNAs in eutherian mammals. These RNAs are thought to help orchestrate the epigenetic transcriptional states of the two X-chromosomes in females as well as of the single X-chromosome in males. To balance X-linked gene expression between the sexes, females undergo transcriptional silencing of most genes on one of the two X-chromosomes in a process termed X-chromosome inactivation. While one X-chromosome is inactivated, the other X-chromosome remains active. Moreover, with a few notable exceptions, the originally established epigenetic transcriptional profiles of the two X-chromosomes is maintained as such through many rounds of cell division, essentially for the life of the organism. The stable and divergent transcriptional fates of the two X-chromosomes, despite residing in a shared nucleoplasm, make X-inactivation a paradigm of epigenetic transcriptional regulation. Originally proposed in 1961 by Mary Lyon, the X-inactivation hypothesis has been validated through much experimentation. In the last 25 years, the discovery and functional characterization has firmly established X-linked long noncoding RNAs as key players in choreographing X-chromosome inactivation.