Next Article in Journal
BRCA1/2 Mutations in Vietnamese Patients with Hereditary Breast and Ovarian Cancer Syndrome
Previous Article in Journal
CRL4Cdt2 Ubiquitin Ligase, A Genome Caretaker Controlled by Cdt2 Binding to PCNA and DNA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Analysis Reveals Subdivision of Black Rats (Rattus rattus) in India, Origin of the Worldwide Species Spread

1
Laboratory of Molecular and Conservation Genetics, Department of Zoology, Govt. Vidarbha Institute of Science and Humanities, Amravati 444604, India
2
Department of Integrative Biology, University of Guelph, 50 Stone Road East, Guelph, ON N1G 2W1, Canada
3
Department of Ecology and Evolutionary Biology, Cornell University, Corson Hall, Ithaca, NY 14853-2701, USA
4
Department of Zoology, Amolakchand Mahavidyalaya, Godhani Road, Yavatmal 445001, India
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Genes 2022, 13(2), 267; https://doi.org/10.3390/genes13020267
Submission received: 19 December 2021 / Revised: 21 January 2022 / Accepted: 24 January 2022 / Published: 29 January 2022
(This article belongs to the Section Animal Genetics and Genomics)

Abstract

:
In contrast to the detailed and globally extensive studies on the spread of the commensal black rat, Rattus rattus, there has been relatively little work on the phylogeography of the species within India, from where this spread originated. Taking a genomic approach, we typed 27 R. rattus samples from Peninsular India using the genotyping-by-sequencing (GBS) method. Filtering and alignment of the FASTQ files yielded 1499 genome-wide SNPs. Phylogenomic tree reconstruction revealed a distinct subdivision in the R. rattus population, manifested as two clusters corresponding to the east and west coasts of India. We also identified signals of admixture between these two subpopulations, separated by an Fst of 0.20. This striking genomic difference between the east and west coast populations mirrors what has previously been described with mitochondrial DNA sequencing. It is notable that the west coast population of R. rattus has been spread globally, reflecting the origins of commensalism of the species in Western India and the subsequent transport by humans worldwide.

1. Introduction

The demographic history and origin of commensalism in rodents inevitably reflects the history of human settlements, sedentism, and ancient trades. Commensal rodents have been subject to detailed zooarchaeological studies, given this association with human history, and the profound influence of commensal rodents on human economics and health. Among rodents, black rats (R. rattus) and house mice (Mus musculus) are the world’s most widespread and destructive invasive species [1,2]. Besides causing damage to crop and agricultural products every year, R. rattus serves as the carrier of zoonotic diseases, such as plague, typhoid and leptospirosis [3]. Based on mitochondrial DNA (mtDNA) analyses, phylogeographic studies on R. rattus have revealed a long history of commensalism and maritime transport by humans [4,5]. Both archaeological and genetic evidence point to India as the origin of R. rattus and its dispersal outside of the subcontinent by maritime activities operating from its coasts [4,6,7,8]. Considering that the commensalism of R. rattus originated in India [8], geography and human history are key. The Indian Peninsula borders the Indian Ocean along two coastlines, the west and east coasts, and there is the Gangetic Plain that separates the Indian Peninsula from the Himalayas, which has been considered the most fertile region in mainland India since historic times. In the northwest of the Indian subcontinent, the rise of the ~5 kya Indus civilization is considered a major centre that provided an impetus for the origin of the earliest human civilization in South Asia [9,10], and the available evidence suggests that R. rattus first became commensal with the rise of that civilization [4,6,8]. In this way, the beginnings of agriculture, with subsequent transformation to village and city life, most probably paved the way to commensalism for the R. rattus population living in the northwest of the subcontinent. Clearly, India is very important in understanding commensalism in R. rattus, with a need for further studies.
Genetics has a long history of applications in zooarcheological studies, including providing insights into ancient human migrations and colonizations, and the domestication of animals. Previous research on M. musculus exemplifies how commensal rodents can be used as a model biomarker to study ancient human history. Several mtDNA-based studies elucidate the maritime transportation of M. musculus by Vikings [11,12,13]. Genetic evidence in support of the transport of house mice to Australia from the British Isles in the late 18th and 19th centuries by maritime activities has also been obtained [14]. The present study examines the phylogeography and population genomics of R. rattus in India using genome-wide SNP data. So far, most of the phylogeographic studies on R. rattus have focused on mtDNA, and four mitochondrial lineages (I–IV) of the R. rattus complex (RrC) have been reported in mainland Asia [4,15], with commensalism associated with the RrC lineage-I. Despite large sample sizes, most of the studies attempted on R. rattus outside India demonstrate low diversity [16,17,18]. Taking into consideration the power of genomics and the role of India in the dispersal of R. rattus in Eurasia and elsewhere, understanding the genetic diversity of R. rattus within Peninsular India is of the utmost importance. Taking a genome-wide approach, we genotyped R. rattus samples from a wide range of localities in India, with emphasis on the west and east coasts, which are regions that allowed the spread of commensal R. rattus from Peninsular India (Figure 1A). Since prehistoric times, both the west and east coasts witnessed flourishing of human settlements and kingdoms that had maritime links with West and Southeast Asia [19]. The uses of population genomic approaches have started to provide better insights into the complex histories of invasions of species worldwide [20,21,22]. Thus, the main objective of this study was to use the combined power of genome-wide SNPs to decipher the population structure in spatially distributed R. rattus populations across India, and to provide insight into commensalism in the R. rattus population in light of the knowledge of the archaeology and history of India.

