Phylogeography and Origin of Indian Domestic Goats

The evolution of contemporary livestock species: Insights from mitochondrial genome

The domestication of animals marks a pivotal moment in human history, profoundly influencing our demographic and cultural progress. This process has led to significant genetic, behavioral, and physical changes in livestock species compared to their wild ancestors. Understanding the evolutionary history and genetic diversity of livestock species is crucial, and mitochondrial DNA (mtDNA) has emerged as a robust marker for investigating molecular diversity in animals. Its highly conserved gene content across animal species, minimal duplications, absence of introns, and short intergenic regions make mtDNA analysis ideal for such studies. Mitochondrial DNA analysis has uncovered distinct cattle domestication events dating back to 8000 years BC in Southwestern Asia. The sequencing of water buffalo mtDNA in 2004 provided important insights into their domestication history. Caprine mtDNA analysis identified three haplogroups, indicating varied maternal origins. Sheep, domesticated 12,000 years ago, exhibit diverse mtDNA lineages, suggesting multiple domestication events. Ovine mtDNA studies revealed clades A, B, C, and a fourth lineage, group D. The origins of domestic pigs were traced to separate European and Asian events followed by interbreeding. In camels, mtDNA elucidated the phylogeographic structure and genetic differentiation between wild and domesticated species. Horses, domesticated around 3500 BC, show significant mtDNA variability, highlighting their diverse origins. Yaks exhibit unique adaptations for high-altitude environments, with mtDNA analysis providing insights into their adaptation. Chicken mtDNA studies supported a monophyletic origin from Southeast Asia's red jungle fowl, with evidence of multiple origins. This review explores livestock evolution and diversity through mtDNA studies, focusing on cattle, water buffalo, goat, sheep, pig, camel, horse, yak and chicken. It highlights mtDNA's significance in unraveling maternal lineages, genetic diversity, and domestication histories, concluding with insights into its potential application in improving livestock production and reproduction dynamics.

Mitochondrial Genomes of Korean Native Black Goats Reveal Shared Phylogeographic Patterns and Demographic History

