White horses - non-coding sequences drive premature hair greying and predisposition to melanoma

The Grey allele in horses is causing premature hair greying and susceptibility to melanoma. The causal mutation is a 4.6 kb tandem duplication in intron 6 of the Syntaxin 17 gene. A recent study demonstrated that the most common allele at the Grey locus (G3) involves three tandem copies of this sequence, whilst a more rare allele (G2) has two tandem copies and the wild-type allele (G1) only one copy. The G3 allele is causing fast greying and high incidence of skin melanoma, whereas the G2 allele is causing slow greying and no obvious increase in melanoma incidence. Further somatic copy number expansion has been documented in melanoma tissue from Grey horses. Functional studies showed that this intronic sequence acts as a weak melanocyte-specific enhancer that becomes substantially stronger by the copy number expansion. The Grey mutation is associated with upregulated expression of both Syntaxin 17 and the neighbouring NR4A3 gene in Grey horse melanomas. It is still an open question which of these genes is most important for the phenotypic effects or if causality is due to the combined effect of upregulation of both genes. Interestingly, RNAseq data in the Human Protein Atlas give support for a possible role of NR4A3 because it is particularly upregulated in human skin cancer, and it belongs to a cluster of genes associated with skin cancer and melanin biosynthesis. The Grey mutation and its association with melanoma provide a possibility to study the path to tumour development in numerous Grey horses carrying exactly the same predisposing mutation.

Research Note: Phenotypic trends for the multigenerational advanced intercross of the Virginia body weight lines of chickens

Random samples from generation S41 of the Virginia high and low 8-week body weight lines formed the base population for producing a multigenerational reciprocal intercross population. Although genetic mapping from this intercross has been reported, lacking are phenotypic trends across multiple generations. Here, we provide phenotypic information for the parental base population, the F1 reciprocal cross, and subsequent segregating recombinant generations F2 to F17. Heterosis for the selected trait in the F1 was negative for both reciprocal crosses. Phenotypic correlations for the selected trait in the recombinant generations were essentially nil for both males and females as was percent sexual dimorphism and coefficients of variation.

Adaptation to seasonal reproduction and environment-associated factors drive temporal and spatial differentiation in northwest Atlantic herring despite gene flow

Understanding how marine organisms adapt to local environments is crucial for predicting how populations will respond to global climate change. The genomic basis, environmental factors and evolutionary processes involved in local adaptation are however not well understood. Here we use Atlantic herring, an abundant, migratory and widely distributed marine fish with substantial genomic resources, as a model organism to evaluate local adaptation. We examined genomic variation and its correlation with environmental variables across a broad environmental gradient, for 15 spawning aggregations in Atlantic Canada and the United States. We then compared our results with available genomic data of northeast Atlantic populations. We confirmed that population structure lies in a fraction of the genome including likely adaptive genetic variants of functional importance. We discovered 10 highly differentiated genomic regions distributed across four chromosomes. Nine regions show strong association with seasonal reproduction. One region, corresponding to a known inversion on chromosome 12, underlies a latitudinal pattern discriminating populations north and south of a biogeographic transition zone on the Scotian Shelf. Genome-environment associations indicate that winter seawater temperature best correlates with the latitudinal pattern of this inversion. The variation at two so-called 'islands of divergence' related to seasonal reproduction appear to be private to the northwest Atlantic. Populations in the northwest and northeast Atlantic share variation at four of these divergent regions, simultaneously displaying significant diversity in haplotype composition at another four regions, which includes an undescribed structural variant approximately 7.7 Mb long on chromosome 8. Our results suggest that the timing and geographic location of spawning and early development may be under diverse selective pressures related to allelic fitness across environments. Our study highlights the role of genomic architecture, ancestral haplotypes and selection in maintaining adaptive divergence in species with large population sizes and presumably high gene flow.

How Fish Population Genomics Can Promote Sustainable Fisheries: A Road Map

Maintenance of genetic diversity in marine fishes targeted by commercial fishing is a grand challenge for the future. Most of these species are abundant and therefore important for marine ecosystems and food security. Here, we present a road map of how population genomics can promote sustainable fisheries. In these species, the development of reference genomes and whole genome sequencing is key, because genetic differentiation at neutral loci is usually low due to large population sizes and gene flow. First, baseline allele frequencies representing genetically differentiated populations within species must be established. These can then be used to accurately determine the composition of mixed samples, forming the basis for population demographic analysis to inform sustainably set fish quotas. SNP-chip analysis is a cost-effective method for determining baseline allele frequencies and for population identification in mixed samples. Finally, we describe how genetic marker analysis can transform stock identification and management.

Conservation of Affinity Rather Than Sequence Underlies a Dynamic Evolution of the Motif-Mediated p53/MDM2 Interaction in Ray-Finned Fishes

The transcription factor and cell cycle regulator p53 is marked for degradation by the ubiquitin ligase MDM2. The interaction between these 2 proteins is mediated by a conserved binding motif in the disordered p53 transactivation domain (p53TAD) and the folded SWIB domain in MDM2. The conserved motif in p53TAD from zebrafish displays a 20-fold weaker interaction with MDM2, compared to the interaction in human and chicken. To investigate this apparent difference, we tracked the molecular evolution of the p53TAD/MDM2 interaction among ray-finned fishes (Actinopterygii), the largest vertebrate clade. Intriguingly, phylogenetic analyses, ancestral sequence reconstructions, and binding experiments showed that different loss-of-affinity changes in the canonical binding motif within p53TAD have occurred repeatedly and convergently in different fish lineages, resulting in relatively low extant affinities (KD = 0.5 to 5 μM). However, for 11 different fish p53TAD/MDM2 interactions, nonconserved regions flanking the canonical motif increased the affinity 4- to 73-fold to be on par with the human interaction. Our findings suggest that compensating changes at conserved and nonconserved positions within the motif, as well as in flanking regions of low conservation, underlie a stabilizing selection of "functional affinity" in the p53TAD/MDM2 interaction. Such interplay complicates bioinformatic prediction of binding and calls for experimental validation. Motif-mediated protein-protein interactions involving short binding motifs and folded interaction domains are very common across multicellular life. It is likely that the evolution of affinity in motif-mediated interactions often involves an interplay between specific interactions made by conserved motif residues and nonspecific interactions by nonconserved disordered regions.

Low Mutation Load in a Supergene Underpinning Alternative Male Mating Strategies in Ruff (Calidris pugnax)

