Genome‐wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26–0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise F_ST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome‐wide smoothed F_ST, 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty‐four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.     

Assessing introgression of Sahelian zebu genes into native Bos taurus breeds in Burkina Faso

A total of 350 samples were analyzed to estimate zebu gene proportions into two different taurine cattle breeds of Burkina Faso (Lobi and N’Dama) using 38 microsatellites and various statistical methodologies. West African and East African zebu samples were sequentially used as reference parental populations. Furthermore, N’Dama cattle from Congo, the composite South African Bonsmara cattle breed and a pool of European cattle were used successively as second parental populations. Independently of the methodology applied: (a) the use of West African zebu samples gave higher admixture coefficients than the East African zebu; (b) the higher zebu proportions were estimated when the European cattle was used as parental population 2; and (c) the use of the N’Dama population from Congo as parental population 2 gave the more consistent zebu proportion estimates for both the Lobi and the N’Dama breeds. In any case, the zebu admixture proportions estimated were not negligible and were always higher in the N’Dama cattle than in the Lobi cattle of Burkina Faso. This suggested that the introgression of Sahelian zebu genes into the taurine cattle of Southern West Africa can follow a complex pattern that can depend on local agro-ecological features. The current research pointed out that the estimation of admixture coefficients is highly dependent on both the assumptions underlying the methodologies applied and the selection of parental populations. Our analyses suggest that either too high or nil genetic identity between the parental and the expectedly derived populations must be avoided.     

A polymorphism in randomly amplified DNA that differentiates the Y chromosomes of Bos indicus and Bos taurus

A small number of west African Bos taurus cattle breeds, including the N'Dama, constitute a valuable genetic resource by virtue of their ability to remain productive under trypanosomiasis challenge. However, introgression of Bos indicus genes into the trypanotolerant breeds, particularly by introduction of zebu bulls, is a threat to this resource. This work describes the characterization and cloning of a bovine randomly amplified polymorphic DNA (RAPD) that is generated in polymerase chain reaction (PCR) with the 10 base primer ILO1065 from Bos indicus male templates, but not from B. taurus male templates or female templates of either type. Male-specific sequences with homology to the RAPD also occur in B. taurus breeds. This suggests that the polymorphism may be due to base substitution(s) in an ILO1065 priming site, or insertion/deletion events either affecting priming sites or occurring between sites on the cattle Y chromosome. We have shown that cattle, whether of B. indicus or B. taurus phenotype, which possess a typically B. indicus metaphase Y chromosome on the basis of QFQ banding, have a B. indicus ILO1065-generated genotype. The ILO1065-primed RAPD can be used in a simple dot blot assay as a probe of RAPD-PCR products, to provide a convenient, reliable and effective means of detecting introgression of zebu genes in B. taurus cattle populations.     

Complete sequence of the yak (Bos grunniens) mitochondrial genome and its evolutionary relationship with other ruminants

The yak (Bos grunniens) is the most important domesticated species in the Qinhai-Tibetan Plateau. In present study, the complete sequence of the yak mitochondrial genome was determined. Sequence analysis revealed that there are no differences with cattle in the yak mitochondrial genome organization. Interestingly, within the D-loop, the conserved sequence blocks are less conserved than surrounding regions. Neighbor-Joining (NJ) trees based on single genes, gene sets and concatenated genes of mitochondrial genome were constructed. The analysis identified the yak as a sister group of a cattle/zebu clade. Based on substitutions in 22 tRNA genes, 12S rRNA gene and 16S rRNA gene, the dating of divergence between yak and cattle/zebu, and yak and water buffalo, was proposed to have occurred 4.38-5.32 and 10.54-13.85 million years before present, respectively. This is consistent with the paleontologyical data. Yak and sheep/goat divergent dating predicts that their divergence occurred at 13.14-27.99 million years before the present day.     

Comprehensive Evolutionary and Expression Analysis of FCS-Like Zinc finger Gene Family Yields Insights into Their Origin,Expansion and Divergence

