Search

As human complex diseases are influenced by the interaction between genetics and the environment, identifying gene–environment interactions (⁠ G × E ⁠) is crucial for understanding disease mechanisms and predicting risk. Developing robust quantitative tools for G × E analysis can enhance the study of complex diseases. However, many existing methods that explore G × E focus on the interplay between an environmental factor and genetic variants, exclusively for common or rare variants. In this study, we developed MAGEIT_RAN and MAGEIT_FIX to identify interactions between an environmental factor and a set of genetic markers, including both rare and common variants, based on the MinQue for Summary statistics. The genetic main effects in MAGEIT_RAN and MAGEIT_FIX are modeled as random and fixed effects, respectively. Simulation studies showed that both tests had type I error under control, with MAGEIT_RAN being the most powerful test. Applying MAGEIT to a genome-wide analysis of gene–alcohol interactions on hypertension and seated systolic blood pressure in the Multiethnic Study of Atherosclerosis revealed genes like EIF2AK2 , CCNDBP1 , and EPB42 influencing blood pressure through alcohol interaction. Pathway analysis identified 1 apoptosis and survival pathway involving PKR and 2 signal transduction pathways associated with hypertension and alcohol intake, demonstrating MAGEIT_RAN's ability to detect biologically relevant gene–environment interactions. [ABSTRACT FROM AUTHOR]

Background: Nitric oxide (NO) is pivotal in regulating the activity of NBS-LRR specific R genes, crucial components of the plant’s immune system. It is noteworthy that previous research has not included a genome-wide analysis of NO-responsive NBS-LRR genes in plants. Results: The current study examined 29 NO-induced NBS-LRR genes from Arabidopsis thaliana, along with two monocots (rice and maize) and two dicots (soybean and tomato) using genome-wide analysis tools. These NBS-LRR genes were subjected to comprehensive characterization, including analysis of their physio-chemical properties, phylogenetic relationships, domain and motif identification, exon/intron structures, cis-elements, protein-protein interactions, prediction of S-Nitrosylation sites, and comparison of transcriptomic and qRT-PCR data. Results showed the diverse distribution of NBS-LRR genes across chromosomes, and variations in amino acid number, exons/introns, molecular weight, and theoretical isoelectric point, and they were found in various cellular locations like the plasma membrane, cytoplasm, and nucleus. These genes predominantly harbor the NB-ARC superfamily, LRR, LRR_8, and TIR domains, as also confirmed by motif analysis. Additionally, they feature species-specific PLN00113 superfamily and RX-CC_like domain in dicots and monocots, respectively, both responsive to defense against pathogen attacks. The NO-induced NBS-LRR genes of Arabidopsis reveal the presence of cis-elements responsive to phytohormones, light, stress, and growth, suggesting a wide range of responses mediated by NO. Protein-protein interactions, coupled with the prediction of S-Nitrosylation sites, offer valuable insights into the regulatory role of NO at the protein level within each respective species. Conclusion: These above findings aimed to provide a thorough understanding of the impact of NO on NBS-LRR genes and their relationships with key plant species. [ABSTRACT FROM AUTHOR]

Background: Type 2 diabetes is one of the leading causes of chronic kidney failure, which increases globally and represents a significant threat to public health. People from the Middle East represent one of the largest immigrant groups in Europe today. Despite poor glucose regulation and high risk for early-onset insulin-deficient type 2 diabetes, they have better kidney function and lower rates of all-cause and cardiovascular-specific mortality compared with people of European ancestry. Here, we assessed the genetic basis of estimated glomerular filtration rate (eGFR) and other metabolic traits in people of Iraqi ancestry living in southern Sweden.Methods: Genome-wide association study (GWAS) analyses were performed in 1201 Iraqi-born residents of the city of Malmö for eGFR and ten other metabolic traits using linear mixed-models to account for family structure.Results: The strongest association signal was detected for eGFR in CST9 (rs13037490; P value = 2.4 × 10-13), a locus previously associated with cystatin C-based eGFR; importantly, the effect (major) allele here contrasts the effect (minor) allele in other populations, suggesting favorable selection at this locus. Additional novel genome-wide significant loci for eGFR (ERBB4), fasting glucose (CAMTA1, NDUFA10, TRIO, WWC1, TRAPPC9, SH3GL2, ABCC11), quantitative insulin-sensitivity check index (METTL16), and HbA1C (CAMTA1, ME1, PAK1, RORA) were identified.Conclusions: The genetic effects discovered here may help explain why people from the Middle East have better kidney function than those of European descent. Genetic predisposition to preserved kidney function may also underlie the observed survival benefits in Middle Eastern immigrants with type 2 diabetes. [ABSTRACT FROM AUTHOR]

