Phylogenetic analysis of the tribe Bovini using microsatellites

Characterization of olfactory receptor repertoires provides insights into the high-altitude adaptation of the yak based on the chromosome-level genome

Olfaction in vertebrates plays pivotal parts in many aspects, such as localizing prey or food, mating behavior, avoiding predators, and social communication. Yak (Bos grunniens) is the only Bos species that can thrive in high-altitude areas. In view of the critical role of olfactory receptors (ORs) in the specific recognition of diverse stimuli, investigating the evolutionary dynamics of ORs in the yak means a lot. In this study, we used the chromosome-level genome of the yak to identify the ORs genes and discussed the effects of high altitude on the yak's olfaction by comparing the yak with other low-altitude living Bos species (Bos frontalis (gayal), Bos gaurus (gaur), Bos indicus (zebu) and Bos taurus (cattle)). The yak had 400 OR genes, including 264 functional genes, 16 partial genes and 120 OR pseudo genes. There were 387 OR genes mapped to yak 31 chromosomes, and chromosomes 13 and 8 had the most OR genes and functional OR genes. Among these five Bos species, yak had the least number of OR gene subfamilies, OR genes and functional OR genes, while the total number of OR genes in gayal (n = 784) was almost twice as many as that of yak, indicating that the olfaction of yak may be less developed. In addition, the phylogenetic relationships of the functional Bos OR genes were illustrated, which comprised 79 families and 466 subfamilies distributed in two classes (Class I and Class II). There were 76 OR gene subfamilies shared by these five Bos species and 17 OR gene subfamilies were unique to the yak. The potential odor specificity of 44 yak OR genes was identified through the similarity to human OR protein sequences. Remarkably, yak lacks β-ionone and Isovaleric acid(IVA)-related ORs, which may be related to the decline of high-altitude herbaceous plant diversity and underdeveloped yak sweat glands. The conserved motifs of OR genes were highly conserved in Bos species. These results provided a solid foundation for further studies on the molecular mechanisms of the yak's adaptation to the high-altitude environment in olfaction.

Evaluation of Genetic Diversity and Structure of Turkish Water Buffalo Population by Using 20 Microsatellite Markers

The present study was aimed to investigate the genetic diversity among 17 Turkish water buffalo populations. A total of 837 individuals from 17 provincial populations were genotyped, using 20 microsatellites markers. The microsatellite markers analyzed were highly polymorphic with a mean number of alleles of (7.28) ranging from 6 (ILSTS005) to 17 (ETH003). The mean observed and expected heterozygosity values across all polymorphic loci in all studied buffalo populations were 0.61 and 0.70, respectively. Observed heterozygosity varied from 0.55 (Bursa (BUR)) to 0.70 (Muş (MUS)). It was lower than expected heterozygosity in most of the populations indicating a deviation from Hardy-Weinberg equilibrium. The overall value for the polymorphic information content of noted microsatellite loci was 0.655, indicating their suitability for genetic diversity analysis in buffalo. The mean F_IS value was 0.091 and all loci were observed significantly deviated from Hardy-Weinberg Equilibrium (HWE), most likely based on non-random breeding. The 17 buffalo populations were genetically less diverse as indicated by a small mean F_ST value (0.032 ± 0.018). The analysis of molecular variance (AMOVA) analysis indicated that about 2% of the total genetic diversity was clarified by population distinctions and 88 percent corresponded to differences among individuals. The information produced by this study can be used to establish a base of national conservation and breeding strategy of water buffalo population in Turkey.

Distribution of the variant at AKIRIN2: c.*188G > A in Chinese cattle

Marbling score (MS), is an economically important trait in cattle. Previous results showed that a SNP (c.*188G > A) of akirin 2 (AKIRIN2) gene was associated with MS in Japanese Black cattle and Korean cattle. However, the distribution of the genotypic frequency of the single nucleotide polymorphism (SNP) has not been explored in Chinese cattle. In this study, we used polymerase chain reaction (PCR) and DNA sequencing to detect the variation in 1296 individuals from 39 Chinese cattle breeds, one semi-wild bovine species (Dulong) and three introduced breeds (Angus, Holstein and Brahman). Our study found the frequency of the A allele at this locus roughly diminished from north to south in Chinese cattle, and we detected statistically significant differences between Angus and Brahman (p < 0.05), Dulong and another two breeds (Angus and Holstein; p < 0.01) using Chi-Square Independence Test. Our results reflected the variation of AKIRIN2: c.*188G > A in Chinese cattle, which would help us better understand Chinese cattle genetic resources and provide reference for further research.

The complete mitochondrial genome of Indian gaur, Bos gaurus and its phylogenetic implications

The gaur is the largest extant cattle species and distributed across South and Southeast Asia. Around 85% of its current global population resides in India, however there has been a gradual decrease in the gaur population over the last two decades due to various anthropogenic activities. Mitochondrial genome is considered as an important tool for species identification and monitoring the populations of conservation concern and therefore it becomes an obligation to sequence the mitochondrial genome of Indian gaur. We report here for the first time 16,345 bp mitochondrial genome of four Indian gaur sequenced using two different approaches. Mitochondrial genome consisted of 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a control region. Among the 37 genes, 28 were positioned on the H-strand and 9 were positioned on the L-strand. The overall base composition appeared to be 33.5% A, 27.2% T, 25.9% C and 13.4% G, which yielded a higher AT content. The phylogenetic analysis using complete mitochondrial genome sequences unambiguously suggested that gaur is the maternal ancestor of domestic mithun. Moreover, it also clearly distinguished the three sub species of B. gaurus i.e. B. gaurus gaurus, B. gaurus readei and B. gaurus hubbacki. Among the three sub species, B. gaurus gaurus was genetically closer to B. gaurus readei as compared to B. gaurus hubbacki. The findings of our study provide an insight into the genetic structure and evolutionary history of Indian gaur.

