L1 insertion within SPEF2 gene is associated with increased litter size in the Finnish Yorkshire population

The prolactin receptor: A cross-species comparison of gene structure, transcriptional regulation, tissue-specificity, and genetic variation

The conserved and multifaceted functions of prolactin (PRL) are coordinated through varied distribution and expression of its cell-surface receptor (PRLR) across a range of tissues and physiological states. The resultant heterogeneous expression of PRLR mRNA and protein across different organs and cell types supports a wide range of PRL-regulated processes including reproduction, lactation, development, and homeostasis. Genetic variation within the PRLR gene also accounts for several phenotypes impacting agricultural production and human pathology. The goal of this review is to highlight the many elements that control differential expression of the PRLR across tissues, and the various phenotypes that exist across species due to variation in the PRLR gene.

Review: Balancing Selection for Deleterious Alleles in Livestock

Harmful alleles can be under balancing selection due to an interplay of artificial selection for the variant in heterozygotes and purifying selection against the variant in homozygotes. These pleiotropic variants can remain at moderate to high frequency expressing an advantage for favorable traits in heterozygotes, while harmful in homozygotes. The impact on the population and selection strength depends on the consequence of the variant both in heterozygotes and homozygotes. The deleterious phenotype expressed in homozygotes can range from early lethality to a slightly lower fitness in the population. In this review, we explore a range of causative variants under balancing selection including loss-of-function variation (i.e., frameshift, stop-gained variants) and regulatory variation (affecting gene expression). We report that harmful alleles often affect orthologous genes in different species, often influencing analogous traits. The recent discoveries are mainly driven by the increasing genomic and phenotypic resources in livestock populations. However, the low frequency and sometimes subtle effects in homozygotes prevent accurate mapping of such pleiotropic variants, which requires novel strategies to discover. After discovery, the selection strategy for deleterious variants under balancing selection is under debate, as variants can contribute to the heterosis effect in crossbred animals in various livestock species, compensating for the loss in purebred animals. Nevertheless, gene-assisted selection is a useful tool to decrease the frequency of the harmful allele in the population, if desired. Together, this review marks various deleterious variants under balancing selection and describing the functional consequences at the molecular, phenotypic, and population level, providing a resource for further study.

Search for Associations of FSHR, INHA, INHAB, PRL, TNP2 and SPEF2 Genes Polymorphisms with Semen Quality in Russian Holstein Bulls (Pilot Study)

The aim of the study was to search for new mutations in the previously studied gene loci of follicle-stimulating hormone receptor (FSHR), inhibin α (INHA), inhibin β A (INHAB), prolactin (PRL), transition protein 2 (TNP2), and sperm flagella 2 (SPEF2) by sequencing, as well as the search for associations of previously identified mutations at these loci with fresh semen quality in Russian Holstein bulls. Phenotypic data from 189 bulls was collected. Data was analyzed for most bulls for three years of semen collection. The maximum value of each semen quality indicator (doublet ejaculate volume, sperm concentration, progressive motility and total number of spermatozoa) were selected. SNPs were identified in the FSHR, INHA, INHAB, TNP2, SPEF2 genes. The PRL gene did not have polymorphism. Significant (p < 0.05) associations of polymorphisms in the FSHR gene with double ejaculate volume, concentration and total number of spermatozoa were identified. Polymorphism in the INHA gene was significantly associated (p < 0.05) with sperm concentration. Polymorphism in the INHAB gene was significantly associated (p < 0.05) with doublet ejaculate volume and total number of spermatozoa. Polymorphisms in the TNP2 and SPEF2 genes did not have significant associations with semen quality. The SNPs studied in our pilot work may be considered as candidate genetic markers in the selection of bulls.

Genome-Wide Detection of Structural Variations Reveals New Regions Associated with Domestication in Small Ruminants

