Gene expression across mammalian organ development.

The evolution of gene expression in mammalian organ development remains largely uncharacterized. Here we report the transcriptomes of seven organs (cerebrum, cerebellum, heart, kidney, liver, ovary and testis) across developmental time points from early organogenesis to adulthood for human, rhesus macaque, mouse, rat, rabbit, opossum and chicken. Comparisons of gene expression patterns identified correspondences of developmental stages across species, and differences in the timing of key events during the development of the gonads. We found that the breadth of gene expression and the extent of purifying selection gradually decrease during development, whereas the amount of positive selection and expression of new genes increase. We identified differences in the temporal trajectories of expression of individual genes across species, with brain tissues showing the smallest percentage of trajectory changes, and the liver and testis showing the largest. Our work provides a resource of developmental transcriptomes of seven organs across seven species, and comparative analyses that characterize the development and evolution of mammalian organs.

A single introduction of wild rabbits triggered the biological invasion of Australia.

Biological invasions are a major cause of environmental and economic disruption. While ecological factors are key determinants of their success, the role of genetics has been more challenging to demonstrate. The colonization of Australia by the European rabbit is one of the most iconic and devastating biological invasions in recorded history. Here, we show that despite numerous introductions over a 70-y period, this invasion was triggered by a single release of a few animals that spread thousands of kilometers across the continent. We found genetic support for historical accounts that these were English rabbits imported in 1859 by a settler named Thomas Austin and traced the origin of the invasive population back to his birthplace in England. We also find evidence of additional introductions that established local populations but have not spread geographically. Combining genomic and historical data we show that, contrary to the earlier introductions, which consisted mostly of domestic animals, the invasive rabbits had wild ancestry. In New Zealand and Tasmania, rabbits also became a pest several decades after being introduced. We argue that the common denominator of these invasions was the arrival of a new genotype that was better adapted to the natural environment. These findings demonstrate how the genetic composition of invasive individuals can determine the success of an introduction and provide a mechanism by which multiple introductions can be required for a biological invasion.

The role of historical biogeography in shaping colour morph diversity in the common wall lizard.

The maintenance of polymorphisms often depends on multiple selective forces, but less is known on the role of stochastic or historical processes in maintaining variation. The common wall lizard (Podarcis muralis) is a colour polymorphic species in which local colour morph frequencies are thought to be modulated by natural and sexual selection. Here, we used genome-wide single-nucleotide polymorphism data to investigate the relationships between morph composition and population biogeography at a regional scale, by comparing morph composition with patterns of genetic variation of 54 populations sampled across the Pyrenees. We found that genetic divergence was explained by geographic distance but not by environmental features. Differences in morph composition were associated with genetic and environmental differentiation, as well as differences in sex ratio. Thus, variation in colour morph frequencies could have arisen via historical events and/or differences in the permeability to gene flow, possibly shaped by the complex topography and environment. In agreement with this hypothesis, colour morph diversity was positively correlated with genetic diversity and rates of gene flow and inversely correlated with the likelihood of the occurrence of bottlenecks. Concurrently, we did not find conclusive evidence for selection in the two colour loci. As an illustration of these effects, we observed that populations with higher proportions of the rarer yellow and yellow-orange morphs had higher genetic diversity. Our results suggest that processes involving a decay in overall genetic diversity, such as reduced gene flow and/or bottleneck events have an important role in shaping population-specific morph composition via non-selective processes.

Disentangling the contemporary and historical effects of landscape on the population genomic variation of two bird species restricted to the highland forest enclaves of northeastern Brazil.

Investigating the impact of landscape features on patterns of genetic variation is crucial to understand spatially dependent evolutionary processes. Here, we assess the population genomic variation of two bird species (Conopophaga cearae and Sclerurus cearensis) through the Caatinga moist forest enclaves in northeastern Brazil. To infer the evolutionary dynamics of bird populations through the Late Quaternary, we used genome-wide polymorphism data obtained from double-digestion restriction-site-associated DNA sequencing (ddRADseq), and integrated population structure analyses, historical demography models, paleodistribution modeling, and landscape genetics analyses. We found the population differentiation among enclaves to be significantly related to the geographic distance and historical resistance across the rugged landscape. The climate changes at the end of the Pleistocene to the Holocene likely triggered synchronic population decline in all enclaves for both species. Our findings revealed that both geographic distance and historical connectivity through highlands are important factors that can explain the current patterns of genetic variation. Our results further suggest that levels of population differentiation and connectivity cannot be explained purely on the basis of contemporary environmental conditions. By combining historical demographic analyses and niche modeling predictions in a historical framework, we provide strong evidence that climate fluctuations of the Quaternary promoted population differentiation and a high degree of temporal synchrony among population size changes in both species.