Plant evolution is characterized by frequent genome duplication events. Expansion of habitat resulted in the origin of many novel genes and genome duplication events which in turn resulted in the expansion of many regulatory gene families. The plant-specific FCS-Like Zinc finger (FLZ) gene family is characterized by the presence of a FCS-Like Zinc finger (FLZ) domain which mediates the protein-protein interaction. In this study, we identified that the expansion of FLZ gene family size in different species is correlated with ancestral and lineage-specific whole genome duplication events. The subsequent gene loss found to have a greater role in determining the size of this gene family in many species. However, genomic block duplications played the significant role in the expansion of FLZ gene family in some species. Comparison of Arabidopsis thaliana and Oryza sativa FLZ gene family revealed monocot and dicot specific evolutionary trends. The FLZ genes were found to be under high purifying selection. The spatiotemporal expression analyses of Arabidopsis thaliana FLZ gene family revealed that majority of the members are highly expressed in reproductive organs. FLZ genes were also found to be highly expressed during vegetative-to-reproductive phase transition which is correlated with the proposed role of this gene family in sugar signaling. The comparison of sequence, structural and expression features of duplicated genes identified lineage-specific redundancy and divergence. This extensive evolutionary analysis and expression analysis of Arabidopsis thaliana FLZ genes will pave the way for further functional analysis of FLZ genes.     

Analysis of copy number variants in the cattle genome

Copy number variation (CNV) is likely to be an important component of heritable variation in livestock. To characterise CNVs in cattle, we performed a genome wide survey to determine the number, location and gene content of these genomic features. A tiling oligonucleotide array with ~385,000 probes was used for comparative genomic hybridisation of both taurine and zebu cattle. Using a conservative set of calling criteria, a total of 51 CNV were detected that collectively spanned approximately half of one percent of the bovine genome. The size of the average CNV within each animal ranged from 213 kb up to 335 kb. Half of the CNV were detected in a single animal only, whilst the remainder was independently identified in multiple individuals. Analysis was performed to determine the gene content for each CNV region. This revealed that the majority of CNV (82%) spanned at least one gene, with a number of CNV containing genes which are known to control aspects of phenotypic variation in cattle. Whilst additional studies are required to determine the impact of individual CNV, this study confirmed them as an important class of genomic variation in cattle.     

Diversity and phylogeny of mitochondrial DNA isolated from mithun Bos frontalis located in Bhutan

We sequenced the 16S rRNA gene in mitochondrial DNA to characterize mithun located in Bhutan and to increase our understanding of its origin. We compared mithun with yak, European cattle, Bhutanese zebu and Indian zebu. Sequencing revealed low nucleotide diversity within the mithun population and their phylogenetic proximity to gaur. A close relationship between Bhutanese mithun and gaur was confirmed by an additional comparison with wild gaur specimens from three locations in Bhutan. Direct domestication of mithun from gaur was supported, while maternal contribution from the cattle lineage during domestication was not supported.     

Hair follicle measurements in some tropical and temperate breeds of cattle in Nigeria

Anatomical studies of the skin of mature White Fulani, N'Dama, Muturu and the Friesian cattle showed that there are significant regional and breeds differences in hair follicle density and depth. The hair follicle distribution per cm² of skin was higher on the neck, ventral and the upper limb sites (2000 – 2496) than on the dorsal, midside and eye margins (1344 – 1920). Low densities were recorded on the horn margin, perineum and the lower limb sites (576 – 960). This regional distribution was consistent in all breeds except the White Fulani in which follicle densities were significantly lower for most body sites. Hair follicle depth was greater in the Friesians (1.63 ± 0.68 mm) than in the N'Dama (0.98 ± 0.13 mm), the White Fulani (0.91 ± 0.31 mm) and the Muturu (0.84 ± 0.25 mm) cattle. Among body sites, follicle depth was significantly greater on the horn margin, lower and upper limbs and the scrotum (1.16 – 1.48 mm) than on other body sites (0.86 – 0.98 mm). A significant negative correlation was found between hair follicle density and follicle depth.     

Genome-wide analysis reveals the ancient and recent admixture history of East African Shorthorn Zebu from Western Kenya