During domestication processes, changes in selective pressures induce multiple phenotypical, physiological, and behavioral changes in target species. The rise of next-generation sequencing has provided a chance to study the genetics bases of these changes, most of the time based on single nucleotide polymorphisms (SNPs). However, several studies have highlighted the impact of structural variations (SVs) on individual fitness, particularly in domestic species. We aimed at unraveling the role of SVs during the domestication and later improvement of small ruminants by analyzing whole-genome sequences of 40 domestic sheep and 11 of their close wild relatives (Ovis orientalis), and 40 goats and 18 of their close wild relatives (Capra aegagrus). Using a combination of detection tools, we called 45,796 SVs in Ovis and 15,047 SVs in Capra genomes, including insertions, deletions, inversions, copy number variations, and chromosomal translocations. Most of these SVs were previously unreported in small ruminants. 69 and 45 SVs in sheep and goats, respectively, were in genomic regions with neighboring SNPs highly differentiated between wilds and domestics (i.e., putatively related to domestication). Among them, 25 and 20 SVs were close to or overlapping with genes related to physiological and morpho-anatomical traits linked with productivity (e.g., size, meat or milk quality, wool color), reproduction, or immunity. Finally, several of the SVs differentiated between wilds and domestics would not have been detected by screening only the differentiation of SNPs surrounding them, highlighting the complementarity of SVs and SNPs based approaches to detect signatures of selection.

SINE Insertion in the Intron of Pig GHR May Decrease Its Expression by Acting as a Repressor

Simple Summary

GH/IGF axis genes play a central role in the regulation of skeletal accretion during development and growth, and thus represent candidate genes for growth traits. Retrotransposon insertion polymorphisms are major contributors to structural variations. They tend to generate large effect mutations resulting in variations in target gene activity and phenotype due to the fact that they carry functional elements, such as enhancers, insulators, or promoters. In the present study, RIPs in four GH/IGF axis genes (GH, GHR, IGF1, and IGF1R) were investigated by comparative genomics and PCR. Four RIPs in the GHR gene and one RIP in the IGF1 gene were identified. Further analysis revealed that one RIP in the first intron of GHR might play a role in the regulation of GHR expression by acting as a repressor. These findings contribute to the understanding of the role of RIPs in the genetic variation of GH/IGF axis genes and phenotypic variation in pigs.

Abstract

The genetic diversity of the GH/IGF axis genes and their association with the variation of gene expression and phenotypic traits, principally represented by SNPs, have been extensively reported. Nevertheless, the impact of retrotransposon insertion polymorphisms (RIPs) on the GH/IGF axis gene activity has not been reported. In the present study, bioinformatic prediction and PCR verification were performed to screen RIPs in four GH/IGF axis genes (GH, GHR, IGF1 and IGF1R). In total, five RIPs, including one SINE RIP in intron 3 of IGF1, one L1 RIP in intron 7 of GHR, and three SINE RIPs in intron 1, intron 5 and intron 9 of GHR, were confirmed by PCR, displaying polymorphisms in diverse breeds. Dual luciferase reporter assay revealed that the SINE insertion in intron 1 of GHR significantly repressed the GHR promoter activity in PK15, Hela, C2C12 and 3T3-L1 cells. Furthermore, qPCR results confirmed that this SINE insertion was associated with a decreased expression of GHR in the leg muscle and longissimus dorsi, indicating that it may act as a repressor involved in the regulation of GHR expression. In summary, our data revealed that RIPs contribute to the genetic variation of GH/IGF axis genes, whereby one SINE RIP in the intron 1 of GHR may decrease the expression of GHR by acting as a repressor.

Genetic Diversity and Population Structures in Chinese Miniature Pigs Revealed by SINE Retrotransposon Insertion Polymorphisms, a New Type of Genetic Markers

Simple Summary

Our previous studies suggested that the short interspersed nuclear element (SINE) retrotransposon insertion polymorphisms (RIPs), as a new type of molecular marker developed very recently, are ideal molecular markers and have the potential to be used for population genetic analysis and molecular breeding in pigs and possibly it can be extended to other livestock animals as well. However, no report is available for the application of SINE RIPs in population genetic analysis in livestock, including pigs. Here, we evaluated 30 SINE RIPs in several indigenous Chinese miniature pig breeds, including three subpopulations of Bama pigs (BM-cov, BM-clo, and BM-inb). BM-cov is a subpopulation conserved in the national conservation farm, and BM-clo is a closed population maintained over 30 years with only 2 boars and 14 sows imported from its original area, while BM-inb herd is an 18 generation continuous inbreeding line based on the BM-clo population. To our knowledge, it is the first time to report the genetic diversity, breed differentiation, and population structures for these populations by using SINE RIPs, and which suggests the feasibility of SINE RIPs in pig genetic analysis.