2. Materials and Methods

2.1. Sample Collection and DNA Extraction

Thirty-four R. rattus were live trapped from different locations in India (Figure 1A and Table S1). After removing the tail tip, all sampled R. rattus individuals were released. The same individuals were sequenced earlier by us for Cyt b and D-loop sequencing [8]. The study was approved by the Animal Ethical Committee of the Govt. Vidarbha Institute of Science and Humanities, Amravati, India. All tail tips were preserved in ethanol until subject to genomic DNA extraction using Qiagen Blood and Tissue DNA extraction kit (Qiagen, Germantown, MD, USA). Sampling of rats was carried out in rural communities so as to avoid the possibility of admixture by recent arrival in the cities through transport by trains and buses. Three house mice, M. m. castaneus, were also sampled to be used as an outgroup.

2.2. DNA Sequencing and SNP Genotyping

A genotyping-by-sequencing (GBS) approach was used to prepare the libraries and genotype genome-wide SNPs for the 34 R. rattus and 3 M. musculus samples. GBS was performed at the Cornell Genomic Diversity Facility following the protocol described by [21,22,23]. After removing adaptors and filtering the FASTQ file, the R. rattus sequences were aligned with the reference genome for SNP calling in TASSEL-GBS pipeline version 5 [24]. A minimum base call of five reads and a maximum locus missing data of 20% were defined. We used the whole genome of Rattus norvegicus (GCA_000317375) as a reference genome for SNP calling. After SNP calling, the metadata were exported in a VCF file format for downstream analyses using TASSEL 5 Standalone version [24]. Seven samples were omitted due to poor typing quality. Thus, a resulting vcf file comprising 27 R. rattus and 3 M. m. castaneus samples was processed for further analysis.

2.3. Phylogenomic Tree Construction

Initially, we used a distance matrix to build a simple neighbour-joining phylogenetic tree within the TASSEL program. To further revisit the phylogenetic relationship, we constructed a maximum likelihood (ML) tree using RAxML v8.2.X [25] and further modified the tree in Figtree 1.4.2 [26].

2.4. Genetic Subdivision

The genetic subdivision of R. rattus in India was explored by Bayesian clustering using the program STRUCTURE 2.3.2 [27,28]. The number of clusters (K) was identified by running the program for 1 million MCMC chains, following a burn-in of 1 million clusters. A batch run approach was used with 5 iterations for each cluster ranging from K = 1 to 4 incorporating admixture with the allele frequency-correlated model. The most probable number of clusters was selected using the second-order rate of change of log probability of the data between successive values of K (∆K) [29] in the program STRUCTURE HARVESTER v0.6.8 [30]. The whole analysis was repeated three times to confirm the accuracy of the analysis. ARLEQUIN 3.5.1 [31] was used to further examine genetic differentiation using FST.

3. Results

3.1. Genomic Variation and Diversity Indices

