A method to analyze gene expression profiles from hippocampal neurons electrophysiologically recorded in vivo.

Hippocampal pyramidal neurons exhibit diverse spike patterns and gene expression profiles. However, their relationships with single neurons are not fully understood. In this study, we designed an electrophysiology-based experimental procedure to identify gene expression profiles using RNA sequencing of single hippocampal pyramidal neurons whose spike patterns were recorded in living mice. This technique involves a sequence of experiments consisting of in vivo juxtacellular recording and labeling, brain slicing, cell collection, and transcriptome analysis. We demonstrated that the expression levels of a subset of genes in individual hippocampal pyramidal neurons were significantly correlated with their spike burstiness, submillisecond-level spike rise times or spike rates, directly measured by in vivo electrophysiological recordings. Because this methodological approach can be applied across a wide range of brain regions, it is expected to contribute to studies on various neuronal heterogeneities to understand how physiological spike patterns are associated with gene expression profiles.

Neuronal DSCAM regulates the peri-synaptic localization of GLAST in Bergmann glia for functional synapse formation.

In the central nervous system, astrocytes enable appropriate synapse function through glutamate clearance from the synaptic cleft; however, it remains unclear how astrocytic glutamate transporters function at peri-synaptic contact. Here, we report that Down syndrome cell adhesion molecule (DSCAM) in Purkinje cells controls synapse formation and function in the developing cerebellum. Dscam-mutant mice show defects in CF synapse translocation as is observed in loss of function mutations in the astrocytic glutamate transporter GLAST expressed in Bergmann glia. These mice show impaired glutamate clearance and the delocalization of GLAST away from the cleft of parallel fibre (PF) synapse. GLAST complexes with the extracellular domain of DSCAM. Riluzole, as an activator of GLAST-mediated uptake, rescues the proximal impairment in CF synapse formation in Purkinje cell-selective Dscam-deficient mice. DSCAM is required for motor learning, but not gross motor coordination. In conclusion, the intercellular association of synaptic and astrocyte proteins is important for synapse formation and function in neural transmission.

A predisposed motor bias shapes individuality in vocal learning.

The development of individuality during learned behavior is a common trait observed across animal species; however, the underlying biological mechanisms remain understood. Similar to human speech, songbirds develop individually unique songs with species-specific traits through vocal learning. In this study, we investigate the developmental and molecular mechanisms underlying individuality in vocal learning by utilizing F1 hybrid songbirds (Taeniopygia guttata cross with Taeniopygia bichenovii), taking an integrating approach combining experimentally controlled systematic song tutoring, unbiased discriminant analysis of song features, and single-cell transcriptomics. When tutoring with songs from both parental species, F1 hybrid individuals exhibit evident diversity in their acquired songs. Approximately 30% of F1 hybrids selectively learn either song of the two parental species, while others develop merged songs that combine traits from both species. Vocal acoustic biases during vocal babbling initially appear as individual differences in songs among F1 juveniles and are maintained through the sensitive period of song vocal learning. These vocal acoustic biases emerge independently of the initial auditory experience of hearing the biological father's and passive tutored songs. We identify individual differences in transcriptional signatures in a subset of cell types, including the glutamatergic neurons projecting from the cortical vocal output nucleus to the hypoglossal nuclei, which are associated with variations of vocal acoustic features. These findings suggest that a genetically predisposed vocal motor bias serves as the initial origin of individual variation in vocal learning, influencing learning constraints and preferences.

The admixed brushtail possum genome reveals invasion history in New Zealand and novel imprinted genes.

Combining genome assembly with population and functional genomics can provide valuable insights to development and evolution, as well as tools for species management. Here, we present a chromosome-level genome assembly of the common brushtail possum (Trichosurus vulpecula), a model marsupial threatened in parts of their native range in Australia, but also a major introduced pest in New Zealand. Functional genomics reveals post-natal activation of chemosensory and metabolic genes, reflecting unique adaptations to altricial birth and delayed weaning, a hallmark of marsupial development. Nuclear and mitochondrial analyses trace New Zealand possums to distinct Australian subspecies, which have subsequently hybridised. This admixture allowed phasing of parental alleles genome-wide, ultimately revealing at least four genes with imprinted, parent-specific expression not yet detected in other species (MLH1, EPM2AIP1, UBP1 and GPX7). We find that reprogramming of possum germline imprints, and the wider epigenome, is similar to eutherian mammals except onset occurs after birth. Together, this work is useful for genetic-based control and conservation of possums, and contributes to understanding of the evolution of novel mammalian epigenetic traits.

