Background: Receptor-like kinases (RLKs) play key roles during development and in responses to the environment. Despite the relevance of the RLK family and the completion of the tomato genome sequencing, the tomato RLK family has not yet been characterized, and a framework for functional predictions of the members of the family is lacking. Results: To generate a complete list of all the members of the tomato RLK family, we performed a phylogenetic analysis using the Arabidopsis family as a template. A total of 647 RLKs were identified in the tomato genome, which were organized into the same subfamily clades as Arabidopsis RLKs. Only eight of 58 RLK subfamilies exhibited specific expansion/reduction compared to their Arabidopsis counterparts. We also characterized the LRRII-RLK family by phylogeny, genomic analysis, expression profile and interaction with the virulence factor from begomoviruses, the nuclear shuttle protein (NSP). The LRRII subfamily members from tomato and Arabidopsis were highly conserved in both sequence and structure. Nevertheless, the majority of the orthologous pairs did not display similar conservation in the gene expression profile, indicating that these orthologs may have diverged in function after speciation. Based on the fact that members of the Arabidopsis LRRII subfamily (AtNIK1, AtNIK2 and AtNIK3) interact with the begomovirus nuclear shuttle protein (NSP), we examined whether the tomato orthologs of NIK, BAK1 and NsAK genes interact with NSP of Tomato Yellow Spot Virus (ToYSV). The tomato orthologs of NSP interactors, SlNIKs and SlNsAK, interacted specifically with NSP in yeast and displayed an expression pattern consistent with the pattern of geminivirus infection. In addition to suggesting a functional analogy between these phylogenetically classified orthologs, these results expand our previous observation that NSP-NIK interactions are neither virus-specific nor host-specific. Conclusions: The tomato RLK superfamily is made-up of 647 proteins that form a monophyletic tree with the Arabidopsis RLKs and is divided into 58 subfamilies. Few subfamilies have undergone expansion/reduction, and only six proteins were lineage-specific. Therefore, the tomato RLK family shares functional and structural conservation with Arabidopsis. For the LRRII-RLK members SlNIK1 and SlNIK3, we observed functions analogous to those of their Arabidopsis counterparts with respect to protein-protein interactions and similar expression profiles, which predominated in tissues that support high efficiency of begomovirus infection. Therefore, NIK-mediated antiviral signaling is also likely to operate in tomato, suggesting that tomato NIKs may be good targets for engineering resistance against tomato-infecting begomoviruses.

Reproductive traits such as number of stillborn piglets (SB) and number of teats (NT) have been evaluated in many genome-wide association studies (GWAS). Most of these GWAS were performed under the assumption that these traits were normally distributed. However, both SB and NT are discrete (e.g. count) variables. Therefore, it is necessary to test for better fit of other appropriate statistical models based on discrete distributions. In addition, although many GWAS have been performed, the biological meaning of the identified candidate genes, as well as their functional relationships still need to be better understood. Here, we performed and tested a Bayesian treatment of a GWAS model assuming a Poisson distribution for SB and NT in a commercial pig line. To explore the biological role of the genes that underlie SB and NT and identify the most likely candidate genes, we used the most significant single nucleotide polymorphisms (SNPs), to collect related genes and generated gene-transcription factor (TF) networks. Comparisons of the Poisson and Gaussian distributions showed that the Poisson model was appropriate for SB, while the Gaussian was appropriate for NT. The fitted GWAS models indicated 18 and 65 significant SNPs with one and nine quantitative trait locus (QTL) regions within which 18 and 57 related genes were identified for SB and NT, respectively. Based on the related TF, we selected the most representative TF for each trait and constructed a gene-TF network of gene-gene interactions and identified new candidate genes. Our comparative analyses showed that the Poisson model presented the best fit for SB. Thus, to increase the accuracy of GWAS, counting models should be considered for this kind of trait. We identified multiple candidate genes (e.g. PTP4A2, NPHP1, and CYP24A1 for SB and YLPM1, SYNDIG1L, TGFB3, and VRTN for NT) and TF (e.g. NF-κB and KLF4 for SB and SOX9 and ELF5 for NT), which were consistent with known newborn survival traits (e.g. congenital heart disease in fetuses and kidney diseases and diabetes in the mother) and mammary gland biology (e.g. mammary gland development and body length).