Abstract

RIPs have been developed as effective genetic markers and popularly applied for genetic analysis in plants, but few reports are available for domestic animals. Here, we established 30 new molecular markers based on the SINE RIPs, and applied them for population genetic analysis in seven Chinese miniature pigs. The data revealed that the closed herd (BM-clo), inbreeding herd (BM-inb) of Bama miniature pigs were distinctly different from the BM-cov herds in the conservation farm, and other miniature pigs (Wuzhishan, Congjiang Xiang, Tibetan, and Mingguang small ear). These later five miniature pig breeds can further be classified into two clades based on a phylogenetic tree: one included BM-cov and Wuzhishan, the other included Congjiang Xiang, Tibetan, and Mingguang small ear, which was well-supported by structure analysis. The polymorphic information contents estimated by using SINE RIPs are lower than the predictions based on microsatellites. Overall, the genetic distances and breed-relationships between these populations revealed by 30 SINE RIPs generally agree with their evolutions and geographic distributions. We demonstrated the potential of SINE RIPs as new genetic markers for genetic monitoring and population structure analysis in pigs, which can even be extended to other livestock animals.

Genome-Wide Association Study and Transcriptome Differential Expression Analysis of the Feather Rate in Shouguang Chickens

The feather rate phenotype in chicks, including early-feathering and late-feathering phenotypes, are widely used as a sexing system in the poultry industry. The objective of this study was to obtain candidate genes associated with the feather rate in Shouguang chickens. In the present study, we collected 56 blood samples and 12 hair follicle samples of flight feathers from female Shouguang chickens. Then we identified the chromosome region associated with the feather rate by genome-wide association analysis (GWAS). We also performed RNA sequencing and analyzed differentially expressed genes between the early-feathering and late-feathering phenotypes using HISAT2, StringTie, and DESeq2. We identified a genomic region of 10.0–13.0 Mb of chromosome Z, which is statistically associated with the feather rate of Shouguang chickens at one-day old. After RNA sequencing analysis, 342 differentially expressed known genes between the early-feathering (EF) and late-feathering (LF) phenotypes were screened out, which were involved in epithelial cell differentiation, intermediate filament organization, protein serine kinase activity, peptidyl-serine phosphorylation, retinoic acid binding, and so on. The sperm flagellar 2 gene (SPEF2) and prolactin receptor (PRLR) gene were the only two overlapping genes between the results of GWAS and differential expression analysis, which implies that SPEF2 and PRLR are possible candidate genes for the formation of the chicken feathering phenotype in the present study. Our findings help to elucidate the molecular mechanism of the feather rate in chicks.

Biallelic mutations in Sperm flagellum 2 cause human multiple morphological abnormalities of the sperm flagella (MMAF) phenotype

Male patients with multiple morphological abnormalities of the sperm flagella (MMAF) are infertile and exhibit absent, short, coiled, bent and/or irregular sperm flagella. Mutations in the SPEF2 gene reduce sperm motility and cause sperm tail defects in animal models and humans. In the present study, we performed a genetic analysis on an MMAF patient and identified novel biallelic mutations in the SPEF2 gene. The biallelic mutations were confirmed by Sanger sequencing and in silico analysis revealed that, these variations were deleterious. The expression of truncated SPEF2 protein was reduced significantly in the patient's spermatozoa. The spermatozoa harbored biallelic mutations and showed severe ultrastructural defects in the axoneme and mitochondrial sheath. Our data suggest that biallelic mutations in SPEF2 can cause severe sperm flagellum defects, thus providing a novel candidate genetic pathogen for the human MMAF phenotype.

Retrotransposons evolution and impact on lncRNA and protein coding genes in pigs

BACKGROUND

