Ecological adaptation is of major relevance to speciation and sustainable population management, but the underlying genetic factors are typically hard to study in natural populations due to genetic differentiation caused by natural selection being confounded with genetic drift in subdivided populations. Here, we use whole genome population sequencing of Atlantic and Baltic herring to reveal the underlying genetic architecture at an unprecedented detailed resolution for both adaptation to a new niche environment and timing of reproduction. We identify almost 500 independent loci associated with a recent niche expansion from marine (Atlantic Ocean) to brackish waters (Baltic Sea), and more than 100 independent loci showing genetic differentiation between spring- and autumn-spawning populations irrespective of geographic origin. Our results show that both coding and non-coding changes contribute to adaptation. Haplotype blocks, often spanning multiple genes and maintained by selection, are associated with genetic differentiation.

Chicken body weight is an economically important trait and great genetic progress has been accomplished in genetic selective for body weight. To identify genes and chromosome regions associated with body weight, we performed a genome-wide association study using the chicken 60 k SNP panel in a chicken F2 resource population derived from the cross between Silky Fowl and White Plymouth Rock. A total of 26 SNP effects involving 9 different SNP markers reached 5% Bonferroni genome-wide significance. A chicken chromosome 4 (GGA4) region approximately 8.6 Mb in length (71.6–80.2 Mb) had a large number of significant SNP effects for late growth during weeks 7–12. The LIM domain-binding factor 2 (LDB2) gene in this region had the strongest association with body weight for weeks 7–12 and with average daily gain for weeks 6–12. This GGA4 region was previously reported to contain body weight QTL. GGA1 and GGA18 had three SNP effects on body weight with genome-wide significance. Some of the SNP effects with the significance of “suggestive linkage” overlapped with previously reported results.

Rose-comb, a classical monogenic trait of chickens, is characterized by a drastically altered comb morphology compared to the single-combed wild-type. Here we show that Rose-comb is caused by a 7.4 Mb inversion on chromosome 7 and that a second Rose-comb allele arose by unequal crossing over between a Rose-comb and wild-type chromosome. The comb phenotype is caused by the relocalization of the MNR2 homeodomain protein gene leading to transient ectopic expression of MNR2 during comb development. We also provide a molecular explanation for the first example of epistatic interaction reported by Bateson and Punnett 104 years ago, namely that walnut-comb is caused by the combined effects of the Rose-comb and Pea-comb alleles. Transient ectopic expression of MNR2 and SOX5 (causing the Pea-comb phenotype) occurs in the same population of mesenchymal cells and with at least partially overlapping expression in individual cells in the comb primordium. Rose-comb has pleiotropic effects, as homozygosity in males has been associated with poor sperm motility. We postulate that this is caused by the disruption of the CCDC108 gene located at one of the inversion breakpoints. CCDC108 is a poorly characterized protein, but it contains a MSP (major sperm protein) domain and is expressed in testis. The study illustrates several characteristic features of the genetic diversity present in domestic animals, including the evolution of alleles by two or more consecutive mutations and the fact that structural changes have contributed to fast phenotypic evolution. Comb morphology is a trait that shows considerable variability among domestic chickens. The Rose-comb mutation causes a drastically altered shape of the comb, whereas the Pea-comb mutation leads to a considerable reduction in the size of the comb. The combined effect of Rose-comb and Pea-comb is a comb shaped like a walnut, and the phenotype is consequently named walnut-comb. Both Pea-comb and Rose-comb are caused by structural changes in the genome leading to altered expression of important transcription factors. In a previous study we showed that Pea-comb is caused by misexpression of SOX5 during the development of the comb. In this study we report that Rose-comb is caused by a large inversion on chicken chromosome 7. The inversion moves the MNR2 gene to a new genomic location. This leads to misexpression of MNR2 during comb development, similar to the defect causing Pea-comb. Roosters that are homozygous for the Rose-comb inversion show poor sperm motility, and our results suggest that this is caused by the disruption of the CCDC108 gene that is located at one of the inversion breakpoints. CCDC108 is well conserved between chickens and humans, and this study establishes CCDC108 as a candidate gene for sperm motility disorders in humans.