In recent years, there has been increased interest in the study of the molecular processes that affect semen traits. In this study, our aim was to identify quantitative trait loci (QTL) regions associated with four semen traits (motility, progressive motility, number of sperm cells per ejaculate and total morphological defects) in two commercial pig lines (L1: Large White type and L2: Landrace type). Since the number of animals with both phenotypes and genotypes was relatively small in our dataset, we conducted a weighted single-step genome-wide association study, which also allows unequal variances for single nucleotide polymorphisms. In addition, our aim was also to identify candidate genes within QTL regions that explained the highest proportions of genetic variance. Subsequently, we performed gene network analyses to investigate the biological processes shared by genes that were identified for the same semen traits across lines. We identified QTL regions that explained up to 10.8% of the genetic variance of the semen traits on 12 chromosomes in L1 and 11 chromosomes in L2. Sixteen QTL regions in L1 and six QTL regions in L2 were associated with two or more traits within the population. Candidate genes SCN8A, PTGS2, PLA2G4A, DNAI2, IQCG and LOC102167830 were identified in L1 and NME5, AZIN2, SPATA7, METTL3 and HPGDS in L2. No regions overlapped between these two lines. However, the gene network analysis for progressive motility revealed two genes in L1 (PLA2G4A and PTGS2) and one gene in L2 (HPGDS) that were involved in two biological processes i.e. eicosanoid biosynthesis and arachidonic acid metabolism. PTGS2 and HPGDS were also involved in the cyclooxygenase pathway. We identified several QTL regions associated with semen traits in two pig lines, which confirms the assumption of a complex genetic determinism for these traits. A large part of the genetic variance of the semen traits under study was explained by different genes in the two evaluated lines. Nevertheless, the gene network analysis revealed candidate genes that are involved in shared biological pathways that occur in mammalian testes, in both lines.

A complete approach for genome-wide selection (GWS) involves reliable statistical genetics models and methods. Reports on this topic are common for additive genetic models but not for additive-dominance models. The objective of this paper was (i) to compare the performance of 10 additive-dominance predictive models (including current models and proposed modifications), fitted using Bayesian, Lasso and Ridge regression approaches; and (ii) to decompose genomic heritability and accuracy in terms of three quantitative genetic information sources, namely, linkage disequilibrium (LD), co-segregation (CS) and pedigree relationships or family structure (PR). The simulation study considered two broad sense heritability levels (0.30 and 0.50, associated with narrow sense heritabilities of 0.20 and 0.35, respectively) and two genetic architectures for traits (the first consisting of small gene effects and the second consisting of a mixed inheritance model with five major genes). G-REML/G-BLUP and a modified Bayesian/Lasso (called BayesA*B* or t-BLASSO) method performed best in the prediction of genomic breeding as well as the total genotypic values of individuals in all four scenarios (two heritabilities x two genetic architectures). The BayesA*B*-type method showed a better ability to recover the dominance variance/additive variance ratio. Decomposition of genomic heritability and accuracy revealed the following descending importance order of information: LD, CS and PR not captured by markers, the last two being very close. Amongst the 10 models/methods evaluated, the G-BLUP, BAYESA*B* (−2,8) and BAYESA*B* (4,6) methods presented the best results and were found to be adequate for accurately predicting genomic breeding and total genotypic values as well as for estimating additive and dominance in additive-dominance genomic models.

In the era of genome-wide selection (GWS), genotype-by-environment (G×E) interactions can be studied using genomic information, thus enabling the estimation of SNP marker effects and the prediction of genomic estimated breeding values (GEBV) for young candidates for selection in different environments. Although G×E studies in pigs are scarce, the use of artificial insemination has enabled the distribution of genetic material from sires across multiple environments. Given the relevance of reproductive traits, such as the total number born (TNB) and the variation in environmental conditions encountered by commercial dams, understanding G×E interactions can be essential for choosing the best sires for different environments. The present work proposes a two-step reaction norm approach for G×E analysis using genomic information. The first step provided estimates of environmental effects (herd-year-season, HYS), and the second step provided estimates of the intercept and slope for the TNB across different HYS levels, obtained from the first step, using a random regression model. In both steps, pedigree (A) and genomic (G) relationship matrices were considered. The genetic parameters (variance components, h² and genetic correlations) were very similar when estimated using the A and G relationship matrices. The reaction norm graphs showed considerable differences in environmental sensitivity between sires, indicating a reranking of sires in terms of genetic merit across the HYS levels. Based on the G matrix analysis, SNP by environment interactions were observed. For some SNP, the effects increased at increasing HYS levels, while for others, the effects decreased at increasing HYS levels or showed no changes between HYS levels. Cross-validation analysis demonstrated better performance of the genomic approach with respect to traditional pedigrees for both the G×E and standard models. The genomic reaction norm model resulted in an accuracy of GEBV for “juvenile” boars varying from 0.14 to 0.44 across different HYS levels, while the accuracy of the standard genomic prediction model, without reaction norms, varied from 0.09 to 0.28. These results show that it is important and feasible to consider G×E interactions in evaluations of sires using genomic prediction models and that genomic information can increase the accuracy of selection across environments. Copyright © 2014. American Society of Animal Science .