A paradox in evolutionary biology is how supergenes can maintain high fitness despite reduced effective population size, the suppression of recombination, and the expected accumulation of mutational load. The ruff supergene involves 2 rare inversion haplotypes (satellite and faeder). These are recessive lethals but with dominant effects on male mating strategies, plumage, and body size. Sequence divergence to the wild-type (independent) haplotype indicates that the inversion could be as old as 4 million years. Here, we have constructed a highly contiguous genome assembly of the inversion region for both the independent and satellite haplotypes. Based on the new data, we estimate that the recombination event(s) creating the satellite haplotype occurred only about 70,000 yr ago. Contrary to expectations for supergenes, we find no substantial expansion of repeats and only a modest mutation load on the satellite and faeder haplotypes despite high sequence divergence to the non-inverted haplotype (1.46%). The essential centromere protein N (CENPN) gene is disrupted by the inversion and is as well conserved on the inversion haplotypes as on the noninversion haplotype. These results suggest that the inversion may be much younger than previously thought. The low mutation load, despite recessive lethality, may be explained by the introgression of the inversion from a now extinct lineage.

Low mutation rate in epaulette sharks is consistent with a slow rate of evolution in sharks

Sharks occupy diverse ecological niches and play critical roles in marine ecosystems, often acting as apex predators. They are considered a slow-evolving lineage and have been suggested to exhibit exceptionally low cancer rates. These two features could be explained by a low nuclear mutation rate. Here, we provide a direct estimate of the nuclear mutation rate in the epaulette shark (Hemiscyllium ocellatum). We generate a high-quality reference genome, and resequence the whole genomes of parents and nine offspring to detect de novo mutations. Using stringent criteria, we estimate a mutation rate of 7×10-10 per base pair, per generation. This represents one of the lowest directly estimated mutation rates for any vertebrate clade, indicating that this basal vertebrate group is indeed a slowly evolving lineage whose ability to restore genetic diversity following a sustained population bottleneck may be hampered by a low mutation rate.

Introgression underlies phylogenetic uncertainty but not parallel plumage evolution in a recent songbird radiation

Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails", a group of five species in the songbird genus Motacilla: one species, Motacilla alba, shows wide intra-specific plumage variation, while the four others form two pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the two species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome, to quantify the extent of discordance across genomic regions, and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investigation of genome-wide patterns of gene tree heterogeneity to help understanding the mechanisms underlying phenotypic evolution.

Community-wide genome sequencing reveals 30 years of Darwin's finch evolution

A fundamental goal in evolutionary biology is to understand the genetic architecture of adaptive traits. Using whole-genome data of 3955 of Darwin's finches on the Galápagos Island of Daphne Major, we identified six loci of large effect that explain 45% of the variation in the highly heritable beak size of Geospiza fortis, a key ecological trait. The major locus is a supergene comprising four genes. Abrupt changes in allele frequencies at the loci accompanied a strong change in beak size caused by natural selection during a drought. A gradual change in Geospiza scandens occurred across 30 years as a result of introgressive hybridization with G. fortis. This study shows how a few loci with large effect on a fitness-related trait contribute to the genetic potential for rapid adaptive radiation.

The RNA-binding protein ZC3H11A interacts with the nuclear poly(A)-binding protein PABPN1 and alters polyadenylation of viral transcripts

Nuclear mRNA metabolism is regulated by multiple proteins, which either directly bind to RNA or form multiprotein complexes. The RNA-binding protein ZC3H11A is involved in nuclear mRNA export, NF-κB signaling, and is essential during mouse embryo development. Furthermore, previous studies have shown that ZC3H11A is important for nuclear-replicating viruses. However, detailed biochemical characterization of the ZC3H11A protein has been lacking. In this study, we established the ZC3H11A protein interactome in human and mouse cells. We demonstrate that the nuclear poly(A)-binding protein PABPN1 interacts specifically with the ZC3H11A protein and controls ZC3H11A localization into nuclear speckles. We report that ZC3H11A specifically interacts with the human adenovirus type 5 (HAdV-5) capsid mRNA in a PABPN1-dependent manner. Notably, ZC3H11A uses the same zinc finger motifs to interact with PABPN1 and viral mRNA. Further, we demonstrate that the lack of ZC3H11A alters the polyadenylation of HAdV-5 capsid mRNA. Taken together, our results suggest that the ZC3H11A protein may act as a novel regulator of polyadenylation of nuclear mRNA.

Ablation of ZC3H11A causes early embryonic lethality and dysregulation of metabolic processes

ZC3H11A (zinc finger CCCH domain-containing protein 11A) is a stress-induced mRNA-binding protein required for efficient growth of nuclear-replicating viruses. The cellular functions of ZC3H11A during embryonic development are unknown. Here, we report the generation and phenotypic characterization of Zc3h11a knockout (KO) mice. Heterozygous null Zc3h11a mice were born at the expected frequency without distinguishable phenotypic differences compared with wild-type mice. In contrast, homozygous null Zc3h11a mice were missing, indicating that Zc3h11a is crucial for embryonic viability and survival. Zc3h11a ^-/- embryos were detected at the expected Mendelian ratios up to late preimplantation stage (E4.5). However, phenotypic characterization at E6.5 revealed degeneration of Zc3h11a ^-/- embryos, indicating developmental defects around the time of implantation. Transcriptomic analyses documented a dysregulation of glycolysis and fatty acid metabolic pathways in Zc3h11a^-/- embryos at E4.5. Proteomic analysis indicated a tight interaction between ZC3H11A and mRNA-export proteins in embryonic stem cells. CLIP-seq analysis demonstrated that ZC3H11A binds a subset of mRNA transcripts that are critical for metabolic regulation of embryonic cells. Furthermore, embryonic stem cells with an induced deletion of Zc3h11a display an impaired differentiation toward epiblast-like cells and impaired mitochondrial membrane potential. Altogether, the results show that ZC3H11A is participating in export and posttranscriptional regulation of selected mRNA transcripts required to maintain metabolic processes in embryonic cells. While ZC3H11A is essential for the viability of the early mouse embryo, inactivation of Zc3h11a expression in adult tissues using a conditional KO did not lead to obvious phenotypic defects.

The genomic basis and environmental correlates of local adaptation in the Atlantic horse mackerel (Trachurus trachurus)

Understanding how populations adapt to their environment is increasingly important to prevent biodiversity loss due to overexploitation and climate change. Here we studied the population structure and genetic basis of local adaptation of Atlantic horse mackerel, a commercially and ecologically important marine fish that has one of the widest distributions in the eastern Atlantic. We analyzed whole-genome sequencing and environmental data of samples collected from the North Sea to North Africa and the western Mediterranean Sea. Our genomic approach indicated low population structure with a major split between the Mediterranean Sea and the Atlantic Ocean and between locations north and south of mid-Portugal. Populations from the North Sea are the most genetically distinct in the Atlantic. We discovered that most population structure patterns are driven by a few highly differentiated putatively adaptive loci. Seven loci discriminate the North Sea, two the Mediterranean Sea, and a large putative inversion (9.9 Mb) on chromosome 21 underlines the north-south divide and distinguishes North Africa. A genome-environment association analysis indicates that mean seawater temperature and temperature range, or factors correlated to them, are likely the main environmental drivers of local adaptation. Our genomic data broadly support the current stock divisions, but highlight areas of potential mixing, which require further investigation. Moreover, we demonstrate that as few as 17 highly informative SNPs can genetically discriminate the North Sea and North African samples from neighboring populations. Our study highlights the importance of both, life history and climate-related selective pressures in shaping population structure patterns in marine fish. It also supports that chromosomal rearrangements play a key role in local adaptation with gene flow. This study provides the basis for more accurate delineation of the horse mackerel stocks and paves the way for improving stock assessments.