Retrotransposons are the major determinants of genome sizes and they have shaped both genes and genomes in mammalian organisms, but their overall activity, diversity, and evolution dynamics, particularly their impact on protein coding and lncRNA genes in pigs remain largely unknown.

RESULTS

In the present study, we performed de novo detection of retrotransposons in pigs by using multiple pipelines, four distinct families of pig-specific L1 s classified into 51 distinct subfamilies and representing four evolution models and three expansion waves of pig-specific SINEs represented by three distinct families were identified. ERVs were classified into 18 families and found two most "modern" subfamilies in the pig genome. The transposition activity of pig L1 was verified by experiment, the sense and antisense promoter activities of young L1 5'UTRs and ERV LTRs and expression profiles of young retrotransposons in multiple tissues and cell lines were also validated. Furthermore, retrotransposons had an extensive impact on lncRNA and protein coding genes at both the genomic and transcriptomic levels. Most protein coding and lncRNA (> 80%) genes contained retrotransposon insertions, and about half of protein coding genes (44.30%) and one-fourth (24.13%) of lncRNA genes contained the youngest retrotransposon insertions. Nearly half of protein coding genes (43.78%) could generate chimeric transcripts with retrotransposons. Significant distribution bias of retrotransposon composition, location, and orientation in lncRNA and protein coding genes, and their transcripts, were observed.

CONCLUSIONS

In the current study, we characterized the classification and evolution profile of retrotransposons in pigs, experimentally proved the transposition activity of the young pig L1 subfamily, characterized the sense and antisense expression profiles and promoter activities of young retrotransposons, and investigated their impact on lncRNA and protein coding genes by defining the mobilome landscapes at the genomic and transcriptomic levels. These findings help provide a better understanding of retrotransposon evolution in mammal and their impact on the genome and transcriptome.

Balancing selection on a recessive lethal deletion with pleiotropic effects on two neighboring genes in the porcine genome

Livestock populations can be used to study recessive defects caused by deleterious alleles. The frequency of deleterious alleles including recessive lethal alleles can stay at high or moderate frequency within a population, especially if recessive lethal alleles exhibit an advantage for favourable traits in heterozygotes. In this study, we report such a recessive lethal deletion of 212kb (del) within the BBS9 gene in a breeding population of pigs. The deletion produces a truncated BBS9 protein expected to cause a complete loss-of-function, and we find a reduction of approximately 20% on the total number of piglets born from carrier by carrier matings. Homozygous del/del animals die mid- to late-gestation, as observed from high increase in numbers of mummified piglets resulting from carrier-by-carrier crosses. The moderate 10.8% carrier frequency (5.4% allele frequency) in this pig population suggests an advantage on a favourable trait in heterozygotes. Indeed, heterozygous carriers exhibit increased growth rate, an important selection trait in pig breeding. Increased growth and appetite together with a lower birth weight for carriers of the BBS9 null allele in pigs is analogous to the phenotype described in human and mouse for (naturally occurring) BBS9 null-mutants. We show that fetal death, however, is induced by reduced expression of the downstream BMPER gene, an essential gene for normal foetal development. In conclusion, this study describes a lethal 212kb deletion with pleiotropic effects on two different genes, one resulting in fetal death in homozygous state (BMPER), and the other increasing growth (BBS9) in heterozygous state. We provide strong evidence for balancing selection resulting in an unexpected high frequency of a lethal allele in the population. This study shows that the large amounts of genomic and phenotypic data routinely generated in modern commercial breeding programs deliver a powerful tool to monitor and control lethal alleles much more efficiently.

We report a large deletion within the BBS9 gene that induces late fetal mortality in homozygous affected animals in a commercial pig population. This late fetal mortality causes the fetus to become encapsulated and desiccated during the remaining time of the pregnancy, a process called mummification. The unusually high carrier frequency for this lethal deletion (10.8%) likely results from its strong positive association with growth rate in heterozygous individuals, an important selection trait in the pig breeding industry. Interestingly, we show that the positive effect on growth is induced by a heterozygous loss-of-function of the BBS9 gene, associated with obesity in human and mouse. However, late fetal mortality is induced by insufficient expression of the BMPER gene located directly downstream of the deletion which affects its regulatory elements required for gene expression. Together, our study shows an unique example of allelic pleiotropy in which one allele (deletion) is responsible for both increased growth and late fetal mortality by affecting two different genes.