The R. rattus samples, originating from 16 locations from the west coast, east coast, and central India, were genotyped to yield 1499 genome-wide SNPs. The average minor allele frequencies were calculated as 0.1750, and the heterozygous proportion was found to be 0.2944.

3.2. Phylogenomics and Genetic Subdivision in R. rattus

With strong bootstrap support, the genome-wide SNPs resulted in an ML tree, which showed the splitting of R. rattus samples into two distinct clusters, corresponding to the east and west coasts of India (Figure 1A,B). Moreover, the R. rattus sampled from central India were grouped with west coast samples. In addition to this, the analysis in STRUCTURE 2.3.2 reconfirms the clustering pattern revealed by the phylogenomic tree (Figure 1C). The highest value of ∆K for K = 2, detected by Evanno’s method, reveals the partitioning of the contemporary R. rattus population in India into two lineages (Figure S1 and Table S2). Additionally, the average genetic differentiation (Fst) calculated for these two clusters was 0.204, indicating strong disparity between the spatially separated R. rattus populations. Of note, genetic signals of admixture were also detected in five R. rattus samples obtained from central India and the east coast (Figure 1C).

4. Discussion

In recent years, rodents, including house mice, M. musculus, and black rats, R. rattus, have been studied extensively as a proxy to understand ancient human settlements and movements. Considering the lack of information and importance of understanding the genomic diversity of R. rattus within India, we used the genotyping-by-sequencing (GBS) approach to study samples obtained from the west and east coasts and central India to understand their genetic relationship. Both our phylogenomic tree and clustering analysis show distinct subdivisions in the R. rattus population, corresponding to the west and east coasts. This genomic differentiation in the R. rattus population within Peninsular India corroborates the substantial mitochondrial diversity and differentiation that we have previously reported in the R. rattus population in India [8]. Both in our previous mitochondrial DNA study and in the current genomic analysis, we found clear divergence between the west coast and the east coast. Although we have not been able to perform accurate dating for either of our analyses, the close correspondence between the genomic and mitochondrial data suggests substantial differentiation between the west coast and east coast populations that would have occurred at a time well prior to the human settlements and the origin of commensalism. This is substantiated by the mitochondrial DNA tree (Figure 2 in reference [8]), which shows a level of differentiation between the west coast and east coast populations of a similar magnitude as between major phylogenetic lineages within R. rattus previously described in [4]. In Peninsular India, the west and east regions fall under two different ecological zones [32] that might have shaped the genome of R. rattus differently.
Considering the first association between R. rattus and humans, the advent of agriculture followed by a village/urban way of life in the northwest of the Indian subcontinent is the most plausible explanation for the origin of commensalism in the wild-living R. rattus population, and, once commensal, humans unintentionally moved them around [10,33]. Mehrgarh, a 7 kya site, was the epicentre of the origin of Indus valley civilization in the northwest of the subcontinent [9,34]. Most probably, it was in this region where commensalism first originated in R. rattus and further spread to other parts of the sub-continent. Many sites, ranging from early (~7–6.5 kya) to late (~5–2.8 kya) Indus, have been excavated in India [35]. Importantly, ten of the samples genotyped in this study originated from the northwest (Table S1: Rr1, Rr2, Rr3, Rr4, Rr5, Rr6, Rr9, Rr12, Rr13, Rr14). At many of the Neolithic Indus valley sites in the northwest, commensal R. rattus bones have been found and recorded [36,37]. In terms of early human settlements, compared to the northwest, the eastern and central Indian regions are less explored. However, a few studies suggest the origin of agriculture in other regions of India, including the Deccan Plateau in central India, as well as in the east (i.e., parts of present-day Bihar and Orissa), to 4.8–3.5 kya [32,38]. In addition, the presence of Neolithic practices of Austro-Asiatic linguistic tribes, such as Munda, Korku, and Santhal, in central and Eastern India is well documented [39].
Of the two lineages of R. rattus found on the east and west coasts, the basal positions of east coast R. rattus populations in both mtDNA [8] and genomic phylogenies indicate them to be ancestral. For mtDNA, there is evidence of greater genetic diversity in the east coast population, and the west coast population may be derived from the east coast population by having spread westwards well before R. rattus became commensal [8]. The signatures of the genetic admixture found in the R. rattus samples of central India and the east coast could be the outcome of ancient hybridization between the lineages, or it might have resulted from the intermixing of two lineages in this region, as a result of human transport. Considering the intermediate position of central India between the east and west coasts, admixture in this region is to be expected, although the R. rattus from central India predominantly have a west coast genetic signature. The contribution of the east coast rat genome to the genome in these central India individuals is very small. Given this characteristic of the central India R. rattus, and the relatively close proximity of the central India and some east coast populations, it is not surprising that there is admixture in both.
The genomic studies reported here confirm the coherence of the R. rattus population on the west coast of India already suggested by mtDNA analysis [8]. It is here that R. rattus became commensal with humans, and from here that R. rattus spread globally [4,8]. The mtDNA diversity for globally distributed R. rattus is low, as would be expected when largely derived from this one source area (the west coast of India) [4,8]. Interestingly, that is also true within that source area, also consistent with the human-mediated spread of commensal R. rattus throughout the west coast from a site of initial commensalism in the Indus valley. This could be the reason for the low mtDNA diversity on the west coast, rather than ancient derivation from the east coast populations. The present study indicates that the globally distributed R. rattus is not only derived from one particular geographically localised mtDNA lineage [8], it is derived from a population on the west coast of India that has distinct nuclear genome characteristics. This is evident because we have sampled all the way along the west coast of India, and these rats show genomic coherence, but are consistently different from those along the east coast (Figure 1). It is striking that the R. rattus on the east coast of India are also commensal, but, as we have said, they are genomically different from those on the west coast, and more genetically diverse in mtDNA (see Figure 5 in reference [8]). Similarly, there is not the same indication of a rapid spread from a single commensal source area. More sampling and genomic analyses are needed to fully understand this difference between the east and west coast. It would also be informative to include studies of ancient DNA from Indus and other Neolithic sites in the subcontinent to incorporate a time dimension into the work. India, as the origin of the spread of R. rattus, has much to offer to our understanding of the phylogeography of this species. It should be noted that this phylogeography and, in particular, knowledge of the source area is key information when trying to elucidate the invasiveness and other properties of the introduced populations of black rats, e.g., as vectors of human disease [40]. Taking this more broadly, it is important to emphasize that genetic and genomic studies of the black rat can help our understanding of various aspects of this most destructive, invasive mammal, and its management, as illustrated by recent studies defining eradication units [41,42,43].

