Terhi Iso-Touru

A 660-Kb Deletion with Antagonistic Effects on Fertility and Milk Production Segregates at High Frequency in Nordic Red Cattle: Additional Evidence for the Common Occurrence of Balancing Selection in Livestock

In dairy cattle, the widespread use of artificial insemination has resulted in increased selection intensity, which has led to spectacular increase in productivity. However, cow fertility has concomitantly severely declined. It is generally assumed that this reduction is primarily due to the negative energy balance of high-producing cows at the peak of lactation. We herein describe the fine-mapping of a major fertility QTL in Nordic Red cattle, and identify a 660-kb deletion encompassing four genes as the causative variant. We show that the deletion is a recessive embryonically lethal mutation. This probably results from the loss of RNASEH2B, which is known to cause embryonic death in mice. Despite its dramatic effect on fertility, 13%, 23% and 32% of the animals carry the deletion in Danish, Swedish and Finnish Red Cattle, respectively. To explain this, we searched for favorable effects on other traits and found that the deletion has strong positive effects on milk yield. This study demonstrates that embryonic lethal mutations account for a non-negligible fraction of the decline in fertility of domestic cattle, and that associated positive effects on milk yield may account for part of the negative genetic correlation. Our study adds to the evidence that structural variants contribute to animal phenotypic variation, and that balancing selection might be more common in livestock species than previously appreciated.

Mitogenomic meta-analysis identifies two phases of migration in the history of eastern Eurasian sheep

Despite much attention, history of sheep (Ovis aries) evolution, including its dating, demographic trajectory and geographic spread, remains controversial. To address these questions, we generated 45 complete and 875 partial mitogenomic sequences, and performed a meta-analysis of these and published ovine mitochondrial DNA sequences (n = 3,229) across Eurasia. We inferred that O. orientalis and O. musimon share the most recent female ancestor with O. aries at approximately 0.790 Ma (95% CI: 0.637–0.934 Ma) during the Middle Pleistocene, substantially predating the domestication event (∼8–11 ka). By reconstructing historical variations in effective population size, we found evidence of a rapid population increase approximately 20–60 ka, immediately before the Last Glacial Maximum. Analyses of lineage expansions showed two sheep migratory waves at approximately 4.5–6.8 ka (lineages A and B: ∼6.4–6.8 ka; C: ∼4.5 ka) across eastern Eurasia, which could have been influenced by prehistoric West–East commercial trade and deliberate mating of domestic and wild sheep, respectively. A continent-scale examination of lineage diversity and approximate Bayesian computation analyses indicated that the Mongolian Plateau region was a secondary center of dispersal, acting as a “transportation hub” in eastern Eurasia: Sheep from the Middle Eastern domestication center were inferred to have migrated through the Caucasus and Central Asia, and arrived in North and Southwest China (lineages A, B, and C) and the Indian subcontinent (lineages B and C) through this region. Our results provide new insights into sheep domestication, particularly with respect to origins and migrations to and from eastern Eurasia.

Mapping of fertility traits in Finnish Ayrshire by genome-wide association analysis

A whole-genome scan using single marker association was used to detect chromosome regions associated with seven female fertility traits in Finnish Ayrshire dairy cattle. The phenotypic data consisted of de-regressed estimated breeding values for 340 bulls which were estimated using a single trait model. Genotypes were obtained with the Illumina BovineSNP50 panel and a total of 35 630 informative, high-quality single nucleotide polymorphism (SNP) markers were used. The association analysis was performed using a mixed-model approach which fitted a fixed effect for each SNP and a random polygenic effect. We detected eleven genome-wide significant associations on eight different chromosomes. With at least chromosome-wise significance after Bonferroni correction, sixteen SNPs on nine chromosomes showed significant associations with one or more fertility traits. The results confirmed quantitative trait loci on three chromosomes (1, 2 and 20) for fertility traits previously reported for the same breed and one on chromosome four previously detected in Holstein cattle.

Fine-mapping QTL for mastitis resistance on BTA9 in three Nordic red cattle breeds.