Intraflagellar transport protein IFT20 is essential for male fertility and spermiogenesis in mice

Intraflagellar transport (IFT) is a conserved mechanism thought to be essential for the assembly and maintenance of cilia and flagella. However, little is known about its role in mammalian sperm flagella formation. To fill this gap, we disrupted the Ift20 gene in male germ cells. Homozygous mutant mice were infertile with significantly reduced sperm counts and motility. In addition, abnormally shaped elongating spermatid heads and bulbous round spermatids were found in the lumen of the seminiferous tubules. Electron microscopy revealed increased cytoplasmic vesicles, fiber-like structures, abnormal accumulation of mitochondria and a decrease in mature lysosomes. The few developed sperm had disrupted axonemes and some retained cytoplasmic lobe components on the flagella. ODF2 and SPAG16L, two sperm flagella proteins failed to be incorporated into sperm tails of the mutant mice, and in the germ cells, both were assembled into complexes with lighter density in the absence of IFT20. Disrupting IFT20 did not significantly change expression levels of IFT88, a component of IFT-B complex, and IFT140, a component of IFT-A complex. Even though the expression level of an autophagy core protein that associates with IFT20, ATG16, was reduced in the testis of the Ift20 mutant mice, expression levels of other major autophagy markers, including LC3 and ubiquitin were not changed. Our studies suggest that IFT20 is essential for male fertility and spermiogenesis in mice, and its major function is to transport cargo proteins for sperm flagella formation. It also appears to be involved in removing excess cytoplasmic components.

Detection of one large insertion/deletion (indel) and two novel SNPs within the <i>SPEF2</i> gene and their associations with male piglet reproduction traits

Abstract. The sperm flagella 2 (SPEF2) gene is essential for normal sperm tail development and male fertility. To fully characterize the structure of the mutation and to further study the function of the pig SPEF2 gene, we explored the insertion/deletion (indel) and novel single-nucleotide polymorphisms (SNPs) within the pig SPEF2 gene, and tested their associations with the testicular sizes in male Large White (LW) and Landrace (LD) pigs from China. Herein, a large insertion located at the SPEF2 gene in chromosome 16 was found, and two alleles of "I" (insertion) and "D" (deletion) were designated. Allele "D" was dominant in all analyzed pigs. Two novel SNPs (namely (NC_010458) g.19642G > A, resulting in AfaI aCRS PCR–PFLP, and g.19886C > G, resulting in EcoRI aCRS PCR–PFLP) were found in LW and LD pigs. Association testing revealed that g.19886C > G was significantly associated with the testis long circumference (TLC) in LW pigs (P < 0.05), suggesting that this SNP would be the DNA marker for the marker-assisted selection (MAS) in reproduction traits. This preliminary result indicates that the pig SPEF2 gene had significant effects on male reproduction traits. These findings could not only extend the spectrum of genetic variations in the pig SPEF2 gene but also contribute to implementing MAS in genetics and breeding in pigs.

Identification of candidate genes for chicken early- and late-feathering

Previous studies suggest that prolactin receptor (Prlr) is a potential causative gene for chicken early- (EF) and late-feathering (LF) phenotypes. In this study, we evaluated candidate genes for this trait and determined the expression of 3 genes, including Prlr, sperm flagellar protein 2 (Spef2), and their fusion gene, in the skins of one-day-old EF and LF chicks using RT-qPCR. Data indicated that Prlr expression in the skin did not show significant difference between EF and LF chicks, suggesting Prlr may not be a suitable candidate gene. In contrast, Spef2 expression in the skin displayed a significant difference between EF and LF chicks (P < 0.01), suggesting that Spef2 may be a good candidate gene for chicken feathering. Moreover, dPrlr/dSpef2, the fusion gene, was also a good candidate gene as it was expressed only in LF chicks. However, the expression of the fusion gene was much lower than that of Prlr Additionally, using strand-specific primers, we found that the fusion gene was transcribed in 2 directions (one from dPrlr promoter, another from dSpef2 promoter), which could result in the formation of a double strand RNA. In conclusion, both Spef2 and the fusion gene are good candidate genes for chicken feathering, but Prlr is not. The research on the function and regulation of the candidate genes will help elucidate the molecular basis of the chicken feathering trait.