Ghost introgression in ricefishes of the genus Adrianichthys in an ancient Wallacean lake.

Because speciation might have been promoted by ancient introgression from an extinct lineage, it is important to detect the existence of 'ghost introgression' in focal taxa and examine its contribution to their diversification. In this study, we examined possible ghost introgression and its contributions to the diversification of ricefishes of the genus Adrianichthys in Lake Poso, an ancient lake on Sulawesi Island, in which some extinctions are known to have occurred. Population-genomic analysis revealed that two extant Adrianichthys species, A. oophorus and A. poptae are reproductively isolated from each other. Comparisons of demographic models demonstrated that introgression from a ghost population, which diverged from the common ancestor of A. oophorus and A. poptae, is essential for reconstructing the demographic history of Adrianichthys. The best model estimated that the divergence of the ghost population greatly predated the divergence between A. oophorus and A. poptae, and that the ghost population secondarily contacted the two extant species within Lake Poso more recently. Genome scans and simulations detected a greatly divergent locus, which cannot be explained without ghost introgression. This locus was also completely segregated between A. oophorus and A. poptae. These findings suggest that variants that came from a ghost population have contributed to the divergence between A. oophorus and A. poptae, but the large time-lag between their divergence and ghost introgression indicates that the contribution of introgression may be restricted.

Single-cell RNA sequencing of intestinal immune cells in neonatal necrotizing enterocolitis.

PURPOSE: Necrotizing enterocolitis (NEC) causes fatal intestinal necrosis in neonates, but its etiology is unknown. We analyzed the intestinal immune response to NEC. METHODS: Using single-cell RNA sequencing (scRNA-seq), we analyzed the gene expression profiles of intestinal immune cells from four neonates with intestinal perforation (two with NEC and two without NEC). Target mononuclear cells were extracted from the lamina propria of the resected intestines. RESULTS: In all four cases, major immune cells, such as T cells (15.1-47.7%), B cells (3.1-19.0%), monocytes (16.5-31.2%), macrophages (1.6-17.4%), dendritic cells (2.4-12.2%), and natural killer cells (7.5-12.8%), were present in similar proportions to those in the neonatal cord blood. Gene set enrichment analysis showed that the MTOR, TNF-α, and MYC signaling pathways were enriched in T cells of the NEC patients, suggesting upregulated immune responses related to inflammation and cell proliferation. In addition, all four cases exhibited a bias toward cell-mediated inflammation, based on the predominance of T helper 1 cells. CONCLUSION: Intestinal immunity in NEC subjects exhibited stronger inflammatory responses compared to non-NEC subjects. Further scRNA-seq and cellular analysis may improve our understanding of the pathogenesis of NEC.

Gene expression profiling of the masticatory muscle tendons and Achilles tendons under tensile strain in the Japanese macaque Macaca fuscata.

Both Achilles and masticatory muscle tendons are large load-bearing structures, and excessive mechanical loading leads to hypertrophic changes in these tendons. In the maxillofacial region, hyperplasia of the masticatory muscle tendons and aponeurosis affect muscle extensibility resulting in limited mouth opening. Although gene expression profiles of Achilles and patellar tendons under mechanical strain are well investigated in rodents, the gene expression profile of the masticatory muscle tendons remains unexplored. Herein, we examined the gene expression pattern of masticatory muscle tendons and compared it with that of Achilles tendons under tensile strain conditions in the Japanese macaque Macaca fuscata. Primary tenocytes isolated from the masticatory muscle tendons (temporal tendon and masseter aponeurosis) and Achilles tendons were mechanically loaded using the tensile force and gene expression was analyzed using the next-generation sequencing. In tendons exposed to tensile strain, we identified 1076 differentially expressed genes with a false discovery rate (FDR) < 10-10. To identify genes that are differentially expressed in temporal tendon and masseter aponeurosis, an FDR of < 10-10 was used, whereas the FDR for Achilles tendons was set at > 0.05. Results showed that 147 genes are differentially expressed between temporal tendons and masseter aponeurosis, out of which, 125 human orthologs were identified using the Ensemble database. Eight of these orthologs were related to tendons and among them the expression of the glycoprotein nmb and sphingosine kinase 1 was increased in temporal tendons and masseter aponeurosis following exposure to tensile strain. Moreover, the expression of tubulin beta 3 class III, which promotes cell cycle progression, and septin 9, which promotes cytoskeletal rearrangements, were decreased in stretched Achilles tendon cells and their expression was increased in stretched masseter aponeurosis and temporal tendon cells. In conclusion, cyclic strain differentially affects gene expression in Achilles tendons and tendons of the masticatory muscles.