The Atlantic herring is one of the most abundant vertebrates on earth but its nucleotide diversity is moderate (π = 0.3%), only three-fold higher than in human. Here, we present a pedigree-based estimation of the mutation rate in this species. Based on whole-genome sequencing of four parents and 12 offspring, the estimated mutation rate is 2.0 × 10^-9 per base per generation. We observed a high degree of parental mosaicism indicating that a large fraction of these de novo mutations occurred during early germ cell development. The estimated mutation rate – the lowest among vertebrates analyzed to date – partially explains the discrepancy between the rather low nucleotide diversity in herring and its huge census population size. But a species like the herring will never reach its expected nucleotide diversity because of fluctuations in population size over the millions of years it takes to build up high nucleotide diversity.

DNA cytosine methylation is an important epigenetic modification that has significant effects on a variety of biological processes in animals. Avian species hold a crucial position in evolutionary history. In this study, we used whole-genome bisulfite sequencing (MethylC-seq) to generate single base methylation profiles of lungs in two genetically distinct and highly inbred chicken lines (Fayoumi and Leghorn) that differ in genetic resistance to multiple pathogens, and we explored the potential regulatory role of DNA methylation associated with immune response differences between the two chicken lines. The MethylC-seq was used to generate single base DNA methylation profiles of Fayoumi and Leghorn birds. In addition, transcriptome profiling using RNA–seq from the same chickens and tissues were obtained to interrogate how DNA methylation regulates gene transcription on a genome-wide scale. The general DNA methylation pattern across different regions of genes was conserved compared to other species except for hyper-methylation of repeat elements, which was not observed in chicken. The methylation level of miRNA and pseudogene promoters was high, which indicates that silencing of these genes may be partially due to promoter hyper-methylation. Interestingly, the promoter regions of more recently evolved genes tended to be more highly methylated, whereas the gene body regions of evolutionarily conserved genes were more highly methylated than those of more recently evolved genes. Immune-related GO (Gene Ontology) terms were significantly enriched from genes within the differentially methylated regions (DMR) between Fayoumi and Leghorn, which implicates DNA methylation as one of the regulatory mechanisms modulating immune response differences between these lines. This study establishes a single-base resolution DNA methylation profile of chicken lung and suggests a regulatory role of DNA methylation in controlling gene expression and maintaining genome transcription stability. Furthermore, profiling the DNA methylomes of two genetic lines that differ in disease resistance provides a unique opportunity to investigate the potential role of DNA methylation in host disease resistance. Our study provides a foundation for future studies on epigenetic modulation of host immune response to pathogens in chickens.

The dwarf phenotype characterizes the smallest of rabbit breeds and is governed largely by the effects of a single dwarfing allele with an incompletely dominant effect on growth. Dwarf rabbits typically weigh under 1 kg and have altered craniofacial morphology. The dwarf allele is recessive lethal and dwarf homozygotes die within a few days of birth. The dwarf phenotype is expressed in heterozygous individuals and rabbits from dwarf breeds homozygous for the wild-type allele are normal, although smaller when compared to other breeds. Here, we show that the dwarf allele constitutes a ∼12.1 kb deletion overlapping the promoter region and first three exons of the HMGA2 gene leading to inactivation of this gene. HMGA2 has been frequently associated with variation in body size across species. Homozygotes for null alleles are viable in mice but not in rabbits and probably not in humans. RNA-sequencing analysis of rabbit embryos showed that very few genes (4–29 genes) were differentially expressed among the three HMGA2/dwarf genotypes, suggesting that dwarfism and inviability in rabbits are caused by modest changes in gene expression. Our results show that HMGA2 is critical for normal expression of IGF2BP2, which encodes an RNA-binding protein. Finally, we report a catalog of regions of elevated genetic differentiation between dwarf and normal-size rabbits, including LCORL-NCAPG, STC2, HOXD cluster, and IGF2BP2. Levels and patterns of genetic diversity at the LCORL-NCAPG locus further suggest that small size in dwarf breeds was enhanced by crosses with wild rabbits. Overall, our results imply that small size in dwarf rabbits results from a large effect, loss-of-function (LOF) mutation in HMGA2 combined with polygenic selection.