Key message: Genome-wide identification, classification, functional characterization and expression analysis of Auxin Responsive Factor (ARF) gene family in wheat reveal their attributes and role during leaf rust infection. Auxins are important plant growth regulators that also impact plant-pathogen interaction. Auxin responsive factors (ARF) are plant specific transcription factors that control responses to auxins. Whole genome investigation of ARF gene family is limited in allohexaploid wheat (Triticum aestivum L.). Comprehensive study of this gene family was carried out by employing the currently available reference genome sequence of wheat. In total, 27 ARF genes were identified and located on the wheat genome as well as were positioned on wheat chromosome arms. Additionally, examination of the predicted genes unveiled a decent degree of relatedness within and among the phylogenetic clades. Leaf rust, caused by the obligate biotrophic fungal pathogen Puccinia triticina, is responsible for drastic loss of wheat crop worldwide reducing grain yield by 10–90%. Expression profiling of ARF genes in retort to leaf rust infection indicated their differential regulation during this plant-pathogen interaction. Highest expression of ARF genes were observed at 12 hpi that was maintained up to 72 hpi during incompatible interaction, whereas the high expression levels receded at 48 hpi during compatible interactions. Few of the identified ARF genes were likely to be post-transcriptionally regulated by microRNAs. Many light and stress responsive elements were detected in the promoter regions of ARF genes. Microsynteny analysis showed the conservation of ARF genes within the members of the Poaceae family. This study provides fundamental details for understanding the different types of ARF genes in wheat and there putative roles during leaf rust–wheat interaction. [ABSTRACT FROM AUTHOR]

Purpose: To synthesise data from genome-wide studies reporting molecular signature of eutopic endometrium through the phases of the menstrual cycle in endometriosis. Methods: Extraction of data from publications reporting genetic signatures characterising endometrium associated with endometriosis. The nomenclature of extracted differentially expressed transcripts and proteins was adopted according to the HUGO Gene Nomenclature Committee (HGNC). Loci were further sorted according to the different phases of the menstrual cycle, i.e. menstrual (M), proliferative (P), secretory (S), early-secretory (ES), mid-secretory (MS), late-secretory (LS), and not specified (N/S) if the endometrial dating was not available. Enrichment analysis was performed using the DAVID bioinformatics tool. Results: Altered molecular changes were reported by 21 studies, including 13 performed at the transcriptomic, 6 at proteomic, and 2 at epigenomic level. Extracted data resulted in a catalogue of total 670 genetic causes with available 591 official gene symbols, i.e. M = 3, P = 188, S = 81, ES = 82, MS = 173, LS = 36, and N/S = 28. Enriched pathways included oestrogen signalling pathway, extracellular matrix organization, and endothelial cell chemotaxis. Our study revealed that knowledge of endometrium biology in endometriosis is fragmented due to heterogeneity of published data. However, 15 genes reported as dysregulated by at least two studies within the same phase and 33 significantly enriched GO-BP terms/KEGG pathways associated with different phases of the menstrual cycle were identified. Conclusions: A multi-omics insight into molecular patterns underlying endometriosis could contribute towards identification of endometrial pathological mechanisms that impact fertility capacities of women with endometriosis. [ABSTRACT FROM AUTHOR]

The objectives of this study were to identify genomic regions by Bayesian methods (BayesA, BayesB, or BayesN) that fit fixed-length haplotypes or SNPs using GenSel. Covariates for haplo-type alleles of five lengths (125, 250, 500 kb, 1 or 2 Mb) were generated, and rare haplotypes were removed at three thresholds (1, 5, or 10%). Subsequently, we performed gene network analyses to investigate the biological processes shared by genes that were identified for the same across traits. Genotypes at 41,034 SNPs that were common on OvineSNP50 panel were phased for 751 Scottish Blackface (SBF) lambs. This is the first study to quantify the proportion of genetic variance using haplotypes across the whole genome in an SBF population. The genetic variance explained of haplotype-based GWAS was higher than that of SNP-based GWAS in across traits studied. In this population, fitting 500kb haplotypes with a 1% frequency threshold resulted in the highest proportion of genetic variance explained for nematode resistance and fitting 2Mb haplotypes with a 10% frequency threshold improved genetic variance explained for body weight comparable to fitting SNPs by BayesB. Candidate genes, including CXCR4, STAT4, CCL1, CCL2, CCL3, CCL5, CCL8, CCL16, CCL18, CARMIL2, and HSPA14 were identified for nematode resistance and ADH5, PPP3CA, and FABP4 for body weight traits. Network analysis provided annotation results linking to all identified candidate genes. This study supported previous results from GWAS of nematode resistance and body weight and revealed additional regions in the ovine genome associated with these economically important traits. These results suggest that network analysis can provide new information regarding biological mechanisms and genes leading to complex phenotypes, like nematode resistance and body weight of lamb. [ABSTRACT FROM AUTHOR]