Cattle microsatellite markers successfully established diversity status of Arunachali yak (only registered yak breed of India)

Yak diversity of the country has remained predominantly unexplored for a long time. Among the 169 registered livestock breeds of India, the sole representation from yak genetic resources is the Arunachali yak. This study for the first time investigated genetic diversity status of Arunachali yak using 26 bovine microsatellite markers. All the markers recommended for cattle except one (ILSTS05) amplified with yak genome. Allelic genotype pattern overlapped between yak and cattle across 25 microsatellite loci and a total of 233 alleles were detected in yak. The number of observed alleles across loci ranged from 3–16 with an average of 9.32±0.70. Observed heterozygosity (0.552±0.04) was less than the expected heterozygosity (0.648±0.035) pointing towards heterozygote deficiency in the population. In addition, positive value of FIS index (0.143±0.043) suggested considerable inbreeding. There was no indication of a recent bottleneck event in this population based on heterozygosity excess tests as well as mode-shift analysis. In summary, bovine microsatellite markers proved to be a valuable tool for characterization of Indian yak population. Arunachali yak represents an interesting gene pool with moderate level of diversity. Inbreeding in population calls for sincere efforts to formulate breeding policy so that this precious germplasm is conserved with substantial genetic diversity.

Effect of complement component 5 polymorphisms on mastitis resistance in Egyptian buffalo and cattle

Mastitis is one of the costliest diseases affecting the world's dairy industry. The important contribution of complement Component 5 (C5) to phagocytosis, which plays a major role in the defence of the bovine mammary gland against infection, makes this component of innate immunity a potential contributor in defending udder against mastitis. The objectives of this study were to sequence and analyse the whole coding region of the C5 gene in Egyptian buffalo and cattle, to detect any nucleotide variations (polymorphisms) and to investigate their associations with milk somatic cell score (SCS) as an indicator of mastitis in dairy animals. We sequenced a buffalo C5 cDNA fragment of 5336 bp (KP221293) and a cattle C5 cDNA fragment of 5303 bp (KP221294), which included the whole coding region and 3-UTR. Buffalo and cattle C5 cDNA shared sequence identity of 99%. The predicted complement C5 proteins consist of 1677 amino acid residues in both animals, one amino acid less than in humans and three amino acids more than in mouse C5 protein. Comparing cDNA sequences of different animals revealed nine novel SNPs in buffalo and seven SNPs in cattle, with two of them being novel. The association analysis revealed that five SNPs in buffalo are highly associated with SCS; indicating the contribution of complement C5 variants in buffalo mastitis resistance. No significant associations were detected between C5 variants and SCS in cattle. This is the first report about C5 variants in buffalo and its association with SCS.

New Insights on Water Buffalo Genomic Diversity and Post-Domestication Migration Routes From Medium Density SNP Chip Data

The domestic water buffalo is native to the Asian continent but through historical migrations and recent importations, nowadays has a worldwide distribution. The two types of water buffalo, i.e., river and swamp, display distinct morphological and behavioral traits, different karyotypes and also have different purposes and geographical distributions. River buffaloes from Pakistan, Iran, Turkey, Egypt, Romania, Bulgaria, Italy, Mozambique, Brazil and Colombia, and swamp buffaloes from China, Thailand, Philippines, Indonesia and Brazil were genotyped with a species-specific medium-density 90K SNP panel. We estimated the levels of molecular diversity and described population structure, which revealed historical relationships between populations and migration events. Three distinct gene pools were identified in pure river as well as in pure swamp buffalo populations. Genomic admixture was seen in the Philippines and in Brazil, resulting from importations of animals for breed improvement. Our results were largely consistent with previous archeological, historical and molecular-based evidence for two independent domestication events for river- and swamp-type buffaloes, which occurred in the Indo-Pakistani region and close to the China/Indochina border, respectively. Based on a geographical analysis of the distribution of diversity, our evidence also indicated that the water buffalo spread out of the domestication centers followed two major divergent migration directions: river buffaloes migrated west from the Indian sub-continent while swamp buffaloes migrated from northern Indochina via an east-south-eastern route. These data suggest that the current distribution of water buffalo diversity has been shaped by the combined effects of multiple migration events occurred at different stages of the post-domestication history of the species.

Genetic variation in N- and C-terminal regions of bovine DNAJA1 heat shock protein gene in African, Asian and American cattle

DNAJA1 or heat shock protein 40 (Hsp40) is associated with heat adaptation in various organisms. We amplified and sequenced a total of 1,142 bp of bovine Hsp40 gene representing the critical N-terminal (NTR) and C-terminal (CTR) regions in representative samples of African, Asian and American cattle breeds. Eleven and 9 different haplotypes were observed in the NTR in Asian and African breeds respectively while in American Brangus, only two mutations were observed resulting in two haplotypes. The CTR appears to be highly conserved between cattle and yak. In-silico functional analysis with PANTHER predicted putative deleterious functional impact of c.161 T>A; p. V54Q while alignment of bovine and human NTR-J domains revealed that p.Q19H, p.E20Q and p. E21X mutations occurred in helix 2 and p.V54Q missense mutation occurred in helix 3 respectively. The 124 bp insertion found in the yak DNAJA1 ortholog may have significant functional relevance warranting further investigation. Our results suggest that these genetic differences may be concomitant with population genetic history and possible functional consequences for climate adaptation in bovidae.