This study explores the phylogeny of Korean native black goats through analysis of their complete mitochondrial DNA. The National Institute of Animal Science has gathered genetic material on purebred goats from isolated regions such as Tongyeong, Dangjin, and Jangsu, and is actively breeding them on a national level. These populations, however, are small and exhibit high inbreeding rates, highlighting the urgent need to preserve genetic diversity. The haplotype diversity within this native group is 0.659, with 39 haplotypes identified. By contrast, including international breeds in the analysis increases the overall haplotype diversity to 0.925 with 203 haplotypes identified, highlighting the limited genetic diversity among native black goats. For phylogenetic assessment, a neighbor-joining tree and median-joining network were constructed using identified haplogroups (A, B, C, D, G, and F) from prior studies. The results pinpoint the native black goats as closely related to, but distinct from, Haplogroup A with a bootstrap value of 98, establishing them as a separate clade (A'). This supports the notion of a shared ancestry with various global populations. This research provides essential data on the origins and evolutionary history of Korean native black goats, supporting conservation and breeding efforts aimed at enhancing genetic diversity.

Mitochondrial DNA D-Loop Polymorphisms among the Galla Goats Reveals Multiple Maternal Origins with Implication on the Functional Diversity of the HSP70 Gene

Despite much attention given to the history of goat evolution in Kenya, information on the origin, demographic history, dispersal route, and genetic diversity of Galla goats remains unclear. Here, we examined the genetic background, diversity, demographic history, and population genetic variation of Galla goats using mtDNA D-loop and HSP70 single-nucleotide polymorphism markers. The results revealed 90 segregating sites and 68 haplotypes in a 600-bp mtDNA D-loop sequence. The overall mean mitochondrial haplotype diversity was 0.993. The haplotype diversities ranged between 0.8939 ± 0.0777 and 1.0000 ± 0.0221 in all populations supporting high genetic diversity. Mitochondrial phylogenetic analysis revealed three Galla goat haplogroups (A, G, and D), supporting multiple maternal ancestries, of which haplogroup A was the most predominant. Analysis of molecular variance (AMOVA) showed considerable variation within populations at 94.39%, evidence of high genetic diversity. Bimodal mismatch distribution patterns were observed while most populations recorded negative results for Tajima and Fu's Fs neutrality tests supporting population expansion. Genetic variation among populations was also confirmed using HSP70 gene fragment sequences, where six polymorphic sites which defined 21 haplotypes were discovered. Analysis of molecular variance revealed a significant FST index value of 0.134 and a high FIS index value of 0.746, an indication of inbreeding. This information will pave the way for conservation strategies and informed breeding to improve Galla or other goat breeds for climate-smart agriculture.

Design and validation of high-density SNP array of goats and population stratification of Indian goat breeds

Goats are the supporting pillars of rural economy contributing significantly to meat and milk production in India. It is a species targeted for fulfilling the interdependent goals of poverty reduction and creation of employment for supporting the rural income. The increased demand for goat products necessitates their genetic characterization and improvement to augment the production of native breeds. Bi-allelic, genome wide, densely placed single nucleotide polymorphism (SNP) markers are most suitable for this purpose. This paper describes the design and validation of an Affymetrix Axiom-based high-density (HD) SNP chip for goats. The array was designed using a panel of 225 samples from 15 diverse goat breeds of India. In total, more than 38 million high quality SNPs were subjected to stringent filtering and 626,975 SNPs were finally tiled on the array. The average coverage of SNPs in our chip is one SNP per four kilobase (kb), providing a denser coverage of the goat genome than previously available arrays. The HD chip (Axiom_Cahi) was validated by genotyping 443 samples from 26 indigenous goat breeds/populations. The results revealed 95.83% markers to be highly informative and polymorphic in Indian goats. Multivariate analysis indicated population structuring, as 15 breeds could be segregated using the designed array. Phylogenetic analysis suggested stratification of breeds by geographic proximity. This HD SNP chip for goats is a valuable resource for genomic selection, genome wide association as well as population genetic studies in goats.

Explicating the genetic diversity and population structure of Saanen × Beetal goats using pedigree analysis

Pedigree analysis is required to assess the genetic diversity and population structure of a close breeding population in order to effectively manage the breeding program and keep inbreeding under acceptable limits. Saanen × Beetal is a crossbred population of goats, reared at the National Dairy Research Institute (NDRI), Karnal, for the last five decades. This germplasm has been acclimated to a tropical climate and has a higher milk potential and prolificacy. The objective of this study was to elucidate the genetic diversity, population structure, and inbreeding in the flock of the Saanen × Beetal goats. The data were collected from the Animal Genetics and Breeding Division of ICAR-NDRI, Karnal, for 2603 animals from the year 1971 to 2021. Animals born between 2014 and 2017 were considered as a reference cohort. Results revealed that the average generation interval was 3.44 years for the complete pedigree. The average inbreeding coefficient and the average relatedness were 4.20% and 6.87%, respectively, for the complete pedigree and 10.78% and 10.80% for the reference population. Higher inbreeding coefficient and average relatedness in the reference cohort demonstrated the impact of the enclosed gene pool and demands immediate intervention for managing diversity in the closed nucleus under study. Ancestors contributing 50% of the gene pool were 8 and 3 for the complete pedigree and reference cohort, respectively, which illustrates the fact that very few ancestors were responsible for genetic diversity in the flock, which results in the decline of effective population size. Effective numbers of founders (fe), ancestors (fa), and founder genome equivalents (fg) were 15, 7, and 3.11, respectively. The (fe/fa) ratio in the reference population was 2.14, indicating the occurrence of the bottleneck effect in the flock. We observed that inbreeding was non-significant for all reproductive traits except for age at first service and age at first kidding. To lessen inbreeding and augment genetic diversity in the flock, the stratified breeding plan needs to be followed, where mate selection would be based on relatedness. Furthermore, the introduction of unrelated Saanen and Beetal crosses will help alleviate the inbreeding accumulation.

A. B. D. A, Abbas K, Hussain T, Ramasamy S, Manomohan V, Tapsoba ASR, Pichler R, Babar ME, Periasamy K. Indigenous cattle of Sri Lanka: Genetic and phylogeographic relationship with Zebu of Indus Valley and South Indian origin

The present study reports the population structure, genetic admixture and phylogeography of cattle breeds of Sri Lanka viz. Batu Harak, Thawalam and White cattle. Moderately high level of genetic diversity was observed in all the three Sri Lankan zebu cattle breeds. Estimates of inbreeding for Thawalam and White cattle breeds were relatively high with 6.1% and 7.2% respectively. Genetic differentiation of Sri Lankan Zebu (Batu Harak and White cattle) was lowest with Red Sindhi among Indus Valley Zebu while it was lowest with Hallikar among the South Indian cattle. Global F statistics showed 6.5% differences among all the investigated Zebu cattle breeds and 1.9% differences among Sri Lankan Zebu breeds. The Sri Lankan Zebu cattle breeds showed strong genetic relationships with Hallikar cattle, an ancient breed considered to be ancestor for most of the Mysore type draught cattle breeds of South India. Genetic admixture analysis revealed high levels of breed purity in Lanka White cattle with >97% Zebu ancestry. However, significant taurine admixture was observed in Batu Harak and Thawalam cattle. Two major Zebu haplogroups, I1 and I2 were observed in Sri Lankan Zebu with the former predominating the later in all the three breeds. A total of 112 haplotypes were observed in the studied breeds, of which 50 haplotypes were found in Sri Lankan Zebu cattle. Mismatch analysis revealed unimodal distribution in all the three breeds indicating population expansion. The sum of squared deviations (SSD) and raggedness index were non-significant in both the lineages of all the three breeds except for I1 lineage of Thawalam cattle (P<0.01) and I2 lineage of Batu Harak cattle (P<0.05). The results of neutrality tests revealed negative Tajima's D values for both the lineages of Batu Harak (P>0.05) and White cattle (P>0.05) indicating an excess of low frequency polymorphisms and demographic expansion. Genetic dilution of native Zebu cattle germplasm observed in the study is a cause for concern. Hence, it is imperative that national breeding organizations consider establishing conservation units for the three native cattle breeds to maintain breed purity and initiate genetic improvement programs.

Peeping into Mitochondrial Diversity of Andaman Goats: Unveils Possibility of Maritime Transport with Diversified Geographic Signaling

Andaman and Nicobar Islands, a part of South-East Asia, is enriched with the presence of native breeds of livestock (cattle, pig, goat) and poultry. There are two native goat breeds, viz., Andaman local goat and Teressa goat in Andaman and Nicobar Islands. However, to date, the origin and genetic makeup of these two breeds have not been detailed. Therefore, the present study describes the genetic makeup of Andaman goats through analysis of mitochondrial D-loop sequence for sequence polymorphism, phylogeographical signaling and population expansion events. The genetic diversity of the Teressa goat was less compared to the Andaman local goat due to its sole presence on Teressa Island. Out of 38 well-defined haplotypes of Andaman goats, the majority of haplotypes belonged to haplogroup A followed by haplogroup B and haplogroup D. The result of mismatch distribution and neutrality tests indicated no population expansion event of haplogroup A and B. Finally, based on poor geographical signaling, we hypothesize that Andaman goats have been imported to these Islands either through multidirectional diffusion or unidirectional diffusion. We justify our hypothesis of multidirectional diffusion on the basis of observation of the haplotype and nucleotide diversity of Andaman goats. Simultaneously, the probability of unidirectional diffusion of goats in these islands from the Indian subcontinent in different spells of domestication events through maritime routes cannot be ignored.

Human and Animal Fascioliasis: Origins and Worldwide Evolving Scenario

Fascioliasis is a plant- and waterborne zoonotic parasitic disease caused by two trematode species: (i) Fasciola hepatica in Europe, Asia, Africa, the Americas, and Oceania and (ii) F. gigantica, which is restricted to Africa and Asia. Fasciolid liver flukes infect mainly herbivores as ruminants, equids, and camelids but also omnivore mammals as humans and swine and are transmitted by freshwater Lymnaeidae snail vectors. Two phases may be distinguished in fasciolid evolution. The long predomestication period includes the F. gigantica origin in east-southern Africa around the mid-Miocene, the F. hepatica origin in the Near-Middle East of Asia around the latest Miocene to Early Pliocene, and their subsequent local spread. The short postdomestication period includes the worldwide spread by human-guided movements of animals in the last 12,000 years and the more recent transoceanic anthropogenic introductions of F. hepatica into the Americas and Oceania and of F. gigantica into several large islands of the Pacific with ships transporting livestock in the last 500 years. The routes and chronology of the spreading waves followed by both fasciolids into the five continents are redefined on the basis of recently generated knowledge of human-guided movements of domesticated hosts. No local, zonal, or regional situation showing disagreement with historical records was found, although in a few world zones the available knowledge is still insufficient. The anthropogenically accelerated evolution of fasciolids allows us to call them "peridomestic endoparasites." The multidisciplinary implications for crucial aspects of the disease should therefore lead the present baseline update to be taken into account in future research studies.

Mitochondrial DNA D-loop sequence analysis reveals high variation and multiple maternal origins of indigenous Tanzanian goat populations

The Small East African (SEA) goat are widely distributed in different agro-ecological zones of Tanzania. We report the genetic diversity, maternal origin, and phylogenetic relationship among the 12 Tanzanian indigenous goat populations, namely Fipa, Songwe, Tanga, Pwani, Iringa, Newala, Lindi, Gogo, Pare, Maasai, Sukuma, and Ujiji, based on the mitochondrial DNA (mtDNA) D-loop. High haplotype (H d = 0.9619-0.9945) and nucleotide (π = 0.0120-0.0162) diversities were observed from a total of 389 haplotypes. The majority of the haplotypes (n = 334) belonged to Haplogroup A which was consistent with the global scenario on the genetic pattern of maternal origin of all goat breeds in the world. Haplogroup G comprised of 45 haplotypes drawn from all populations except the Ujiji goat population while Haplogroup B with 10 haplotypes was dominated by Ujiji goats (41%). Tanzanian goats shared four haplotypes with the Kenyan goats and two with goats from South Africa, Namibia, and Mozambique. There was no sharing of haplotypes observed between individuals from Tanzanian goat populations with individuals from North or West Africa. The indigenous goats in Tanzania have high genetic diversity defined by 389 haplotypes and multiple maternal origins of haplogroup A, B, and G. There is a lot of intermixing and high genetic variation within populations which represent an abundant resource for selective breeding in the different agro-ecological regions of the country.

Microsatellite markers based genetic diversity and differentiation of Balangir goat population of Odisha

This study was executed to investigate the genetic diversity of Balangir, a local goat population from Odisha, and its differentiation from the geographically closest registered goat breeds, Ganjam and Black Bengal by utilizing 22 polymorphic microsatellite markers. The genotypic status of individuals at each locus was identified by an automated DNA sequencer and allelic data was analyzed for genetic diversity parameters. It resulted in the acquisition of original and comparable information confirming the genetic distinctness of the Balangir goat population. The observed number of alleles varied between 4 (RM4) and 16 (OMHC1) with 8.23±0.64 alleles per locus in Balangir goats. A moderate level of observed heterozygosity (0.57±0.05) indicated sufficiently existing genetic diversity in this goat population which varied between 0.09 (OarJMP29) and 1 (ILSTS22). Expected heterozygosity (0.64±0.04) oscillated between 0.18 (ILSTS34) and 0.89 (OMHC1). Within the population, heterozygosity deficiency estimate (FIS) was significantly positive (0.13±0.05) and 73% of the investigated loci contributed towards the homozygote excess. Consequently, significant deviation from Hardy–Weinberg equilibrium was observed at 12 of 22 analyzed loci. Assessment of Balangir goat population for mutation drift equilibrium indicated that the population did not suffer a reduction in effective population size in the last few generations. Four different approaches utilized to study genetic relationships (F statistics, analysis of molecular variance, phylogenetic relationship genetic distance, and individual assignment) indicated that the Balangir goat population was genetically differentiated from the two registered goat breeds. The study highlighted that Balangir is discrete Indian goat germplasm with sufficient genetic variability and needs to be registered as a goat breed of India.

Transcriptomic diversity in longissimus thoracis muscles of Barbari and Changthangi goat breeds of India

The present study is an attempt to examine the differential expression of genes in longissimus thoracis muscles between meat and wool type Indian goat breeds. Barbari goat is considered the best meat breed while Changthangi is famous for its fine fibre quality. RNA sequencing data was generated from four biological replicates of longissimus thoracis muscles of Barbari and Changthangi goats. A clear demarcation could be observed between the breeds in terms of expression of genes associated with lipid metabolism (FASN, SCD, THRSP, DGAT2 and FABP3). Most significant genes with high connectivity identified by gene co-expression network analysis were associated with triacylglycerol biosynthesis pathway in Barbari goat. Highly interactive genes identified in Changthangi goat were mainly associated with muscle fibre type. This study provides an insight into the differential expression of genes in longissimus thoracis muscles between Barbari and Changthangi goats that are adapted to and reared in different agro-climatic regions.

Genetic diversity and population structure of Jamunapari goat in India using pedigree analysis

Pedigree records of 6821 Jamunapari goats of India were collected from 1980 to 2011 and used to evaluate the population structure and genetic diversity in this flock. Animals born between 2009 and 2011 represented the current reference population. The average pedigree completeness index (PCI) and numbers of equivalent complete generations (EqG) were estimated for the entire (PCI = 0.18, EqG = 2.24) and reference (PCI = 0.31, EqG = 3.45) populations. The average generation interval was 3.33 years. The average inbreeding coefficient and the average relatedness were 0.46 and 1.06%, respectively, for the entire population and 0.77 and 3.87% for the reference population. The rate of inbreeding was 0.06% per generation. The effective population size (N_e), estimated from increases in inbreeding coefficients between the first and third equivalent complete generations, was 52.65, but periodic introductions of unrelated breeding males resulted in average inbreeding levels in the reference population that were lower than those predicted from the estimate of N_e. Effective numbers of founders (f_e), ancestors (f_a), founder genomes equivalents (f_g), and non-founder genomes (f_ng) were 51, 39, 25.8, and 48.2, respectively. The f_e/f_a ratio in the reference population was 1.31 and indicated that occasional bottlenecks had occurred in the population. The 14 most influential ancestors contributed 50% of the genetic variability in the reference population, with a maximum individual contribution of 9.25%. Approximately 1.9% of the initial heterozygosity had been lost from the population, indicating that substantial genetic diversity still exists in this flock.

Global Goat! Is the Expanding Goat Population an Important Reservoir of Cryptosporidium?

Goats are a primary or additional income source for many families in resource-poor areas. Although often considered inferior to other livestock, the resilience of goats and their ability to thrive in a range of environments means that that they are of particular value. Furthermore, goats emit less methane than other livestock species. In these same areas, it is well-documented that cryptosporidiosis has a substantial impact on infant morbidity and mortality, as well as reducing child growth and development. As Cryptosporidium also causes diarrheal disease in goats, the question arises whether goats may represent a reservoir of infection to humans. Epidemiological studies regarding the potential for transmission of Cryptosporidium between goats and humans have largely concluded that Cryptosporidium species infecting goats are not zoonotic. However, these studies are mostly from developed countries, where goat husbandry is smaller, management routines differ greatly from those of developing countries, contact between goats and their owners is more limited, and cryptosporidiosis has less impact on human health. In this article, background information on goat husbandry in different countries is provided, along with information on Cryptosporidium prevalence among goats, at both the species and sub-species levels, and the potential for zoonotic transmission. The intention is to indicate data gaps that should be filled and to increase awareness of the role of goats as providers for low-income families, often living in areas where cryptosporidiosis is endemic and where appropriate baseline interventions could have a positive impact, regardless of species of goat or parasite.

Indonesian native goats (Capra hircus) reveal highest genetic frequency of mitochondrial DNA haplogroup B in the world

The objective of this study was to determine mtDNA sequences of the 481 bp HV1 region from two Indonesian native goat breeds, Kacang and Marica, to confirm the phylogeographic distribution of caprine haplogroup B in Southeast Asia. Based on these sequences, 12 haplotypes were observed and categorized into the predominant haplogroup B and minority haplogroup A, indicating that Indonesian native goats present the highest frequency (0.950) of the haplogroup B in the world. These results strongly emphasize previous observations of the haplogroup B frequencies tending to increase southeastward in Southeast Asia. Additionally, this suggests that goats primary bred into Southeast Asia might predominantly carry haplogroup B. This could essentially contribute to the understanding of the origin, propagation route and/or introgression history of Southeast Asian goats.

Molecular and phenometric characterization of Bhakarwali goat breed of India

Bhakarwali is recently registered as 34th goat breed of India. It is distributed in the hilly tracts of Jammu and Kashmir. The breed is distinct with superior qualities such as high temperature resistance and milk, meat and fiber productivity under the low input system. It’s characterization at phenotypic level was carried out by surveying the breeding tract and at genotypic level by microsatellite markers. Information on body traits, performance traits and managemental practices were collected by interviewing the goat keepers. All the microsatellite loci selected for diversity analysis were polymorphic and a total of 190 alleles were identified across the 23 microsatellite loci. OMHC1 depicted the highest number of alleles (15) while ILSTS065 had the lowest (2) with 8.26±0.663 mean number of alleles per locus. Expected number of alleles ranged from 1.065 (ILSTS044) to 6.755 (OMHC1) with a mean value of 3.613±0.367 alleles per locus. The observed heterozygosity ranged from 0.063 (ILSTS044) to 0.915 (OMHC1) with a mean of 0.629±0.045. Corresponding values of expected heterozygosity varied between 0.061 (ILSTS044) to 0.852 (OMHC1) with a mean of 0.639±0.043. Heterozygote deficiency was negligible as average FIS value was only 0.002±0.033. Bottleneck was examined using all the three mutations models and was found to be absent. Normal L-shaped curve indicated lack of mode shift in the population. This is the first-hand report on current diversity status of Bhakarwali goat and is expected to be useful in planning conservation and in facilitating their effective use in future breeding programs.

Genetic diversity analysis among Indian goat breeds based on mitochondrial DNA

India ranks second in goat population with 34 genetically recognized and registered breeds. Information on their diversity and origin and ancestry is little known. Hence, the mtDNA based genetic diversity analysis of Indian goats; targeting mitochondrial HVR1 region from 21 Indian breeds belonging to different geographical regions was undertaken. A total of 124 haplotypes were identified and haplotype diversity estimate ranged from 0.67 to 1.0 with an average value of 0.99. The average nucleotide diversity was minimum (0.02) in Kanniadu and maximum in Surti breeds. Analysis of molecular variance revealed 5.16% variation among the breeds and 94.84% within breeds indicating weak phylogeographic structure. Neighbor-joining tree analysis revealed that the maximum number of individuals of Indian goats fall under A and few in B and C lineages. Principal component analyses of the Indian goat breeds revealed that Kanniadu goats clustered distantly from rest of the breeds of the country. Mantel test revealed a significant correlation between FST and geographical distance (r=0.29) among the goat breeds. The mismatch distribution analysis of the Indian goat breeds revealed bimodal distribution patterns. The analysis revealed that Kanniadu is highly distinct from the rest of the breeds.

Maternal Origins and Haplotype Diversity of Seven Russian Goat Populations Based on the D-loop Sequence Variability

Simple Summary

Russia has diverse specifically selected and multipurpose goat resources. However, the origin of the local goats is still enigmatic. In this study, we sequenced and analyzed mitochondrial DNA (mtDNA) fragments of seven Russian local goat populations to provide the first insight into their maternal lineage.

Abstract

The territory of modern Russia lies on the crossroads of East and West and covers various geographical environments where diverse groups of local goats originated. In this work, we present the first study on the maternal origin of Russian local goats, including Altai Mountain (n = 9), Dagestan Downy (n = 18), Dagestan Local (n = 12), Dagestan Milk (n = 15), Karachaev (n = 21), Orenburg (n = 10), and Soviet Mohair (n = 7) breeds, based on 715 bp D-loop mitochondrial DNA (mtDNA) sequences. Saanen goats (n = 5) were used for comparison. Our findings reveal a high haplotype (HD = 0.843–1.000) and nucleotide diversity (π = 0.0112–0.0261). A total of 59 haplotypes were determined in the Russian goat breeds, in which all differed from the haplotypes of the Saanen goats. The haplotypes identified in Altai Mountain, Orenburg, Soviet Mohair, and Saanen goats were breed specific. Most haplotypes (56 of 59) were clustered together with samples belonging to haplogroup A, which was in accordance with the global genetic pattern of maternal origin seen in most goats worldwide. The haplotypes that were grouped together with rare haplogroups D and G were found in the Altai Mountain breed and haplogroup C was detected in the Soviet Mohair breed. Thus, our findings revealed that local goats might have been brought to Russia via various migration routes. In addition, haplotype sharing was found in aboriginal goat populations from overlapping regions, which might be useful information for their official recognition status.

Maternal origins, population structure and demographic history of ten Chinese indigenous goat breeds from Yunnan

Yunnan as a frontier zone that connects China with South and Southeast Asia, has 11 well-recognized goat breeds. However, the knowledge about maternal origins, population structure and demographic history of Chinese indigenous goats from Yunnan is limited. In this study, we analysed a 481-bp fragment of first hypervariable segment (HVSI) of the mitochondrial DNA (mtDNA) control region sequences of 749 individuals from 10 Yunnan indigenous goat breeds, of which 556 sequences were newly determined. There were 110 polymorphic sites that defined 158 haplotypes among all sequences. The haplotype and nucleotide diversity of these breeds ranged from 0.782 ± 0.079 to 0.982 ± 0.015 and from 0.028 ± 0.003 to 0.043 ± 0.005, respectively. Phylogenetic analysis identified two lineages A and B, of which the lineage A had higher frequency (68.1%) and distributed in all Yunnan breeds. We combined previously reported sequences with our sequences belonging to the lineage B and detected two subclades B1 and B2, in which the B1 subclade shared individuals from Eastern Asia, Southeast Asia and Southern Asia. Given higher level of diversity and more unique haplotypes, the B2 subclade probably originated from Southwestern China. The haplotype network, analysis of molecular variance (AMOVA) and a Mantel test revealed no significant phylogeographic structuring among Yunnan goat breeds. This can be explained by high gene flow and genetic admixture among these breeds from different geographic regions in Yunnan. Additionally, both the lineages A and B reflected different demographic histories. This study will provide a scientific basis for the conservation and utilization of Yunnan indigenous goats.

Comprehensive analysis of mitochondrial DNA based genetic diversity in Indian goats

Thirty four distinct breeds and many non-descript populations represent the caprine diversity of India. Genetic characterization of breeds is an essential element in designing breeding strategies and preserving genetic diversity. Considering the popularity of mitochondrial DNA for phylogeographical studies, this study involved an extensive analysis of population structure and genetic diversity of 28 defined breeds and 5 lesser known populations representing all four major agro-climatic zones of India using mitochondrial DNA markers. Analysis of hypervariable region 1 of mtDNA control region in 443 goats together with 22 reference sequences, delineated 341 distinct haplotypes belonging to four maternal haplogroups; A, B, C and D, with haplogroup A representing 90% of the individuals. The haplotype and nucleotide diversity indices of Indian goats were 0.998 ± 0.001 and 0.028 ± 0.001, respectively indicating abundant genetic variability. Estimates of population demographic parameters from mismatch analysis suggested a relatively good fit to the model of either spatial or demographic expansion of Indian goats. AMOVA analysis and topology of MJ network suggested lack of phylogeographic structure in domestic goats, which can be attributed to unstructured animal breeding, dwindling pastures and nomadic pastoralism. Genetic differentiation between goats from different agro-ecological regions was in accordance with their geographical propinquity.

The Eurasian Steppe is an important goat propagation route: A phylogeographic analysis using mitochondrial DNA and Y-chromosome sequences of Kazakhstani goats

Goats (Capra hircus) were domesticated in the Fertile Crescent and propagated all over the world. The Silk Road through the Eurasian Steppe belt is a possible propagation route for domestic goats to Central Asia. Kazakhstan is in close geographical proximity to domestication centers and covers the majority of the Eurasian Steppe belt. In this study, we examined the genetic diversity and phylogeographic structure of Kazakhstani goats. The mtDNA sequences of 141 Kazakhstani goats were categorized into haplogroups A, C, and D, of which haplogroup A was predominant (97%), whereas haplogroups C and D were detected at low frequencies (1.4% each). The Kazakhstani haplotypes C were thzen categorized into Asian mtDNA type. Sequence analysis of the SRY gene on the Y-chromosome in 67 male Kazakhstani goats revealed two haplotypes: Y1A (64%) and Y2A (36%). Analysis of the distribution of mtDNA haplogroups and SRY haplotypes from Eurasia and Africa demonstrated genetic similarity among animals from Kazakhstan, Mongolia, and Northwest China located on the Eurasian Steppe belt. These phylogeographic results suggested that the Eurasian Steppe belt was an important propagation route for goats to Central Asia.

Genetic diversity and maternal lineages of south Indian goats

In India, goats are considered to be one of the important livestock species that reinforce the rural economy. Even though India has 23 well-recognized goat breeds, the knowledge about their genetic diversity and domestication history is limited. In this study, we have analyzed the genetic diversity of 104 goats representing five different south Indian breeds using mtDNA D-loop region. The haplotype diversity of the breeds ranged from 0.9524 ± 0.0403 (Malabari) to 0.9921 ± 0.0154 (Kanni aadu). Analysis of molecular variance showed only 1.05% variation among breeds. On other hand, the variation within breed was remarkably high (98.95%) which suggested the weak phylogeographic structure of south Indian goats. The phylogenetic analysis revealed three haplogroups representing maternal lineages namely A, B and D. The analysis of 466 Indian goat sequences showed an additional lineage C. As reported in the previous studies, a major fraction of analyzed goats fell into haplogroup A. Our study confirms the presence of three maternal lineages for south Indian domestic goats.

Mitochondrial DNA variation reveals maternal origins and demographic dynamics of Ethiopian indigenous goats

The Horn of Africa forms one of the two main historical entry points of domestics into the continent and Ethiopia is particularly important in this regard. Through the analysis of mitochondrial DNA (mtDNA) d-loop region in 309 individuals from 13 populations, we reveal the maternal genetic variation and demographic dynamics of Ethiopian indigenous goats. A total of 174 variable sites that generated 231 haplotypes were observed. They defined two haplogroups that were present in all the 13 study populations. Reference haplotypes from the six globally defined goat mtDNA haplogroups show the two haplogroups present in Ethiopia to be A and G, the former being the most predominant. Although both haplogroups are characterized by an increase in effective population sizes (Ne) predating domestication, they also have experienced a decline in Ne at different time periods, suggesting different demographic histories. We observed seven haplotypes, six were directly linked to the central haplotypes of the two haplogroups and one was central to haplogroup G. The seven haplotypes were common between Ethiopia, Kenya, Egypt, and Saudi Arabia populations, suggesting common maternal history and the introduction of goats into East Africa via Egypt and the Arabian Peninsula, respectively. While providing new mtDNA data from a historically important region, our results suggest extensive intermixing of goats mediated by human socio-cultural and economic interactions. These have led to the coexistence of the two haplogroups in different geographic regions in Ethiopia resulting in a large caprine genetic diversity that can be exploited for genetic improvement.

Genetic origin of goat populations in Oman revealed by mitochondrial DNA analysis

The Sultanate of Oman has a complex mosaic of livestock species and production systems, but the genetic diversity, demographic history or origins of these Omani animals has not been expensively studied. Goats might constitute one of the most abundant and important domestic livestock species since the Neolithic transition. Here, we examined the genetic diversity, origin, population structure and demographic history of Omani goats. Specifically, we analyzed a 525-bp fragment of the first hypervariable region of the mitochondrial DNA (mtDNA) control region from 69 Omani individuals and compared this fragment with 17 mtDNA sequences from Somalia and Yemen as well as 18 wild goat species and 1,198 previously published goat sequences from neighboring countries. The studied goat breeds show substantial diversity. The haplotype and nucleotide diversities of Omani goats were found equal to 0.983 ± 0.006 and 0.0284 ± 0.014, respectively. The phylogenetic analyses allowed us to classify Omani goats into three mtDNA haplogroups (A, B and G): haplogroup A was found to be predominant and widely distributed and accounted for 80% of all samples, and haplogroups B and G exhibited low frequencies. Phylogenetic comparisons with wild goats revealed that five of the native Omani goat populations originate from Capra aegagrus. Furthermore, most comparisons of pairwise population F_ST values within and between these five Omani goat breeds as well as between Omani goats and nine populations from nearby countries were not significant. These results suggest strong gene flow among goat populations caused by the extensive transport of goats and the frequent movements of human populations in ancient Arabia. The findings improve our understanding of the migration routes of modern goats from their region of domestication into southeastern Arabia and thereby shed light on human migratory and commercial networks during historical times.

Cytochrome b conservation between six camel breeds reared in Egypt

This study was aimed to assess cytochrome b conservation in six breeds of camels reared in Egypt and to compare its sequence with those of other livestock species. The 208-bp fragments from camel mtDNA cyto b were amplified using PCR for 54 camels belonging to 6 camel breeds reared in Egypt. The alignment of camel cyto b sequences showed the presence of two polymorphic sites resulting in four haplotypes and their nucleotide sequences were submitted to GenBank under the accession numbers: KX909894-KX909897. The genetic distances between tested camel breeds were zero between Baladi, Fallahi and Maghrabi breeds whereas they were at low value between the other three breeds: Mowaled, Sodany and Somali. Neighbor-joining showed 4 branches; one of them include most of the tested animals and another one contains 2 Somali animals which is considered a specific haplotype for this breed. The other two branches are mixed between Sodani and Mowaled breeds. Neighbor-joining tree was constructed between cyto b sequences of our tested camels and their sequences from livestock species include Camelus dromedaries, Camelus bactrianus, Ovis aries, Capra hircus, Bubalus bubalis, Bos Taurus and Sus scrofa. The result confirmed that our camel breeds belong to Camelus dromedaries and are clearly separated from other species. It is concluded that cyto b sequence is highly conserved among all camel breeds reared in Egypt which belong to Camelus dromedaries in addition to the advantage of cyto b in differentiation between different livestock sources which enables it to widely use for the adulteration detection in mixed meat.

Polymorphisms, differentiation, and phylogeny of 10 Tibetan goat populations inferred from mitochondrial D-loop sequences

China has abundant population of Tibetan goats, but available information regarding genetic structure and phylogenetic status remains limited. Here, 130 mitochondrial D-loop sequences of individuals from 10 Tibetan goat populations located in distinct plateau areas were analyzed. Eighty-six haplotypes were defined, among which 97.7% were group-specific haplotypes. Haplotype and nucleotide diversity indices were 0.990 ± 0.003 and 0.0145 ± 0.0013, respectively. The pairwise Wright's F-statistics ranged from -0.028 to 0.385, and over half of them were greater than 0.05, indicating apparent genetic differentiation among the populations. AMOVA analysis (F_ST = 0.0858) manifested that the genetic structure has become weak. Phylogenetic trees revealed four haplogroups (A, B, C, and D), suggesting that Tibetan goats had four origins. Mismatch distribution analyses and neutrality tests indicated that at least one population expansion event occurred during the demographic history of Tibetan goat. These results will provide a more complete understanding of Tibetan goat genetic resources.

Genetic Variation of Goat Interferon Regulatory Factor 3 Gene and Its Implication in Goat Evolution

The immune systems are fundamentally vital for evolution and survival of species; as such, selection patterns in innate immune loci are of special interest in molecular evolutionary research. The interferon regulatory factor (IRF) gene family control many different aspects of the innate and adaptive immune responses in vertebrates. Among these, IRF3 is known to take active part in very many biological processes. We assembled and evaluated 1356 base pairs of the IRF3 gene coding region in domesticated goats from Africa (Nigeria, Ethiopia and South Africa) and Asia (Iran and China) and the wild goat (Capra aegagrus). Five segregating sites with θ value of 0.0009 for this gene demonstrated a low diversity across the goats’ populations. Fu and Li tests were significantly positive but Tajima’s D test was significantly negative, suggesting its deviation from neutrality. Neighbor joining tree of IRF3 gene in domesticated goats, wild goat and sheep showed that all domesticated goats have a closer relationship than with the wild goat and sheep. Maximum likelihood tree of the gene showed that different domesticated goats share a common ancestor and suggest single origin. Four unique haplotypes were observed across all the sequences, of which, one was particularly common to African goats (MOCH-K14-0425, Poitou and WAD). In assessing the evolution mode of the gene, we found that the codon model d_N/d_S ratio for all goats was greater than one. Phylogenetic Analysis by Maximum Likelihood (PAML) gave a ω₀ (d_N/d_S) value of 0.067 with LnL value of -6900.3 for the first Model (M1) while ω2 = 1.667 in model M2 with LnL value of -6900.3 with positive selection inferred in 3 codon sites. Mechanistic empirical combination (MEC) model for evaluating adaptive selection pressure on particular codons also confirmed adaptive selection pressure in three codons (207, 358 and 408) in IRF3 gene. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat IRF3 led us to conclude that the gene evolution may have been influenced by domestication processes in goats.

Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability

Background

The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA).