Penicillium expansum is a necrotrophic plant pathogen that has a broad range of fruit hosts and is responsible for blue mold rot. This rot has a significant impact on the fruit industry, causing economic losses during storage, transport, and sale. Polygalacturonase (PGs) enzymes are essential for pathogenic fungi to break down pectin, as they hydrolyze the polygalacturonic acid chain along the oxygen bridge. This enzyme has been identified in decayed tissue and is linked to multiple soft rot diseases. This study encompasses a thorough examination of the PGs gene family within the genome of P. expansum , utilizing a genome-wide analysis approach. As a result, a total of seven PG genes were successfully found. These genes were phylogenetically divided into four distinct clades. The gene structural characteristics of each group exhibited a fundamental level of conservation. Furthermore, the gene architectures and motif compositions exhibited robust evidence in favor of the credibility of the phylogenetic relationship. Motif analysis reveals that most of the motifs predicted had Glycosyl hydrolases family 28 conserved domains (PF00295). Analysis of RNA-seq data and the PG genes identified by RT-qPCR were mostly upregulated. Hence, comprehending the characteristics of PePGs holds significant scientific importance and will offer comprehensive insights into the PG superfamily in P. expansum and will make a valuable contribution to the ongoing validation of functional genes. • Genome-wide analysis reveals seven PG gene with GH28 domain in Penicillium expansum genome. • The seven PG genes are subgroup into four families and distributed on seven chromosomes. • Phylogenetic analysis reveals that P. expansum PGs are closely related to P. italicum. • RT-qPCR analysis confirms that most PG genes in P. expansum genome were upregulated. [ABSTRACT FROM AUTHOR]

Atherosclerosis develops with age and is partially controlled by genetics. Research to date has identified common variants with small effects on atherosclerosis related traits. We aimed to use family-based genome-wide linkage analysis to identify chromosomal regions potentially harboring rare variants with larger effects for atherosclerosis related traits. Participants included 2205 individuals from the Long Life Family Study (LLFS), which recruited families with exceptional longevity from Boston, New York, Pittsburgh, and Denmark. Participants underwent B-mode ultrasonography of the carotid arteries to measure intima-media thickness (IMT), inter-adventitial diameter (IAD), and plaque presence and severity. We conducted residual heritability and genome-wide linkage analyses adjusted for age, age2, sex, and field center using pedigree-based maximum-likelihood methods in SOLAR. All carotid traits were significantly heritable with a range of 0.68 for IAD to 0.38 for IMT. We identified three chromosomal regions with linkage to IAD (3q13; max LOD 5.3), plaque severity (17q22-q23, max LOD 3.2), and plaque presence (17q24, max LOD 3.1). No common allelic variants within these linkage peaks were associated with the carotid artery traits. We identified three chromosomal regions with evidence of linkage to carotid artery diameter and atherosclerotic plaque in exceptionally long-lived families. Since common allelic variants within our linkage peaks did not account for our findings, future follow-up resequencing of these regions in LLFS families should help advance our understanding of atherosclerosis, CVD, and healthy vascular aging. Image 1 • Measures of carotid atherosclerosis assessed via ultrasound were significantly heritable. • We identified 3 discrete chromosomal regions with significant linkage to carotid atherosclerosis. • These regions contain genes with previous evidence of association with atherosclerotic traits. • Common SNPs in these regions were not associated with atherosclerosis; thus, these families may harbor rare genetic variants. [ABSTRACT FROM AUTHOR]