A long-standing hybrid population between Pacific and Atlantic herring in a subarctic fjord of Norway

Atlantic herring (Clupea harengus) and Pacific herring (C. pallasii) are sister species that split from a common ancestor about 2 million years ago. Balsfjord, a subarctic fjord in Northern Norway, harbors an outpost population of Pacific herring within the range of the Atlantic herring. We used whole genome sequencing to show that gene flow from Atlantic herring into the Balsfjord population has generated a stable hybrid population that has persisted for thousands of generations. The Atlantic herring ancestry in Balsfjord was estimated in the range 25-26%. The old age and large proportion of introgressed regions suggest there are no obvious genetic incompatibilities between species. Introgressed regions were widespread in the genome and large, with some in excess of 1 Mb, and they were overrepresented in low-recombination regions. We show that the distribution of introgressed material is non-random; introgressed sequence blocks in different individuals are shared more often than expected by chance. Furthermore, introgressed regions tend to show elevated divergence (FST) between Atlantic and Pacific herring. Together, our results suggest that introgression of genetic material has facilitated adaptation in the Balsfjord population. The Balsfjord population provides a rare example of a stable interspecies hybrid population that has persisted over thousands of years.

A frame-shift mutation in COMTD1 is associated with impaired pheomelanin pigmentation in chicken

The biochemical pathway regulating the synthesis of yellow/red pheomelanin is less well characterized than the synthesis of black/brown eumelanin. Inhibitor of gold (IG phenotype) is a plumage colour variant in chicken that provides an opportunity to further explore this pathway since the recessive allele (IG) at this locus is associated with a defect in the production of pheomelanin. IG/IG homozygotes display a marked dilution of red pheomelanin pigmentation, whilst black pigmentation (eumelanin) is only slightly affected. Here we show that a 2-base pair insertion (frame-shift mutation) in the 5th exon of the Catechol-O-methyltransferase containing domain 1 gene (COMTD1), expected to cause a complete or partial loss-of-function of the COMTD1 enzyme, shows complete concordance with the IG phenotype within and across breeds. We show that the COMTD1 protein is localized to mitochondria in pigment cells. Knockout of Comtd1 in a mouse melanocytic cell line results in a reduction in pheomelanin metabolites and significant alterations in metabolites of glutamate/glutathione, riboflavin, and the tricarboxylic acid cycle. Furthermore, COMTD1 overexpression enhanced cellular proliferation following chemical-induced transfection, a potential inducer of oxidative stress. These observations suggest that COMTD1 plays a protective role for melanocytes against oxidative stress and that this supports their ability to produce pheomelanin.

ZC3H11A loss of function enhances NF-κB signaling through defective IκBα protein expression

ZC3H11A is a cellular protein associated with the transcription export (TREX) complex that is induced during heat-shock. Several nuclear-replicating viruses exploit the mRNA export mechanism of ZC3H11A protein for their efficient replication. Here we show that ZC3H11A protein plays a role in regulation of NF-κB signal transduction. Depletion of ZC3H11A resulted in enhanced NF-κB mediated signaling, with upregulation of numerous innate immune related mRNAs, including IL-6 and a large group of interferon-stimulated genes. IL-6 upregulation in the absence of the ZC3H11A protein correlated with an increased NF-κB transcription factor binding to the IL-6 promoter and decreased IL-6 mRNA decay. The enhanced NF-κB signaling pathway in ZC3H11A deficient cells correlated with a defect in IκBα inhibitory mRNA and protein accumulation. Upon ZC3H11A depletion The IκBα mRNA was retained in the cell nucleus resulting in failure to maintain normal levels of the cytoplasmic IκBα mRNA and protein that is essential for its inhibitory feedback loop on NF-κB activity. These findings indicate towards a previously unknown mechanism of ZC3H11A in regulating the NF-κB pathway at the level of IkBα mRNA export.

Spatiotemporal variations in retrovirus-host interactions among Darwin's finches

Endogenous retroviruses (ERVs) are inherited remnants of retroviruses that colonized host germline over millions of years, providing a sampling of retroviral diversity across time. Here, we utilize the strength of Darwin’s finches, a system synonymous with evolutionary studies, for investigating ERV history, revealing recent retrovirus-host interactions in natural populations. By mapping ERV variation across all species of Darwin’s finches and comparing with outgroup species, we highlight geographical and historical patterns of retrovirus-host occurrence, utilizing the system for evaluating the extent and timing of retroviral activity in hosts undergoing adaptive radiation and colonization of new environments. We find shared ERVs among all samples indicating retrovirus-host associations pre-dating host speciation, as well as considerable ERV variation across populations of the entire Darwin’s finches’ radiation. Unexpected ERV variation in finch species on different islands suggests historical changes in gene flow and selection. Non-random distribution of ERVs along and between chromosomes, and across finch species, suggests association between ERV accumulation and the rapid speciation of Darwin’s finches.

Endogenous retroviruses (ERV) are inherited remains of retroviruses that have colonized host genomes during evolution. Here the authors observe considerable species-specific ERV variation among Darwin’s finches, reflecting historic retrovirus-host interactions.

A baseline for the genetic stock identification of Atlantic herring, Clupea harengus, in ICES Divisions 6.a, 7.b-c

Atlantic herring in International Council for Exploration of the Sea (ICES) Divisions 6.a, 7.b-c comprises at least three populations, distinguished by temporal and spatial differences in spawning, which have until recently been managed as two stocks defined by geographical delineators. Outside of spawning the populations form mixed aggregations, which are the subject of acoustic surveys. The inability to distinguish the populations has prevented the development of separate survey indices and separate stock assessments. A panel of 45 single-nucleotide polymorphisms, derived from whole-genome sequencing, were used to genotype 3480 baseline spawning samples (2014-2021). A temporally stable baseline comprising 2316 herring from populations known to inhabit Division 6.a was used to develop a genetic assignment method, with a self-assignment accuracy greater than 90%. The long-term temporal stability of the assignment model was validated by assigning archive (2003-2004) baseline samples (270 individuals) with a high level of accuracy. Assignment of non-baseline samples (1514 individuals) from Divisions 6.a, 7.b-c indicated previously unrecognized levels of mixing of populations outside of the spawning season. The genetic markers and assignment models presented constitute a 'toolbox' that can be used for the assignment of herring caught in mixed survey and commercial catches in Division 6.a into their population of origin with a high level of accuracy.