Results

We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12–10 ka ago.

Conclusions

Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2342-2) contains supplementary material, which is available to authorized users.

Toll-like receptor responses to Peste des petits ruminants virus in goats and water buffalo

Ovine rinderpest or goat plague is an economically important and contagious viral disease of sheep and goats, caused by the Peste des petits ruminants virus (PPRV). Differences in susceptibility to goat plague among different breeds and water buffalo exist. The host innate immune system discriminates between pathogen associated molecular patterns and self antigens through surveillance receptors known as Toll like receptors (TLR). We investigated the role of TLR and cytokines in differential susceptibility of goat breeds and water buffalo to PPRV. We examined the replication of PPRV in peripheral blood mononuclear cells (PBMC) of Indian domestic goats and water buffalo and demonstrated that the levels of TLR3 and TLR7 and downstream signalling molecules correlation with susceptibility vs resistance. Naturally susceptible goat breeds, Barbari and Tellichery, had dampened innate immune responses to PPRV and increased viral loads with lower basal expression levels of TLR 3/7. Upon stimulation of PBMC with synthetic TLR3 and TLR7 agonists or PPRV, the levels of proinflammatory cytokines were found to be significantly higher while immunosuppressive interleukin (IL) 10 levels were lower in PPRV resistant Kanni and Salem Black breeds and water buffalo at transcriptional level, correlating with reduced viralloads in infected PBMC. Water buffalo produced higher levels of interferon (IFN) α in comparison with goats at transcriptional and translational levels. Pre-treatment of Vero cells with human IFNα resulted in reduction of PPRV replication, confirming the role of IFNα in limiting PPRV replication. Treatment with IRS66, a TLR7 antagonist, resulted in the reduction of IFNα levels, with increased PPRV replication confirming the role of TLR7. Single nucleotide polymorphism analysis of TLR7 of these goat breeds did not show any marked nucleotide differences that might account for susceptibility vs resistance to PPRV. Analyzing other host genetic factors might provide further insights on susceptibility to PPRV and genetic polymorphisms in the host.