Cytochrome P450s (CYPs) constitute the largest group of enzymes in plants and are involved in a variety of processes related to growth and protection. However, the CYP gene superfamily in pear (Pyrus bretschneideri) and their characteristics is unclear. Through a comprehensive genome-wide analysis, this article identified a total of 74 CYP genes in the P. bretschneideri genome, which were categorized into fourteen families. Motif analysis reveals that most of the ten motifs predicted were with the p450 conserved domain. The majority of the CYP genes have exon arrangements. Furthermore, promoter analysis unveiled a multitude of cis-acting elements associated with diverse responsiveness including hormones, light responsive, anoxic specific inducibility and anaerobic induction. Analysis of the transcriptome data reveal that about 80% of the pear CYPs genes were upregulated and they were positively correlated with the antioxidant's parameters such as total flavonoids and total phenol content as well as ABTS and DPPH radicals. RT-qPCR analysis confirmed that the CYP genes could be regulated in pear. Collectively, our results reveal comprehensive insights into the CYP superfamily in pear and make a valuable contribution to the ongoing process of functional validation. • The enhanced Sporidiobolus pararoseus Y16 induced several CYPs genes in pear. • Seventy-four CYP genes were identified in Pyrus. bretschneideri genome and grouped into 14 families. • Pear CYPs genes were upregulated and correlated positively with the antioxidant's parameters. • The gene structure is composed of mainly two exon organization and with p450 conserved domain. [ABSTRACT FROM AUTHOR]

Idiosyncratic drug-induced hepatotoxicity is a formidable challenge for rational drug discovery and development, as well as the science of personalized medicine. There is evidence that hereditary factors, in part, contribute to drug toxicity. This expert analysis and review offer the insights gained, and the challenges ahead, for genome-wide association studies (GWASs) of idiosyncratic drug-induced hepatotoxicity. Published articles on genome-wide and subsequent replication studies were systematically searched in the PubMed electronic database. We found that the genetic risk variants that were identified genome-wide, and replication confirmed, are mainly related to polymorphisms in the human leukocyte antigen (HLA) region that include HLA-DQB1*06:02 for amoxicillin-clavulanate, HLA-B*57:01 for flucloxacillin, HLA-DRB1*15:01 for lumiracoxib, and HLA-DRB1*07:01 for lapatinib and ximelagatran-induced hepatotoxicity. Additionally, polymorphisms in ST6 β-galactosamide α-2, 6-sialyltranferase-1 ( ST6GAL1), which plays a role in systemic inflammatory response, and variants in intron of family with sequence similarity-65 member-B ( FAM65B) that play roles in liver inflammation displayed association with flucloxacillin and antituberculosis drug-induced hepatotoxicity, respectively. Taken together, these GWAS findings offer molecular leads on the central role that the immune system plays in idiosyncratic drug-induced hepatotoxicity. We conclude the expert review with a brief discussion of the salient challenges ahead. These include, for example, the need for discursive discovery paradigms that incorporate alternating GWASs and candidate gene studies, as well as the study of the environtome, the entire complement of environmental factors, including science and innovation policies that enact on global society and the human host, and by extension, on susceptibility for idiosyncratic drug-induced hepatotoxicity. [ABSTRACT FROM AUTHOR]

Precursor lesions that progress into colorectal cancer (CRC) could be largely classified into sessile serrated lesions (SSLs), traditional serrated adenoma (TSA), and tubular adenoma (TA). We aimed to determine whether high expression of trefoil factor 1 (TFF1) is closely associated with serrated lesions, particularly SSLs. The samples were divided into the first (12 SSLs, 5 TSAs, and 15 TAs) and second cohorts (15 SSLs, 9 TSAs, and 15 TAs). First, we investigated TFF1 expression in isolated gland samples using array-based and reverse-transcription PCR. Second, we performed immunohistochemical analysis of TFF1 expression in paraffin-embedded tissues obtained from SSL, TSA, TA, and hyperplastic polyp (HP) samples. In addition, we compared TFF1 mRNA levels between SSLs and HPs. TFF1 expression was significantly higher in SSLs than in TSA and TA in both cohorts. Additionally, immunohistochemical staining of TFF1 in the HP, SSL, TSA, and TA samples revealed significant differences in the immunohistochemical scores of TFF1 among the four types of lesions (higher expression in SSLs than in the other three lesions). Finally, there were significant differences in TFF1 mRNA expression levels between SSLs and HPs in paraffin-embedded tissues. However, there was considerable overlap in the immunohistochemical scores and expression levels of TFF1 transcripts between SSLs and HPs. The current findings may help elucidate the molecular mechanisms involved in serrated lesion development. In addition, we suggest that despite the limited practical application, upregulation of TFF1 transcripts may help differentiate SSLs from other lesions. • We revealed higher expression of trefoil factor 1 (TFF1) in SSLs than in traditional serrated adenomas or tubular adenomas, suggesting that TTF1 promotes the development of SSLs. • Using paraffin-embedded samples, we found higher TFF1 mRNA levels in SSLs than in hyperplastic polyps. • We suggest that upregulation of TTF1 expression is closely associated with the development of serrated lesions, particularly SSLs. [ABSTRACT FROM AUTHOR]

