Whole-Genome Analysis Deciphers Population Structure and Genetic Introgression Among Bovine Species.

It is known that throughout history and presently, taurine (Bos taurus) and indicine/zebu (Bos indicus) cattle were crossed with other bovine species (e.g., gayal, gaur, banteng, yak, wisent, and bison). Information on the role of interspecific hybridization to facilitate faster adaptation of the newly arrived domestic species to new environments is poorly known. Herein, we collected 266 samples of bovine species of the taurine, zebu, yak, and gaur from West Europe, Southwest China, Indian subcontinent, and Southeast Asia to conduct the principal component analysis (PCA), admixture, gene flow, and selection signature analyses by using SNPs distributed across the bovine autosomes. The results showed that the genetic relationships between the zebu, yak, and gaur mirrored their geographical origins. Three ancestral components of the European taurine, East Asian taurine, and Indian zebu were found in domestic cattle, and the bidirectional genetic introgression between the Diqing cattle and Zhongdian yak was also detected. Simultaneously, the introgressed genes from the Zhongdian yak to the Diqing cattle were mainly enriched with immune-related pathways, and the ENPEP, FLT1, and PIK3CA genes related to the adaptation of high-altitude hypoxia were detected. Additionally, we found the genetic components of the Zhongdian yak had introgressed into Tibetan cattle. The 30 selected genes were detected in Tibetan cattle, which were significantly enriched in the chemokine signaling pathway. Interestingly, some genes (CDC42, SLC39A2, and EPAS1) associated with hypoxia response were discovered, in which CDC42 and SLC39A2 played important roles in angiogenesis and erythropoiesis, and heart function, respectively. This result showed that genetic introgression was one of the important ways for the environmental adaptation of domestic cattle.

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

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.

VARV B22R homologue as phylogenetic marker gene for Capripoxvirus classification and divergence time dating.

Sheeppox disease is associated with significant losses in sheep production world over. The sheep pox virus, the goatpox virus, and the lumpy skin disease virus cannot be distinguished by conventional serological tests. Identification of these pathogens needs molecular methods. In this study, seven genes viz. EEV maturation protein-F12L, Virion protein-D3R, RNA polymerase subunit-A5R, Virion core protein-A10L, EEV glycoprotein-A33R, VARV B22R homologue, and Kelch like protein-A55R that cover the start, middle, and end of the genome were selected. These genes were amplified from Roumanian-Fanar vaccine strain and Jaipur virulent strain, cloned, and sequenced. On analysis with the available database sequences, VARV B22R homologue was identified as a marker for phylogenetic reconstruction for classifying the sheeppox viruses of the ungulates. Further, divergence time dating with VARV B22R gene accurately predicted the sheeppox disease outbreak involving Jaipur virulent strain.

Genome sequencing of an Indian peste des petits ruminants virus isolate, Izatnagar/94, and its implications for virus diversity, divergence and phylogeography.

Peste des petits ruminants is an important transboundary disease infecting small ruminants. Genome or gene sequence analysis enriches our knowledge about the evolution and transboundary nature of the causative agent of this disease, peste des petits ruminants virus (PPRV). Although analysis using whole genome sequences of pathogens leads to more precise phylogenetic relationships, when compared to individual genes or partial sequences, there is still a need to identify specific genes/genomic regions that can provide evolutionary assessments consistent with those predicted with full-length genome sequences. Here the virulent Izatnagar/94 PPRV isolate was assembled and compared to all available complete genome sequences (currently in the NCBI database) to estimate nucleotide diversity and to deduce evolutionary relationships between genes/genomic regions and the full length genomes. Our aim was to identify the preferred candidate gene for use as a phylogenetic marker, as well as to predict divergence time and explore PPRV phylogeography. Among all the PPRV genes, the H gene was identified to be the most diverse with the highest evolutionary relationship with the full genome sequences. Hence it is considered as the most preferred candidate gene for phylogenetic study with 93% identity set as a nucleotide cutoff. A whole genome nucleotide sequence cutoff value of 94% permitted specific differentiation of PPRV lineages. All the isolates examined in the study were found to have a most recent common ancestor in the late 19th or in the early 20th century with high posterior probability values. The Bayesian skyline plot revealed a decrease in genetic diversity among lineage IV isolates since the start of the vaccination program and the network analysis localized the ancestry of PPRV to Africa.

Mitochondrial DNA analyses revealed low genetic diversity in the endangered Indian wild ass Equus hemionus khur.

The Indian wild ass Equus hemionus khur, belonging to ass-like equid branch, inhabits the dry and arid desert of the Little Rann of Kutch, Gujarat. The E. h. khur is the sole survivor of Asiatic wild ass species/subspecies in South Asia. To provide first ever insights into the genetic diversity, phylogeny, and demography of the endangered Indian wild ass, we sampled 52 free-ranging individuals from the Little Rann of Kutch by using a non-invasive methodology. The sequencing of 230 bp in cytochrome b (Cyt b) and displacement loop (D-loop) region revealed that current ∼4000 extant population of Indian wild ass harbours low genetic diversity. Phylogenetic analyses confirmed that E. h. khur, E. h. onager, and E. h. kulan belong to a single strict monophyletic clade. Therefore, we suggest the delimitation of the five E. hemionus subspecies in vogue to a single species E. hemionus. The application of molecular clock confirmed that the Asiatic wild ass had undergone diversification 0.65 Million years ago. Demographic measurements assessed using a Bayesian skyline plot demonstrated decline in the maternal effective population size of the Indian wild ass during different periods; these periods coincided with the origin and rise of the Indus civilization in the northwest of the Indian subcontinent during the Neolithic. In conclusion, maintaining high genetic diversity in the existing isolated population of 4000 Indian wild asses inhabiting the wild ass sanctuary is important compared with subspecies preservation alone.