Rapid adaptive radiation of Darwin's finches depends on ancestral genetic modules

Recent adaptive radiations are models for investigating mechanisms contributing to the evolution of biodiversity. An unresolved question is the relative importance of new mutations, ancestral variants, and introgressive hybridization for phenotypic evolution and speciation. Here, we address this issue using Darwin's finches and investigate the genomic architecture underlying their phenotypic diversity. Admixture mapping for beak and body size in the small, medium, and large ground finches revealed 28 loci showing strong genetic differentiation. These loci represent ancestral haplotype blocks with origins predating speciation events during the Darwin's finch radiation. Genes expressed in the developing beak are overrepresented in these genomic regions. Ancestral haplotypes constitute genetic modules for selection and act as key determinants of the unusual phenotypic diversity of Darwin's finches. Such ancestral haplotype blocks can be critical for how species adapt to environmental variability and change.

Hmga2 deficiency is associated with allometric growth retardation, infertility, and behavioral abnormalities in mice

The high mobility group AT-hook 2 (HMGA2) protein works as an architectural regulator by binding AT-rich DNA sequences to induce conformational changes affecting transcription. Genomic deletions disrupting HMGA2 coding sequences and flanking noncoding sequences cause dwarfism in mice and rabbits. Here, CRISPR/Cas9 was used in mice to generate an Hmga2 null allele that specifically disrupts only the coding sequence. The loss of one or both alleles of Hmga2 resulted in reduced body size of 20% and 60%, respectively, compared to wild-type littermates as well as an allometric reduction in skull length in Hmga2-/- mice. Both male and female Hmga2-/- mice are infertile, whereas Hmga2+/- mice are fertile. Examination of reproductive tissues of Hmga2-/- males revealed a significantly reduced size of testis, epididymis, and seminal vesicle compared to controls, and 70% of knock-out males showed externalized penis, but no cryptorchidism was observed. Sperm analyses revealed severe oligospermia in mutant males and slightly decreased sperm viability, increased DNA damage but normal sperm chromatin compaction. Testis histology surprisingly revealed a normal seminiferous epithelium, despite the significant reduction in testis size. In addition, Hmga2-/- mice showed a significantly reduced exploratory behavior. In summary, the phenotypic effects in mouse using targeted mutagenesis confirmed that Hmga2 is affecting prenatal and postnatal growth regulation, male reproductive tissue development, and presents the first indication that Hmga2 function is required for normal mouse behavior. No specific effect, despite an allometric reduction, on craniofacial development was noted in contrast to previous reports of an altered craniofacial development in mice and rabbits carrying deletions of both coding and noncoding sequences at the 5' part of Hmga2.

In Situ Measurements of Thermal Ion Temperature in the Martian Ionosphere

In situ measurements of ionospheric and thermospheric temperatures are experimentally challenging because orbiting spacecraft typically travel supersonically with respect to the cold gas and plasma. We present O 2 + temperatures in Mars' ionosphere derived from data measured by the SupraThermal And Thermal Ion Composition instrument onboard the Mars Atmosphere and Volatile EvolutioN spacecraft. We focus on data obtained during nine special orbit maneuvers known as Deep Dips, during which MAVEN lowered its periapsis altitude from the nominal 150 to 120 km for 1 week in order to sample the ionospheric main peak and approach the homopause. We use two independent techniques to calculate ion temperatures from the measured energy and angular widths of the supersonic ram ion beam. After correcting for background and instrument response, we are able to measure ion temperatures as low as 100 K with associated uncertainties as low as 10%. It is theoretically expected that ion temperatures will converge to the neutral temperature at altitudes below the exobase region (∼180-200 km) due to strong collisional coupling; however, no evidence of the expected thermalization is observed. We have eliminated several possible explanations for the observed temperature difference between ions and neutrals, including Coulomb collisions with electrons, Joule heating, and heating caused by interactions with the spacecraft. The source of the energy maintaining the high ion temperatures remains unclear, suggesting that a fundamental piece of physics is missing from existing models of the Martian ionosphere.

The feather pattern autosomal barring in chicken is strongly associated with segregation at the MC1R locus

Color patterns within individual feathers are common in birds but little is known about the genetic mechanisms causing such patterns. Here, we investigate the genetic basis for autosomal barring in chicken, a horizontal striping pattern on individual feathers. Using an informative backcross, we demonstrate that the MC1R locus is strongly associated with this phenotype. A deletion at SOX10, underlying the dark brown phenotype on its own, affects the manifestation of the barring pattern. The coding variant L133Q in MC1R is the most likely causal mutation for autosomal barring in this pedigree. Furthermore, a genetic screen across six different breeds showing different patterning phenotypes revealed that the most striking shared characteristics among these breeds were that they all carried the MC1R alleles Birchen or brown. Our data suggest that the presence of activating MC1R mutations enhancing pigment synthesis is an important mechanism underlying pigmentation patterns on individual feathers in chicken. We propose that MC1R and its antagonist ASIP play a critical role for determining within-feather pigmentation patterns in birds by acting as activator and inhibitor possibly in a Turing reaction-diffusion model.

Porcine ZBED6 regulates growth of skeletal muscle and internal organs via multiple targets

ZBED6 (zinc finger BED domain containing protein 6) is a transcription factor unique to placental mammals and its interaction with the IGF2 (insulin-like growth factor 2) locus plays a prominent role in the regulation of postnatal skeletal muscle growth. Here, we generated lean Bama miniature pigs by generating ZBED6-knockout (ZBED6^−/−) and investigated the mechanism underlying ZBED6 in growth of muscle and internal organs of placental mammals. ZBED6^−/− pigs show markedly higher lean mass, lean mass rate, larger muscle fiber area and heavier internal organs (heart and liver) than wild-type (WT) pigs. The striking phenotypic changes of ZBED6^-/- pigs coincided with remarkable upregulation of IGF2 mRNA and protein expression across three tissues (gastrocnemius muscle, longissimus dorsi, heart). Despite a significant increase in liver weight, ZBED6^-/- pigs show comparable levels of IGF2 expression to those of WT controls. A mechanistic study revealed that elevated methylation in the liver abrogates ZBED6 binding at the IGF2 locus, explaining the unaltered hepatic IGF2 expression in ZBED6^-/- pigs. These results indicate that a ZBED6-IGF2-independent regulatory pathway exists in the liver. Transcriptome analysis and ChIP-PCR revealed new ZBED6 target genes other than IGF2, including cyclin dependent kinase inhibitor 1A (CDKN1A) and tsukushi, small leucine rich proteoglycan (TSKU), that regulates growth of muscle and liver, respectively.