5. Conclusions

India is crucial as the source area of the global spread of the black rat, R. rattus. Here we conduct the first population genomic study of Indian R. rattus, and demonstrate a similar genetic subdivision into east and west coast populations as found with mt DNA. Thus, the R. rattus worldwide range expansion originated from a west coast Indian population with defined nuclear genome characteristics as well as its well-known mt DNA characteristics. This understanding is valuable not only for future genome-based phylogeographic studies but also with respect genetic approaches to managing R. rattus as a devastating pest worldwide.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13020267/s1: Figure S1: a STRUCTURE HARVESTER output illustrating the values of delta K for the most probable number of clusters; Table S1: details of the sampling locations across the west coast, east coast and central India; Table S2: estimates of the most probable number of genetically distinct groups (K) indicated by STRUCTURE for the R. rattus population using the method of Evanno et al., (2005); VCF File: Dataset_R. rattus.

Author Contributions

M.B., A.A. collected the samples, S.F., M.B. and A.A. performed GBS and population genetic analyses, M.B. and J.B.S. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by DBT-CREST award (BT/IN/DBT CREST Awards/18/MMB/2012–2013) awarded to Mumtaz Baig for his postdoctoral research.

Institutional Review Board Statement

This study was conducted using already collected samples published previously and did not require any ethical approval or permits.

Informed Consent Statement

Not applicable.

Data Availability Statement

The VCF file is provided as Electronic Supplementary Material.

Acknowledgments