Genetic differentiations between randomly and selectively bred pig populations in Yunnan, China.

To assess the genetic diversity between randomly and selectively bred populations, we sequenced 438 bp of the mitochondrial DNA control region from 102 pigs. These samples represented four native pig breeds, one nucleus and one conservation herd from Yunnan, China. Twenty haplotypes with sixteen polymorphic sites were identified. The number of haplotypes in the nucleus herd of Saba pig and the conservation herd of Banna miniature pig were restricted to three and one, respectively, while the randomly bred pig populations exhibited over six haplotypes. Notably, haplotype diversity in randomly bred populations was significantly greater than the selectively bred populations (h=0.732 vs. 0.425 and 0, exact test, P<=0.0036). These findings demonstrate that selective breeding generated low genetic diversity compared to randomly bred pig breeds. A timely intervention and well programmed breeding approach would stop further genetic diversity reduction in the nucleus and conservation herds of native pig breeds. Otherwise, selective breeding would dramatically reduce genetic diversity in only several years, indicating that sharp contradictions exist between breeding, conservation and genetic diversity. Genetic relationships are discussed based on net genetic distances among pig populations.

Genetic affinities within a large global collection of pathogenic Leptospira: implications for strain identification and molecular epidemiology.

Leptospirosis is an important zoonosis with widespread human health implications. The non-availability of accurate identification methods for the individualization of different Leptospira for outbreak investigations poses bountiful problems in the disease control arena. We harnessed fluorescent amplified fragment length polymorphism analysis (FAFLP) for Leptospira and investigated its utility in establishing genetic relationships among 271 isolates in the context of species level assignments of our global collection of isolates and strains obtained from a diverse array of hosts. In addition, this method was compared to an in-house multilocus sequence typing (MLST) method based on polymorphisms in three housekeeping genes, the rrs locus and two envelope proteins. Phylogenetic relationships were deduced based on bifurcating Neighbor-joining trees as well as median joining network analyses integrating both the FAFLP data and MLST based haplotypes. The phylogenetic relationships were also reproduced through Bayesian analysis of the multilocus sequence polymorphisms. We found FAFLP to be an important method for outbreak investigation and for clustering of isolates based on their geographical descent rather than by genome species types. The FAFLP method was, however, not able to convey much taxonomical utility sufficient to replace the highly tedious serotyping procedures in vogue. MLST, on the other hand, was found to be highly robust and efficient in identifying ancestral relationships and segregating the outbreak associated strains or otherwise according to their genome species status and, therefore, could unambiguously be applied for investigating phylogenetics of Leptospira in the context of taxonomy as well as gene flow. For instance, MLST was more efficient, as compared to FAFLP method, in clustering strains from the Andaman island of India, with their counterparts from mainland India and Sri Lanka, implying that such strains share genetic relationships and that leptospiral strains might be frequently circulating between the islands and the mainland.

Zebu cattle are an exclusive legacy of the South Asia neolithic.

Animal domestication was a major step forward in human prehistory, contributing to the emergence of more complex societies. At the time of the Neolithic transition, zebu cattle (Bos indicus) were probably the most abundant and important domestic livestock species in Southern Asia. Although archaeological evidence points toward the domestication of zebu cattle within the Indian subcontinent, the exact geographic origins and phylogenetic history of zebu cattle remains uncertain. Here, we report evidence from 844 zebu mitochondrial DNA (mtDNA) sequences surveyed from 19 Asiatic countries comprising 8 regional groups, which identify 2 distinct mitochondrial haplogroups, termed I1 and I2. The marked increase in nucleotide diversity (P < 0.001) for both the I1 and I2 haplogroups within the northern part of the Indian subcontinent is consistent with an origin for all domestic zebu in this area. For haplogroup I1, genetic diversity was highest within the Indus Valley among the three hypothesized domestication centers (Indus Valley, Ganges, and South India). These data support the Indus Valley as the most likely center of origin for the I1 haplogroup and a primary center of zebu domestication. However, for the I2 haplogroup, a complex pattern of diversity is detected, preventing the unambiguous pinpointing of the exact place of origin for this zebu maternal lineage. Our findings are discussed with respect to the archaeological record for zebu domestication within the Indian subcontinent.

Multiple maternal origins of chickens: out of the Asian jungles.

Domestic chickens have long been important to human societies for food, religion, entertainment, and decorative uses, yet the origins and phylogeography of chickens through Eurasia remain uncertain. Here, we assessed their origins and phylogeographic history by analyzing the mitochondrial DNA hypervariable segment I (HVS-I) for 834 domestic chickens (Gallus gallus domesticus) across Eurasia as well as 66 wild red jungle fowls (Gallus gallus) from Southeast Asia and China. Phylogenetic analyses revealed nine highly divergent mtDNA clades (A-I) in which seven clades contained both the red jungle fowls and domestic chickens. There was no breed-specific clade in the chickens. The clades A, B, and E are distributed ubiquitously in Eurasia, while the other clades were restricted to South and Southeast Asia. Clade C was mainly distributed in Japan and Southeast China, while clades F and G were exclusive to Yunnan, China. The geographic distribution of clade D was closely related to the distribution of the pastime of cock fighting. Statistical tests detect population expansion within each subclade. These distinct distribution patterns and expansion signatures suggest that different clades may originate from different regions, such as Yunnan, South and Southwest China and/or surrounding areas (i.e., Vietnam, Burma, and Thailand), and the Indian subcontinent, respectively, which support the theory of multiple origins in South and Southeast Asia.