Cognitive genomics of learning delay and low level of social performance monitoring in macaque.

Cognitive skills and the underlying neural architecture are under the influence of genetics. Cognitive genomics research explores the triadic relationship between genes, brain, and cognition, with its major strategy being genotype-driven. Here we show that an inverse strategy is feasible to identify novel candidate genes for particular neuro-cognitive phenotypes in macaques. Two monkeys, originally involved in separate psychological studies, exhibited learning delay and low levels of social performance monitoring. In one monkey, mirror neurons were fewer compared to controls and mu suppression was absent in the frontal cortex. The other monkey showed heightened visual responsiveness in both frontal cortex and dopamine-rich midbrain, with a lack of inter-areal synchronization. Exome analyses revealed that the two monkeys were most likely cousins and shared variants in MAP2, APOC1, and potentially HTR2C. This phenotype-driven strategy in cognitive genomics provides a useful means to clarify the genetic basis of phenotypic variation and develop macaque models of neuropsychiatric disorders.

Nicotinic acetylcholine receptors in a songbird brain.

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that mediate fast synaptic transmission and cell signaling, which contribute to learning, memory, and the execution of motor skills. Birdsong is a complex learned motor skill in songbirds. Although the existence of 15 nAChR subunits has been predicted in the avian genome, their expression patterns and potential contributions to song learning and production have not been comprehensively investigated. Here, we cloned all the 15 nAChR subunits (ChrnA1-10, B2-4, D, and G) from the zebra finch brain and investigated the mRNA expression patterns in the neural pathways responsible for the learning and production of birdsong during a critical period of song learning. Although there were no detectable hybridization signals for ChrnA1, A6, A9, and A10, the other 11 nAChR subunits were uniquely expressed in one or more major subdivisions in the song nuclei of the songbird brain. Of these 11 subunits, ChrnA3-5, A7, and B2 were differentially regulated in the song nuclei compared with the surrounding anatomically related regions. ChrnA5 was upregulated during the critical period of song learning in the lateral magnocellular nucleus of the anterior nidopallium. Furthermore, single-cell RNA sequencing revealed ChrnA7 and B2 to be the major subunits expressed in neurons of the vocal motor nuclei HVC and robust nucleus of the arcopallium, indicating the potential existence of ChrnA7-homomeric and ChrnB2-heteromeric nAChRs in limited cell populations. These results suggest that relatively limited types of nAChR subunits provide functional contributions to song learning and production in songbirds.

Uterus-specific transcriptional regulation underlies eggshell pigment production in Japanese quail.

The precursor of heme, protoporphyrin IX (PPIX), accumulates abundantly in the uteri of birds, such as Japanese quail, Coturnix japonica, which has brown-speckled eggshells; however, the molecular basis of PPIX production in the uterus remains largely unknown. Here, we investigated the cause of low PPIX production in a classical Japanese quail mutant exhibiting white eggshells by comparing its gene expression in the uterus with that of the wild type using transcriptome analysis. We also performed genetic linkage analysis to identify the causative genomic region of the white eggshell phenotype. We found that 11 genes, including 5'-aminolevulinate synthase 1 (ALAS1) and hephaestin-like 1 (HEPHL1), were specifically upregulated in the wild-type uterus and downregulated in the mutant. We mapped the 172 kb candidate genomic region on chromosome 6, which contains several genes, including a part of the paired-like homeodomain 3 (PITX3), which encodes a transcription factor. ALAS1, HEPHL1, and PITX3 were expressed in the apical cells of the luminal epithelium and lamina propria cells of the uterine mucosa of the wild-type quail, while their expression levels were downregulated in the cells of the mutant quail. Biochemical analysis using uterine homogenates indicated that the restricted availability of 5'-aminolevulinic acid is the main cause of low PPIX production. These results suggest that uterus-specific transcriptional regulation of heme-biosynthesis-related genes is an evolutionarily acquired mechanism of eggshell pigment production in Japanese quail. Based on these findings, we discussed the molecular basis of PPIX production in the uteri of Japanese quails.

In vitro monitoring of HTR2A-positive neurons derived from human-induced pluripotent stem cells.