Mitochondrial DNA diversity of Anatolian indigenous domestic goats

Anatolia has been an important region for civilizations and agricultural revolution as a major domestication centre for livestock species. Goats (Capra hircus) were among the earliest domesticated animals in this region. In this study, genetic diversity of Anatolian goat breeds was characterized by comparison of mitochondrial DNA hypervariable region 1. A total of 295 individuals, including 99 Anatolian Black goats, 96 Angora goats and 100 Kilis goats, were used. Haplogroup A was found to be the dominant haplogroup in all three breeds. The highest haplogroup diversity, including haplogroups A, B2, C and G, was observed in the Anatolian Black breed. Haplogroup D was only observed in Kilis and Angora goats. Haplogroup G was found in Angora and Anatolian Black breeds. The Anatolian goat breeds had high genetic diversity values and a weak phylogeographical structure. The nucleotide diversity values were found to be higher than those in previously studied goat breeds. The fact that Anatolia is a domestication centre and its geographical position as a junction of trade routes may have caused the higher genetic diversity of Anatolian goat breeds.

Retrospective and prospective perspectives on zoonotic brucellosis

Members of the genus Brucella are pathogenic bacteria exceedingly well adapted to their hosts. The bacterium is transmitted by direct contact within the same host species or accidentally to secondary hosts, such as humans. Human brucellosis is strongly linked to the management of domesticated animals and ingestion of their products. Since the domestication of ungulates and dogs in the Fertile Crescent and Asia in 12000 and 33000 ya, respectively, a steady supply of well adapted emergent Brucella pathogens causing zoonotic disease has been provided. Likewise, anthropogenic modification of wild life may have also impacted host susceptibility and Brucella selection. Domestication and human influence on wild life animals are not neutral phenomena. Consequently, Brucella organisms have followed their hosts’ fate and have been selected under conditions that favor high transmission rate. The “arm race” between Brucella and their preferred hosts has been driven by genetic adaptation of the bacterium confronted with the evolving immune defenses of the host. Management conditions, such as clustering, selection, culling, and vaccination of Brucella preferred hosts have profound influences in the outcome of brucellosis and in the selection of Brucella organisms. Countries that have controlled brucellosis systematically used reliable smooth live vaccines, consistent immunization protocols, adequate diagnostic tests, broad vaccination coverage and sustained removal of the infected animals. To ignore and misuse tools and strategies already available for the control of brucellosis may promote the emergence of new Brucella variants. The unrestricted use of low-efficacy vaccines may promote a “false sense of security” and works towards selection of Brucella with higher virulence and transmission potential.

