Origin of agouti-melanistic polymorphism in Wild Black Rats (Rattus rattus) inferred from Mc1r gene sequences

Region-Specific Genetic Diversity of Black Rats (Rattus rattus Complex) in Southeast and East Asia Shaped by Rapid Population Expansion Events

Among the six mitochondrial DNA lineages of the black rat (Rattus rattus Complex; RrC), lineages II and IV are widespread in Southeast and East Asia. This study explored their demographic history using 17 new sequences from the Miyako Islands in the Ryukyu archipelago, together with 178 publicly available cytochrome b sequences. We defined six and two haplotype groups showing rapid population expansion signals in Lineages II and IV, respectively. The six haplotype groups of Lineage II were represented by haplotypes from 1) Myanmar/Bangladesh/Northeast India, 2) Laos, 3) Thailand, 4) Indonesia/Philippines, 5) Vietnam/southern China, and 6) the Ryukyu archipelago. These expansion times were estimated using time-dependent evolutionary rates to be 115,300 years ago (ya), 128,500 ya, 9600 ya, 10,600 ya, 7200 ya, and 1400 ya, respectively, although all had large confidence intervals. The two groups of Lineage IV were recovered from the mainland and islands of Southeast Asia with predicted expansion times of 197,000 ya and 5800 ya, respectively. These results suggest that climatic fluctuations during the last 200,000 years of the Quaternary, affected the population dynamics in subtropical areas at different times. Furthermore, the results of the younger rapid expansion events of RrC suggest the possibility of agricultural advancement and dispersal of Neolithic farmers to different areas within the mainland and islands of Southeast Asia during the Holocene. A subset of rats from the Miyako Islands were found to have the same lineage IV haplotypes as those in Southeast Asia, suggesting a recent introduction of these new lineages.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

A 1-bp deletion in Mc1r in a Norway rat (Rattus norvegicus) from Sado Island, Japan gives rise to a yellowish color variant: an insight into mammalian MC1R variants

The melanocortin-1 receptor gene (MC1R) controls production of the pigments eumelanin and pheomelanin. Changes in MC1R lead to variation in coat color in mammals, which can range from entirely black (melanism) to yellowish. In this study, we report a case of a wild-caught Norway rat (Rattus norvegicus) from Sado Island, Japan with a yellowish coat color. Upon sequencing the whole coding region of the Mc1r gene (954 bp), we found a 1-bp deletion at site 337 (c.337del), indicative of a frameshift mutation, which was characterized as a severe loss-of-function or null mutation. A spectrophotometer was used to measure coat color, revealing that the rat had a distinctly lighter coat, based on lightness score, than mice with homozygous similar loss-of-function mutations. This implies that loss-of-function mutations can yield different phenotypes in murine rodents. The loss-of-function-mutant rat exhibited a contrasting coat pattern consisting of darker and lighter colors along its dorsal and ventral sides, respectively. Similar patterns have been observed in homozygous MC1R-deficient mutants in other mammals, implying that the countershading pattern can still be expressed despite the absence of MC1R in the melanocyte.

Mouse Hair Significantly Lightened Through Replacement of the Cysteine Residue in the N-Terminal Domain of Mc1r Using the CRISPR/Cas9 System

A loss-of-function mutation in the melanocortin 1 receptor gene (MC1R), which switches off the eumelanin production, causes yellowish coat color variants in mammals. In a wild population of sables (Martes zibellina) in Hokkaido, Japan, the mutation responsible for a bright yellow coat color variant was inferred to be a cysteine replacement at codon 35 of the N-terminal extracellular domain of the Mc1r receptor. In the present study, we validated these findings by applying genome editing on Mc1r in mouse strains C3H/HeJ and C57BL/6N, altering the codon for cysteine (Cys33Phe). The resulting single amino acid substitution (Cys33Phe) and unintentionally generated frameshift mutations yielded a color variant exhibiting substantially brighter body color, indicating that the Cys35 replacement produced sufficient MC1R loss of function to confirm that this mutation is responsible for producing the Hokkaido sable yellow color variant. Notably, the yellowish mutant mouse phenotype exhibited brown coloration in subapical hair on the dorsal side in both the C3H/HeJ and C57BL/6N strains, despite the inability of the latter to produce the agouti signaling protein (Asip). This darker hair and body coloration was not apparent in the Hokkaido sable variant, implying the presence of an additional genetic system shaping yellowish hair variability.

Multiple Conserved Elements Structuring Inverted Repeats in the Mammalian Coat Color-Related Gene Asip