The serotonin 5-HT2A receptor (5-HT2AR) has been receiving increasing attention because its genetic variants have been associated with a variety of neurological diseases. To elucidate the pathogenesis of the neurological diseases associated with 5-HT2AR gene (HTR2A) variants, we have previously established a protocol to induce HTR2A-expressing neurons from human-induced pluripotent stem cells (hiPSCs). Here, we investigated the maturation stages and electrophysiological properties of HTR2A-positive neurons induced from hiPSCs and constructed an HTR2A promoter-specific reporter lentivirus to label the neurons. We found that neuronal maturity increased over time and that HTR2A expression was induced at the late stage of neuronal maturation. Furthermore, we demonstrated successful labelling of the HTR2A-positive neurons, which had fluorescence and generated repetitive action potentials in response to depolarizing currents and an inward current during the application of TCB-2, a selective agonist of 5-HT2ARs, respectively. These results indicated that our in vitro model mimicked the in vivo dynamics of 5-HT2AR. Therefore, in vitro monitoring of the function of HTR2A-positive neurons induced from hiPSCs could help elucidate the pathophysiological mechanisms of neurological diseases associated with genetic variations of the HTR2A gene.

Population history and genomic admixture of sea snakes of the genus Laticauda in the West Pacific.

Speciation in the open ocean has long been studied, but it remains largely elusive what factors promote or inhibit speciation in such an open environment. Marine amniotes, which evolved from terrestrial ancestors, provide valuable opportunities for studying speciation in the ocean because of their evident aquatic origins. Sea snakes are phylogenetically related to terrestrial elapid snakes and consist of two monophyletic groups (Hydrophiini and Laticaudini). These two groups migrated from land to water almost at the same time, but species diversities are remarkably different: there are approx. 60 species in 16 genera described for hydrophiins, whereas only eight species in the genus Laticauda are described for laticaudins. Here, we provide a high-quality reference genome assembly of a laticaudin L. colubrina with a scaffold N50 value of 40 Mbp, and focused on laticaudins to consider why they have seldom speciated. We performed whole-genome shotgun sequencing of several species of laticaudins sampled in their southmost (Vanuatu) and northmost (Ryukyu) habitats. Demographic histories of Vanuatu and Ryukyu populations suggest that populations of broadly distributed major species are geographically structured. Each species is genetically clearly distinguished, but there is a considerable amount of gene flow between two sibling species distributed sympatrically in Vanuatu. In addition, inter-species genomic admixture is ubiquitously observed among laticaudins even between phylogenetically distant species. Broad distribution of major species combined with such genetic mixability might have prevented laticaudins from genetic isolation and speciation.

Redundant type II cadherins define neuroepithelial cell states for cytoarchitectonic robustness.

Individual cell shape and integrity must precisely be orchestrated during morphogenesis. Here, we determine function of type II cadherins, Cdh6, Cdh8, and Cdh11, whose expression combinatorially demarcates the mouse neural plate/tube. While CRISPR/Cas9-based single type II cadherin mutants show no obvious phenotype, Cdh6/8 double knockout (DKO) mice develop intermingled forebrain/midbrain compartments as these two cadherins' expression opposes at the nascent boundary. Cdh6/8/11 triple, Cdh6/8 or Cdh8/11 DKO mice further cause exencephaly just within the cranial region where mutated cadherins' expression merges. In the Cdh8/11 DKO midbrain, we observe less-constricted apical actin meshwork, ventrally-directed spreading, and occasional hyperproliferation among dorsal neuroepithelial cells as origins for exencephaly. These results provide rigid evidence that, by conferring distinct adhesive codes to each cell, redundant type II cadherins serve essential and shared roles in compartmentalization and neurulation, both of which proceed under the robust control of the number, positioning, constriction, and fluidity of neuroepithelial cells.

Transcriptome analysis revealed misregulated gene expression in blastoderms of interspecific chicken and Japanese quail F1 hybrids.