A QTL affecting clinical mastitis and/or somatic cell score (SCS) has been reported previously on chromosome 9 from studies in 16 families from the Swedish Red and White (SRB), Finnish Ayrshire (FA) and Danish Red (DR) breeds. In order to refine the QTL location, 67 markers were genotyped over the whole chromosome in the 16 original families and 18 additional half-sib families. This enabled linkage disequilibrium information to be used in the analysis. Data were analysed by an approach that combines information from linkage and linkage disequilibrium, which allowed the QTL affecting clinical mastitis to be mapped to a small interval (<1 cM) between the markers BM4208 and INRA084. This QTL showed a pleiotropic effect on SCS in the DR and SRB breeds. Haplotypes associated with variations in mastitis resistance were identified. The haplotypes were predictive in the general population and can be used in marker-assisted selection. Pleiotropic effects of the mastitis QTL were studied for three milk production traits and eight udder conformation traits. This QTL was also associated with yield traits in DR but not in FA or SRB. No QTL were found for udder conformation traits on chromosome 9.

Genome-wide association analysis of milk yield traits in Nordic Red Cattle using imputed whole genome sequence variants.

BackgroundThe Nordic Red Cattle consisting of three different populations from Finland, Sweden and Denmark are under a joint breeding value estimation system. The long history of recording of production and health traits offers a great opportunity to study production traits and identify causal variants behind them. In this study, we used whole genome sequence level data from 4280 progeny tested Nordic Red Cattle bulls to scan the genome for loci affecting milk, fat and protein yields.ResultsUsing a genome-wise significance threshold, regions on Bos taurus chromosomes 5, 14, 23, 25 and 26 were associated with fat yield. Regions on chromosomes 5, 14, 16, 19, 20 and 25 were associated with milk yield and chromosomes 5, 14 and 25 had regions associated with protein yield. Significantly associated variations were found in 227 genes for fat yield, 72 genes for milk yield and 30 genes for protein yield. Ingenuity Pathway Analysis was used to identify networks connecting these genes displaying significant hits. When compared to previously mapped genomic regions associated with fertility, significantly associated variations were found in 5 genes common for fat yield and fertility, thus linking these two traits via biological networks.ConclusionThis is the first time when whole genome sequence data is utilized to study genomic regions affecting milk production in the Nordic Red Cattle population. Sequence level data offers the possibility to study quantitative traits in detail but still cannot unambiguously reveal which of the associated variations is causative. Linkage disequilibrium creates difficulties to pinpoint the causative genes and variations. One solution to overcome these difficulties is the identification of the functional gene networks and pathways to reveal important interacting genes as candidates for the observed effects. This information on target genomic regions may be exploited to improve genomic prediction.

Mitochondrial DNA and Y-chromosomal diversity in ancient populations of domestic sheep (Ovis aries) in Finland: comparison with contemporary sheep breeds

BackgroundSeveral molecular and population genetic studies have focused on the native sheep breeds of Finland. In this work, we investigated their ancestral sheep populations from Iron Age, Medieval and Post-Medieval periods by sequencing a partial mitochondrial DNA D-loop and the 5’-promoter region of the SRY gene. We compared the maternal (mitochondrial DNA haplotypes) and paternal (SNP oY 1) genetic diversity of ancient sheep in Finland with modern domestic sheep populations in Europe and Asia to study temporal changes in genetic variation and affinities between ancient and modern populations.ResultsA 523-bp mitochondrial DNA sequence was successfully amplified for 26 of 36 sheep ancient samples i.e. five, seven and 14 samples representative of Iron Age, Medieval and Post-Medieval sheep, respectively. Genetic diversity was analyzed within the cohorts. This ancient dataset was compared with present-day data consisting of 94 animals from 10 contemporary European breeds and with GenBank DNA sequence data to carry out a haplotype sharing analysis. Among the 18 ancient mitochondrial DNA haplotypes identified, 14 were present in the modern breeds. Ancient haplotypes were assigned to the highly divergent ovine haplogroups A and B, haplogroup B being the major lineage within the cohorts. Only two haplotypes were detected in the Iron Age samples, while the genetic diversity of the Medieval and Post-Medieval cohorts was higher. For three of the ancient DNA samples, Y-chromosome SRY gene sequences were amplified indicating that they originated from rams. The SRY gene of these three ancient ram samples contained SNP G-oY 1, which is frequent in modern north-European sheep breeds.ConclusionsOur study did not reveal any sign of major population replacement of native sheep in Finland since the Iron Age. Variations in the availability of archaeological remains may explain differences in genetic diversity estimates and patterns within the cohorts rather than demographic events that occurred in the past. Our ancient DNA results fit well with the genetic context of domestic sheep as determined by analyses of modern north-European sheep breeds.

Divergent evolution in the cytoplasmic domains of PRLR and GHR genes in Artiodactyla