The Kenyan East African zebu cattle are valuable and widely used genetic resources. Previous studies using microsatellite loci revealed the complex history of these populations with the presence of taurine and zebu genetic backgrounds. Here, we estimate at genome-wide level the genetic composition and population structure of the East African Shorthorn Zebu (EASZ) of western Kenya. A total of 548 EASZ from 20 sub-locations were genotyped using the Illumina BovineSNP50 v. 1 beadchip. STRUCTURE analysis reveals admixture with Asian zebu, African and European taurine cattle. The EASZ were separated into three categories: substantial (⩾12.5%), moderate (1.56%<X<12.5%) and non-introgressed (⩽1.56%) according to the European taurine genetic proportion. The non-European taurine introgressed animals (n=425) show an unfluctuating zebu and taurine ancestry of 0.84±0.009 s.d. and 0.16±0.009 s.d., respectively, with significant differences in African taurine (AT) and Asian zebu backgrounds across chromosomes (P<0.0001). In contrast, no such differences are observed for the European taurine ancestry (P=0.1357). Excluding European introgressed animals, low and nonsignificant genetic differentiation and isolation by distance are observed among sub-locations (F_st=0.0033, P=0.09; r=0.155, P=0.07). Following a short population expansion, a major reduction in effective population size (N_e) is observed from approximately 240 years ago to present time. Our results support ancient zebu × AT admixture in the EASZ population, subsequently shaped by selection and/or genetic drift, followed by a more recent exotic European cattle introgression.     

Microsatellite DNA Variation and the Evolution, Domestication and Phylogeography of Taurine and Zebu Cattle (Bos Taurus and Bos Indicus)

Genetic variation at 20 microsatellite loci was surveyed to determine the evolutionary relationships and molecular biogeography of 20 different cattle populations from Africa, Europe and Asia. Phylogenetic reconstruction and multivariate analysis highlighted a marked distinction between humpless (taurine) and humped (zebu) cattle, providing strong support for a separate origin for domesticated zebu cattle. A molecular clock calculation using bison (Bison sp.) as an outgroup gave an estimated divergence time between the two subspecies of 610,000-850,000 years. Substantial differences in the distribution of alleles at 10 of these loci were observed between zebu and taurine cattle. These markers subsequently proved very useful for investigations of gene flow and admixture in African populations. When these data were considered in conjunction with previous mitochondrial and Y chromosomal studies, a distinctive male-mediated pattern of zebu genetic introgression was revealed. The introgression of zebu-specific alleles in African cattle afforded a high resolution perspective on the hybrid nature of African cattle populations and also suggested that certain West African populations of valuable disease-tolerant taurine cattle are under threat of genetic absorption by migrating zebu herds.     

中国部分地方黄牛品种线粒体D-loop区的遗传变异与进化分析

测定了13个黄牛品种125个个体的线粒体D-loop区段的全序列,包括12个中国地方黄牛品种的123个个体和德国黄牛2个个体,并进行了分析。结果显示,共检测到93个变异位点,57个单倍型,平均核苷酸差异(average number ofnucleotide differences,k)为22.708,核苷酸多样度(nucleotide diversity,π)为0.0251±0.00479,单倍型多样度(haplotypediversity,Hd)为0.888±0.026,表明我国黄牛品种遗传多样性非常丰富。构建的Neighbor-Joining进化树显示这13个品种主要分成两大类型:普通牛和瘤牛;新发现的特殊类型Ⅲ只有一个西藏阿沛甲咂牛的个体,它与牦牛D-loop序列最相近,证明西藏地区的黄牛与牦牛之间存在基因渗入现象。普通牛和瘤牛在日喀则驼峰牛中占的比例分别是64.3%和35.7%,在阿沛甲咂牛中占的比例分别是50.0%和50.0%,证明了西藏的黄牛也有瘤牛类型。云南牛品种的单倍型非常丰富证明了云南在中国黄牛起源上的重要地位;在27个中国黄牛品种中(本研究11个品种以及GenBank上的16个品种)找到了中国瘤牛的核心单倍型i1,并且对它进行了讨论。同时证明了西藏瘤牛独立于中国瘤牛核心类群的特殊性。     

