Genetic Structure and Diversity of Black Francolin in Uttarakhand, Western Himalaya, India
Total Page:16
File Type:pdf, Size:1020Kb
Journal of Wildlife and Biodiversity 4(1): 29-39 (2020) -(http://jwb.araku.ac.ir/) Research Article DOI: 10.22120/jwb.2019.115300.1093 Genetic structure and diversity of Black Francolin in Uttarakhand, Western Himalaya, India populations (2.99%). Overall the populations of 1* 2 Priyanka Negi , Atul Kathait , Tripti black francolin were genetically variable with Negi3, Pramesh Lakhera1 high adaptive potential in Uttarakhand, Western 1* Himalaya. Department of Zoology, Eternal University Baru sahib, India Keywords: Conservation, heterozygosity, 2 School of Life Science, ApeejaySatya University, homozygosity, genetic diversity, microsatellite Gurgaon, India markers, polymorphic sites 3School of Environment and Natural Resources, Doon University, Dehradun, India Introduction *email: [email protected] The increased pressure of habitat fragmentation, Received: 04 October 2019 / Revised: 12 November 2019 / Accepted: 02 November 2019 / Published online: 05 January overhunting and induced anthropogenic activity, 2020. Ministry of Sciences, Research and Technology, Arak put the number of species in threat. The University, Iran. Himalayan geographical region is facing the Abstract same. The protection of species requires a The present study evaluates the genetic diversity thorough understanding of species biology as well as of the level and distribution of genetic of black francolin (Francolinus francolinus diversity (Neel andEllstrand 2003). The latter is asiae) in Uttarakhand on the basis of an important factor for the adaptation and microsatellite loci. For this purpose, we survival of the natural populations in the examined five populations from three changing environment (Lande 1988; Landeand geographical zones of Uttarakhand, Western Shannon 1996). However current conservation Himalaya. Microsatellite markers were practice focused on captive breeding, which polymorphic with the number of alleles per reduces genetic diversity and this can severely locus ranging from 4-21, effective number of affect the capability of a population to cope with alleles per locus from 1.34 to 4.93, the the fluctuations in the environment (Landeand Polymorphic Information Content (PIC) value Shannon 1996, Frankham 2002, 2010). ranged from 0.22 to 0.85. The averaged The Himalayan ecosystem is a kind of sky observed heterozygosity across all loci was island, where the distribution of species restricts to their geographical boundaries. This provides Ho=0.320.12 and averaged expected uniqueness in genetic diversity among the heterozygosity He=0.510.06 respectively. The species. Galliformes are widely distributed in genetic structure showed that there were two the entire Himalayan range. The black francolin genetically distinct clusters. The Lesser is one of the key species of genus francolin with Himalayan and Himalayan foothill population the presence of five of species and six sub forming a single cluster and population of the species (Forcinaet al. 2012). Black francolin Tarai region forming another cluster. The (Francolinus francolinus asiae) widely pairwise FST results showed a sizeable genetic distributed from Kashmir to West Bengal along difference between the population of higher and with the foothills of central Nepal, Bihar, and lower altitudes. The AMOVA showed that thence through Punjab, Uttar Pradesh, Madhya higher levels of variation were observed among Pradesh to the Chilka lake in Odisha (Ali and Ripley 1983, McGowan and Kirwan 2015). individuals within populations (64.36%) and Although the Bird Life International (2016) lower differentiation observed among listed black francolin is stable at global level, but 30 | Journal of Wildlife and Biodiversity 4(1): 29-39 (2020) there were some reports indicate that its Material and Methods population is declining in many parts of its native range (Behbash et al. 2010, Riaz et al. Study Area 2011). The Uttarakhand (28044’ to 31028’N and 77035’ Despite habitat preference, feeding habit and to 81001’ E) encompasses 53,483 Km2 and is population abundance (Negiand and Lakhera situated in the Western Himalayan region of 2016, 2019), there is no significant study on India. It is India ‘s youngest mountain state and genetic diversity of the black francolin in India. according to biogeographical classification by the Wildlife Institute of India, the state comes In Indian Himalayan region, there is a strong into the Biogeographic Zones 2B Western need for conservation programs based on Himalaya and 7B Shiwaliks (Rodgers et al. genetic data. In the present work, our aim is to 2000). Uttarakhand is divided into five distinct carry out the first in-depth genetic study of black topographic regions i) snow cover Trans francolin in the Garhwal Himalayan region to Himalaya, ii) Sub Himalaya, iii) Lesser get useful information and offer logical Himalaya, iv) Bhabar Himalayan foothills and management recommendations for the purpose v) Tarai. The black francolin is distributed from of sustainable use and long-term protection of 100m to 2200m asl. So, to cover all the species within an adaptive conservation representative topographic zones, five sites were framework. In this work, the genetic makeup of selected. Two sites were selected from Lesser the black francolin resident in Western Himalaya (Rudraprayag and Jaunsar), one from Himalaya was characterized using microsatellite Bhabar Himalayan Foothills (Kotdwar) and two (short tandem repeats, STR) DNA markers. from Tarai region (Haridwar and Udham Singh Nagar) (Fig. 1). A total 56 samples were collected from the period of 2013 to 2015. Figure 1. Different sampling sites of black francolin from Uttarakhand, India 31 | Journal of Wildlife and Biodiversity 4(1): 29-39 (2020) DNA extraction and genotyping 64D4, MCW104, MCW 145). The important A combination of the standard phenol- details of the markers are given in Table1. PCR chloroform method was used for the extraction was performed in 25l reaction mixture, of DNA from tissues (Sambrook and Russell containing 50-100ng DNA, 100ng of each 2001) while a modified phenol-chloroform primer, 2.5mM of each dNTPs, 10X Taq assay method (Barbanera et al. 2005) and Qiagen buffer and 3U Taq polymerase. DNA isolation kit were also used for feathers. In The amplification conditions were: 4 min initial view to examine the quality and concentration of denaturation at 94°C, followed by 35 cycles of extracted DNA, 2 μl of each DNA sample was denaturation at 94°C for 40 s, annealing at a loaded into a 0.8% agarose gel containing specific temperature (Table 1) for 50 s and ethidium bromide. The concentration of DNA extension at 72°C for 1 min with a final was visualized under UV transilluminator. A extension at 72°C for 5 min. Subsequently, the total of eight polymorphic cross-amplified PCR products were separated by capillary microsatellite markers was standardized for the electrophoresis on analyzer ABI 3500XL genetic analysis of black francolin (Negi and (Applied Biosystems) and the allelic sizes of the Lakhera, 2018). All individuals were genotyped fragments were estimated based on with eight dinucleotide microsatellite loci fluorescently labeled forward primers (FAM, (MCW81, MCW330, MCW183, MCW11, and HEX) using the GENEMAPPER 3.7 MCW206, (Applied Biosystems, USA) and Peak analyzer 1.0. Table 1. Description of microsatellite loci used in the study Annealing Fluorescent Amplicon Loci 5’-3’ Sequence temperature Label size 0 TA ( C) MCW F GTTGCTGAGAGCCTGGTGCAG FAM 112-135bp 600C 81 RACACCCTGTATGTGGAATTACTTCTC MCW F CTTGACAGTGATGCATTAAATG HEX 221-249bp 500C 206 R ACATCTAGAATTGACTGTTCAC F GATCTCACCAGTATGAGCTGC P6A1 FAM 190-234bp 520C R TCTCACACTGTAACACAGTGC MCW FGTTATATTTAATGTCCACTTGTCAATGATG HEX 96-120 bp 600C 11 R TAAACCACTTCACATGGAGCCT MCW FTGCTGGACCTCATCAGTCTGACAG HEX 256-300bp 550C 330 RCAAACAAAATGTTCTCATAGAGTTCCTGC MCW F GATCCCAGTGTCGAGTATCCGA FAM 296-326bp 550C 183 R CTGAGATTTACTGGAGCCTGCC MCW FCAATTTAACTTTATTCTCCAAATTTGGCT HEX 180-210bp 520C 145 RGAGTAAACACAATGGCAACGGAAA MCW FCTTTTTAGCACAACTCAAGCTGTGAG FAM 210-225bp 550C 104 R CACAGACTTGCACAGCTGTGACC 32 | Journal of Wildlife and Biodiversity 4(1): 29-39 (2020) Data analysis Carlo (MCMCs), with 10 independent runs for Genetic diversity and population differentiation: each K, using the admixture model with Allele frequencies and genetic diversity correlated allele frequencies. The most reliable measures were calculated using GenAlEX 6.5 value of K was chosen according to Evannoet al. (Peakall and Smouse 2012) and FSTAT 2.9.4 2005, implemented in STRUCTURE (Goudet 1995). These measures included the HARVESTER (Earl and VonHoldt 2012). number of alleles (Na), number of private alleles Analysis of Bottleneck: For analysis of the (Ne), expected heterozygosity or Gene Diversity recent detection of bottlenecks in different (HE) and observed heterozygosity (HO). demographic population software GENEPOP 4.1 (Rousset 2008) was used to infer BOTTLENECK 1.202 (Piryet al. 1999) was possible departure from the Hardy-Weinberg used. In the present study, a one-tailed Wilcoxon Equilibrium (HWE) at each locus and evidence signed-rank test of heterozygosity of excess of linkage disequilibrium. Using allelic (Luikartet al. 1998) detected under the Stepwise frequencies, Polymorphic Information Content Mutation Model (SMM), the Infinite Allele (PIC), a measure of marker informativeness was Model (IAM) and the Two-Phase Model (TPM). calculated with CERVUS 3.0 (Kalinowskiet al. 2007). Genetic differences between populations Results were analyzed with pairwise FST (Weir and Genetic Diversity: A total