Relocation to avoid costs: A hypothesis on red carotenoid-based signals based on recent CYP2J19 gene expression data.

In many vertebrates, the enzymatic oxidation of dietary yellow carotenoids generates red keto-carotenoids giving color to ornaments. The oxidase CYP2J19 is here a key effector. Its purported intracellular location suggests a shared biochemical pathway between trait expression and cell functioning. This might guarantee the reliability of red colorations as individual quality signals independent of production costs. We hypothesize that the ornament type (feathers vs. bare parts) and production costs (probably CYP2J19 activity compromising vital functions) could have promoted tissue-specific gene relocation. We review current avian tissue-specific CYP2J19 expression data. Among the ten red-billed species showing CYP2J19 bill expression, only one showed strong hepatic expression. Moreover, a phylogenetically-controlled analysis of 25 red-colored species shows that those producing red bare parts are less likely to have strong hepatic CYP2J19 expression than species with only red plumages. Thus, both production costs and shared pathways might have contributed to the evolution of red signals.

Molecular parallelisms between pigmentation in the avian iris and the integument of ectothermic vertebrates.

Birds exhibit striking variation in eye color that arises from interactions between specialized pigment cells named chromatophores. The types of chromatophores present in the avian iris are lacking from the integument of birds or mammals, but are remarkably similar to those found in the skin of ectothermic vertebrates. To investigate molecular mechanisms associated with eye coloration in birds, we took advantage of a Mendelian mutation found in domestic pigeons that alters the deposition of yellow pterin pigments in the iris. Using a combination of genome-wide association analysis and linkage information in pedigrees, we mapped variation in eye coloration in pigeons to a small genomic region of ~8.5kb. This interval contained a single gene, SLC2A11B, which has been previously implicated in skin pigmentation and chromatophore differentiation in fish. Loss of yellow pigmentation is likely caused by a point mutation that introduces a premature STOP codon and leads to lower expression of SLC2A11B through nonsense-mediated mRNA decay. There were no substantial changes in overall gene expression profiles between both iris types as well as in genes directly associated with pterin metabolism and/or chromatophore differentiation. Our findings demonstrate that SLC2A11B is required for the expression of pterin-based pigmentation in the avian iris. They further highlight common molecular mechanisms underlying the production of coloration in the iris of birds and skin of ectothermic vertebrates.

A loss-of-function mutation in RORB disrupts saltatorial locomotion in rabbits.

Saltatorial locomotion is a type of hopping gait that in mammals can be found in rabbits, hares, kangaroos, and some species of rodents. The molecular mechanisms that control and fine-tune the formation of this type of gait are unknown. Here, we take advantage of one strain of domesticated rabbits, the sauteur d'Alfort, that exhibits an abnormal locomotion behavior defined by the loss of the typical jumping that characterizes wild-type rabbits. Strikingly, individuals from this strain frequently adopt a bipedal gait using their front legs. Using a combination of experimental crosses and whole genome sequencing, we show that a single locus containing the RAR related orphan receptor B gene (RORB) explains the atypical gait of these rabbits. We found that a splice-site mutation in an evolutionary conserved site of RORB results in several aberrant transcript isoforms incorporating intronic sequence. This mutation leads to a drastic reduction of RORB-positive neurons in the spinal cord, as well as defects in differentiation of populations of spinal cord interneurons. Our results show that RORB function is required for the performance of saltatorial locomotion in rabbits.

Regulatory changes in pterin and carotenoid genes underlie balanced color polymorphisms in the wall lizard.

Reptiles use pterin and carotenoid pigments to produce yellow, orange, and red colors. These conspicuous colors serve a diversity of signaling functions, but their molecular basis remains unresolved. Here, we show that the genomes of sympatric color morphs of the European common wall lizard (Podarcis muralis), which differ in orange and yellow pigmentation and in their ecology and behavior, are virtually undifferentiated. Genetic differences are restricted to two small regulatory regions near genes associated with pterin [sepiapterin reductase (SPR)] and carotenoid [beta-carotene oxygenase 2 (BCO2)] metabolism, demonstrating that a core gene in the housekeeping pathway of pterin biosynthesis has been coopted for bright coloration in reptiles and indicating that these loci exert pleiotropic effects on other aspects of physiology. Pigmentation differences are explained by extremely divergent alleles, and haplotype analysis revealed abundant transspecific allele sharing with other lacertids exhibiting color polymorphisms. The evolution of these conspicuous color ornaments is the result of ancient genetic variation and cross-species hybridization.