Hybrid incompatibility, such as sterility and inviability, prevents gene flow between closely-related populations as a reproductive isolation barrier. F1 hybrids between chickens and Japanese quail (hereafter, referred to as quail), exhibit a high frequency of developmental arrest at the preprimitive streak stage. To investigate the molecular basis of the developmental arrest at the preprimitive streak stage in chicken-quail F1 hybrid embryos, we investigated chromosomal abnormalities in the hybrid embryos using molecular cytogenetic analysis. In addition, we quantified gene expression in parental species and chicken- and quail-derived allele-specific expression in the hybrids at the early blastoderm and preprimitive streak stages by mRNA sequencing. Subsequently, we compared the directions of change in gene expression, including upregulation, downregulation, or no change, from the early blastoderm stage to the preprimitive streak stage between parental species and their hybrids. Chromosome analysis revealed that the cells of the hybrid embryos contained a fifty-fifty mixture of parental chromosomes, and numerical chromosomal abnormalities were hardly observed in the hybrid cells. Gene expression analysis revealed that a part of the genes that were upregulated from the early blastoderm stage to the preprimitive streak stage in both parental species exhibited no upregulation of both chicken- and quail-derived alleles in the hybrids. GO term enrichment analysis revealed that these misregulated genes are involved in various biological processes, including ribosome-mediated protein synthesis and cell proliferation. Furthermore, the misregulated genes included genes involved in early embryonic development, such as primitive streak formation and gastrulation. These results suggest that numerical chromosomal abnormalities due to a segregation failure does not cause the lethality of chicken-quail hybrid embryos, and that the downregulated expression of the genes that are involved in various biological processes, including translation and primitive streak formation, mainly causes the developmental arrest at the preprimitive streak stage in the hybrids.

AUTS2 Regulation of Synapses for Proper Synaptic Inputs and Social Communication.

Impairments in synapse development are thought to cause numerous psychiatric disorders. Autism susceptibility candidate 2 (AUTS2) gene has been associated with various psychiatric disorders, such as autism and intellectual disabilities. Although roles for AUTS2 in neuronal migration and neuritogenesis have been reported, its involvement in synapse regulation remains unclear. In this study, we found that excitatory synapses were specifically increased in the Auts2-deficient primary cultured neurons as well as Auts2 mutant forebrains. Electrophysiological recordings and immunostaining showed increases in excitatory synaptic inputs as well as c-fos expression in Auts2 mutant brains, suggesting that an altered balance of excitatory and inhibitory inputs enhances brain excitability. Auts2 mutant mice exhibited autistic-like behaviors including impairments in social interaction and altered vocal communication. Together, these findings suggest that AUTS2 regulates excitatory synapse number to coordinate E/I balance in the brain, whose impairment may underlie the pathology of psychiatric disorders in individuals with AUTS2 mutations.

Towards HCP-Style macaque connectomes: 24-Channel 3T multi-array coil, MRI sequences and preprocessing.

Macaque monkeys are an important animal model where invasive investigations can lead to a better understanding of the cortical organization of primates including humans. However, the tools and methods for noninvasive image acquisition (e.g. MRI RF coils and pulse sequence protocols) and image data preprocessing have lagged behind those developed for humans. To resolve the structural and functional characteristics of the smaller macaque brain, high spatial, temporal, and angular resolutions combined with high signal-to-noise ratio are required to ensure good image quality. To address these challenges, we developed a macaque 24-channel receive coil for 3-T MRI with parallel imaging capabilities. This coil enables adaptation of the Human Connectome Project (HCP) image acquisition protocols to the in-vivo macaque brain. In addition, we adapted HCP preprocessing methods to the macaque brain, including spatial minimal preprocessing of structural, functional MRI (fMRI), and diffusion MRI (dMRI). The coil provides the necessary high signal-to-noise ratio and high efficiency in data acquisition, allowing four- and five-fold accelerations for dMRI and fMRI. Automated FreeSurfer segmentation of cortex, reconstruction of cortical surface, removal of artefacts and nuisance signals in fMRI, and distortion correction of dMRI all performed well, and the overall quality of basic neurobiological measures was comparable with those for the HCP. Analyses of functional connectivity in fMRI revealed high sensitivity as compared with those from publicly shared datasets. Tractography-based connectivity estimates correlated with tracer connectivity similarly to that achieved using ex-vivo dMRI. The resulting HCP-style in vivo macaque MRI data show considerable promise for analyzing cortical architecture and functional and structural connectivity using advanced methods that have previously only been available in studies of the human brain.

MacaquePose: A Novel "In the Wild" Macaque Monkey Pose Dataset for Markerless Motion Capture.