BackgroundProlactin receptor (PRLR) and growth hormone receptor (GHR) belong to the large superfamily of class 1 cytokine receptors. Both of them have been identified as candidate genes affecting key quantitative traits, like growth and reproduction in livestock. We have previously studied the molecular anatomy of the cytoplasmic domain of GHR in different cattle breeds and artiodactyl species. In this study we have analysed the corresponding cytoplasmic signalling region of PRLR.ResultsWe sequenced PRLR gene exon 10, coding for the major part of the cytoplasmic domain, from cattle, American bison, European bison, yak, sheep, pig and wild boar individuals. We found different patterns of variation in the two receptors within and between ruminants and pigs. Pigs and bison species have no variation within GHR exon 10, but show high haplotype diversity for the PRLR exon 10. In cattle, PRLR shows lower diversity than GHR. The Bovinae PRLR haplotype network fits better the known phylogenetic relationships between the species than that of the GHR, where differences within cattle breeds are larger than between the different species in the subfamily. By comparison with the wild boar haplotypes, a high number of subsequent nonsynonymous substitutions seem to have accumulated in the pig PRLR exon 10 after domestication.ConclusionBoth genes affect a multitude of traits that have been targets of selection after domestication. The genes seem to have responded differently to different selection pressures imposed by human artificial selection. The results suggest possible effects of selective sweeps in GHR before domestication in the pig lineage or species divergence in the Bison lineage. The PRLR results may be explained by strong directional selection in pigs or functional switching.

Temporal Fluctuation in North East Baltic Sea Region Cattle Population Revealed by Mitochondrial and Y-Chromosomal DNA Analyses

Background Ancient DNA analysis offers a way to detect changes in populations over time. To date, most studies of ancient cattle have focused on their domestication in prehistory, while only a limited number of studies have analysed later periods. Conversely, the genetic structure of modern cattle populations is well known given the undertaking of several molecular and population genetic studies. Results Bones and teeth from ancient cattle populations from the North-East Baltic Sea region dated to the Prehistoric (Late Bronze and Iron Age, 5 samples), Medieval (14), and Post-Medieval (26) periods were investigated by sequencing 667 base pairs (bp) from the mitochondrial DNA (mtDNA) and 155 bp of intron 19 in the Y-chromosomal UTY gene. Comparison of maternal (mtDNA haplotypes) genetic diversity in ancient cattle (45 samples) with modern cattle populations in Europe and Asia (2094 samples) revealed 30 ancient mtDNA haplotypes, 24 of which were shared with modern breeds, while 6 were unique to the ancient samples. Of seven Y-chromosomal sequences determined from ancient samples, six were Y2 and one Y1 haplotype. Combined data including Swedish samples from the same periods (64 samples) was compared with the occurrence of Y-chromosomal haplotypes in modern cattle (1614 samples). Conclusions The diversity of haplogroups was highest in the Prehistoric samples, where many haplotypes were unique. The Medieval and Post-Medieval samples also show a high diversity with new haplotypes. Some of these haplotypes have become frequent in modern breeds in the Nordic Countries and North-Western Russia while other haplotypes have remained in only a few local breeds or seem to have been lost. A temporal shift in Y-chromosomal haplotypes from Y2 to Y1 was detected that corresponds with the appearance of new mtDNA haplotypes in the Medieval and Post-Medieval period. This suggests a replacement of the Prehistoric mtDNA and Y chromosomal haplotypes by new types of cattle.

Fine mapping of quantitative trait loci for mastitis resistance on bovine chromosome 11

Quantitative trait loci (QTL) affecting clinical mastitis (CM) and somatic cell score (SCS) were mapped on bovine chromosome 11. The mapping population consisted of 14 grandsire families belonging to three Nordic red cattle breeds: Finnish Ayrshire (FA), Swedish Red and White (SRB) and Danish Red. The families had previously been shown to segregate for udder health QTL. A total of 524 progeny tested bulls were included in the analysis. A linkage map including 33 microsatellite and five SNP markers was constructed. We performed combined linkage disequilibrium and linkage analysis (LDLA) using the whole data set. Further analyses were performed for FA and SRB separately to study the origin of the identified QTL/haplotype and to examine if it was common in both populations. Finally, different two-trait models were fitted. These postulated either a pleiotropic QTL affecting both traits; two linked QTL, each affecting one trait; or one QTL affecting a single trait. A QTL affecting CM was fine-mapped. In FA, a haplotype having a strong association with a high negative effect on mastitis resistance was identified. The mapping precision of an earlier detected SCS-QTL was not improved by the LDLA analysis because of lack of linkage disequilibrium between the markers used and the QTL in the region.