Identification of Reproduction-Related Gene Polymorphisms Using Whole Transcriptome Sequencing in the Large White Pig Population

Recent developments in high-throughput sequencing techniques have enabled large-scale analysis of genetic variations and gene expression in different tissues and species, but gene expression patterns and genetic variations in livestock are not well-characterized. In this study, we have used high-throughput transcriptomic sequencing of the Finnish Large White to identify gene expression patterns and coding polymorphisms within the breed in the testis and oviduct. The main objective of this study was to identify polymorphisms within genes that are highly and specifically expressed in male and/or female reproductive organs. The differential expression (DE) analysis underlined 1234 genes highly expressed in the testis and 1501 in the oviduct. Furthermore, we used a novel in-house R-package hoardeR for the identification of novel genes and their orthologs, which underlined 55 additional DE genes based on orthologs in the human, cow, and sheep. Identification of polymorphisms in the dataset resulted in a total of 29,973 variants, of which 10,704 were known coding variants. Fifty-seven nonsynonymous SNPs were present among genes with high expression in the testis and 67 were present in the oviduct, underlining possible influential genes for reproduction traits. Seven genes (PGR, FRAS1, TCF4, ADAT1, SPAG6, PIWIL2, and DNAH8) with polymorphisms were highlighted as reproduction-related based on their biological function. The expression and SNPs of these genes were confirmed using RT-PCR and Sanger sequencing. The identified nonsynonymous mutations within genes highly expressed in the testis or oviduct provide a list of candidate genes for reproduction traits within the pig population and enable identification of biomarkers for sow and boar fertility.

Genome-wide association study and ancestral origins of the slick-hair coat in tropically adapted cattle

The slick hair coat (SLICK) is a dominantly inherited trait typically associated with tropically adapted cattle that are from Criollo descent through Spanish colonization of cattle into the New World. The trait is of interest relative to climate change, due to its association with improved thermo-tolerance and subsequent increased productivity. Previous studies localized the SLICK locus to a 4 cM region on chromosome (BTA) 20 and identified signatures of selection in this region derived from Senepol cattle. The current study compares three slick-haired Criollo-derived breeds including Senepol, Carora, and Romosinuano and three additional slick-haired cross-bred lineages to non-slick ancestral breeds. Genome-wide association (GWA), haplotype analysis, signatures of selection, runs of homozygosity (ROH), and identity by state (IBS) calculations were used to identify a 0.8 Mb (37.7-38.5 Mb) consensus region for the SLICK locus on BTA20 in which contains SKP2 and SPEF2 as possible candidate genes. Three specific haplotype patterns are identified in slick individuals, all with zero frequency in non-slick individuals. Admixture analysis identified common genetic patterns between the three slick breeds at the SLICK locus. Principal component analysis (PCA) and admixture results show Senepol and Romosinuano sharing a higher degree of genetic similarity to one another with a much lesser degree of similarity to Carora. Variation in GWA, haplotype analysis, and IBS calculations with accompanying population structure information supports potentially two mutations, one common to Senepol and Romosinuano and another in Carora, effecting genes contained within our refined location for the SLICK locus.

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.

We report the identification of a large deletion encompassing four genes and the demonstration of its negative effect on fertility in Nordic Red dairy cattle. We show that this deletion is recessively lethal (homozygous embryos die) and therefore, when carrier cows are mated to carrier bulls, there is a high risk of embryonic mortality. As a result, chances of insemination failure are higher for such matings. Surprisingly, despite its negative effect, the deletion is frequent in Nordic Red cattle. We show that this high frequency may be a consequence of the fact that the deletion is associated with increased milk production and therefore selected for. Due to increased levels of inbreeding resulting from the widespread use of artificial insemination, such recessive lethal alleles may account for a non-negligible fraction of the reduction in fertility observed in cattle.