Genome-wide characterization of perfect microsatellites in yak (Bos grunniens)

Microsatellites or simple sequence repeats (SSRs) constitute a significant portion of genomes and play an important role in gene function and genome organization. The availability of a complete genome sequence for yak (Bos grunniens) has made it possible to carry out genome-wide analysis of microsatellites in this species. We analyzed the abundance and density of perfect SSRs in the yak genome. We found a total of 723,172 SSRs with 1-6 bp nucleotide motifs, indicating that about 0.47 % of the yak whole genome sequence (2.66 Gb) comprises perfect SSRs, the average length of which was 17.34 bp/Mb. The average frequency and density of perfect SSRs was 272.18 loci/Mb and 4719.25 bp/Mb, respectively. The proportion of the six classes of perfect SSRs was not evenly distributed in the yak genome. Mononucleotide repeats (44.04 %) with a total number of 318,435 and a average length of 14.71 bp appeared to be the most abundant SSRs class, while the percentages of dinucleotide, trinucleotide, pentanucleotide, tetranucleotide and hexanucleotide repeats was 24.11 %, 15.80 %, 9.50 %, 6.40 % and 0.15 %, respectively. Different repeat classes of SSRs varied in their repeat number with the highest being 1206. Our results suggest that 15 motifs comprised the predominant categories with a frequency above 1 loci/Mb: A, AC, AT, AG, AGC, AAC, AAT, ACC, ATTT, GTTT, AATG, CTTT, ATGG, AACTG and ATCTG.

Revisiting AFLP fingerprinting for an unbiased assessment of genetic structure and differentiation of taurine and zebu cattle

Background

Descendants from the extinct aurochs (Bos primigenius), taurine (Bos taurus) and zebu cattle (Bos indicus) were domesticated 10,000 years ago in Southwestern and Southern Asia, respectively, and colonized the world undergoing complex events of admixture and selection. Molecular data, in particular genome-wide single nucleotide polymorphism (SNP) markers, can complement historic and archaeological records to elucidate these past events. However, SNP ascertainment in cattle has been optimized for taurine breeds, imposing limitations to the study of diversity in zebu cattle. As amplified fragment length polymorphism (AFLP) markers are discovered and genotyped as the samples are assayed, this type of marker is free of ascertainment bias. In order to obtain unbiased assessments of genetic differentiation and structure in taurine and zebu cattle, we analyzed a dataset of 135 AFLP markers in 1,593 samples from 13 zebu and 58 taurine breeds, representing nine continental areas.

Results

We found a geographical pattern of expected heterozygosity in European taurine breeds decreasing with the distance from the domestication centre, arguing against a large-scale introgression from European or African aurochs. Zebu cattle were found to be at least as diverse as taurine cattle. Western African zebu cattle were found to have diverged more from Indian zebu than South American zebu. Model-based clustering and ancestry informative markers analyses suggested that this is due to taurine introgression. Although a large part of South American zebu cattle also descend from taurine cows, we did not detect significant levels of taurine ancestry in these breeds, probably because of systematic backcrossing with zebu bulls. Furthermore, limited zebu introgression was found in Podolian taurine breeds in Italy.

Conclusions

The assessment of cattle diversity reported here contributes an unbiased global view to genetic differentiation and structure of taurine and zebu cattle populations, which is essential for an effective conservation of the bovine genetic resources.

Molecular characterization and expression profile of ghrelin gene during different reproductive phases in buffalo (Bubalus bubalis)

Ghrelin, a novel motilin-related endogenous ligand for growth hormone secretagouge receptor, is implicated in various biological functions, including regulation of female reproduction. But the presence of ghrelin and its role in reproductive functions in buffalo, a species with poor reproductive efficiency, is not known. In the present study full-length ghrelin cDNA was isolated from bubaline abomasum, which encodes the entire prepropeptide of 116 amino acids. The deduced amino acid sequence of ghrelin of buffalo showed >95% and 31% identity with that of ruminants (cattle, sheep, and goat) and humans, respectively. Analysis of synonymous and nonsynonymous nucleotide substitutions in the coding region of ghrelin indicated that these sequences of different species have been under purifying selection. The 3995-bp amplicon of ghrelin gene consisting of 4 exons and 3 introns was cloned with genomic DNA from buffalo. Further, ghrelin expression was determined by quantitative real-time PCR, in situ hybridization, and immunohistochemistry in bubaline endometrial tissues at different stages of the estrous cycle and early pregnancy. Our results indicated the persistent expression of ghrelin mRNA and protein in the endometrium during stage I (day 3-5), stage II (day 6-15), and stage III (day 16-21) of the estrous cycle and also during early (~day 30-40) pregnancy. Immunohistochemistry and quantitative real-time PCR experiments indicated the relatively higher expression of ghrelin in the endometrium during stage II (day 6-15) of the estrous cycle and early pregnancy than during stage I (day 3-5) and stage III (day 16-21) of the estrous cycle, but no statistically significant difference in ghrelin expression was observed among stages. To conclude, the results of the present study indicate the persistent expression of ghrelin in the uterine endometrium throughout the estrous cycle and in early pregnancy which might be helpful in determining its role in buffalo reproduction.

Inferring demography from runs of homozygosity in whole-genome sequence, with correction for sequence errors