Genetic structure of Chinese indigenous goats and the special geographical structure in the Southwest China as a geographic barrier driving the fragmentation of a large population

Background

China has numerous native domestic goat breeds, however, extensive studies are focused on the genetic diversity within the fewer breeds and limited regions, the population demograogic history and origin of Chinese goats are still unclear. The roles of geographical structure have not been analyzed in Chinese goat domestic process. In this study, the genetic relationships of Chinese indigenous goat populations were evaluated using 30 microsatellite markers.

Methodology/Principal Findings

Forty Chinese indigenous populations containing 2078 goats were sampled from different geographic regions of China. Moderate genetic diversity at the population level (H_S of 0.644) and high population diversity at the species level (H_T value of 0.737) were estimated. Significant moderate population differentiation was detected (F_ST value of 0.129). Significant excess homozygosity (F_IS of 0.105) and recent population bottlenecks were detected in thirty-six populations. Neighbour-joining tree, principal components analysis and Bayesian clusters all revealed that Chinese goat populations could be subdivided into at least four genetic clusters: Southwest China, South China, Northwest China and East China. It was observed that the genetic diversity of Northern China goats was highest among these clusters. The results here suggested that the goat populations in Southwest China might be the earliest domestic goats in China.

Conclusions/Significance

Our results suggested that the current genetic structure of Chinese goats were resulted from the special geographical structure, especially in the Western China, and the Western goat populations had been separated by the geographic structure (Hengduan Mountains and Qinling Mountains-Huaihe River Line) into two clusters: the Southwest and Northwest. It also indicated that the current genetic structure was caused by the geographical origin mainly, in close accordance with the human’s migration history throughout China. This study provides a fundamental genetic profile for the conservation of these populations and better to understand the domestication process and origin of Chinese goats.