Genetic Basis of De Novo Appearance of Carotenoid Ornamentation in Bare Parts of Canaries.

Unlike wild and domestic canaries (Serinus canaria), or any of the three dozen species of finches in genus Serinus, the domestic urucum breed of canaries exhibits bright red bills and legs. This novel trait offers a unique opportunity to understand the mechanisms of bare-part coloration in birds. To identify the mutation producing the colorful phenotype, we resequenced the genome of urucum canaries and performed a range of analyses to search for genotype-to-phenotype associations across the genome. We identified a nonsynonymous mutation in the gene BCO2 (beta-carotene oxygenase 2, also known as BCDO2), an enzyme involved in the cleavage and breakdown of full-length carotenoids into short apocarotenoids. Protein structural models and in vitro functional assays indicate that the urucum mutation abrogates the carotenoid-cleavage activity of BCO2. Consistent with the predicted loss of carotenoid-cleavage activity, urucum canaries tended to have increased levels of full-length carotenoid pigments in bill tissue and reduced levels of carotenoid-cleavage products (apocarotenoids) in retinal tissue compared with other breeds of canaries. We hypothesize that carotenoid-based bare-part coloration might be readily gained, modified, or lost through simple switches in the enzymatic activity or regulation of BCO2 and this gene may be an important mediator in the evolution of bare-part coloration among bird species.

Pterin-based pigmentation in animals.

Pterins are one of the major sources of bright coloration in animals. They are produced endogenously, participate in vital physiological processes and serve a variety of signalling functions. Despite their ubiquity in nature, pterin-based pigmentation has received little attention when compared to other major pigment classes. Here, we summarize major aspects relating to pterin pigmentation in animals, from its long history of research to recent genomic studies on the molecular mechanisms underlying its evolution. We argue that pterins have intermediate characteristics (endogenously produced, typically bright) between two well-studied pigment types, melanins (endogenously produced, typically cryptic) and carotenoids (dietary uptake, typically bright), providing unique opportunities to address general questions about the biology of coloration, from the mechanisms that determine how different types of pigmentation evolve to discussions on honest signalling hypotheses. Crucial gaps persist in our knowledge on the molecular basis underlying the production and deposition of pterins. We thus highlight the need for functional studies on systems amenable for laboratory manipulation, but also on systems that exhibit natural variation in pterin pigmentation. The wealth of potential model species, coupled with recent technological and analytical advances, make this a promising time to advance research on pterin-based pigmentation in animals.

Changes in brain architecture are consistent with altered fear processing in domestic rabbits.

The most characteristic feature of domestic animals is their change in behavior associated with selection for tameness. Here we show, using high-resolution brain magnetic resonance imaging in wild and domestic rabbits, that domestication reduced amygdala volume and enlarged medial prefrontal cortex volume, supporting that areas driving fear have lost volume while areas modulating negative affect have gained volume during domestication. In contrast to the localized gray matter alterations, white matter anisotropy was reduced in the corona radiata, corpus callosum, and the subcortical white matter. This suggests a compromised white matter structural integrity in projection and association fibers affecting both afferent and efferent neural flow, consistent with reduced neural processing. We propose that compared with their wild ancestors, domestic rabbits are less fearful and have an attenuated flight response because of these changes in brain architecture.

Testing the carotenoid-based sexual signalling mechanism by altering CYP2J19 gene expression and colour in a bird species.

Ornaments can evolve to reveal individual quality when their production/maintenance costs make them reliable as 'signals' or if their expression level is intrinsically linked to condition by some unfalsifiable mechanism (indices). The latter has been mostly associated with traits constrained by body size. In red ketocarotenoid-based colorations, that link could, instead, be established with cell respiration at the inner mitochondrial membrane (IMM). The production mechanism could be independent of resource (yellow carotenoids) availability, thus discarding costs linked to allocation trade-offs. A gene coding for a ketolase enzyme (CYP2J19) responsible for converting dietary yellow carotenoids to red ketocarotenoids has recently been described. We treated male zebra finches with an antioxidant designed to penetrate the IMM (mitoTEMPO) and a thyroid hormone (triiodothyronine) with known hypermetabolic effects. Among hormone controls, MitoTEMPO downregulated CYP2J19 in the bill (a red ketocarotenoid-based ornament), supporting the mitochondrial involvement in ketolase function. Both treatments interacted when increasing hormone dosage, indicating that mitochondria and thyroid metabolisms could simultaneously regulate coloration. Moreover, CYP2J19 expression was positively correlated to redness but also to yellow carotenoid levels in the blood. However, treatment effects were not annulated when controlling for blood carotenoid variability, which suggests that costs linked to resource availability could be minor.