Video-based markerless motion capture permits quantification of an animal's pose and motion, with a high spatiotemporal resolution in a naturalistic context, and is a powerful tool for analyzing the relationship between the animal's behaviors and its brain functions. Macaque monkeys are excellent non-human primate models, especially for studying neuroscience. Due to the lack of a dataset allowing training of a deep neural network for the macaque's markerless motion capture in the naturalistic context, it has been challenging to apply this technology for macaques-based studies. In this study, we created MacaquePose, a novel open dataset with manually labeled body part positions (keypoints) for macaques in naturalistic scenes, consisting of >13,000 images. We also validated the application of the dataset by training and evaluating an artificial neural network with the dataset. The results indicated that the keypoint estimation performance of the trained network was close to that of a human-level. The dataset will be instrumental to train/test the neural networks for markerless motion capture of the macaques and developments of the algorithms for the networks, contributing establishment of an innovative platform for behavior analysis for non-human primates for neuroscience and medicine, as well as other fields using macaques as a model organism.

Loss of olfaction in sea snakes provides new perspectives on the aquatic adaptation of amniotes.

Marine amniotes, a polyphyletic group, provide an excellent opportunity for studying convergent evolution. Their sense of smell tends to degenerate, but this process has not been explored by comparing fully aquatic species with their amphibious relatives in an evolutionary context. Here, we sequenced the genomes of fully aquatic and amphibious sea snakes and identified repertoires of chemosensory receptor genes involved in olfaction. Snakes possess large numbers of the olfactory receptor (OR) genes and the type-2 vomeronasal receptor (V2R) genes, and expression profiling in the olfactory tissues suggests that snakes use the ORs in the main olfactory system (MOS) and the V2Rs in the vomeronasal system (VNS). The number of OR genes has decreased in sea snakes, and fully aquatic species lost MOS which is responsible for detecting airborne odours. By contrast, sea snakes including fully aquatic species retain a number of V2R genes and a well-developed VNS for smelling underwater. This study suggests that the sense of smell also degenerated in sea snakes, particularly in fully aquatic species, but their residual olfactory capability is distinct from that of other fully aquatic amniotes. Amphibious species show an intermediate status between terrestrial and fully aquatic snakes, implying their importance in understanding the process of aquatic adaptation.

Considerable Synteny and Sequence Similarity of Primate Chromosomal Region VIIq31.

Human chromosome 7 has been the focus of many behavioral, genetic, and medical studies because it carries genes related to cancer and neurodevelopment. We examined the evolution of the chromosome 7 homologs, and the 7q31 region in particular, using chromosome painting analyses and 3 paint probes derived from (i) the whole of chimpanzee chromosome VII (wcVII), (ii) human 7q31 (h7q31), and (iii) the chimpanzee homolog VIIq31 (cVIIq31). The wcVII probe was used instead of the whole human chromosome 7 because the chimpanzee contains additional C-bands and revealed large areas of synteny conservation as well as fragmentation across 20 primate species. Analyses focusing specifically on the 7q31 homolog and vicinity revealed considerable conservation across lineages with 2 exceptions. First, the probes verified an insertion of repetitive sequence at VIIq22 in chimpanzees and bonobos and also detected the sequence in most subtelomeres of the African apes. Second, a paracentric inversion with a breakpoint in the cVIIq31 block was found in the common marmoset, confirming earlier studies. Subsequent in silico comparative genome analysis of 17 primate species revealed that VIIq31.1 is more significantly conserved at the sequence level than other regions of chromosome VII, which indicates that its components are likely responsible for critical shared traits across the order, including conditions necessary for proper human development and wellbeing.

Nonsense mutation in PMEL is associated with yellowish plumage colour phenotype in Japanese quail.

The L strain of Japanese quail exhibits a plumage phenotype that is light yellowish in colour. In this study, we identified a nonsense mutation in the premelanosome protein (PMEL) gene showing complete concordance with the yellowish plumage within a pedigree as well as across strains by genetic linkage analysis of an F2 intercross population using approximately 2,000 single nucleotide polymorphisms (SNPs) that were detected by double digest restriction site-associated DNA sequencing (ddRAD-seq). The yellowish plumage was inherited in an autosomal recessive manner, and the causative mutation was located within an 810-kb genomic region of the LGE22C19W28_E50C23 linkage group (LGE22). This region contained the PMEL gene that is required for the normal melanosome morphogenesis and eumelanin deposition. A nonsense mutation that leads to a marked truncation of the deduced protein was found in PMEL of the mutant. The gene expression level of PMEL decreased substantially in the mutant. Genotypes at the site of the nonsense mutation were fully concordant with plumage colour phenotypes in 196 F2 offspring. The nonsense mutation was not found in several quail strains with non-yellowish plumage. Thus, the yellowish plumage may be caused by the reduced eumelanin content in feathers because of the loss of PMEL function.