In the agouti signaling gene protein (Asip) of the house mouse (Mus musculus), inverted repeat (IR) arrays are known to exist in a non-coding region adjacent to the ventral-specific promoter region and the accompanying two exons (exons 1A and 1A'), which are around 100 kb upstream from the amino acid coding regions of exons 2, 3, and 4. To determine the gene structure of mammalian Asip and to elucidate trends in its evolution, non-coding sequences of six rodent (mouse, rat, Chinese hamster, squirrel, guinea pig, and naked mole rat) and three non-rodent (rabbit, human, and cow) species were retrieved from databases and compared. Our homology search analyses revealed the presence of three to five highly conserved non-coding elements (CNE). These CNEs were found to form IRs in rodents and lagomorphs. Combinations of IRs were further shown to build symmetric, long IR arrays. Intra- and inter-specific comparisons of the sequences of three universal CNEs showed homogeneity between CNE pairs within species. This implies that certain evolutionary constraints maintained the IR structure in the rodent and rabbit species.

Potential Causative Mutation for Melanism in Rats Identified in the Agouti Signaling Protein Gene (Asip) of the Rattus rattus Species Complex on Okinawa Island, Japan

The occurrence of black fur, or melanism, in many mammalian species is known to be linked to DNA sequence variation in the agouti signaling protein (Asip) gene, which is a major determinant of eumelanin and pheomelanin pigments in coat color. We investigated 38 agouti (i.e., banded wildtype) and four melanistic Rattus rattus species complex (RrC) lineage II specimens from Okinawa Island, Ryukyu Islands, Japan, for genetic variation in three exons and associated flanking regions in the Asip gene. On Okinawa, a predicted loss-of-function mutation caused by a cysteine to serine amino acid change at p.124C>S (c.370T>A) in the highly conserved functional domain of Asip was found in melanistic rats, but was absent in agouti specimens, suggesting that the p.124C>S mutation is responsible for the observed melanism. Phylogeographic analysis found that Asip sequences from Okinawan RrC lineage II, including both agouti and melanistic specimens, differed from: 1) both agouti and melanistic RrC lineage I from Otaru, Hokkaido, Japan, and 2) agouti RrC lineages I and II from South Australia. This suggests the possibility of in-situ mutation of the Asip gene, either within the RrC lineage II population on Okinawa or in an unsampled RrC lineage II population with biogeographic links to Okinawa, although incomplete lineage sorting could not be ruled out.

Molecular phylogenetic characterization of common murine rodents from Manipur, Northeast India

The Indian subcontinent and Southeast Asia are hotspots of murine biodiversity, but no species from the Arakan Mountain system that demarcates the border between the two areas has been subjected to molecular phylogenetic analyses. We examined the mitochondrial cytochrome b gene sequences in six murine species (the Rattus rattus species complex, R. norvegicus, R. nitidus, Berylmys manipulus, Niviventer sp. and Mus musculus) from Manipur, which is located at the western foot of the mountain range. The sequences of B. manipulus and Niviventer sp. examined here were distinct from available congeneric sequences in the databases, with sequence divergences of 10-15%. Substantial degrees of intrapopulation divergence were detected in R. nitidus and the R. rattus species complex from Manipur, implying ancient habitation of the species in this region, while the recent introduction by modern and prehistoric human activities was suggested for R. norvegicus and M. musculus, respectively. In the nuclear gene Mc1r, also analyzed here, the R. rattus species complex from Manipur was shown to possess allelic sequences related to those from the Indian subcontinent in addition to those from East Asia. These results not only fill gaps in the phylogenetic knowledge of each taxon examined but also provide valuable insight to better understand the biogeographic importance of the Arakan Mountain system in generating the species and genetic diversity of murine rodents.

Genetic characterization of Okinawan black rats showing coat color polymorphisms of white spotting and melanism

We examined pelage color variation in wild populations of black rats (the Rattus rattus species complex) in the Yambaru forest area, northern Okinawa Island, Ryukyu Archipelago, Japan. Our field study revealed that 8.7% (38/438) and 0.2% (4/2500) of rats exhibited two types of coat color: white spotting and melanism, respectively. Using 34 representative animals, the phylogeography of the population was inferred using a nuclear gene marker, i.e., sequences (954 bp) of the melanocortin-1 receptor (Mc1r) gene responsible for the melanistic form in black rats. Four sequences from Okinawa were characterized as R. tanezumi, the Asian strain of black rat. Notably, neither of the phenotypic characters of white spotting or melanism was associated with the Mc1r haplotypes. Analysis of mitochondrial cytochrome b (Cytb) sequences (1140 bp) revealed that four haplotypes recovered from Okinawa clustered with the clade of R. tanezumi and differed by one or more bases from haplotypes at other localities in Japan and Asian countries. Thus, both variants may have arisen in the native rat population of Okinawa without interaction with the lineage of R. rattus, which exhibits a worldwide distribution and displays such coat color variants. The Yambaru population of black rats has thus experienced its own evolutionary history in allopatry for a substantial period of time (e.g., 10,000 years), which has preserved valuable genetic polymorphisms and will be useful for assessing the ecological consequences of genetic variation in natural populations.

Multiple geographic origins of commensalism and complex dispersal history of Black Rats

The Black Rat (Rattus rattus) spread out of Asia to become one of the world's worst agricultural and urban pests, and a reservoir or vector of numerous zoonotic diseases, including the devastating plague. Despite the global scale and inestimable cost of their impacts on both human livelihoods and natural ecosystems, little is known of the global genetic diversity of Black Rats, the timing and directions of their historical dispersals, and the risks associated with contemporary movements. We surveyed mitochondrial DNA of Black Rats collected across their global range as a first step towards obtaining an historical genetic perspective on this socioeconomically important group of rodents. We found a strong phylogeographic pattern with well-differentiated lineages of Black Rats native to South Asia, the Himalayan region, southern Indochina, and northern Indochina to East Asia, and a diversification that probably commenced in the early Middle Pleistocene. We also identified two other currently recognised species of Rattus as potential derivatives of a paraphyletic R. rattus. Three of the four phylogenetic lineage units within R. rattus show clear genetic signatures of major population expansion in prehistoric times, and the distribution of particular haplogroups mirrors archaeologically and historically documented patterns of human dispersal and trade. Commensalism clearly arose multiple times in R. rattus and in widely separated geographic regions, and this may account for apparent regionalism in their associated pathogens. Our findings represent an important step towards deeper understanding the complex and influential relationship that has developed between Black Rats and humans, and invite a thorough re-examination of host-pathogen associations among Black Rats.