High-density lipoprotein receptor SCARB1 is required for carotenoid coloration in birds.

Yellow, orange, and red coloration is a fundamental aspect of avian diversity and serves as an important signal in mate choice and aggressive interactions. This coloration is often produced through the deposition of diet-derived carotenoid pigments, yet the mechanisms of carotenoid uptake and transport are not well-understood. The white recessive breed of the common canary (Serinus canaria), which carries an autosomal recessive mutation that renders its plumage pure white, provides a unique opportunity to investigate mechanisms of carotenoid coloration. We carried out detailed genomic and biochemical analyses comparing the white recessive with yellow and red breeds of canaries. Biochemical analysis revealed that carotenoids are absent or at very low concentrations in feathers and several tissues of white recessive canaries, consistent with a genetic defect in carotenoid uptake. Using a combination of genetic mapping approaches, we show that the white recessive allele is due to a splice donor site mutation in the scavenger receptor B1 (SCARB1; also known as SR-B1) gene. This mutation results in abnormal splicing, with the most abundant transcript lacking exon 4. Through functional assays, we further demonstrate that wild-type SCARB1 promotes cellular uptake of carotenoids but that this function is lost in the predominant mutant isoform in white recessive canaries. Our results indicate that SCARB1 is an essential mediator of the expression of carotenoid-based coloration in birds, and suggest a potential link between visual displays and lipid metabolism.

Signatures of Selection on Standing Genetic Variation Underlie Athletic and Navigational Performance in Racing Pigeons.

Racing pigeons have been selectively bred to find their way home quickly over what are often extremely long distances. This breed is of substantial commercial value and is also an excellent avian model to gain empirical insights into the evolution of traits associated with flying performance and spatial orientation. Here, we investigate the molecular basis of the superior athletic and navigational capabilities of racing pigeons using whole-genome and RNA sequencing data. We inferred multiple signatures of positive selection distributed across the genome of racing pigeons. The strongest signature overlapped the CASK gene, a gene implicated in the formation of neuromuscular junctions. However, no diagnostic alleles were found between racing pigeons and other breeds, and only a small proportion of highly differentiated variants were exclusively detected in racing pigeons. We can thus conclude that very few individual genetic changes, if any, are either strictly necessary or sufficient for superior athletics and navigation. Gene expression analysis between racing and nonracing breeds revealed modest differences in muscle (213) and brain (29). These transcripts, however, showed only slightly elevated levels of genetic differentiation between the two groups, suggesting that most differential expression is not causative but likely a consequence of alterations in regulatory networks. Our results show that the unique suite of traits that enable fast flight, long endurance, and accurate navigation in racing pigeons, do not result from few loci acting as master switches but likely from a polygenic architecture that leveraged standing genetic variation available at the onset of the breed formation.

Maximum likelihood approach suggests positive selection in platelet integrin αIIbß3 in mammalian species.

Platelet integrin αIIbß3 is crucial for platelet aggregation. Although structural and functional characteristics of this protein have been extensively studied, the evolutionary pattern studies of this protein complex in mammals are scarce. Here, we addressed this question using maximum likelihood approaches to identify codons that are evolving under positive selection. Likelihood of positive selection was estimated using CODEML implemented in PAML software applied to integrin αIIbß3 derived from nucleotide sequences of 10 different mammalian species. Four codons in mature αIIb-subunit (corresponding to residues 150, 184, 193, and 370) and three codons in mature ß3-subunit (corresponding to residues 129, 440, and 444) showed signs of positive selection with posterior probabilities over 95%. The different amino acids observed for each of the positively selected residues detected showed different physicochemical properties. These results open new research avenues to understand the physiological importance of specific residues and should allow for a better understanding of the function and the different interactions of each residue within the mature protein.

A non-coding region near Follistatin controls head colour polymorphism in the Gouldian finch.