Extensive variation and sub-structuring in lineage A mtDNA in Indian sheep: genetic evidence for domestication of sheep in India

Previous studies on mitochondrial DNA analysis of sheep from different regions of the world have revealed the presence of two major- A and B, and three minor- C, D and E maternal lineages. Lineage A is more frequent in Asia and lineage B is more abundant in regions other than Asia. We have analyzed mitochondrial DNA sequences of 330 sheep from 12 different breeds of India. Neighbor-joining analysis revealed lineage A, B and C in Indian sheep. Surprisingly, multidimensional scaling plot based on F_ST values of control region of mtDNA sequences showed significant breed differentiation in contrast to poor geographical structuring reported earlier in this species. The breed differentiation in Indian sheep was essentially due to variable contribution of two major lineages to different breeds, and sub- structuring of lineage A, possibly the latter resulting from genetic drift. Nucleotide diversity of this lineage was higher in Indian sheep (0.014 ± 0.007) as compared to that of sheep from other regions of the world (0.009 ± 0.005 to 0.01 ± 0.005). Reduced median network analysis of control region and cytochrome b gene sequences of Indian sheep when analyzed along with available published sequences of sheep from other regions of the world showed that several haplotypes of lineage A were exclusive to Indian sheep. Given the high nucleotide diversity in Indian sheep and the poor sharing of lineage A haplotypes between Indian and non-Indian sheep, we propose that lineage A sheep has also been domesticated in the east of Near East, possibly in Indian sub-continent. Finally, our data provide support that lineage B and additional lineage A haplotypes of sheep might have been introduced to Indian sub-continent from Near East, probably by ancient sea trade route.