Whole-genome sequence is potentially the richest source of genetic data for inferring ancestral demography. However, full sequence also presents significant challenges to fully utilize such large data sets and to ensure that sequencing errors do not introduce bias into the inferred demography. Using whole-genome sequence data from two Holstein cattle, we demonstrate a new method to correct for bias caused by hidden errors and then infer stepwise changes in ancestral demography up to present. There was a strong upward bias in estimates of recent effective population size (Ne) if the correction method was not applied to the data, both for our method and the Li and Durbin (Inference of human population history from individual whole-genome sequences. Nature 475:493-496) pairwise sequentially Markovian coalescent method. To infer demography, we use an analytical predictor of multiloci linkage disequilibrium (LD) based on a simple coalescent model that allows for changes in Ne. The LD statistic summarizes the distribution of runs of homozygosity for any given demography. We infer a best fit demography as one that predicts a match with the observed distribution of runs of homozygosity in the corrected sequence data. We use multiloci LD because it potentially holds more information about ancestral demography than pairwise LD. The inferred demography indicates a strong reduction in the Ne around 170,000 years ago, possibly related to the divergence of African and European Bos taurus cattle. This is followed by a further reduction coinciding with the period of cattle domestication, with Ne of between 3,500 and 6,000. The most recent reduction of Ne to approximately 100 in the Holstein breed agrees well with estimates from pedigrees. Our approach can be applied to whole-genome sequence from any diploid species and can be scaled up to use sequence from multiple individuals.

Splice variants and promoter methylation status of the Bovine Vasa Homology (Bvh) gene may be involved in bull spermatogenesis

Background

Vasa is a member of the DEAD-box protein family that plays an indispensable role in mammalian spermatogenesis, particularly during meiosis. Bovine vasa homology (Bvh) of Bos taurus has been reported, however, its function in bovine testicular tissue remains obscure. This study aimed to reveal the functions of Bvh and to determine whether Bvh is a candidate gene in the regulation of spermatogenesis in bovine, and to illustrate whether its transcription is regulated by alternative splicing and DNA methylation.

Results

Here we report the molecular characterization, alternative splicing pattern, expression and promoter methylation status of Bvh. The full-length coding region of Bvh was 2190 bp, which encodes a 729 amino acid (aa) protein containing nine consensus regions of the DEAD box protein family. Bvh is expressed only in the ovary and testis of adult cattle. Two splice variants were identified and termed Bvh-V4 (2112 bp and 703 aa) and Bvh-V45 (2040 bp and 679 aa). In male cattle, full-length Bvh (Bvh-FL), Bvh-V4 and Bvh-V45 are exclusively expressed in the testes in the ratio of 2.2:1.6:1, respectively. Real-time PCR revealed significantly reduced mRNA expression of Bvh-FL, Bvh-V4 and Bvh-V45 in testes of cattle-yak hybrids, with meiotic arrest compared with cattle and yaks with normal spermatogenesis (P < 0.01). The promoter methylation level of Bvh in the testes of cattle-yak hybrids was significantly greater than in cattle and yaks (P < 0.01).

Conclusion

In the present study, Bvh was isolated and characterized. These data suggest that Bvh functions in bovine spermatogenesis, and that transcription of the gene in testes were regulated by alternative splice and promoter methylation.

Genetic variation and differentiation of bison (Bison bison) subspecies and cattle (Bos taurus) breeds and subspecies

The genetic relationship of American plains bison (Bison bison bison) and wood bison (Bison bison athabascae) was quantified and compared with that among breeds and subspecies of cattle. Plains bison from 9 herds (N = 136), wood bison from 3 herds (N = 65), taurine cattle (Bos taurus taurus) from 14 breeds (N = 244), and indicine cattle (Bos taurus indicus) from 2 breeds (N = 53) were genotyped for 29 polymorphic microsatellite loci. Bayesian cluster analyses indicate 3 groups, 2 of which are plains bison and 1 of which is wood bison with some admixture, and genetic distances do not show plains bison and wood bison as distinct groups. Differentiation of wood bison and plains bison is also significantly less than that of cattle breeds and subspecies. These and other genetic data and historical interbreeding of bison do not support recognition of extant plains bison and wood bison as phylogenetically distinct subspecies.

Genetic structure and diversity of coffee (Coffea) across Africa and the Indian Ocean islands revealed using microsatellites

BACKGROUND AND AIMS

The coffee genus (Coffea) comprises 124 species, and is indigenous to the Old World Tropics. Due to its immense economic importance, Coffea has been the focus of numerous genetic diversity studies, but despite this effort it remains insufficiently studied. In this study the genetic diversity and genetic structure of Coffea across Africa and the Indian Ocean islands is investigated.

METHODS

Genetic data were produced using 13 polymorphic nuclear microsatellite markers (simple sequence repeats, SSRs), including seven expressed sequence tag-SSRs, and the data were analysed using model- and non-model-based methods. The study includes a total of 728 individuals from 60 species.

KEY RESULTS

Across Africa and the Indian Ocean islands Coffea comprises a closely related group of species with an overall pattern of genotypes running from west to east. Genetic structure was identified in accordance with pre-determined geographical regions and phylogenetic groups. There is a good relationship between morpho-taxonomic species delimitations and genetic units. Genetic diversity in African and Indian Ocean Coffea is high in terms of number of alleles detected, and Madagascar appears to represent a place of significant diversification in terms of allelic richness and species diversity.

CONCLUSIONS

Cross-species SSR transferability in African and Indian Ocean islands Coffea was very efficient. On the basis of the number of private alleles, diversification in East Africa and the Indian Ocean islands appears to be more recent than in West and West-Central Africa, although this general trend is complicated in Africa by the position of species belonging to lineages connecting the main geographical regions. The general pattern of phylogeography is not in agreement with an overall east to west (Mascarene, Madagascar, East Africa, West Africa) increase in genome size, the high proportion of shared alleles between the four regions or the high numbers of exclusive shared alleles between pairs or triplets of regions.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000 ), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000 ) was aligned with the cattle MHC RH map (cR12000 ) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.