Powerful array-based single-nucleotide polymorphism-typing platforms have recently heralded a new era in which genome-wide studies are conducted with increasing frequency. A genetic polymorphism associated with population pharmacokinetics (PK) is typically analyzed using nonlinear mixed-effect models (NLMM). Applying NLMM to large-scale data, such as those generated by genome-wide studies, raises several issues related to the assumption of random effects as follows: (i) computation time: it takes a long time to compute the marginal likelihood; (ii) convergence of iterative calculation: an adaptive Gauss-Hermite quadrature is generally used to estimate NLMM; however, iterative calculations may not converge in complex models; and (iii) random-effects misspecification leads to slightly inflated type-I error rates. As an alternative effective approach to resolving these issues, in this article, we propose a generalized estimating equation (GEE) approach for analyzing population PK data. In general, GEE analysis does not account for interindividual variability in PK parameters; therefore, the usual GEE estimators cannot be interpreted straightforwardly, and their validities have not been justified. Here, we propose valid inference methods for using GEE even under conditions of interindividual variability and provide theoretical justifications of the proposed GEE estimators for population PK data. In numerical evaluations by simulations, the proposed GEE approach exhibited high computational speed and stability relative to the NLMM approach. Furthermore, the NLMM analysis was sensitive to the misspecification of the random-effects distribution, and the proposed GEE inference is valid for any distributional form. We provided an illustration by using data from a genome-wide pharmacogenomic study of an anticancer drug. Copyright © 2013 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

Background: Congenital heart defects (CHD) is the most common cause of death from a congenital structure abnormality in newborns and is often associated with fetal loss. There are many types of CHD. Human genetic studies have identified genes that are responsible for the inheritance of a particular type of CHD and for some types of CHD previously thought to be sporadic. However, occasionally different members of the same family might have anatomically distinct defects - for instance, one member with atrial septal defect, one with tetralogy of Fallot, and one with ventricular septal defect. Our objective is to identify susceptibility loci for CHD in families affected by distinct defects. The occurrence of these apparently discordant clinical phenotypes within one family might hint at a genetic framework common to most types of CHD. Results: We performed a genome-wide linkage analysis using MOD score analysis in families with diverse CHD. Significant linkage was obtained in two regions, at chromosome 15 (15q26.3, Pempirical= 0.0004) and at chromosome 18 (18q21.2, Pempirical = 0.0005). Conclusions: In these two novel regions four candidate genes are located: SELS, SNRPA1, and PCSK6 on 15q26.3, and TCF4 on 18q21.2. The new loci reported here have not previously been described in connection with CHD. Although further studies in other cohorts are needed to confirm these findings, the results presented here together with recent insight into how the heart normally develops will improve the understanding of CHD. [ABSTRACT FROM AUTHOR]

Background: Malaria is caused by Plasmodium parasites, which are transmitted via the bites of infected Anopheline mosquitoes. Midgut invasion is a major bottleneck for Plasmodium development inside the mosquito vectors. Malaria parasites in the midgut are surrounded by a hostile environment rich in digestive enzymes, while a rapidly responding immune system recognizes Plasmodium ookinetes and recruits killing factors from the midgut and surrounding tissues, dramatically reducing the population of invading ookinetes before they can successfully traverse the midgut epithelium. Understanding molecular details of the parasite-vector interactions requires precise measurement of nascent protein synthesis in the mosquito during Plasmodium infection. Current expression profiling primarily monitors alterations in steady-state levels of mRNA, but does not address the equally critical issue of whether the proteins encoded by the mRNAs are actually synthesized. Results: In this study, we used sucrose density gradient centrifugation to isolate actively translating Anopheles gambiae mRNAs based upon their association with polyribosomes (polysomes). The proportion of individual gene transcripts associated with polysomes, which is determined by RNA deep sequencing, reflects mRNA translational status. This approach led to identification of 1017 mosquito transcripts that were primarily regulated at the translational level after ingestion of Plasmodium falciparum-infected blood. Caspar, a negative regulator of the NFkappaB transcription factor Rel2, appears to be substantially activated at the translational levels during Plasmodium infection. In addition, transcripts of Dcr1, Dcr2 and Drosha, which are involved in small RNA biosynthesis, exhibited enhanced associations with polysomes after P. falciparum challenge. This observation suggests that mosquito microRNAs may play an important role in reactions against Plasmodium invasion. Conclusions: We analyzed both total cellular mRNAs and mRNAs that are associated with polysomes to simultaneously monitor transcriptomes and nascent protein synthesis in the mosquito. This approach provides more accurate information regarding the rate of protein synthesis, and identifies some mosquito factors that might have gone unrecognized because expression of these proteins is regulated mainly at the translational level rather than at the transcriptional level after mosquitoes ingest a Plasmodium-infected blood meal. [ABSTRACT FROM AUTHOR]