The lean meat rate is an important economic trait for the swine industry and it is determined by muscle growth and development. A single base change in intron 3 of the insulin-like growth factor 2 (IGF2) gene increases meat production in pigs by disrupting a binding site for zinc finger BED domain containing protein 6 (ZBED6). Chinese indigenous pig breeds carrying the homozygous IGF2 wild-type allele produce low lean meat. We thus generate a lean pig model in Chinese Bama pig by knocking out ZBED6. In this model, we demonstrate that ZBED6 KO increases muscle and internal organ growth through ZBED6-IGF2 axis and other target genes. These results not only open new strategies for lean meat breeding in Chinese indigenous pigs, but also provide new insights to the global function of ZBED6 in organ growth and development.

A multispecies BCO2 beak color polymorphism in the Darwin's finch radiation

Carotenoid-based polymorphisms are widespread in populations of birds, fish, and reptiles,¹ but generally little is known about the factors affecting their maintenance in populations.² We report a combined field and molecular-genetic investigation of a nestling beak color polymorphism in Darwin's finches. Beaks are pink or yellow, and yellow is recessive.³ Here we show that the polymorphism arose in the Galápagos half a million years ago through a mutation associated with regulatory change in the BCO2 gene and is shared by 14 descendant species. The polymorphism is probably a balanced polymorphism, maintained by ecological selection associated with survival and diet. In cactus finches, the frequency of the yellow genotype is correlated with cactus fruit abundance and greater hatching success and may be altered by introgressive hybridization. Polymorphisms that are hidden as adults, as here, may be far more common than is currently recognized, and contribute to diversification in ways that are yet to be discovered.

Cis-acting mutation affecting GJA5 transcription is underlying the Melanotic within-feather pigmentation pattern in chickens

Melanotic (Ml) is a mutation in chickens that extends black (eumelanin) pigmentation in normally brown or red (pheomelanin) areas, thus affecting multiple within-feather patterns [J. W. Moore, J. R. Smyth Jr, J. Hered. 62, 215-219 (1971)]. In the present study, linkage mapping using a back-cross between Dark Cornish (Ml/Ml) and Partridge Plymouth Rock (ml+/ml+ ) chickens assigned Ml to an 820-kb region on chromosome 1. Identity-by-descent mapping, via whole-genome sequencing and diagnostic tests using a diverse set of chickens, refined the localization to the genomic region harboring GJA5 encoding gap-junction protein 5 (alias connexin 40) previously associated with pigmentation patterns in zebrafish. An insertion/deletion polymorphism located in the vicinity of the GJA5 promoter region was identified as the candidate causal mutation. Four different GJA5 transcripts were found to be expressed in feather follicles and at least two showed differential expression between genotypes. The results showed that Melanotic constitutes a cis-acting regulatory mutation affecting GJA5 expression. A recent study established the melanocortin-1 receptor (MC1R) locus and the interaction between the MC1R receptor and its antagonist agouti-signaling protein as the primary mechanism underlying variation in within-feather pigmentation patterns in chickens. The present study advances understanding the mechanisms underlying variation in plumage color in birds because it demonstrates that the activity of connexin 40/GJA5 can modulate the periodic pigmentation patterns within individual feathers.

ZBED6 counteracts high-fat diet-induced glucose intolerance by maintaining beta cell area and reducing excess mitochondrial activation

AIMS/HYPOTHESIS

ZBED6 (zinc finger, BED-type containing 6) is known to regulate muscle mass by suppression of Igf2 gene transcription. In insulin-producing cell lines, ZBED6 maintains proliferative capacity at the expense of differentiation and beta cell function. The aim was to study the impact of Zbed6 knockout on beta cell function and glucose tolerance in C57BL/6 mice.

METHODS

Beta cell area and proliferation were determined in Zbed6 knockout mice using immunohistochemical analysis. Muscle and fat distribution were assessed using micro-computed tomography. Islet gene expression was assessed by RNA sequencing. Effects of a high-fat diet were analysed by glucose tolerance and insulin tolerance tests. ZBED6 was overexpressed in EndoC-βH1 cells and human islet cells using an adenoviral vector. Beta cell cell-cycle analysis, insulin release and mitochondrial function were studied in vitro using propidium iodide staining and flow cytometry, ELISA, the Seahorse technique, and the fluorescent probes JC-1 and MitoSox.

RESULTS

Islets from Zbed6 knockout mice showed lowered expression of the cell cycle gene Pttg1, decreased beta cell proliferation and decreased beta cell area, which occurred independently from ZBED6 effects on Igf2 gene expression. Zbed6 knockout mice, but not wild-type mice, developed glucose intolerance when given a high-fat diet. The high-fat diet Zbed6 knockout islets displayed upregulated expression of oxidative phosphorylation genes and genes associated with beta cell differentiation. In vitro, ZBED6 overexpression resulted in increased EndoC-βH1 cell proliferation and a reduced glucose-stimulated insulin release in human islets. ZBED6 also reduced mitochondrial JC-1 J-aggregate formation, mitochondrial oxygen consumption rates (OCR) and mitochondrial reactive oxygen species (ROS) production, both at basal and palmitate + high glucose-stimulated conditions. ZBED6-induced inhibition of OCR was not rescued by IGF2 addition. ZBED6 reduced levels of the mitochondrial regulator PPAR-γ related coactivator 1 protein (PRC) and bound its promoter/enhancer region. Knockdown of PRC resulted in a lowered OCR.

CONCLUSIONS/INTERPRETATION

It is concluded that ZBED6 is required for normal beta cell replication and also limits excessive beta cell mitochondrial activation in response to an increased functional demand. ZBED6 may act, at least in part, by restricting PRC-mediated mitochondrial activation/ROS production, which may lead to protection against beta cell dysfunction and glucose intolerance in vivo.

ZBED6 regulates Igf2 expression partially through its regulation of miR483 expression

The expression of Igf2 in mammals shows a complex regulation involving multiple promoters and epigenetic mechanisms. We previously identified a novel regulatory mechanism based on the interaction between the transcriptional factor ZBED6 and Igf2 intron. Disruption of the ZBED6-Igf2 interaction leads to a dramatic up-regulation of IGF2 expression postnatally. In the current study we characterize an additional layer of regulation involving miR483 encoded by another Igf2 intron. We found a highly significant up-regulation of miR483 expression when the ZBED6-Igf2 axis is disrupted in transgenic mice. Furthermore, CRISPR/Cas9 mediated knock-out of miR483 in C2C12 myoblast cells, both wild-type and cells with disrupted ZBED6-Igf2 axis (Igf2dGGCT), resulted in down-regulation of Igf2 expression and a reduced proliferation rate. This was further validated using miR483 mimics and inhibitors. RNA-seq analysis revealed a significant enrichment of genes involved in the PI3K-Akt signaling pathway among genes down-regulated in miR483-/- cells, including Igf2 down-regulation. The opposite pattern was observed in Igf2dGGCT cells, where Igf2 is up-regulated. Our data suggest a positive feedback between miR483 and Igf2 promoter activity, strongly affecting how ZBED6 controls Igf2 expression in various cell types.