Understanding and predicting complex traits: knowledge from cattle

The genetic architecture of complex traits in cattle includes very large numbers of loci affecting any given trait. Most of these loci have small effects but occasionally there are loci with moderate-to-large effects segregating due to recent selection for the mutant allele. Genomic markers capture most but not all of the additive genetic variance for traits, probably because there are causal mutations with low allele frequency and therefore in incomplete linkage disequilibrium with the markers. The prediction of genetic value from genomic markers can achieve high accuracy by using statistical models that include all markers and assuming that marker effects are random variables drawn from a specified prior distribution. Recent effective population size is in the order of 100 within cattle breeds and ≈ 2500 animals with genotypes and phenotypes are sufficient to predict the genetic value of animals with an accuracy of 0.65. Recent effective population size for humans is much larger, in the order of 10,000-15,000, and more than 145,000 records would be required to reach a similar accuracy for people. However, our calculations assume that genomic markers capture all the genetic variance. This may be possible in the future as causal polymorphisms are genotyped using genome sequence data.

Evolution of endogenous sequences of banana streak virus: what can we learn from banana (Musa sp.) evolution?

Endogenous plant pararetroviruses (EPRVs) are viral sequences of the family Caulimoviridae integrated into the nuclear genome of numerous plant species. The ability of some endogenous sequences of Banana streak viruses (eBSVs) in the genome of banana (Musa sp.) to induce infections just like the virus itself was recently demonstrated (P. Gayral et al., J. Virol. 83:6697-6710, 2008). Although eBSVs probably arose from accidental events, infectious eBSVs constitute an extreme case of parasitism, as well as a newly described strategy for vertical virus transmission in plants. We investigated the early evolutionary stages of infectious eBSV for two distinct BSV species-GF (BSGFV) and Imové (BSImV)-through the study of their distribution, insertion polymorphism, and structure evolution among selected banana genotypes representative of the diversity of 60 wild Musa species and genotypes. To do so, the historical frame of host evolution was analyzed by inferring banana phylogeny from two chloroplast regions-matK and trnL-trnF-as well as from the nuclear genome, using 19 microsatellite loci. We demonstrated that both BSV species integrated recently in banana evolution, circa 640,000 years ago. The two infectious eBSVs were subjected to different selective pressures and showed distinct levels of rearrangement within their final structure. In addition, the molecular phylogenies of integrated and nonintegrated BSVs enabled us to establish the phylogenetic origins of eBSGFV and eBSImV.

Water buffalo genome science comes of age

The water buffalo is vital to the lives of small farmers and to the economy of many countries worldwide. Not only are they draught animals, but they are also a source of meat, horns, skin and particularly the rich and precious milk that may be converted to creams, butter, yogurt and many cheeses. Genome analysis of water buffalo has advanced significantly in recent years. This review focuses on currently available genome resources in water buffalo in terms of cytogenetic characterization, whole genome mapping and next generation sequencing. No doubt, these resources indicate that genome science comes of age in the species and will provide knowledge and technologies to help optimize production potential, reproduction efficiency, product quality, nutritional value and resistance to diseases. As water buffalo and domestic cattle, both members of the Bovidae family, are closely related, the vast amount of cattle genetic/genomic resources might serve as shortcuts for the buffalo community to further advance genome science and biotechnologies in the species.

Assessment of cattle genetic introgression into domestic yak populations using mitochondrial and microsatellite DNA markers

Hybridization between yak Poephagus grunniens and taurine Bos taurus or indicine B. indicus cattle has been widely practiced throughout the yak geographical range, and gene flow is expected to have occurred between these species. To assess the impact of cattle admixture on domestic yak, we examined 1076 domestic yak from 29 populations collected in China, Bhutan, Nepal, India, Pakistan, Kyrgyzstan, Mongolia and Russia using mitochondrial DNA and 17 autosomal microsatellite loci. A cattle diagnostic marker-based analysis reveals cattle-specific mtDNA and/or autosomal microsatellite allele introgression in 127 yak individuals from 22 populations. The mean level of cattle admixture across the populations, calculated using allelic information at 17 autosomal microsatellite loci, remains relatively low (mY(cattle) = 2.66 +/- 0.53% and Q(cattle) = 0.69 +/- 2.58%), although it varies a lot across populations as well as among individuals within population. Although the level of cattle admixture shows a clear geographical structure, with higher levels of admixture in the Qinghai-Tibetan Plateau and Mongolian and Russian regions, and lower levels in the Himalayan and Pamir Plateau region, our results indicate that the level of cattle admixture is not significantly correlated with the altitude across geographical regions as well as within geographical region. Although yak-cattle hybridization is primarily driven to produce F(1) hybrids, our results show that the subsequent gene flow between yak and cattle took place and has affected contemporary genetic make-up of domestic yak. To protect yak genetic integrity, hybridization between yak and cattle should be tightly controlled.

Survival of the thriftiest: restricted nurture reveals the thrifty nature of a growth gene in Bos indicus