Domestication process of the goat revealed by an analysis of the nearly complete mitochondrial protein-encoding genes

Goats (Capra hircus) are one of the oldest domesticated species, and they are kept all over the world as an essential resource for meat, milk, and fiber. Although recent archeological and molecular biological studies suggested that they originated in West Asia, their domestication processes such as the timing of population expansion and the dynamics of their selection pressures are little known. With the aim of addressing these issues, the nearly complete mitochondrial protein-encoding genes were determined from East, Southeast, and South Asian populations. Our coalescent time estimations suggest that the timing of their major population expansions was in the Late Pleistocene and significantly predates the beginning of their domestication in the Neolithic era (≈10,000 years ago). The ω (ratio of non-synonymous rate/synonymous substitution rate) for each lineage was also estimated. We found that the ω of the globally distributed haplogroup A which is inherited by more than 90% of goats examined, turned out to be extremely low, suggesting that they are under severe selection pressure probably due to their large population size. Conversely, the ω of the Asian-specific haplogroup B inherited by about 5% of goats was relatively high. Although recent molecular studies suggest that domestication of animals may tend to relax selective constraints, the opposite pattern observed in our goat mitochondrial genome data indicates the process of domestication is more complex than may be presently appreciated and cannot be explained only by a simple relaxation model.