Discrete colour morphs coexisting within a single population are common in nature. In a broad range of organisms, sympatric colour morphs often display major differences in other traits, including morphology, physiology or behaviour. Despite the repeated occurrence of this phenomenon, our understanding of the genetics that underlie multi-trait differences and the factors that promote the long-term maintenance of phenotypic variability within a freely interbreeding population are incomplete. Here, we investigated the genetic basis of red and black head colour in the Gouldian finch (Erythrura gouldiae), a classic polymorphic system in which naturally occurring colour morphs also display differences in aggressivity and reproductive success. We show that the candidate locus is a small (approx. 70 kb) non-coding region mapping to the Z chromosome near the Follistatin (FST) gene. Unlike recent findings in other systems where phenotypic morphs are explained by large inversions containing hundreds of genes (so-called supergenes), we did not identify any structural rearrangements between the two haplotypes using linked-read sequencing technology. Nucleotide divergence between the red and black alleles was high when compared to the remainder of the Z chromosome, consistent with their maintenance as balanced polymorphisms over several million years. Our results illustrate how pleiotropic phenotypes can arise from simple genetic variation, probably regulatory in nature.

A genomic map of clinal variation across the European rabbit hybrid zone.

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole-genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein-protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome-wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.

The genomic architecture of population divergence between subspecies of the European rabbit.

The analysis of introgression of genomic regions between divergent populations provides an excellent opportunity to determine the genetic basis of reproductive isolation during the early stages of speciation. However, hybridization and subsequent gene flow must be relatively common in order to localize individual loci that resist introgression. In this study, we used next-generation sequencing to study genome-wide patterns of genetic differentiation between two hybridizing subspecies of rabbits (Oryctolagus cuniculus algirus and O. c. cuniculus) that are known to undergo high rates of gene exchange. Our primary objective was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. On the basis of 326,000 polymorphisms, we found low to moderate overall levels of differentiation between subspecies, and fewer than 200 genomic regions dispersed throughout the genome showing high differentiation consistent with a signature of reduced gene flow. Most differentiated regions were smaller than 200 Kb and contained very few genes. Remarkably, 30 regions were each found to contain a single gene, facilitating the identification of candidate genes underlying reproductive isolation. This gene-level resolution yielded several insights into the genetic basis and architecture of reproductive isolation in rabbits. Regions of high differentiation were enriched on the X-chromosome and near centromeres. Genes lying within differentiated regions were often associated with transcription and epigenetic activities, including chromatin organization, regulation of transcription, and DNA binding. Overall, our results from a naturally hybridizing system share important commonalities with hybrid incompatibility genes identified using laboratory crosses in mice and flies, highlighting general mechanisms underlying the maintenance of reproductive barriers.

A prominent role of KRAB-ZNF transcription factors in mammalian speciation?

The mechanisms of speciation have been one of the most debated topics in evolutionary biology. Among all reproductive barriers, postzygotic reproductive isolation is perhaps the one that has attracted the most attention from geneticists. Despite remarkable advances in the identification of loci involved in Drosophila speciation, little is known about the genes, functions, and biochemical interactions of the molecules underlying hybrid sterility and inviability in mammals. Here, we discuss the main evolutionary and molecular features that make transcription factors (TFs), especially the family of zinc finger proteins with a Krüppel-associated box domain (KRAB-ZNF), strong candidates to play an important role in postzygotic reproductive isolation. Motivated by the recent identification of the gene encoding PR domain zinc finger protein 9 (Prdm9; a KRAB-ZNF gene) as the first hybrid sterility gene identified in mammals, we further propose integrative approaches to study KRAB-ZNF genes with the main goal of characterizing the molecular pathways and interactions involved in hybrid incompatibilities.

Transcriptome-wide patterns of divergence during allopatric evolution.

Recent studies have revealed repeated patterns of genomic divergence associated with species formation. Such patterns suggest that natural selection tends to target a set of available genes, but is also indicative that closely related taxa share evolutionary constraints that limit genetic variability. Studying patterns of genomic divergence among populations within the same species may shed light on the underlying evolutionary processes. Here, we examine transcriptome-wide divergence and polymorphism in the marine copepod Tigriopus californicus, a species where allopatric evolution has led to replicate sets of populations with varying degrees of divergence and hybrid incompatibility. Our analyses suggest that relatively small effective population sizes have resulted in an exponential decline of shared polymorphisms during population divergence and also facilitated the fixation of slightly deleterious mutations within allopatric populations. Five interpopulation comparisons at three different stages of divergence show that nonsynonymous mutations tend to accumulate in a specific set of proteins. These include proteins with central roles in cellular metabolism, such as those encoded in mtDNA, but also include an additional set of proteins that repeatedly show signatures of positive selection during allopatric divergence. Although our results are consistent with a contribution of nonadaptive processes, such as genetic drift and gene expression levels, generating repeatable patterns of genomic divergence in closely related taxa, they also indicate that adaptive evolution targeting a specific set of genes contributes to this pattern. Our results yield insights into the predictability of evolution at the gene level.