Growth hormone (GH) is a part of the somatotropic axis that controls metabolism, growth, development and aging in a wide range of animals. Mutations that reduce GH signaling have been associated with extended life spans and increased longevity in ways similar to what is observed in dietary restriction (DR) models. However, the mechanism by which DR works is not well understood. Here, we show that DR works as a factor in the evolution of the genetic make-up of domestic cattle. In a series of 6864 bovines of seven Bos indicus and tropically adapted Bos taurus breeds, the frequency of a short, wild-type allele of the promoter region of the bovine GH gene, G1 allele, varied from 2.7 to 17.7%. The frequency of the long, domestic G2 allele increased from 88 to 95% along 20 calf crops of commercial Bos indicus cattle of the Nelore breed undergoing selection for increasing post-weaning weight gain with ad libitum nutrient intake. Under DR, however, the G1 allele sustained growth better than the G2 allele, as observed in a series of feeding tests. The G2 allele was even detrimental or abiotropic, as it caused rapid body decay under DR. We observed a reflection symmetry of GH allele substitution effects on body weight under different dietary schemes. The G2 allele is featured as the "demanding allele", because it is optimally fitted to ad libitum nutrient intake. The G1 allele is featured as the "thrifty allele" because it is optimally fitted to DR. Our results show that dietary regimens need not extend lifespan or increase longevity in the sense of age-specific fitness. Instead, adaptation to any particular dietary regimen is just as much a consequence of selection as its cause; dietary regimens work as do any selection force, optimizing genotypic fitness to nutritional conditions.

Phylogenetic reconstruction and the identification of ancient polymorphism in the Bovini tribe (Bovidae, Bovinae)

Background

The Bovinae subfamily incorporates an array of antelope, buffalo and cattle species. All of the members of this subfamily have diverged recently. Not surprisingly, a number of phylogenetic studies from molecular and morphological data have resulted in ambiguous trees and relationships amongst species, especially for Yak and Bison species. A partial phylogenetic reconstruction of 13 extant members of the Bovini tribe (Bovidae, Bovinae) from 15 complete or partially sequenced autosomal genes is presented.

Results

We identified 3 distinct lineages after the Bovini split from the Boselaphini and Tragelaphini tribes, which has lead to the (1) Buffalo clade (Bubalus and Syncerus species) and a more recent divergence leading to the (2) Banteng, Gaur and Mithan and (3) Domestic cattle clades. A fourth lineage may also exist that leads to Bison and Yak. However, there was some ambiguity as to whether this was a divergence from the Banteng/Gaur/Mithan or the Domestic cattle clade. From an analysis of approximately 30,000 sites that were amplified in all species 133 sites were identified with ambiguous inheritance, in that all trees implied more than one mutation at the same site. Closer examination of these sites has identified that they are the result of ancient polymorphisms that have subsequently undergone lineage sorting in the Bovini tribe, of which 53 have remained polymorphic since Bos and Bison species last shared a common ancestor with Bubalus between 5–8 million years ago (MYA).

Conclusion

Uncertainty arises in our phylogenetic reconstructions because many species in the Bovini diverged over a short period of time. It appears that a number of sites with ambiguous inheritance have been maintained in subsequent populations by chance (lineage sorting) and that they have contributed to an association between Yak and Domestic cattle and an unreliable phylogenetic reconstruction for the Bison/Yak clade. Interestingly, a number of these aberrant sites are in coding sections of the genome and their identification may have important implications for studying the neutral rate of mutation at nonsynonymous sites. The presence of these sites could help account for the apparent contradiction between levels of polymorphism and effective population size in domesticated cattle.

Sequence variation and molecular evolution of hormone-sensitive lipase genes in species of bovidae

The partial sequences of exon I of hormone-sensitive lipase (HSL) genes in yak (Bos grunniens), cattle (Bos taurus), zebu (Bos indicus), and buffalo (Bubalus bubalis) were analyzed. Comparisons of these sequences and the deduced amino acid sequences with the homologous HSL gene and protein sequences in other mammalian species including pig (Sus scrofa), human (Homo sapiens), mouse (Mus musculus), and rat (Rattus sp.) retrieved from the GenBank were carried out and finally a phylogenetic tree was constructed using the partial DNA sequences of the HSL genes in all species. The results showed that the homologies of the partial exon I sequences of the HSL genes between yak and cattle, zebu, buffalo, pig, human, mouse, and rat were as high as 99.8%, 99.6%, 97.4%, 90.6%, 88.4%, 83.5%, and 82.3%, respectively. This was accompanied by highly homologous amino acid sequences of the HSLs: 100%, 100%, 98.2%, 94.0%, 92.2%, 89.8%, and 89.8% identity, respectively. There are more transitions, less transversions, and no insertion or deletion in variable nucleotides of the HSL genes between the yak and other species. The majority of the variable mutations was synonymous and was found most frequently at the third codon, followed by the first and second codons, a finding that was in accordance with the neutralism hypothesis for molecular evolution. In the phylogenetic tree, the cattle and zebu were clustered together first, followed by the yak, buffalo, pig, human, mouse, and rat. This was in agreement with taxonomy suggesting that the partial sequences of exon I of the HSL genes were useful in constructing the phylogenetic tree of mammalian species. Among the four species of Bovidae, genetic differentiation in the HSL genes between yak and buffalo is equivalent to that between buffalo and cattle and between buffalo and zebu. Furthermore, the genetic distances in the HSL genes are much smaller between yak, cattle, and zebu than those between each of the three species and the buffalo. Therefore, it is reasonable to consider yak as an independent species of the genus Bos.

Evolution of microsatellite loci in the adaptive radiation of Hawaiian honeycreepers