Biomedical researchers usually test the null hypothesis that there is no difference of the population mean of pharmacokinetics (PK) parameters between genotypes by the Kruskal–Wallis test. Although a monotone increasing pattern with a number of alleles is expected for PK-related genes, the Kruskal–Wallis test does not consider a monotonic response pattern. For detecting such patterns in clinical and toxicological trials, a maximum contrast method has been proposed. We show how that method can be used with pharmacogenomics data to a develop test of association. Further, using simulation studies, we compare the power of the modified maximum contrast method to those of the maximum contrast method and the Kruskal–Wallis test. On the basis of the results of those studies, we suggest rules of thumb for which statistics to use in a given situation. An application of all three methods to an actual genome-wide pharmacogenomics study illustrates the practical relevance of our discussion.The Pharmacogenomics Journal (2009) 9, 137–146; doi:10.1038/tpj.2008.17; published online 23 December 2008 [ABSTRACT FROM AUTHOR]

We recently reported a joint analysis of genome-wide association (GWA) data on 958 sporadic amyotrophic lateral sclerosis (ALS) cases and 932 controls from Ireland and the publicly available data sets from the United States and the Netherlands. The strongest pooled association was rs10260404 in the dipeptidyl-peptidase 6 (DPP6) gene. Here, we sought confirmation of joint analysis signals in both an expanded Irish and a Polish ALS cohort. Among 287 522 autosomal single-nucleotide polymorphisms (SNPs), 27 were commonly associated on joint analysis of the Irish, US and Dutch GWAs. These 27 SNPs were genotyped in an expanded Irish cohort (312 patients with SALS; 259 controls) and an additional Polish cohort (218 patients; 356 controls). Eleven SNPs, including rs10260404, reached a final P-value below 0.05 in the Irish cohort. In the Polish cohort, only one SNP, rs6299711, showed nominal association with ALS. Pooling of data for 1267 patients with ALS and 1336 control subjects did not identify any association reaching Bonferroni significance (P<1.74 × 10−7). The present strategy did not reveal any consistently associated SNP across four populations. The result for DPP6 is surprising, as it has been replicated elsewhere. We discuss the possible interpretations and implications of these findings for future ALS GWA studies both within and between populations.European Journal of Human Genetics (2009) 17, 213–218; doi:10.1038/ejhg.2008.194; published online 5 November 2008 [ABSTRACT FROM AUTHOR]

Background: Service sire has a considerable impact on reproductive success in dairy cattle. Most gene mapping studies for bull fertility have focused on additive effects, while non-additive effects have been largely ignored. The main goal of this study was to assess the relevance of non-additive effects on Sire Conception Rate (SCR) in Holstein dairy cattle. The analysis included 7.5 k Holstein bulls with both SCR records and 57.8 k single nucleotide polymorphism (SNP) markers spanning the entire genome. Results: The importance of non-additive effects was evaluated using an efficient two-step mixed model-based approach. Four genomic regions located on chromosomes BTA8, BTA9, BTA13 and BTA17 showed marked dominance and/or recessive effects. Most of these regions harbor genes, such as ADAM28, DNAJA1, TBC1D20, SPO11, PIWIL3 and TMEM119, that are directly implicated in testis development, male germ line maintenance, and sperm maturation. Conclusions: This study provides further evidence for the relevance of non-additive effects in fitness-related traits, such as male fertility. In addition, these findings may point out new strategies for improving service sire fertility in dairy cattle via marker-assisted selection. [ABSTRACT FROM AUTHOR]