Polymerase chain reaction analysis of mitochondrial DNA polymorphism in N'Dama and Zebu cattle

A study of polymorphisms of mitochondrial DNA (mtDNA) of West African N'Dama (Bos taurus) and East African Zebu (B. indicus) cattle was carried out to obtain information on maternal phylogenetic relationships between these breeds. A relatively large sample size was made possible by using polymerase chain reaction (PCR) amplification of DNA prepared from small blood samples to generate fragments of two known polymorphic mtDNA regions, one within the gene encoding subunit 5 of NADH dehydrogenase and one encompassing the entire D-loop. This approach allowed us to achieve a higher resolution restriction analysis on mtDNA from more animals than would have been possible by conventional methods. PCR-amplified mtDNA of 58 animals from five populations was examined at 26 restriction sites by 16 enzymes. In this way 154 nucleotides of mtDNA were scanned for polymorphism. Six polymorphic sites were located by this means, five of which were within the D-loop and one of which was within the NADH dehydrogenase 5 gene. None of the polymorphisms observed could be con sidered typical of breed or type.     

Genomic signatures for drylands adaptation at gene-rich regions in African zebu cattle

BackgroundIndigenous Sudanese cattle are mainly indicine/zebu (humped) type. They thrive in the harshest dryland environments characterised by high temperatures, long seasonal dry periods, nutritional shortages, and vector disease challenges. Here, we sequenced 60 indigenous Sudanese cattle from six indigenous breeds and analysed the data using three genomic scan approaches to unravel cattle adaptation to the African dryland region.ResultsWe identified a set of gene-rich selective sweep regions, detected mostly on chromosomes 5, 7 and 19, shared across African and Gir zebu. These include genes involved in immune response, body size and conformation, and heat stress response. We also identified selective sweep regions unique to Sudanese zebu. Of these, a 250 kb selective sweep on chromosome 16 spans seven genes, including PLCH2, PEX10, PRKCZ, and SKI, which are involved in alternative adaptive metabolic strategies of insulin signalling, glucose homeostasis, and fat metabolism.ConclusionsOur results suggest that environmental adaptation may involve recent and ancient selection at gene-rich regions, which might be under a common regulatory genetic control, in zebu cattle.     

A Non-Coding Genomic Duplication at the HMX1 Locus Is Associated with Crop Ears in Highland Cattle

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.     

Genome Duplication and Gene Loss Affect the Evolution of Heat Shock Transcription Factor Genes in Legumes

Whole-genome duplication events (polyploidy events) and gene loss events have played important roles in the evolution of legumes. Here we show that the vast majority of Hsf gene duplications resulted from whole genome duplication events rather than tandem duplication, and significant differences in gene retention exist between species. By searching for intraspecies gene colinearity (microsynteny) and dating the age distributions of duplicated genes, we found that genome duplications accounted for 42 of 46 Hsf-containing segments in Glycine max, while paired segments were rarely identified in Lotus japonicas, Medicago truncatula and Cajanus cajan. However, by comparing interspecies microsynteny, we determined that the great majority of Hsf-containing segments in Lotus japonicas, Medicago truncatula and Cajanus cajan show extensive conservation with the duplicated regions of Glycine max. These segments formed 17 groups of orthologous segments. These results suggest that these regions shared ancient genome duplication with Hsf genes in Glycine max, but more than half of the copies of these genes were lost. On the other hand, the Glycine max Hsf gene family retained approximately 75% and 84% of duplicated genes produced from the ancient genome duplication and recent Glycine-specific genome duplication, respectively. Continuous purifying selection has played a key role in the maintenance of Hsf genes in Glycine max. Expression analysis of the Hsf genes in Lotus japonicus revealed their putative involvement in multiple tissue-/developmental stages and responses to various abiotic stimuli. This study traces the evolution of Hsf genes in legume species and demonstrates that the rates of gene gain and loss are far from equilibrium in different species.