Sincere thanks to Nicky Xavier, Pradipta Sarangi, Poonam Thakur for sampling, to Soraia Barbosa for advice in library preparation, and to Avinash Karn for assistance in TASSEL analysis. S.F. and M.B. are grateful to Elizabeth Boulding for computational facilities. Thank you also to the reviewers for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Doherty, T.S.; Glen, A.S.; Nimmo, D.G.; Ritchie, E.G.; Dickman, C.R. Invasive predators and global biodiversity loss. Proc. Natl. Acad. Sci. USA 2016, 113, 11261–11265. [Google Scholar] [CrossRef] [Green Version]
  2. Russell, J.C.; Blackburn, T.M. Invasive alien species: Denialism, disagreement, definitions, and dialogue. Trends Ecol. Evol. 2017, 32, 312–314. [Google Scholar] [CrossRef] [Green Version]
  3. Meerburg, B.G.; Singleton, G.R.; Kijlstra, A. Rodent-borne diseases and their risks for public health. Crit. Rev. Microbiol. 2009, 35, 221–270. [Google Scholar] [CrossRef]
  4. Aplin, K.P.; Suzuki, H.; Chinen, A.A.; Chesser, R.T.; Have, J.T.; Donnellan, S.C.; Austin, J.; Frost, A.; Gonzalez, J.P.; Herbreteau, V.; et al. Multiple geographic origins of commensalism and complex dispersal history of black rats. PLoS ONE 2011, 6, e26357. [Google Scholar] [CrossRef] [Green Version]
  5. Lack, J.B.; Greene, D.U.; Conroy, C.J.; Hamilton, M.J.; Braun, J.K.; Mares, M.A.; van Den Bussche, R.A. Invasion facilitates hybridization with introgression in the Rattus rattus species complex. Mol. Ecol. 2012, 21, 3545–3561. [Google Scholar] [CrossRef]
  6. Niethammer, J. Zur Taxonomie und Ausbreitungsgeschichte der Hausratte (Rattus rattus). Zool. Anz. 1975, 194, 405–415. [Google Scholar]
  7. Dobney, K.; Ervynck, A. Sedentism or urbanism? On the origin of the commensal black rat (Rattus rattus). In Bones and the Man. Studies in Honour of Don Brothwell; O’Connor, T., Ed.; Oxbow: London, UK, 2002; pp. 95–109. [Google Scholar]
  8. Baig, M.; Khan, S.; Eager, H.; Atkulwar, A.; Searle, J.B. Phylogeography of the black rat Rattus rattus in India and the implications for its dispersal history in Eurasia. Biol. Invasions 2019, 21, 417–433. [Google Scholar] [CrossRef]
  9. Allchin, B.; Allchin, F.R. Origins of a Civilization: The Prehistory and Early Archaeology of South Asia; Viking: New Delhi, India, 1997; p. 287. [Google Scholar]
  10. Petraglia, M.D.; Allchin, B. (Eds.) The Evolution and History of Human Populations in South Asia: Inter-Disciplinary Studies in Archaeology, Biological Anthropology, Linguistics and Genetics; Springer: Dordrecht, The Netherlands, 2007. [Google Scholar]
  11. Searle, J.B.; Jones, C.S.; Gündüz, İ.; Scascitelli, M.; Jones, E.P.; Herman, J.S.; Rambau, R.V.; Noble, L.R.; Berry, R.J.; Giménez, M.D.; et al. Of mice and (Viking?) men: Phylogeography of British and Irish house mice. Proc. R. Soc. B Biol. Sci. 2009, 276, 201–207. [Google Scholar] [CrossRef] [Green Version]
  12. Jones, E.P.; Kooij, J.V.; Solheim, R.; Searle, J.B. Norwegian house mice (Mus musculus musculus/domesticus): Distributions, routes of colonization and patterns of hybridization. Mol. Ecol. 2010, 23, 5252–5264. [Google Scholar] [CrossRef]
  13. Jones, E.P.; Jensen, J.; Magnussen, E.; Gregersen, N.; Hansen, H.S.; Searle, J.B. A molecular characterization of the charismatic Faroe house mouse. Biol. J. Linn. Soc. 2011, 101, 471–482. [Google Scholar] [CrossRef] [Green Version]
  14. Gabriel, S.I.; Stevens, M.I.; Mathias, M.L.; Searle, J.B. Of Mice and ‘Convicts’: Origin of the Australian house mouse, Mus musculus. PLoS ONE 2011, 6, e28622. [Google Scholar] [CrossRef] [Green Version]
  15. Aplin, K.P.; Brown, P.R.; Jacob, J.; Krebs, C.; Singleton, G.R. Field Methods for Rodent Studies in Asia and the Indo-Pacific; Australian Centre for International Agricultural Research: Canberra, Australia, 2003.
  16. Colangelo, P.; Abiadh, A.; Aloise, G.; Capizzi, D.; Vasa, E.; Annesi, F.; Castiglia, R. Mitochondrial phylogeography of the black rat supports a single invasion of the western Mediterranean basin. Biol. Invasions 2015, 17, 1859–1868. [Google Scholar] [CrossRef]
  17. Russell, J.C.; Gleeson, D.M.; le Corre, M. The origin of Rattus rattus on the Iles Eparses, Western Indian Ocean. J. Biogeogr. 2011, 38, 1834–1836. [Google Scholar] [CrossRef]
  18. Raman, R.; Sharma, T. Similarity in karyotypes of Rattus rattus with 38 chromosomes from India and other parts of the world. Experientia 1972, 28, 1375–1377. [Google Scholar] [CrossRef]
  19. Thapar, R. The History of Early India from the Origins to AD 1300; Penguin: New Delhi, India, 2015. [Google Scholar]
  20. Puckett, E.E.; Munshi-South, J. Brown rat demography reveals pre-commensal structure in eastern Asia before expansion into Southeast Asia. Genome Res. 2019, 29, 762–770. [Google Scholar] [CrossRef] [Green Version]
  21. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A robust, simple Genotyping-by-Sequencing (GBS) approach for high diversity species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef] [Green Version]
  22. White, T.A.; Perkins, S.E.; Heckel, G.; Searle, J.B. Adaptive evolution during an ongoing range expansion: The invasive bank vole (Myodes glareolus) in Ireland. Mol. Ecol. 2013, 22, 2971–2985. [Google Scholar] [CrossRef]
  23. Glaubitz, J.C.; Casstevens, T.M.; Lu, F.; Harriman, J.; Elshire, R.J.; Sun, Q.; Buckler, E.S. TASSEL-GBS: A high capacity Genotyping by Sequencing analysis pipeline. PLoS ONE 2014, 9, e90346. [Google Scholar] [CrossRef]
  24. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  25. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  26. Rambaut, A. FigTree Software. 2018. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 7 July 2020).
  27. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef] [PubMed]
  28. Falush, D.; Stephens, M.; Pritchard, J.K. Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics 2003, 164, 1567–1587. [Google Scholar] [CrossRef]
  29. Evanno, G.; Regnault, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE. A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Earl, D.A.; von Holdt, B.M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  31. Excoffier, L.; Lischer, H.E. Arlequin (version 3.5.1.2): An integrated software package for population genetics data analysis. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef] [PubMed]
  32. Fuller, D.Q. Non-human genetics, agricultural origins and historical linguistics in South Asia. In The Evolution and History of Human Populations in South Asia: Inter-Disciplinary Studies in Archaeology, Biological Anthropology, Linguistics and Genetics; Petraglia, M.D., Allchin, B., Eds.; Springer: Dordrecht, The Netherlands, 2007; pp. 393–443. [Google Scholar]
  33. Kovács, Z.S. Dispersal history of an invasive rodent in Hungary—Subfossil finds of Rattus rattus. Acta Zool. Acad. Sci. Hun. 2012, 58, 379–394. [Google Scholar]
  34. Jarrige, J.F. The early architectural traditions of Greater Indus as seen from Mehrgarh, Baluchistan. Stud. Hist. Art. 1993, 31, 25–33. [Google Scholar]
  35. Rao, L.S.; Sahu, N.B.; Sahu, P.; Shastry, U.A.; Diwan, S. New light on the excavation of Harappan settlement at Bhirrana. Puratattva 2005, 35, 67–75. [Google Scholar]
  36. Kane, V.S. Animal remains from Rojdi. In Harappan Civilization and Rojdi; Possehl, G.L., Raval, M.R., Eds.; Oxford-IBH: New Delhi, India, 1989; pp. 182–184. [Google Scholar]
  37. Thomas, P.K.; Joglekar, P.P.; Matsushima, Y.; Pawankar, S.J.; Deshpande, A. Subsistence based on animals in the Harappan culture of Gujarat, India. Anthropozoologica 1997, 25, 767–776. [Google Scholar]
  38. Mitre, V. Origin and history of agriculture in the Indian sub-continent. J. Hum. Evol. 1977, 7, 31–36. [Google Scholar] [CrossRef]
  39. Bellwood, P. The origins and spread of agriculture in the Indo-Pacific region: Gradualism, diffusion or revolution and colonization. In The Origins and Spread of Agriculture and Pastoralism in Eurasia; Harris, D.R., Ed.; UCL Press: London, UK, 1996; pp. 465–498. [Google Scholar]
  40. Searle, J.B. The genetics of mammalian invasions: A review. Wildl. Res. 2008, 35, 185–192. [Google Scholar] [CrossRef]
  41. Iannucci, A.; Natali, C.; Ciofi, C.; Sposimo, P.; Fratini, S. Population genetic analysis of invasive black rats: Defining eradication units in the Tuscan Archipelago National Park. Hystrix. Ital. J. Mamm. 2018, 29, 67–73. [Google Scholar] [CrossRef]
  42. Fratini, S.; Natali, C.; Zanet, S.; Iannucci, A.; Capizzi, D.; Sinibaldi, I.; Sposimo, P.; Ciofi, C. Assessment of rodenticide resistance, eradication units, and pathogen prevalence in black rat populations from a Mediterranean biodiversity hotspot (Pontine Archipelago). Biol. Invasions 2020, 22, 1379–1395. [Google Scholar] [CrossRef]
  43. Sjodin, B.M.F.; Irvine, R.L.; Ford, A.T.; Howald, G.R.; Russello, M.A. Rattus population genomics across the Haida Gwaii archipelago provides a framework for guiding invasive species management. Evol. Appl. 2020, 13, 889–904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (A) Sampling locations of R. rattus across India (yellow circles depict west coast and central India locations, while pink circles indicate east coast locations; stars indicate locations highlighted in (C)). (B) SNP-based maximum likelihood phylogenomic tree constructed using 27 R. rattus samples. Clade K1 includes samples from the west coast and central India while clade K2 includes samples from the east coast. Similarly, clade K3 comprises 3 M. m. castaneus samples used as the outgroup. (C) STRUCTURE plot illustrating the population structure in R. rattus, sampled from the west coast, east coast and central India.
