- Open Access
Population genomics reveal rapid genetic differentiation in a recently invasive population of Rattus norvegicus
Frontiers in Zoology volume 18, Article number: 6 (2021)
Invasive species bring a serious effect on local biodiversity, ecosystems, and even human health and safety. Although the genetic signatures of historical range expansions have been explored in an array of species, the genetic consequences of contemporary range expansions have received little attention, especially in mammal species. In this study, we used whole-genome sequencing to explore the rapid genetic change and introduction history of a newly invasive brown rat (Rattus norvegicus) population which invaded Xinjiang Province, China in the late 1970s.
Bayesian clustering analysis, principal components analysis, and phylogenetic analysis all showed clear genetic differentiation between newly introduced and native rat populations. Reduced genetic diversity and high linkage disequilibrium suggested a severe population bottleneck in this colonization event. Results of TreeMix analyses revealed that the introduced rats were derived from an adjacent population in geographic region (Northwest China). Demographic analysis indicated that a severe bottleneck occurred in XJ population after the split off from the source population, and the divergence of XJ population might have started before the invasion of XJ. Moreover, we detected 42 protein-coding genes with allele frequency shifts throughout the genome for XJ rats and they were mainly associated with lipid metabolism and immunity, which could be seen as a prelude to future selection analyses in the novel environment of XJ.
This study presents the first genomic evidence on genetic differentiation which developed rapidly, and deepens the understanding of invasion history and evolutionary processes of this newly introduced rat population. This would add to our understanding of how invasive species become established and aid strategies aimed at the management of this notorious pest that have spread around the world with humans.
Invasive species are not only a major threat for native biodiversity (such as species decline or extinction) and ecosystems [1, 2], but also cause considerable annual damage to agriculture, property, human health and safety, and natural resources [3, 4]. It is vital to understand the dynamics of invasion processes, such as the rapid evolution and introduction history of these species, as well as their dispersal mechanisms [5,6,7,8,9,10]. Clarifying the likely spread of invasive pest is crucial for accurate risk assessment and for optimized management strategies . Genomic approaches can help in developing this understanding and promise to provide higher resolution than previous genetic studies to explore the population structure, demographic history, and molecular evolution of invasive populations [12,13,14,15].
The brown rat (Rattus norvegicus) is one of the most successful mammalian invaders due to its remarkable migration and adaptation abilities [16, 17]. The species is believed to have originated in either northern Asia [18, 19] or Southeast Asia [20, 21], and emerged ∼1.3 million years ago . Now, it has invaded and spread to nearly every major landmass except Antarctica [17, 22]. In China, the brown rat is also widespread except Tibet [16, 23, 24]. However, based on historical records, the brown rat was not observed in Xinjiang Province (XJ) until the late 1970s [25,26,27]. In the middle of the 1970s, brown rats were first detected on trains from Beijing to XJ, which was opened in the middle of the 1960s, and a few years later, they were first observed on land in the eastern XJ, the Turpan-Hami Basin [25, 27]. Thus, the introduction of XJ brown rats is generally attributed to Beijing-XJ railway transportation [25,26,27]. Despite its recent arrival, it has spread throughout XJ and there is now a large population in the region [26, 28, 29].
To further deepen the understanding of invasion history and evolutionary processes of this recently introduced rat population, we used next-generation sequencing data to explore the extent of whole-genome variation among the introduced XJ population and other native populations of brown rats in China. We explored the molecular phylogenetic relationships of the recently invasive brown rats within a broader phylogenetic framework, and tested whether introduced rats undergo rapid molecular differentiation in this new geographical range during such short period of invasion. We also investigated the source of origin and the demographic history of the introduced XJ population and identified genes with allele frequency shifts throughout the genome for XJ rats that might respond to local selection in this new geographical range.
In the present study, we sequenced and analyzed the whole genomes of 50 brown rats from one invasive region and other native regions across China to an average sequencing depth of ~ 16.5× (Fig. 1; Additional file 1: Table S1). After applying stringent quality control criteria, we identified a total of 11.3 million single nucleotide polymorphisms (SNPs) among all the individuals.
Population structure and genetic relationships among rat populations
To identify population structure and the genetic relationships of different rat populations, we performed a series of classical analyses including Bayesian clustering analysis by ADMIXTURE, principal components analysis (PCA), and phylogenetic assessments using whole-genome SNPs (Fig. 2).
Clustering analyses by ADMIXTURE showed that brown rats were separated into XJ and non-XJ populations when the number of presumed ancestral populations (K) was set to 2 which is the best supported number of clusters (Fig. 2a; Additional file 1: Table S2). PCA showed that the introduced XJ rats were widely divergent from the native ones. XJ individuals clustered together in a single area and the first principal component (PC1) clearly separated XJ from other populations (NW, NE, NC, CC, and SC), which are differentiated into separate clusters along PC2 (Fig. 2b). Phylogenetic reconstruction also identified six major clusters which is consistent with both the results of the clustering analysis and the PCA (Fig. 2c). All the population-level nodes in this phylogenetic tree had a bootstrap support of 100%.
Genetic differentiation and genetic diversity
The genome-wide analysis of FST divergence showed that FST among groups with XJ population were much higher than those with other populations (all P < 2.2E− 16, Wilcoxon rank-sum test; Table 1), which further indicates the striking genetic differentiation between XJ rats and the native ones. Nucleotide diversity (π) and linkage disequilibrium (LD) were used to assess the genetic diversity within rat populations. The π in the introduced XJ population was lower than that in the other rat populations (all P < 2.2E− 16, Wilcoxon rank-sum test; Table 2). LD patterns showed similar trends. By calculating the pairwise LD between polymorphic sites for all regions in each population, we found that LD decayed much slower in the XJ population than in other populations (Fig. 3).
Demography and admixture of XJ rats
We inferred the demographic history of the XJ population using 2D (two population) models in δaδi . Divergence with continuous symmetric migration model yielded the best log likelihood and AIC statistic (Additional file 1: Table S3) between XJ and the northern (including NW, NC, and NE) rats. The size of the populations after the split were 4420 and 9523 for XJ and the northern rats, respectively (Fig. 4a), indicating a bottleneck of XJ population after the divergence. The divergence time between them was 2708 years (Fig. 4a), indicating that the divergence of XJ population might have started before the invasion of XJ. The estimated migration rate was quite low as 4.1 × 10− 5 per generation (Fig. 4a). We further tested for evidence of migration and admixture between populations using the maximum likelihood method implemented in Treemix . The population tree without any migration showed that the XJ clade was most closely related to the native NW population, and these two groups formed part of a large clade of populations located in northern China, which also included the NC and NE populations (Fig. 4b), providing further evidence that the XJ rats were derived from northern China, more specifically, northwest China. There was no evidence for migration involving XJ when migration tracks were allowed in TreeMix (Additional file 2: Fig. S1). Furthermore, the F3 test showed no evidence for introgression (Additional file 2: Fig. S2). These results from Treemix were largely consistent with results from δaδi that indicated a very low level of introgression between XJ and other populations.
F ST outliers analysis
We performed a genome-wide FST scan to identify genes with allele-frequency shifts in XJ rats, as a prelude to future selection analyses in the novel environment of XJ. We identified 42 regions with high Z (FST) between XJ and other native rats (Fig. 5). In these outlier regions, we detected 42 protein-coding genes (Additional file 1: Table S4) which were mainly associated with two different functions: lipid metabolism and immunity (Additional file 1: Table S5). These genes were significantly enriched in Gene Ontology (GO) terms involved in arachidonic acid secretion (GO:0050482), arachidonate transport (GO:1903963), fatty acid transport (GO:0015908), lipid catabolic process (GO:0016042), lipid transport (GO:0006869), lipase activity (GO:0016298), negative regulation of leukocyte cell-cell adhesion (GO:1903038), mannose binding (GO:0005537), monosaccharide binding (GO:0048029), etc. They also showed an overrepresentation in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to linoleic acid metabolism (rno00591), fat digestion and absorption (rno04975), arachidonic acid metabolism (rno00590), C-type lectin receptor signaling pathway (rno04625), measles (rno05162), phagosome (rno04145), etc.
In mammals, analyses of genetic differentiation and molecular evolution among species or populations of the same species by genomic approaches have mostly focused on geologic time scales (more than tens of thousands of years) [24, 32,33,34]. However, rapid genetic change often occurs over ecological time scales (e.g., tens of generations or fewer) [35, 36], and multiple studies on non-mammal species have shown that evolutionary changes can occur within dozens of generations [37,38,39,40,41,42]. Our current work presents the first genomic results on rapid genetic differentiation of the recently introduced invasive brown rat. By analyzing whole-genome sequences of 50 individuals from newly introduced and native populations, we determined the phylogenetic placement of the invasive population. Our results demonstrated that the introduced XJ rats have become a recently diverged population over the past decades of invasion. All the introduced rats formed a single well-supported clade in the phylogenetic tree and showed strong genome-wide divergence in genetic structure from native ones in PCA and ADMIXTURE analyses. Such striking genetic differentiation between frontier and native rat populations implies that the invasive population of the brown rat has undergone a rapid genetic change, indicating that genetic differentiation can develop rapidly in brown rats.
With the low nucleotide diversity (π), the effective population of XJ might undergo a more recent reduction than other native ones . This result is consistent with the high linkage disequilibrium (LD), suggesting that the introduced XJ rat might come from a small size of founder population and this colonization event was associated with a severe bottleneck . Demographic analysis revealed that XJ population isolated from the source population ~ 2700 years ago that is much earlier than the invasion of XJ, indicating the divergence of XJ population started before the invasion of XJ. Future studies should broaden the geographical scale in northern China (especially northwest China)  to perform a much more accurate and comprehensive evaluation of the source population of XJ rats. Invasion and post-establishment expansions are often associated with recurrent founder effects and bottlenecks, which points to low genetic diversity and the accumulation of deleterious alleles or mutations [6, 46,47,48,49]. However, invasive populations with reduced genetic diversity may also successfully colonize new environments and expand across wide geographic area [10, 50,51,52]. It is noted that bottlenecks may not necessarily be negative, but help purge deleterious alleles due to genetic drift and non-random mating in small populations . Nevertheless, a loss of genetic variation may limit mutations available for natural selection and is expected to impact negatively on the adaptive capacity of populations [12, 47]. The range expansion of brown rats was a response to relatively recent increases in global trade , although the brown rat has been widespread throughout the world. Additionally, a severe bottleneck occurred about 20 kya in the brown rat . Thus the brown rat may possess the capability of rapid expansion with low effective population size.
Previous studies revealed that geography and environmental heterogeneity shaped genetic structure of brown rats [17, 54]. We also demonstrated that R. norvegicus was differentiated into clades (NW, NC, NE, CC, and SC) corresponding mostly to geographic locations in China. Tree topology inferred from TreeMix showed that the introduced XJ population was closely related to the native NW population, indicating that the introduced rats were derived from NW, a geographic region neighboring to XJ. In this study, our sampling range covered main distribution areas of all the four subspecies of brown rats in China [16, 23]. Tibet province that is large plateau area neighboring to XJ has no distribution of the brown rat [16, 23]. Although some other samples were also provided in the research of Zeng et al. , most of these sample areas were located in South China, and brown rats from these areas (including our sampling area of SC) cluster together and significantly divergent from North ones . Given the large genetic distance of SC and XJ population in our current work, we believe that these unsampled areas in South China could not be the source of XJ population.
In a population with such a recent bottleneck and expansion, disentangling positive selection from neutral forces is challenging. However, we identified genes throughout the genome with significant allele frequency shifts in XJ rats, as a prelude to future selection analyses in the novel environment of XJ; A total of 42 genes were detected in FST outlier regions and they were mainly associated with lipid metabolism and immunity. One unique climatic characteristics of XJ is the dramatic diurnal and seasonal temperature differences . The variation of lipid metabolism genes may regulate the energy metabolism of brown rats , consequently, fitting in with drastic changes in temperature. Rapid evolution in immune genes is a well-known example and presumably occurs because new mutations help organisms to prevail in evolutionary “arms races” with pathogens. Previous studies showed that immune-related genes were under positive selection in population divergence of the three-spined stickleback (Gasterosteus aculeatus)  or the Olympia oyster (Ostrea lurida) , indicating that how to deal with new types of pathogens might be one of the critical issues during range expansion. As the source of a variety of pathogens, immunity plays a great role in keeping health and adapting to new habitats for brown rats. Zeng et al.  recently demonstrated that genes related to immune system have evolved rapidly under positive selection in wild brown rats during their dispersal, indicating that resistance to external pathogens is a key issue for the brown rat in response to novel environmental forces.
This work highlights that rapid genetic change can occur after such a short period of invasion in brown rats, and the successful invasion could come from a small size of founder population. Geography and environmental heterogeneity shaped genetic structure of brown rats, and the adjacent NW population was the source of XJ population. Invasive species can rapidly adapt to the changed conditions by means of genetic changes through the process of evolution [8, 59, 60]. Identifying the underlying genetic causes of invasion success is a key component of research on biological invasion [6, 61, 62]. This study presents the first genomic evidence on genetic differentiation which developed rapidly, and deepens the understanding of invasion history and evolutionary processes of this recently introduced rat population, which would hold clues to their successful invasion and could aid strategies aimed at the management of this notorious pest that have spread around the world with humans.
Samples, study sites, and genome sequencing
A total of 50 adult brown rats were sampled across China including Xinjiang and other 14 provinces, which cover the main distribution area of brown rats in China [16, 23] (Fig. 1; Additional file 1: Table S1). Rats were trapped with snap trap and the nearest trapping locations among collected rats were 5 km apart to avoid sampling closely related individuals. Tails were obtained and stored in ethanol for DNA extraction. The species status of brown rats was confirmed via morphology and mitochondrial cytochrome oxidase subunit I barcode sequences by Sanger sequencing .
Whole-genomic DNA of 38 individuals was extracted from small pieces of tail tissue using a TailGen DNA extraction kit (CWBIO, Beijing, China). The quality and integrity of the extracted DNA was checked by measuring the A260/A280 ratio using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA, USA) and by agarose gel electrophoresis. Library preparation, Illumina sequencing, read alignments, and variant calling were performed by Annoroad Gene Technology (Beijing) Co. Ltd. For each individual, 100 ng genomic DNA was used to construct PCR-based libraries with a 350 base pair (bp) insertion size and sequenced on an Illumina HiSeq X Ten instrument with 150 bp paired-end reads. Our previously generated 12 whole genome sequences of wild-caught brown rat individuals  were reanalyzed in this study. After filtering out raw sequencing reads containing adapters and reads of low quality, the remaining clean reads were mapped to the reference R. norvegicus genome RGSC5.0  using BWA v0.7.12 with default parameters . SAMtools v1.2  was performed to sort reads, and MarkDuplicates in Picard tools v1.13 (http://broadinstitute.github.io/picard/) was used to remove PCR duplicates. Reads mapped to two or more locations were filtered out.
SNP calling and filtering
The Genome Analysis Toolkit (GATK)  HaplotypeCaller protocol was used for SNP calling via local re-assembly of haplotypes for the populations. SNPs were further filtered by applying the following cutoffs with the GATK VariantFiltration protocol: QD < 10.0, FS > 10.0, DP < 4.0, QUAL < 30.0, ReadPosRankSum < − 8.0. Additionally, the sites with a minor allele frequency (MAF) < 0.05 or including more than 10% missing genotypes were filtered out. Then the filtered high-quality SNPs were kept for subsequent analysis.
Population structure analyses
We used ADMIXTURE v1.3.0  to investigate the population structure, with the number of coancestry clusters ranging from 2 to 6. To assess the best value of K, we performed 10-fold cross-validation and determined the K value with the lowest cross-validation error. PCA was performed using EIGENSOFT v6.0.1 . To mitigate the effects of linkage disequilibrium (LD) on genetic structure, we pruned the markers using the “-indep-pairwise 50 5 0.05” option of PLINK . To investigate the relationships within introduced and native populations, a phylogenetic tree from whole-genome SNP data was constructed using SNPhylo . The program was run with 100 bootstrap repetitions, and the genome information of R. rattus  was used as an outgroup.
Genetic differentiation, genetic diversity, and LD analyses
Pairwise genetic differentiation and nucleotide diversity among rat populations were calculated using VCFtools v0.1.16  by means of FST (−weir-fst) and π (−site-pi) with a 100 kb sliding window, respectively. Linkage disequilibrium (LD) decay was calculated with PopLDdecay  with the following parameters: -MaxDist 300 -MAF 0.05 -Miss 0.9.
To investigate alternative divergence scenarios for XJ population, we used the diffusion approximation method of δaδi to analyze two-dimensional joint site frequency spectra (2D-JSFS) . We used the demographic modelling pipeline (dadi_pipeline) of Portik et al.  to conduct all analyses. For all models, we performed consecutive rounds of optimizations following Portik et al. . For each round, we ran multiple replicates and used parameter estimates from the best scoring replicate (highest log-likelihood) to seed searches in the following round. We used the default settings in dadi_pipeline for each round (replicates = 10, 20, 30, 40; maxiter = 3, 5, 10, 15; fold = 3, 2, 2, 1), and optimized parameters using the Nelder-Mead method (optimize_log_fmin). Akaike Information Criteria (AIC) values were used to compare demography models, and the demography model with the lowest AIC was chosen as the best-fitting model. We used the approach described in Mattingsdal et al.  and Choi et al.  to convert the parameter estimates into meaningful biological values using a generation time of 0.5 year and a mutation rate of 2.96 × 10− 9 [24, 32].
Population admixture analysis
Population tree topology was estimated using the maximum likelihood method implemented in TreeMix . TreeMix models the genetic drift at genome-wide polymorphisms to infer relationships between populations. It first estimates a dendrogram of the relationships between sampled populations. Next it compares the covariance structure modeled by this dendrogram to the observed covariance between populations. When populations are more closely related than modeled by a bifurcating tree it suggests that there has been admixture in the history of those populations. TreeMix then adds an edge to the phylogeny, now making it a phylogenetic network. The position and direction of these edges are informative; if an edge originates more basally in the phylogenetic network it indicates that this admixture occurred earlier in time or from a more diverged population. We first inferred the maximum likelihood (ML) tree with the command “-i input -o output.” We then tested trees for one and two migration events (Additional file 2: Fig. S2). The threepop module (F3 test) from the TreeMix package was used to validate the migration events.
F ST outliers analysis
We estimated the FST (XJ/others) values for each window using VCFtools v0.1.16 (Danecek et al., 2011), with a window size of 100 kb and a step size of 50 kb. Next, we Z-transformed FST values and then identified regions with Z (FST) that were greater than 5 standard deviations from the mean [40, 76, 77]. To characterize the molecular functions of the genes contained in these outlier regions, we performed functional enrichment analyses using the clusterProfiler toolkit , where the significance level was set at 0.05 and the P-value was corrected using the Benjamini-Hochberg false discovery rate (FDR).
Availability of data and materials
The genome sequencing data are available on the Genome Sequence Archive database (http://gsa.big.ac.cn/) under accession numbers CRA001635.
Single nucleotide polymorphism
Principal components analysis
Gurevitch J, Padilla DK. Are invasive species a major cause of extinctions? Trends Ecol Evol. 2004;19:470–4.
Clavero M, García-Berthou E. Invasive species are a leading cause of animal extinctions. Trends Ecol Evol. 2005;20:110.
Witmer G. The ecology of vertebrate pests and integrated pest management (IPM). In: Kogan M, Jepson P, editors. Perspectives in ecological theory and integrated pest management. Cambs: Cambridge University Press; 2007. p. 393–410.
Chown SL, Hodgins KA, Griffin PC, Oakeshott JG, Byrne M, Hoffmann AA. Biological invasions, climate change and genomics. Evol Appl. 2015;8:23–46.
Mooney HA, Cleland EE. The evolutionary impact of invasive species. Proc Natl Acad Sci U S A. 2001;98:5446–51.
Sakai AK, Allendorf FW, Holt JS, Lodge DM, Molofsky J, With KA, Baughman S, Cabin RJ, Cohen JE, Ellstrand NC, McCauley DE, O’Neil P, Parker IM, Thompson JN, Weller SG. The population biology of invasive species. Annu Rev Ecol Syst. 2001;32:305–32.
Lee CE. Evolutionary genetics of invasive species. Trends Ecol Evol. 2002;17:386–91.
Whitney KD, Gabler CA. Rapid evolution in introduced species, ‘invasive traits’ and recipient communities: challenges for predicting invasive potential. Divers Distrib. 2008;14:569–80.
Medley KA, Jenkins DG, Hoffman EA. Human-aided and natural dispersal drive gene flow across the range of an invasive mosquito. Mol Ecol. 2015;24:284–95.
Cao LJ, Wei SJ, Hoffmann AA, Wen JB, Chen M. Rapid genetic structuring of populations of the invasive fall webworm in relation to spatial expansion and control campaigns. Divers Distrib. 2016;22:1276–87.
Stepien CA, Brown JE, Neilson ME, Tumeo MA. Genetic diversity of invasive species in the great lakes versus their eurasian source populations: insights for risk analysis. Risk Anal. 2005;25:1043–60.
Puzey J, Vallejo-Marin M. Genomics of invasion: diversity and selection in introduced populations of monkeyflowers (Mimulus guttatus). Mol Ecol. 2014;23:4472–85.
Osada N, Hettiarachchi N, Adeyemi Babarinde I, Saitou N, Blancher A. Whole-genome sequencing of six Mauritian Cynomolgus macaques (Macaca fascicularis) reveals a genome-wide pattern of polymorphisms under extreme population bottleneck. Genome Biol Evol. 2015;7:821–30.
Bourne SD, Hudson J, Holman LE, Rius M. Marine invasion genomics: revealing ecological and evolutionary consequences of biological invasions. In: Oleksiak MF, Rajora OP, editors. Population genomics: marine organisms, population genomics. Berlin: Springer International Publishing; 2018. p. 1–36.
Sharma V, Hecker N, Roscito JG, Foerster L, Langer BE, Hiller M. A genomics approach reveals insights into the importance of gene losses for mammalian adaptations. Nat Commun. 2018;9:1215.
Musser GG, Carleton MD. Superfamily Muroidea. In: Wilson DE, Reeder DM, editors. Mammal species of the world a taxonomic and geographic reference. Baltimore: Johns Hopkins University Press; 2005. p. 894–1531.
Puckett EE, Park J, Combs M, Blum MJ, Bryant JE, Caccone A, Costa F, Deinum EE, Esther A, Himsworth CG, Keightley PD, Ko A, Lundkvist A, McElhinney LM, Morand S, Robins J, Russell J, Strand TM, Suarez O, Yon L, Munshi-South J. Global population divergence and admixture of the brown rat (Rattus norvegicus). Proc Biol Sci. 2016;283:20161762.
Kosoy M, Khlyap L, Cosson JF, Morand S. Aboriginal and invasive rats of genus Rattus as hosts of infectious agents. Vector-Borne Zoonot. 2015;15:3–12.
Puckett EE, Munshi-South J. Brown rat demography reveals pre-commensal structure in eastern Asia before expansion into Southeast Asia. Genome Res. 2019;29:762–70.
Song Y, Lan Z, Kohn MH. Mitochondrial DNA phylogeography of the Norway rat. PLoS One. 2014;9:e88425.
Zeng L, Ming C, Li Y, Su LY, Su YH, Otecko NO, Dalecky A, Donnellan S, Aplin K, Liu XH, Song Y, Zhang ZB, Esmailizadeh A, Sohrabi SS, Nanaei HA, Liu HQ, Wang MS, Ag Atteynine S, Rocamora G, Brescia F, Morand S, Irwin DM, Peng MS, Yao YG, Li HP, Wu DD, Zhang YP. Out of southern East Asia of the brown rat revealed by large-scale genome sequencing. Mol Biol Evol. 2018;35:149–58.
Robinson R. Genetics of the Norway rat. Oxford: Pergamon Press; 1965.
Wang YX. A complete checklist of mammal species and subspecies in China: a taxonomic and geographic reference. Beijing: China Forestry Publishing House; 2003.
Teng H, Zhang Y, Shi C, Mao F, Cai W, Lu L, Zhao F, Sun Z, Zhang J. Population genomics reveals speciation and introgression between Brown Norway rats and their sibling species. Mol Biol Evol. 2017;34:2214–28.
Yu X, Zhang JT, Ye RY, Bao LB, Wu KQ. Notes on rodent composition and its changes both along the railway line and on the trains in Xinjiang--the discovery of Rattus norvegicus in this region. Zool Res. 1980;1:135–8.
Xie Y, Li ZY, William PG, Dianmo L. Invasive species in China—an overview. Biodivers Conserv. 2001;10:1317–41.
Anwar T, Zhang D. Investigation and prediction of the number of Rattus norvegicus population in Xinjiang arid area. J Xingjiang Normal Univ. 2002;21:45–9.
Zhang DM, Jiang T, Mahmut, Enwer, Liu L, Li W, Liao LF. Spreading again along with the railway of Rattus norvegicus and its regulatiag methods in Xinjiang. Endemic Dis Bull. 1997;12:106–8.
Tumur A, Abliz O. The community effect of Rattus norvegicus in Xinjiang arid zone. Life Sci Res. 2002;6:194–8.
Gutenkunst RN, Hernandez RD, Williamson SH, Bustamante CD. Inferring the joint demographic history of multiple populations from multidimensional SNP frequency data. PLoS Genet. 2009;5:e1000695.
Pickrell J, Pritchard J. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012;8:e1002967.
Deinum EE, Halligan DL, Ness RW, Zhang YH, Cong L, Zhang JX, Keightley PD. Recent evolution in Rattus norvegicus is shaped by declining effective population size. Mol Biol Evol. 2015;32:2547–58.
Wang GD, Zhai W, Yang HC, Wang L, Zhong L, Liu YH, Fan RX, Yin TT, Zhu CL, Poyarkov AD, Irwin DM, Hytönen MK, Lohi H, Wu CI, Savolainen P, Zhang YP. Out of southern East Asia: the natural history of domestic dogs across the world. Cell Res. 2016;26:21–33.
Liu Z, Zhang L, Yan Z, Ren Z, Han F, Tan X, Xiang Z, Dong F, Yang Z, Liu G, Wang Z, Zhang J, Que T, Tang C, Li Y, Wang S, Wu J, Li L, Huang C, Roos C, Li M. Genomic mechanisms of physiological and morphological adaptations of limestone langurs to karst habitats. Mol Biol Evol. 2019;37:952–68.
Carroll PS, Hendry AP, Reznick DN, Fox CW. Evolution on ecological time-scales. Funct Ecol. 2007;21:387–93.
Colautti RI, Lau JA. Contemporary evolution during invasion: evidence for differentiation, natural selection, and local adaptation. Mol Ecol. 2015;24:1999–2017.
Terekhanova NV, Logacheva MD, Penin AA, Neretina TV, Barmintseva AE, Bazykin GA, Kondrashov AS, Mugue NS. Fast evolution from precast bricks: genomics of young freshwater populations of threespine stickleback Gasterosteus aculeatus. PLoS Genet. 2014;10:e1004696.
Bernardi G, Azzurro E, Golani D, Miller MR. Genomic signatures of rapid adaptive evolution in the bluespotted cornetfish, a Mediterranean Lessepsian invade. Mol Ecol. 2016;25:3384–96.
Vera M, Díez-del-Molino D, García-Marín JL. Genomic survey provides insights into the evolutionary changes that occurred during European expansion of the invasive mosquitofish (Gambusia holbrooki). Mol Ecol. 2016;25:1089–105.
Willoughby JR, Harder AM, Tennessen JA, Scribner KT. Christie MR rapid genetic adaptation to a novel environment despite a genome-wide reduction in genetic diversity. Mol Ecol. 2018;27:4041–51.
Gillis MK, Walsh MR. Rapid evolution mitigates the ecological consequences of an invasive species (Bythotrephes longimanus) in lakes in Wisconsin. Proc R Soc B. 2017;284:20170814.
Chen Y, Shenkar N, Ni P, Lin Y, Li S, Zhan A. Rapid microevolution during recent range expansion to harsh environments. BMC Evol Biol. 2018;18:187.
Hartl DL, Clark AG, Clark AG. Principles of population genetics. Sunderland: Sinauer associates; 1997.
Li H, Stephan W. Inferring the demographic history and rate of adaptive substitution in drosophila. PLoS Genet. 2006;2:e166.
He J, Kreft H, Gao E, Wang Z, Jiang H. Patterns and drivers of zoogeographical regions of terrestrial vertebrates in China. J Biogeogr. 2017;44:1172–84.
Novak SJ, Mack RN. Genetic bottlenecks in alien plant species: influence of mating systems and introduction dynamics. In: Sax DF, Stachowicz JJ, Gaines SD, editors. Species invasions: insights into ecology, evolution, and biogeography. Sunderland: Sinauer Associates; 2005. p. 201–28.
Kolbe JJ, Glor RE, Schettino LR, Lara AC, Larson A, Losos JB. Multiple sources, admixture, and genetic variation in introduced Anolis lizard populations. Conserv Biol. 2007;21:1612–25.
Russell JC, Abdelkrim J, Fewster RM. Early colonisation population structure of a Norway rat island invasion. Biol Invasions. 2009;11:1557–67.
Williams JL, Fishman L. Genetic evidence for founder effects in the introduced range of houndstongue (Cynoglossum officinale). Biol Invasions. 2014;16:205–16.
Tsutsui ND, Suarez AV, Holway DA, Case TJ. Reduced genetic variation and the success of an invasive species. Proc Natl Acad Sci U S A. 2000;97:5948–53.
Carlson SM, Cunningham CJ, Westley PA. Evolutionary rescue in a changing world. Trends Ecol Evol. 2014;29:521–30.
Grossniklaus U, Kelly WG, Kelly B, Ferguson-Smith AC, Pembrey M, Lindquist S. Transgenerational epigenetic inheritance: how important is it? Nat Rev Genet. 2013;14:228–35.
Facon B, Hufbauer RA, Tayeh A, Loiseau A, Lombaert E, Vitalis R, Guillemaud T, Lundgren JG, Estoup A. Inbreeding depression is purged in the invasive insect Harmonia axyridis. Curr Biol. 2011;21:424–7.
Zhao L, Zhang JX, Zhang YH. Genetic boundary and gene flow between 2 parapatric subspecies of brown rats. Curr Zool. 2020;66:677–88.
Zhong XL, Liu T, Yang SQ. Research on climatic regionalization of the highway in arid area with large temperature difference. Tec Econ Area Commun. 2013;15:13–5.
Gutteridge WE, Coombs GH. Lipid Metabolism. In: Biochemistry of parasitic protozoa. London: Palgrave; 1977. p. 108–20.
Feulner PG, Chain FJ, Panchal M, Eizaguirre C, Kalbe M, Lenz TL, Mundry M, Samonte IE, Stoll M, Milinski M, Reusch TB, Bornberg-Bauer E. Genome-wide patterns of standing genetic variation in a marine population of three-spined sticklebacks. Mol Ecol. 2013;22:635–49.
Silliman K. Population structure, genetic connectivity, and adaptation in the Olympia oyster (Ostrea lurida) along the west coast of North America. Evol Appl. 2019;12:923–39.
Stuart YE, Campbell TS, Hohenlohe PA, Reynolds RG, Revell LJ, Losos JB. Rapid evolution of a native species following invasion by a congener. Science. 2014;346:463–6.
Hoban S, Kelley JL, Lotterhos KE, Antolin MF, Bradburd G, Lowry DB, Poss ML, Reed LK, Storfer A, Whitlock MC. Finding the genomic basis of local adaptation: pitfalls, practical solutions, and future directions. Am Nat. 2016;188:379–97.
Sax DF, Stachowicz JJ, Brown JH, Bruno JF, Dawson MN, Gaines SD, Grosberg RK, Hastings A, Holt RD, Mayfield MM, O'Connor MI, Rice WR. Ecological and evolutionary insights from species invasions. Trends Ecol Evol. 2007;22:465–71.
Kanarek AR, Webb CT. Original article: allee effects, adaptive evolution, and invasion success. Evol Appl. 2010;3:122–35.
Gibbs RA, Weinstock GM, Metzker ML, Muzny DM, Sodergren EJ, Scherer S, Scott G, Steffen D, Worley KC, Burch PE, et al. Genome sequence of the Brown Norway rat yields insights into mammalian evolution. Nature. 2004;428:493–521.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25:1754–60.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA. The genome analysis toolkit: a MapReduce framework for analyzing next–generation DNA sequencing data. Genome Res. 2010;20:1297–303.
Alexander DH, Novembre J, Lange K. Fast model–based estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–64.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stiratification in genome–wide association studies. Nat Genet. 2006;38:904–9.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.
Lee TH, Guo H, Wang X, Kim C, Paterson AH. SNPhylo: a pipeline to construct a phylogenetic tree from huge SNP data. BMC Genomics. 2014;15:162.
Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, McVean G, Durbin R. The variant call format and VCFtools. Bioinformatics. 2011;27:2156–8.
Zhang C, Dong SS, Xu JY, He WM, Yang TL. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics. 2019;35:1786–8.
Portik DM, Leaché AD, Rivera D, Barej MF, Burger M, Hirschfeld M, Rödel MO, Blackburn DC, Fujita MK. Evaluating mechanisms of diversification in a Guineo-Congolian forest frog using demographic model selection. Mol Ecol. 2017;26:5245–63.
Mattingsdal M, Jorde PE, Knutsen H, Jentoft S, Stenseth NC, Sodeland M, Robalo JI, Hansen MM, André C, Blanco GE. Demographic history has shaped the strongly differentiated corkwing wrasse populations in northern Europe. Mol Ecol. 2020;29:160–71.
Choi JY, Purugganan M, Stacy EA. Divergent selection and primary gene flow shape incipient speciation of a riparian tree on Hawaii Island. Mol Biol Evol. 2020;37:695–710.
Rubin CJ, Zody MC, Eriksson J, Meadows JR, Sherwood E, Webster MT, Jiang L, Ingman M, Sharpe T, Ka S, Hallböök F, Besnier F, Carlborg O, Bed'hom B, Tixier-Boichard M, Jensen P, Siegel P, Lindblad-Toh K, Andersson L. Whole-genome resequencing reveals loci under selection during chicken domestication. Nature. 2010;464:587–91.
Axelsson E, Ratnakumar A, Arendt ML, Maqbool K, Webster MT, Perloski M, Liberg O, Arnemo JM, Hedhammar A, Lindblad-Toh K. The genomic signature of dog domestication reveals adaptation to a starch-rich diet. Nature. 2013;495:360–4.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
We are grateful to Weiwei Zhai and Xiangjiang Zhan for their valuable suggestions on population genetics analyses.
This work was supported by the grants from the National Natural Science Foundation of China (Nos. 31672306 and 31872237), the Major Program of National Natural Science Foundation of China (32090022), the Strategic Priority Research Program of the Chinese Academy of Sciences (No. XDB11010400), and the State Key Laboratory of Integrated Management of Pest Insects and Rodents (ChineseIPM1816).
Ethics approval and consent to participate
The sample collection and handling complied with the Institutional Guidelines for Animal Use and Care at the Institute of Zoology, the Chinese Academy of Sciences, China (approval number IOZ12017).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sample information and sequencing statistics. Table S2. Cross-validation errors for different K values. Table S3. Best replicate of each of the optimized demographic models using ∂a∂i. Table S4. 42 genes in outlier regions between XJ and other native populations. Table S5. Functional enrichment of genes with allele frequency shifts throughout the genome for XJ rats.
Tree topology inferred from TreeMix with m = 1 (a) and m = 2 (b). Group IDs correspond to those in Fig. 1. Figure S2. Detection of genetic mixture across all subgroups. Significance of 3 Population Test (Z) represents whether the corresponding subgroup (on Y axis) is of mixed ancestry of other subgroups. Each dot indicates the Z score of a test between the target subgroup and every pair of other subgroups. Positive value suggests a result of unadmixed. All the groups were showed with only positive values, suggesting a relatively unadmixed relationship to other subgroups. Group IDs correspond to those in Fig. 1.
About this article
Cite this article
Chen, Y., Zhao, L., Teng, H. et al. Population genomics reveal rapid genetic differentiation in a recently invasive population of Rattus norvegicus. Front Zool 18, 6 (2021). https://doi.org/10.1186/s12983-021-00387-z
- Population genomics
- Biological invasion
- Demographic history
- Rapid differentiation
- Ancestral range
- Founder effect
- Rattus norvegicus