Brain Transcriptomics of Wild and Domestic Rabbits Suggests That Changes in Dopamine Signaling and Ciliary Function Contributed to Evolution of Tameness

Domestication has resulted in immense phenotypic changes in animals despite their relatively short evolutionary history. The European rabbit is one of the most recently domesticated animals, but exhibits distinct morphological, physiological, and behavioral differences from their wild conspecifics. A previous study revealed that sequence variants with striking allele frequency differences between wild and domestic rabbits were enriched in conserved noncoding regions, in the vicinity of genes involved in nervous system development. This suggests that a large proportion of the genetic changes targeted by selection during domestication might affect gene regulation. Here, we generated RNA-sequencing data for four brain regions (amygdala, hypothalamus, hippocampus, and parietal/temporal cortex) sampled at birth and revealed hundreds of differentially expressed genes (DEGs) between wild and domestic rabbits. DEGs in amygdala were significantly enriched for genes associated with dopaminergic function and all 12 DEGs in this category showed higher expression in domestic rabbits. DEGs in hippocampus were enriched for genes associated with ciliary function, all 21 genes in this category showed lower expression in domestic rabbits. These results indicate an important role of dopamine signaling and ciliary function in the evolution of tameness during rabbit domestication. Our study shows that gene expression in specific pathways has been profoundly altered during domestication, but that the majority of genes showing differential expression in this study have not been the direct targets of selection.

The crest phenotype in domestic chicken is caused by a 197 bp duplication in the intron of HOXC10

The Crest mutation in chicken shows incomplete dominance and causes a spectacular phenotype in which the small feathers normally present on the head are replaced by much larger feathers normally present only in dorsal skin. Using whole-genome sequencing, we show that the crest phenotype is caused by a 197 bp duplication of an evolutionarily conserved sequence located in the intron of HOXC10 on chromosome 33. A diagnostic test showed that the duplication was present in all 54 crested chickens representing eight breeds and absent from all 433 non-crested chickens representing 214 populations. The mutation causes ectopic expression of at least five closely linked HOXC genes, including HOXC10, in cranial skin of crested chickens. The result is consistent with the interpretation that the crest feathers are caused by an altered body region identity. The upregulated HOXC gene expression is expanded to skull tissue of Polish chickens showing a large crest often associated with cerebral hernia, but not in Silkie chickens characterized by a small crest, both homozygous for the duplication. Thus, the 197 bp duplication is required for the development of a large crest and susceptibility to cerebral hernia because only crested chicken show this malformation. However, this mutation is not sufficient to cause herniation because this malformation is not present in breeds with a small crest, like Silkie chickens.

Introduction of an interactive tool (the Dental Trauma Guide) in the undergraduate dental teaching to manage traumatic dental injuries

BACKGROUND/AIMS

Traumatic dental injuries (TDI) are complex problems where lack of proper care may result in serious complications. The need to improve the management of TDI is a frequently addressed concern. Methods of improvement in their diagnosis and management are continuously evolving. The interactive Internet tool, the Dental Trauma Guide (DTG), helps to simplify diagnostic and management dilemmas. However, it is not a freely available tool. The aim of the current study was to assess the knowledge and diagnostic skills of undergraduate dental students with access to the DTG compared with students without such access, in order to validate and promote this tool in dental education.

MATERIALS AND METHODS

Two groups of students were randomly selected where one group of final year dental undergraduate students were exposed to lectures, demonstrations, discussions and tutorials on the management of TDI according to the standard undergraduate curriculum in Sri Lanka. Another test group of 21 students were provided with access to DTG during their training in paediatric dentistry. At the end of the study period, students were assessed on their knowledge of TDI using MCQs (Multiple Choice Questions) and OSCEs (Objective Structured Clinical Examination), based on the DTG.

RESULTS

The students with access to the DTG were more knowledgeable in providing the correct answers to three out of the seven OSCE questions. Evaluation based on the MCQs did not reveal a significant difference (p = .913). However, users of the DTG showed a statistically significant difference with better overall knowledge based on their answers (p = .028). Following this period of evaluation, all of the students were provided with access to the DTG to supplement their learning experience.

CONCLUSION

The Dental Trauma Guide is a useful supplementary tool for undergraduate students to arrive at a correct diagnosis and treatment plan for TDI.

A chromosome-level assembly of blunt snout bream (Megalobrama amblycephala) reveals an expansion of olfactory receptor genes in freshwater fish

The number of olfactory receptor genes (ORs), which are responsible for detecting diverse odor molecules varies extensively among mammals as a result of frequent gene gains and losses that contribute to olfactory specialization. However, how OR expansions/contractions in fish are influenced by habitat and feeding habit and which OR subfamilies are important in each ecological niche is unknown. Here, we report a major OR expansion in a freshwater herbivorous fish, Megalobrama amblycephala, using a highly contiguous, chromosome-level assembly. We evaluate the possible contribution of OR expansion to habitat and feeding specialization by comparing the OR repertoire in 28 phylogenetically and ecologically diverse teleosts. In total, we analyzed > 4,000 ORs including 3,253 intact, 122 truncated, and 913 pseudogenes. The number of intact ORs is highly variable ranging from 20 to 279. We estimate that the most recent common ancestor of Osteichthyes had 62 intact ORs, which declined in most lineages except the freshwater Otophysa clade that has a substantial expansion in subfamily β and ε ORs. Across teleosts, we found a strong association between duplications of β and ε ORs and freshwater habitat. Nearly, all ORs were expressed in the olfactory epithelium (OE) in three tested fish species. Specifically, all the expanded β and ε ORs were highly expressed in OE of M. amblycephala. Together, we provide molecular and functional evidence for how OR repertoires in fish have undergone gain and loss with respect to ecological factors and highlight the role of β and ε OR in freshwater adaptation.

Mutations Upstream of the TBX5 and PITX1 Transcription Factor Genes Are Associated with Feathered Legs in the Domestic Chicken