Previous studies have examined germ-line mutations to infer the processes that generate and maintain variability in microsatellite loci. Few studies, however, have examined patterns to infer processes that act on microsatellite loci over evolutionary time. Here, we examine changes in 8 dinucleotide loci across the adaptive radiation of Hawaiian honeycreepers. The loci were found to be highly variable across the radiation, and we did not detect ascertainment bias with respect to allelic diversity or allele size ranges. In examining patterns at the sequence level, we found that changes in flanking regions, repeat motifs, or repeat interruptions were often shared between closely related species and may be phylogenetically informative. Genetic distance measures based on microsatellites were strongly correlated with those based on mitochondrial DNA (mtDNA) sequences as well as with divergence time up to 3 My. Phylogenetic inferences based on microsatellite genetic distances consistently recovered 2 of the 4 honeycreeper clades observed in a tree based on mtDNA sequences but differed from the mtDNA tree in the relationships among clades. Our results confirm that microsatellite loci may be conserved over evolutionary time, making them useful in population-level studies of species that diverged from the species in which they were characterized as long as 5 Ma. Despite this, we found that their use in phylogenetic inference was limited to closely related honeycreeper species.

A first generation whole genome RH map of the river buffalo with comparison to domestic cattle

Background

The recently constructed river buffalo whole-genome radiation hybrid panel (BBURH₅₀₀₀) has already been used to generate preliminary radiation hybrid (RH) maps for several chromosomes, and buffalo-bovine comparative chromosome maps have been constructed. Here, we present the first-generation whole genome RH map (WG-RH) of the river buffalo generated from cattle-derived markers. The RH maps aligned to bovine genome sequence assembly Btau_4.0, providing valuable comparative mapping information for both species.

Results

A total of 3990 markers were typed on the BBURH₅₀₀₀ panel, of which 3072 were cattle derived SNPs. The remaining 918 were classified as cattle sequence tagged site (STS), including coding genes, ESTs, and microsatellites. Average retention frequency per chromosome was 27.3% calculated with 3093 scorable markers distributed in 43 linkage groups covering all autosomes (24) and the X chromosomes at a LOD ≥ 8. The estimated total length of the WG-RH map is 36,933 cR₅₀₀₀. Fewer than 15% of the markers (472) could not be placed within any linkage group at a LOD score ≥ 8. Linkage group order for each chromosome was determined by incorporation of markers previously assigned by FISH and by alignment with the bovine genome sequence assembly (Btau_4.0).

Conclusion

We obtained radiation hybrid chromosome maps for the entire river buffalo genome based on cattle-derived markers. The alignments of our RH maps to the current bovine genome sequence assembly (Btau_4.0) indicate regions of possible rearrangements between the chromosomes of both species. The river buffalo represents an important agricultural species whose genetic improvement has lagged behind other species due to limited prior genomic characterization. We present the first-generation RH map which provides a more extensive resource for positional candidate cloning of genes associated with complex traits and also for large-scale physical mapping of the river buffalo genome.

Molecular anatomy of the cytoplasmic domain of bovine growth hormone receptor, a quantitative trait locus

Quantitative trait loci (QTL) studies have indicated growth hormone receptor (GHR) as a candidate gene affecting cattle milk yield and composition. In order to characterize genetic variation at GHR in cattle, we studied European and East African breeds with different histories of selection, and Bos grunniens, Ovis aries, Sus scrofa, Bison bison and Rangifer tarandus as references. We sequenced most of the cytoplasmic domain (900 bp of exon 10), 89 bp of exon 8, including the putative causative mutation for the QTL effect, and 390 bp of intron 8 for comparison. In the cytoplasmic domain, seven synonymous and seven non-synonymous single nucleotide polymorphisms (SNP) were identified in cattle. Three non-synonymous SNPs were found in sheep and one synonymous SNP in yak, while other studied species were monomorphic. Three major haplotypes were observed, one unique to African breeds, one unique to European breeds and one shared. Bison and yak haplotypes are derivatives of the European haplotype lineage. Most of the exon 10 non-synonymous cattle SNPs appear at phylogenetically highly conserved sites. The polymorphisms in exon 10 cluster around a ruminant-specific tyrosine residue, suggesting that this site may act as an additional signalling domain of GHR in ruminants. Alternative explanations for the persistent polymorphism include balancing selection, hitch-hiking, pleiotropic or sexually antagonistic fitness effects or relaxed functional constraints.

Genetic variation in the kappa-casein gene (CSN3) of Chinese yak (Bos grunniens) and phylogenetic analysis of CSN3 sequences in the genus Bos

Variants of kappa-casein (CSN3) have been extensively studied in cattle and 13 alleles have been identified at the protein and DNA levels to date. Evolution of some of these alleles and a possible common ancestor remain unclear. Polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) analysis of CSN3 exon IV in domesticated yak revealed a 2-allele polymorphism showing migration patterns different from known cattle variants. The PCR products of both yak CSN3 SSCP alleles were sequenced. All yak had nucleotide sequences coding for Thr in AA position 136 (identical to bovine CSN3*A) and Ala in position 148 (identical to bovine CSN3*B). The sequencing results were confirmed by PCR-RFLP analysis using HindIII and TaqI. A 12-bp insertion in the coding region, representing a repeated nucleotide and AA motif, was found in 1 yak allele. The duplication corresponds to the codons for AA 147 to 150 (Glu-Ala-Ser-Pro) or 148 to 151 (Ala-Ser-Pro-Glu), which are repeated identically. In 21 yak samples genotyped by PCR-SSCP analysis, frequencies for the insertion variant and the short variant were about 68 and 32%, respectively. The loss of the insertion may have led to the ancestral CSN3 allele from which all currently known variants of CSN3 in the genus Bos evolved. This is the first report of polymorphisms in the yak CSN3 gene and may be helpful for future studies on genetic variation within and between yak populations or on associated traits.

Genomic conservation of cattle microsatellite loci in wild gaur (Bos gaurus) and current genetic status of this species in Vietnam

BACKGROUND

