EAT and AAI, along eigenvector 1, which explains ~15% of the total variation. AFT Muturu and N'Dama are close to EAT along the eigenvector 1. Most of the AFH cattle cluster together regardless of their breed memberships, leaving only Ankole, Mursi and Sheko outside the main cluster toward the AFT Muturu and N'Dama. The PCA results also show that Muturu and N'Dama, our representative of AFT population, are separated from the other cattle groups (eigenvector 2, ~2.5% of total variation). Sheko positions close to the AFH, as similarly reported in other studies 5,43 . Genetic clustering analysis using ADMIXTURE 44 corroborates the pattern found in PCA (Fig. 2b and Extended Data Fig. 2). Most of AFH show a similar proportion of taurine ancestry, around 25% on average. Only a few AFH breeds have elevated taurine ancestry: Ankole (53.37 ± 1.49%), Sheko (46.28 ± 2.03%) and Mursi (35.90 ± 2.16%). (Fig. 2b).Genetic distance and diversity. Pairwise F st were calculated to estimate the genetic distances between populations (n = 38) (Extended Data Fig. 3). Taurine (EUT, AST and AFT) show F st values of 0.1568 and 0.3287 on average against AFH and AAI, respectively.Across AFH, pairwise F st between breeds is close to zero, regardless of their phenotypic classification as African Zebu, Sanga or Zenga. Muturu and N'Dama show F st value of 0.1769, 0.1847 and 0.3734 against AFH, EAT and AAI, respectively.The genome-wide autosomal SNPs show reduced levels of heterozygosity in the taurine (0.0021 ± 0.0005/bp) compared to all other populations (0.0048 ± 0.0008/bp). Heterozygosity values of AFH are similarly higher across populations (0.0046 ± 0.0003/bp). AAI shows a higher level of heterozygosity compared to AFH (0.0052 ± 0.0014/bp) (Extended Data Fig. 4). The degree of inbreeding measured by runs of homozygosity (ROH) shows that taurine, including Muturu and N'Dama, have a higher level of inbreeding compared to the other and Ethiopia), the University of Khartoum (Sudan), and the National Biotechnology Development Agency (NABDA) (Nigeria). The following institutions and their personnel provided help for the sampling of the African cattle: ILRI Kapiti Ranch, Ministry of Animal Resources, Fisheries and Range (Sudan), Ol Pejeta Conservancy (Kenya), Institute of Biodiversity (Ethiopia), the Directors of Veterinary Services and the cattle keepers from Ethiopia, Kenya, Uganda and Sudan. ILRI livestock genomics program is supported by the

The East African Shorthorn Zebu (EASZ) cattle are ancient hybrid between Asian zebu × African taurine cattle preferred by local farmers due to their adaptability to the African environment. The genetic controls of these adaptabilities are not clearly understood yet. Here, we genotyped 92 EASZ samples from Kenya (KEASZ) with more than 770,000 SNPs and sequenced the genome of a pool of 10 KEASZ. We observe an even admixed autosomal zebu × taurine genomic structure in the population. A total of 101 and 165 candidate regions of positive selection, based on genome-wide SNP analyses (meta-SS, Rsb, iHS, and ΔAF) and pooled heterozygosity (Hp) full genome sequence analysis, are identified, in which 35 regions are shared between them. A total of 142 functional variants, one novel, have been detected within these regions, in which 30 and 26 were classified as of zebu and African taurine origins, respectively. High density genome-wide SNP analysis of zebu × taurine admixed cattle populations from Uganda and Nigeria show that 25 of these regions are shared between KEASZ and Uganda cattle, and seven regions are shared across the KEASZ, Uganda, and Nigeria cattle. The identification of common candidate regions allows us to fine map 18 regions. These regions intersect with genes and QTL associated with reproduction and environmental stress (e.g., immunity and heat stress) suggesting that the genome of the zebu × taurine admixed cattle has been uniquely selected to maximize hybrid fitness both in terms of reproduction and survivability.

Despite only 8% of cattle being found in Europe, European breeds dominate current genetic resources. This adversely impacts cattle research in other important global cattle breeds, especially those from Africa for which genomic resources are particularly limited, despite their disproportionate importance to the continent’s economies. To mitigate this issue, we have generated assemblies of African breeds, which have been integrated with genomic data for 294 diverse cattle into a graph genome that incorporates global cattle diversity. We illustrate how this more representative reference assembly contains an extra 116.1 Mb (4.2%) of sequence absent from the current Hereford sequence and consequently inaccessible to current studies. We further demonstrate how using this graph genome increases read mapping rates, reduces allelic biases and improves the agreement of structural variant calling with independent optical mapping data. Consequently, we present an improved, more representative, reference assembly that will improve global cattle research.