Feathered leg is a trait in domestic chickens that has undergone intense selection by fancy breeders. Previous studies have shown that two major loci controlling feathered leg are located on chromosomes 13 and 15. Here, we present genetic evidence for the identification of candidate causal mutations at these loci. This was accomplished by combining classical linkage mapping using an experimental cross segregating for feathered leg and high-resolution identical-by-descent mapping using whole-genome sequence data from 167 samples of chicken with or without feathered legs. The first predicted causal mutation is a single-base change located 25 kb upstream of the gene for the forelimb-specific transcription factor TBX5 on chromosome 15. The second is a 17.7-kb deletion located ∼200 kb upstream of the gene for the hindlimb-specific transcription factor PITX1 on chromosome 13. These mutations are predicted to activate TBX5 and repress PITX1 expression, respectively. The study reveals a remarkable convergence in the evolution of the feathered-leg phenotype in domestic chickens and domestic pigeons, as this phenotype is caused by noncoding mutations upstream of the same two genes. Furthermore, the PITX1 causal variants are large overlapping deletions, 17.7 kb in chicken and 44 kb in pigeons. The results of the present study are consistent with the previously proposed model for pigeon that feathered leg is caused by reduced PITX1 expression and ectopic expression of TBX5 in hindlimb buds resulting in a shift of limb identity from hindlimb to more forelimb-like identity.

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.

Asymmetric introgression reveals the genetic architecture of a plumage trait

Genome-wide variation in introgression rates across hybrid zones offers a powerful opportunity for studying population differentiation. One poorly understood pattern of introgression is the geographic displacement of a trait implicated in lineage divergence from genome-wide population boundaries. While difficult to interpret, this pattern can facilitate the dissection of trait genetic architecture because traits become uncoupled from their ancestral genomic background. We studied an example of trait displacement generated by the introgression of head plumage coloration from personata to alba subspecies of the white wagtail. A previous study of their hybrid zone in Siberia revealed that the geographic transition in this sexual signal that mediates assortative mating was offset from other traits and genetic markers. Here we show that head plumage is associated with two small genetic regions. Despite having a simple genetic architecture, head plumage inheritance is consistent with partial dominance and epistasis, which could contribute to its asymmetric introgression.

Hybrid zones are windows into the evolutionary process. Semenov et al. find that the head plumage differences between white wagtail subspecies have a simple genetic basis involving two small genetic regions, in which partially dominant and epistatic interactions help to explain how this sexual signal has become decoupled from other plumage traits.

Ecological adaptation in European eels is based on phenotypic plasticity

The relative role of genetic adaptation and phenotypic plasticity is of fundamental importance in evolutionary ecology [M. J. West-Eberhard, Proc. Natl. Acad. Sci. U.S.A. 102 (suppl. 1), 6543-6549 (2005)]. European eels have a complex life cycle, including transitions between life stages across ecological conditions in the Sargasso Sea, where spawning occurs, and those in brackish and freshwater bodies from northern Europe to northern Africa. Whether continental eel populations consist of locally adapted and genetically distinct populations or comprise a single panmictic population has received conflicting support. Here we use whole-genome sequencing and show that European eels belong to one panmictic population. A complete lack of geographical genetic differentiation is demonstrated. We postulate that this is possible because the most critical life stages-spawning and embryonic development-take place under near-identical conditions in the Sargasso Sea. We further show that within-generation selection, which has recently been proposed as a mechanism for genetic adaptation in eels, can only marginally change allele frequencies between cohorts of eels from different geographic regions. Our results strongly indicate plasticity as the predominant mechanism for how eels respond to diverse environmental conditions during postlarval stages, ultimately solving a long-standing question for a classically enigmatic species.

Ecological adaptation in Atlantic herring is associated with large shifts in allele frequencies at hundreds of loci

Atlantic herring is widespread in North Atlantic and adjacent waters and is one of the most abundant vertebrates on earth. This species is well suited to explore genetic adaptation due to minute genetic differentiation at selectively neutral loci. Here, we report hundreds of loci underlying ecological adaptation to different geographic areas and spawning conditions. Four of these represent megabase inversions confirmed by long read sequencing. The genetic architecture underlying ecological adaptation in herring deviates from expectation under a classical infinitesimal model for complex traits because of large shifts in allele frequencies at hundreds of loci under selection.

Structure and Characterization of Phosphoglucomutase 5 from Atlantic and Baltic Herring-An Inactive Enzyme with Intact Substrate Binding

Phosphoglucomutase 5 (PGM5) in humans is known as a structural muscle protein without enzymatic activity, but detailed understanding of its function is lacking. PGM5 belongs to the alpha-D-phosphohexomutase family and is closely related to the enzymatically active metabolic enzyme PGM1. In the Atlantic herring, Clupea harengus, PGM5 is one of the genes strongly associated with ecological adaptation to the brackish Baltic Sea. We here present the first crystal structures of PGM5, from the Atlantic and Baltic herring, differing by a single substitution Ala330Val. The structure of PGM5 is overall highly similar to structures of PGM1. The structure of the Baltic herring PGM5 in complex with the substrate glucose-1-phosphate shows conserved substrate binding and active site compared to human PGM1, but both PGM5 variants lack phosphoglucomutase activity under the tested conditions. Structure comparison and sequence analysis of PGM5 and PGM1 from fish and mammals suggest that the lacking enzymatic activity of PGM5 is related to differences in active-site loops that are important for flipping of the reaction intermediate. The Ala330Val substitution does not alter structure or biophysical properties of PGM5 but, due to its surface-exposed location, could affect interactions with protein-binding partners.

Apolipoprotein B48 metabolism in chylomicrons and very low-density lipoproteins and its role in triglyceride transport in normo- and hypertriglyceridemic human subjects

BACKGROUND

Renewed interest in triglyceride-rich lipoproteins as causative agents in cardiovascular disease mandates further exploration of the integrated metabolism of chylomicrons and very low-density lipoproteins (VLDL).

METHODS

Novel tracer techniques and an integrated multi-compartmental model were used to determine the kinetics of apoB48- and apoB100-containing particles in the chylomicron and VLDL density intervals in 15 subjects with a wide range of plasma triglyceride levels.

RESULTS

Following a fat-rich meal, apoB48 appeared in the chylomicron, VLDL₁ and VLDL₂ fractions in all subjects. Chylomicrons cleared rapidly from the circulation but apoB48-containing VLDL accumulated, and over the day were 3-fold higher in those with high versus low plasma triglyceride. ApoB48-containing particles were secreted directly into both the chylomicron and VLDL fractions at rates that were similar across the plasma triglyceride range studied. During fat absorption, whilst most triglyceride entered the circulation in chylomicrons, the majority of apoB48 particles were secreted into the VLDL density range.

CONCLUSION

The intestine secretes apoB48-containing particles not only as chylomicrons but also directly into the VLDL₁ and VLDL₂ density ranges both in the basal state and during dietary lipid absorption. Over the day, apoB48-containing particles appear to comprise about 20-25% of circulating VLDL and, especially in those with elevated triglycerides, form part of a slowly cleared 'remnant' particle population, thereby potentially increasing CHD risk. These findings provide a metabolic understanding of the potential consequences for increased CHD risk when slowed lipolysis leads to the accumulation of remnants, especially in individuals with hypertriglyceridemia.