The wild gaur (Bos gaurus) is an endangered wild cattle species. In Vietnam, the total number of wild gaurs is estimated at a maximum of 500 individuals. Inbreeding and genetic drift are current relevant threats to this small population size. Therefore, information about the genetic status of the Vietnamese wild gaur population is essential to develop strategies for conservation and effective long-term management for this species. In the present study, we performed cross-species amplification of 130 bovine microsatellite markers, in order to evaluate the applicability and conservation of cattle microsatellite loci in the wild gaur genome. The genetic diversity of Vietnamese wild gaur was also investigated, based on data collected from the 117 successfully amplified loci.

RESULTS

One hundred-thirty cattle microsatellite markers were tested on a panel of 11 animals. Efficient amplifications were observed for 117 markers (90%) with a total of 264 alleles, and of these, 68 (58.1%) gave polymorphic band patterns. The number of alleles per locus among the polymorphic markers ranged from two to six. Thirteen loci (BM1314, BM2304, BM6017, BMC2228, BMS332, BMS911, CSSM023, ETH123, HAUT14, HEL11, HEL5, ILSTS005 and INRA189) distributed on nine different cattle chromosomes failed to amplify wild gaur genomic DNA. Three cattle Y-chromosome specific microsatellite markers (INRA124, INRA126 and BM861) were also highly specific in wild gaur, only displaying an amplification product in the males. Genotype data collected from the 117 successfully amplified microsatellites were used to assess the genetic diversity of this species in Vietnam. Polymorphic Information Content (PIC) values varied between 0.083 and 0.767 with a mean of 0.252 while observed heterozygosities (Ho) ranged from 0.091 to 0.909 (mean of 0.269). Nei's unbiased mean heterozygosity and the mean allele number across loci were 0.298 and 2.2, respectively.

CONCLUSION

Extensive conservation of cattle microsatellite loci in the wild gaur genome, as shown by our results, indicated a high applicability of bovine microsatellites for genetic characterization and population genetic studies of this species. Moreover, the low genetic diversity observed in Vietnamese wild gaur further underlines the necessity of specific strategies and appropriate management plans to preserve this endangered species from extinction.

Genetic Structure and Differentiation of Three Chinese Indigenous Cattle Populations

Levels of genetic differentiation, gene flow, and genetic structure of three indigenous cattle populations (Luxi, Bohai, and Minnan) and two reference cattle populations (Chinese Holstein and Qinhai yak) in China were estimated using the information from 12 microsatellites, and 141 microsatellite alleles were identified. The mean number of alleles per locus ranged from 2.9005 in yak to 4.9722 in Holstein. The observed heterozygosity ranged from 0.5325 (yak) to 0.7719 (Holstein); 29 private alleles were detected. The global heterozygote deficit across all populations amounted to 58.5% (p < 0.001). The overall significant (p < 0.001) deficit of heterozygotes because of inbreeding within breeds amounted to 43.2%. The five cattle populations were highly differentiated (F (st) = 26.9%, p < 0.001) at all loci. The heterozygote deficit within the population was highest in Luxi cattle and lowest in yak. The average number of effective migrants exchanged per generation was highest (1.149) between Luxi and Holstein, and lowest (0.509) between Luxi and yak. With the application of prior population information, cluster analysis achieved posterior probabilities from 91% to 98% of correctly assigning individuals to populations. Combining the information of cluster analysis, gene flow, and Structure analysis, the five cattle populations belong to three genetic clusters, a taurine (Luxi and Chinese Holstein), a zebu (Bohai and Minnan), and a yak cluster. This indicates that Bohai black is closer to Bos indicus than Luxi cattle. The evolution and development of three indigenous cattle populations are discussed.

Accelerated molecular evolution in Microtus (Rodentia) as assessed via complete mitochondrial genome sequences

Microtus is one of the most taxonomically diverse mammalian genera, including over 60 extant species. These rodents have evolved rapidly, as the genus originated less than 2 million years ago. If these numbers are taken at face value, then an average of 30 microtine speciation events have occurred every million years. One explanation for the rapid rate of cladogenesis in Microtus could be the karyotypic differentiation exhibited across the genus: diploid numbers range from 17 to 64. Despite the striking chromosomal variability within Microtus, phenotypic variation is unremarkable. To determine whether nucleotide substitution rates are also elevated in voles, we sequenced the entire mitochondrial DNA (mtDNA) genome of the Eurasian sibling vole (Microtus rossiaemeridionalis). We compared this genome to another previously sequenced vole mtDNA genome (Microtus kikuchii) and performed pairwise sequence comparisons with the mtDNA genomes of ten additional mammalian genera. We found that microtine mtDNA genomes are evolving more rapidly than any other mammalian lineage we sampled, as gauged by the rate of nucleotide substitution across the entire mtDNA genome as well as at each individual protein-coding gene. Additionally, we compared substitution rates within the cytochrome b gene to seven other rodent genera and found that Microtus mtDNA is evolving fastest. The root cause of accelerated evolution in Microtus remains uncertain, but merits further investigation.

Application of bovine microsatellite markers for genetic diversity analysis of European bison (Bison bonasus)

In this study, the cross-amplification of a commercial multiplex set of 11 cattle (Bos taurus) microsatellites was tested on a panel of 35 European bison (Bison bonasus) individuals. After polymerase chain reaction optimization, all loci cross-amplified successfully in investigated bisons. Number of alleles and observed and expected heterozygosity per locus are in the range of 2-4, 0.086-0.629 and 0.288-0.621 respectively. The availability of a heterologous set of multiplexed microsatellite markers derived from cattle opens an avenue for collecting profound genetic data for efficient conservation management strategies of the European bison.