Caprine Toll-like receptor 8 gene sequence characterization reveals close relationships among ruminant species

TLR8 mediates antiviral immunity by recognizing ssRNA viruses and triggers potent antiviral and antitumor immune responses. In this study, approximately 3.5 Kb nucleotide sequence data of caprine TLR8 gene were generated from one sample each of twelve different Indian goat breeds belonging to different geographical regions. Cloning and characterization of cDNA synthesized from RNA purified from goat spleen revealed TLR8 ORF to be of 3102 nucleotides long coding for 1033 amino acids similar to other ruminant species, that is sheep, buffalo and cattle. The sequence analysis at nucleotide level revealed goat TLR8 to be closer to buffalo sharing 99.6% homology, followed by cattle and sheep. Simple Modular Architecture Research Tool (SMART) used for the structural analysis of goat TLR8 showed the presence of 16 leucine-rich repeats (LRRs) along with single Toll/interleukin-1 receptor (TIR) domain. TIR domain when compared with other livestock species was found to be conserved in ruminant species goat, sheep, cattle and buffalo. The phylogenetic analysis also revealed grouping of all ruminant species together, goat being closer to buffalo followed by cattle and sheep. Total 4 polymorphic sites were observed in TLR8 gene of one specimen goat representing each of 12 different breeds studied, all of which were synonymous and present within the coding region. Of these 4 SNPs, two were in ectodomains, one in TIR domain and one was found to be present in transmembrane domain. PCR-RFLP genotyping of two of the SNPs indicated variations in allele frequencies among different goat breeds. The expression profiling in 13 tissues of goat showed maximum expression of TLR8 gene in kidney followed by spleen, lung and lymph node. Overall, our results indicate conservation of TLR8 gene among the ruminant species and low variation within Indian goat breeds.

Mammalian cortical bone in tension is non-Haversian

Cortical bone, found in the central part of long bones like femur, is known to adapt to local mechanical stresses. This adaptation has been linked exclusively with Haversian remodelling involving bone resorption and formation of secondary osteons. Compared to primary/plexiform bone, the Haversian bone has lower stiffness, fatigue strength and fracture toughness, raising the question why nature prefers an adaptation that is detrimental to bone's primary function of bearing mechanical stresses. Here, we show that in the goat femur, Haversian remodelling occurs only at locations of high compressive stresses. At locations corresponding to high tensile stresses, we observe a microstructure that is non-Haversian. Compared with primary/plexiform bone, this microstructure's mineralisation is significantly higher with a distinctly different spatial pattern. Thus, the Haversian structure is an adaptation only to high compressive stresses rendering its inferior tensile properties irrelevant as the regions with high tensile stresses have a non-Haversian, apparently primary microstructure.