The importance of the ZBED6-IGF2 axis for metabolic regulation in mouse myoblast cells

The transcription factor ZBED6 acts as a repressor of Igf2 and affects directly or indirectly the transcriptional regulation of thousands of genes. Here, we use gene editing in mouse C2C12 myoblasts and show that ZBED6 regulates Igf2 exclusively through its binding site 5'-GGCTCG-3' in intron 1 of Igf2. Deletion of this motif (Igf2^ΔGGCT ) or complete ablation of Zbed6 leads to ~20-fold upregulation of the IGF2 protein. Quantitative proteomics revealed an activation of Ras signaling pathway in both Zbed6^-/- and Igf2^ΔGGCT myoblasts, and a significant enrichment of mitochondrial membrane proteins among proteins showing altered expression in Zbed6^-/- myoblasts. Both Zbed6^-/- and Igf2^ΔGGCT myoblasts showed a faster growth rate and developed myotube hypertrophy. These cells exhibited an increased O₂ consumption rate, due to IGF2 upregulation. Transcriptome analysis revealed ~30% overlap between differentially expressed genes in Zbed6^-/- and Igf2^ΔGGCT myotubes, with an enrichment of upregulated genes involved in muscle development. In contrast, ZBED6-overexpression in myoblasts led to cell apoptosis, cell cycle arrest, reduced mitochondrial activities, and ceased myoblast differentiation. The similarities in growth and differentiation phenotypes observed in Zbed6^-/- and Igf2^ΔGGCT myoblasts demonstrates that ZBED6 affects mitochondrial activity and myogenesis largely through its regulation of IGF2 expression. This study adds new insights how the ZBED6-Igf2 axis affects muscle metabolism.

A chromosome-level assembly of the Atlantic herring genome-detection of a supergene and other signals of selection

The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.

Exploring a Pool-seq-only approach for gaining population genomic insights in nonmodel species

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool-seq data to generate a de novo genome assembly for mining exons, upon which Pool-seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual-based single nucleotide polymorphisms [SNPs]) and from mapping the Pool-seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F ST) between the two introduced populations exceeds that of the naturally sympatric populations (F ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( π ¯  ≈ 0.002 and π ¯  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high-quality reference assembly from a divergent species. We conclude that the Pool-seq-only approach can be suitable for detecting and quantifying genome-wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.

Characterization of the endogenous retrovirus insertion in CYP19A1 associated with henny feathering in chicken

Background

Henny feathering in chickens is determined by a dominant mutation that transforms male-specific plumage to female-like plumage. Previous studies indicated that this phenotype is caused by ectopic expression in skin of CYP19A1 encoding aromatase that converts androgens to estrogen and thereby inhibits the development of male-specific plumage. A long terminal repeat (LTR) from an uncharacterized endogenous retrovirus (ERV) insertion was found in an isoform of the CYP19A1 transcript from henny feathering chicken. However, the complete sequence and the genomic position of the insertion were not determined.

Results

We used publicly available whole genome sequence data to determine the flanking sequences of the ERV, and then PCR amplified the entire insertion and sequenced it using Nanopore long reads and Sanger sequencing. The 7524 bp insertion contains an intact endogenous retrovirus that was not found in chickens representing 31 different breeds not showing henny feathering or in samples of the ancestral red junglefowl. The sequence shows over 99% sequence identity to the avian leukosis virus ev-1 and ev-21 strains, suggesting a recent integration. The ERV 3’LTR, containing a powerful transcriptional enhancer and core promoter with TATA box together with binding sites for EFIII and Ig/EBP inside the CYP19A1 5′ untranslated region, was detected partially in an aromatase transcript, which present a plausible explanation for ectopic expression of aromatase in non-ovarian tissues underlying the henny feathering phenotype.

Conclusions

We demonstrate that the henny feathering allele harbors an insertion of an intact avian leukosis virus at the 5’end of CYP19A1. The presence of this ERV showed complete concordance with the henny feathering phenotype both within a pedigree segregating for this phenotype and across breeds.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0181-4) contains supplementary material, which is available to authorized users.

Identification and validation of genetic variants predictive of gait in standardbred horses

Several horse breeds have been specifically selected for the ability to exhibit alternative patterns of locomotion, or gaits. A premature stop codon in the gene DMRT3 is permissive for “gaitedness” across breeds. However, this mutation is nearly fixed in both American Standardbred trotters and pacers, which perform a diagonal and lateral gait, respectively, during harness racing. This suggests that modifying alleles must influence the preferred gait at racing speeds in these populations. A genome-wide association analysis for the ability to pace was performed in 542 Standardbred horses (n = 176 pacers, n = 366 trotters) with genotype data imputed to ~74,000 single nucleotide polymorphisms (SNPs). Nineteen SNPs on nine chromosomes (ECA1, 2, 6, 9, 17, 19, 23, 25, 31) reached genome-wide significance (p < 1.44 x 10⁻⁶). Variant discovery in regions of interest was carried out via whole-genome sequencing. A set of 303 variants from 22 chromosomes with putative modifying effects on gait was genotyped in 659 Standardbreds (n = 231 pacers, n = 428 trotters) using a high-throughput assay. Random forest classification analysis resulted in an out-of-box error rate of 0.61%. A conditional inference tree algorithm containing seven SNPs predicted status as a pacer or trotter with 99.1% accuracy and subsequently performed with 99.4% accuracy in an independently sampled population of 166 Standardbreds (n = 83 pacers, n = 83 trotters). This highly accurate algorithm could be used by owners/trainers to identify Standardbred horses with the potential to race as pacers or as trotters, according to the genotype identified, prior to initiating training and would enable fine-tuning of breeding programs with designed matings. Additional work is needed to determine both the algorithm’s utility in other gaited breeds and whether any of the predictive SNPs play a physiologically functional role in the tendency to pace or tag true functional alleles.

Certain horse breeds have been developed over generations specifically for the ability to perform alternative patterns of movement, or gaits. Current understanding of the genetic basis for these gaits is limited to one known mutation apparently necessary, but not sufficient, for explaining variability in “gaitedness.” The Standardbred breed includes two distinct groups, trotters, which exhibit a two-beat gait in which the opposite forelimb and hind limb move together, and pacers, which exhibit an alternative two-beat gait where the legs on the same side of the body move together. Our long-term objective is to identify variants underlying the ability of certain Standardbreds to pace. In this study, we were able to identify several regions of the genome highly associated with pacing and, within these regions, a number of specific highly associated variants. Although the biological function of these variants has yet to be determined, we developed a model based on seven variants that was > 99% accurate in predicting whether an individual was a pacer or a trotter in two independent populations. This predictive model can be used by horse owners to make breeding and training decisions related to this economically important trait, and by scientists interested in understanding the biology of coordinated gait development.