Figure 1. (A) Sampling locations of R. rattus across India (yellow circles depict west coast and central India locations, while pink circles indicate east coast locations; stars indicate locations highlighted in (C)). (B) SNP-based maximum likelihood phylogenomic tree constructed using 27 R. rattus samples. Clade K1 includes samples from the west coast and central India while clade K2 includes samples from the east coast. Similarly, clade K3 comprises 3 M. m. castaneus samples used as the outgroup. (C) STRUCTURE plot illustrating the population structure in R. rattus, sampled from the west coast, east coast and central India.
Genes 13 00267 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Baig, M.; Farah, S.; Atkulwar, A.; Searle, J.B. Genomic Analysis Reveals Subdivision of Black Rats (Rattus rattus) in India, Origin of the Worldwide Species Spread. Genes 2022, 13, 267. https://doi.org/10.3390/genes13020267

AMA Style

Baig M, Farah S, Atkulwar A, Searle JB. Genomic Analysis Reveals Subdivision of Black Rats (Rattus rattus) in India, Origin of the Worldwide Species Spread. Genes. 2022; 13(2):267. https://doi.org/10.3390/genes13020267

Chicago/Turabian Style

Baig, Mumtaz, Sameera Farah, Ashwin Atkulwar, and Jeremy B. Searle. 2022. "Genomic Analysis Reveals Subdivision of Black Rats (Rattus rattus) in India, Origin of the Worldwide Species Spread" Genes 13, no. 2: 267. https://doi.org/10.3390/genes13020267

APA Style

Baig, M., Farah, S., Atkulwar, A., & Searle, J. B. (2022). Genomic Analysis Reveals Subdivision of Black Rats (Rattus rattus) in India, Origin of the Worldwide Species Spread. Genes, 13(2), 267. https://doi.org/10.3390/genes13020267

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop