Genomic profiling of pediatric hematologic malignancies and diagnosis of cancer predisposition syndromes: tumor-only versus paired tumor-normal sequencing.

Not available.

PMLB v1.0: an open-source dataset collection for benchmarking machine learning methods.

MOTIVATION: Novel machine learning and statistical modeling studies rely on standardized comparisons to existing methods using well-studied benchmark datasets. Few tools exist that provide rapid access to many of these datasets through a standardized, user-friendly interface that integrates well with popular data science workflows. RESULTS: This release of PMLB (Penn Machine Learning Benchmarks) provides the largest collection of diverse, public benchmark datasets for evaluating new machine learning and data science methods aggregated in one location. v1.0 introduces a number of critical improvements developed following discussions with the open-source community. AVAILABILITY AND IMPLEMENTATION: PMLB is available at https://github.com/EpistasisLab/pmlb. Python and R interfaces for PMLB can be installed through the Python Package Index and Comprehensive R Archive Network, respectively.

Production and characterization of monoclonal antibodies against mandarinfish ranavirus and first identification of pyloric caecum as the major target tissue.

Mandarinfish ranavirus (MRV), also known as a variant of largemouth bass virus (LMBV), is an emerging pathogen in mandarinfish aquaculture. In this study, monoclonal antibodies (mAbs) against MRV were produced and characterized, and 7 mAbs were obtained through Western blotting screening and all 7 mAbs specifically recognized MRV/LMBV but not several piscine iridoviruses as ISKNV, GIV and TFV. By LC MS/MS analysis, the recognized viral proteins by seven mAbs were identified as MRV-pORF47L, MRV-pORF55R, MRV-pORF57L, MRV-pORF77L and MRV-pORF78L, respectively, and all five viral proteins are late expression structural proteins by Western blotting. Based on mAb 1C4, immuno-histochemistry and immuno-histo-fluorescence were performed to re-assess the tissue tropism of MRV. The result showed that abundant reactive signals were observed in infected spleen, kidney as well as intestine and pyloric caecum. Real-time quantitative PCR also demonstrated that spleen as well as pyloric caecum and intestines are the major target tissue upon MRV infection. In infected intestines and pyloric caecum, numerous enlarged, multinucleated cells with intracytoplasmic inclusions were identified as the target cells of MRV, suggesting that MRV serves as a digestive tract pathogen to mandarinfish, which may explain why acute infection of MRV can cause the typical clinicopathology featured by severe ascites.

Evaluating recommender systems for AI-driven biomedical informatics.

MOTIVATION: Many researchers with domain expertise are unable to easily apply machine learning (ML) to their bioinformatics data due to a lack of ML and/or coding expertise. Methods that have been proposed thus far to automate ML mostly require programming experience as well as expert knowledge to tune and apply the algorithms correctly. Here, we study a method of automating biomedical data science using a web-based AI platform to recommend model choices and conduct experiments. We have two goals in mind: first, to make it easy to construct sophisticated models of biomedical processes; and second, to provide a fully automated AI agent that can choose and conduct promising experiments for the user, based on the user's experiments as well as prior knowledge. To validate this framework, we conduct an experiment on 165 classification problems, comparing to state-of-the-art, automated approaches. Finally, we use this tool to develop predictive models of septic shock in critical care patients. RESULTS: We find that matrix factorization-based recommendation systems outperform metalearning methods for automating ML. This result mirrors the results of earlier recommender systems research in other domains. The proposed AI is competitive with state-of-the-art automated ML methods in terms of choosing optimal algorithm configurations for datasets. In our application to prediction of septic shock, the AI-driven analysis produces a competent ML model (AUROC 0.85±0.02) that performs on par with state-of-the-art deep learning results for this task, with much less computational effort. AVAILABILITY AND IMPLEMENTATION: PennAI is available free of charge and open-source. It is distributed under the GNU public license (GPL) version 3. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Immunological ignorance is an enabling feature of the oligo-clonal T cell response to melanoma neoantigens.

The impact of intratumoral heterogeneity (ITH) and the resultant neoantigen landscape on T cell immunity are poorly understood. ITH is a widely recognized feature of solid tumors and poses distinct challenges related to the development of effective therapeutic strategies, including cancer neoantigen vaccines. Here, we performed deep targeted DNA sequencing of multiple metastases from melanoma patients and observed ubiquitous sharing of clonal and subclonal single nucleotide variants (SNVs) encoding putative HLA class I-restricted neoantigen epitopes. However, spontaneous antitumor CD8+ T cell immunity in peripheral blood and tumors was restricted to a few clonal neoantigens featuring an oligo-/monoclonal T cell-receptor (TCR) repertoire. Moreover, in various tumors of the 4 patients examined, no neoantigen-specific TCR clonotypes were identified despite clonal neoantigen expression. Mature dendritic cell (mDC) vaccination with tumor-encoded amino acid-substituted (AAS) peptides revealed diverse neoantigen-specific CD8+ T responses, each composed of multiple TCR clonotypes. Isolation of T cell clones by limiting dilution from tumor-infiltrating lymphocytes (TILs) permitted functional validation regarding neoantigen specificity. Gene transfer of TCRαß heterodimers specific for clonal neoantigens confirmed correct TCR clonotype assignments based on high-throughput TCRBV CDR3 sequencing. Our findings implicate immunological ignorance of clonal neoantigens as the basis for ineffective T cell immunity to melanoma and support the concept that therapeutic vaccination, as an adjunct to checkpoint inhibitor treatment, is required to increase the breadth and diversity of neoantigen-specific CD8+ T cells.

Scaling tree-based automated machine learning to biomedical big data with a feature set selector.

MOTIVATION: Automated machine learning (AutoML) systems are helpful data science assistants designed to scan data for novel features, select appropriate supervised learning models and optimize their parameters. For this purpose, Tree-based Pipeline Optimization Tool (TPOT) was developed using strongly typed genetic programing (GP) to recommend an optimized analysis pipeline for the data scientist's prediction problem. However, like other AutoML systems, TPOT may reach computational resource limits when working on big data such as whole-genome expression data. RESULTS: We introduce two new features implemented in TPOT that helps increase the system's scalability: Feature Set Selector (FSS) and Template. FSS provides the option to specify subsets of the features as separate datasets, assuming the signals come from one or more of these specific data subsets. FSS increases TPOT's efficiency in application on big data by slicing the entire dataset into smaller sets of features and allowing GP to select the best subset in the final pipeline. Template enforces type constraints with strongly typed GP and enables the incorporation of FSS at the beginning of each pipeline. Consequently, FSS and Template help reduce TPOT computation time and may provide more interpretable results. Our simulations show TPOT-FSS significantly outperforms a tuned XGBoost model and standard TPOT implementation. We apply TPOT-FSS to real RNA-Seq data from a study of major depressive disorder. Independent of the previous study that identified significant association with depression severity of two modules, TPOT-FSS corroborates that one of the modules is largely predictive of the clinical diagnosis of each individual. AVAILABILITY AND IMPLEMENTATION: Detailed simulation and analysis code needed to reproduce the results in this study is available at https://github.com/lelaboratoire/tpot-fss. Implementation of the new TPOT operators is available at https://github.com/EpistasisLab/tpot. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Embedding covariate adjustments in tree-based automated machine learning for biomedical big data analyses.

BACKGROUND: A typical task in bioinformatics consists of identifying which features are associated with a target outcome of interest and building a predictive model. Automated machine learning (AutoML) systems such as the Tree-based Pipeline Optimization Tool (TPOT) constitute an appealing approach to this end. However, in biomedical data, there are often baseline characteristics of the subjects in a study or batch effects that need to be adjusted for in order to better isolate the effects of the features of interest on the target. Thus, the ability to perform covariate adjustments becomes particularly important for applications of AutoML to biomedical big data analysis. RESULTS: We developed an approach to adjust for covariates affecting features and/or target in TPOT. Our approach is based on regressing out the covariates in a manner that avoids 'leakage' during the cross-validation training procedure. We describe applications of this approach to toxicogenomics and schizophrenia gene expression data sets. The TPOT extensions discussed in this work are available at https://github.com/EpistasisLab/tpot/tree/v0.11.1-resAdj . CONCLUSIONS: In this work, we address an important need in the context of AutoML, which is particularly crucial for applications to bioinformatics and medical informatics, namely covariate adjustments. To this end we present a substantial extension of TPOT, a genetic programming based AutoML approach. We show the utility of this extension by applications to large toxicogenomics and differential gene expression data. The method is generally applicable in many other scenarios from the biomedical field.

Comparison of protective efficacy between two DNA vaccines encoding DnaK and GroEL against fish nocardiosis.

Fish nocardiosis is a chronic granulomatous bacterial disease mainly caused by three pathogenic bacteria, including Nocardia seriolae, N. asteroids and N. salmonicida. Molecular chaperone DnaK and GroEL were identified to be the common antigens of the three pathogenic Nocardia species in our previous studies. To evaluate the immune protective effect of two DNA vaccines encoding DnaK or GroEL against fish nocardiosis, hybrid snakehead were vaccinated and the immune responses induced by these two vaccines were comparatively analyzed. The results suggested it needed at least 7 d to transport DnaK or GroEL gene from injected muscle to head kidney, spleen and liver and stimulate host's immune system for later protection after immunization by DNA vaccines. Additionally, non-specific immunity parameters (serum lysozyme (LYZ), peroxidase (POD), acid phosphatase (ACP), alkaline phosphatase (AKP) and superoxide dismutase (SOD) activities), specific antibody (IgM) production and immune-related genes (MHCIα, MHCIIα, CD4, CD8α, IL-1ß and TNFα) were used to evaluate the immune responses induced in vaccinated hybrid snakehead. It proved that all the above-mentioned immune activities were significantly enhanced after immunization with these two DNA vaccines. The protective efficacy of pcDNA-DnaK and pcDNA-GroEL DNA vaccines, in terms of relative percentage survival (RPS), were 53.01% and 80.71% respectively. It demonstrated that these two DNA vaccines could increase the survival rate of hybrid snakehead against fish nocardiosis, albeit with variations in immunoprotective effects. Taken together, these results indicated that both pcDNA-DnaK and pcDNA-GroEL DNA vaccines could boost the innate, humoral and cellular immune response in hybrid snakehead and show highly protective efficacy against fish nocardiosis, suggesting that DnaK and GroEL were promising vaccine candidates. These findings will promote the development of DNA vaccines against fish nocardiosis in aquaculture.

Detection of genomic signatures of recent selection in commercial broiler chickens.

BACKGROUND: Identification of the genomic signatures of recent selection may help uncover causal polymorphisms controlling traits relevant to recent decades of selective breeding in livestock. In this study, we aimed at detecting signatures of recent selection in commercial broiler chickens using genotype information from single nucleotide polymorphisms (SNPs). A total of 565 chickens from five commercial purebred lines, including three broiler sire (male) lines and two broiler dam (female) lines, were genotyped using the 60K SNP Illumina iSelect chicken array. To detect genomic signatures of recent selection, we applied two methods based on population comparison, cross-population extended haplotype homozygosity (XP-EHH) and cross-population composite likelihood ratio (XP-CLR), and further analyzed the results to find genomic regions under recent selection in multiple purebred lines. RESULTS: A total of 321 candidate selection regions spanning approximately 1.45 % of the chicken genome in each line were detected by consensus of results of both XP-EHH and XP-CLR methods. To minimize false discovery due to genetic drift, only 42 of the candidate selection regions that were shared by 2 or more purebred lines were considered as high-confidence selection regions in the study. Of these 42 regions, 20 were 50 kb or less while 4 regions were larger than 0.5 Mb. In total, 91 genes could be found in the 42 regions, among which 19 regions contained only 1 or 2 genes, and 9 regions were located at gene deserts. CONCLUSIONS: Our results provide a genome-wide scan of recent selection signatures in five purebred lines of commercial broiler chickens. We found several candidate genes for recent selection in multiple lines, such as SOX6 (Sex Determining Region Y-Box 6) and cTR (Thyroid hormone receptor beta). These genes may have been under recent selection due to their essential roles in growth, development and reproduction in chickens. Furthermore, our results suggest that in some candidate regions, the same or opposite alleles have been under recent selection in multiple lines. Most of the candidate genes in the selection regions are novel, and as such they should be of great interest for future research into the genetic architecture of traits relevant to modern broiler breeding.

Characterization of a novel chicken muscle disorder through differential gene expression and pathway analysis using RNA-sequencing.

BACKGROUND: Improvements in poultry production within the past 50 years have led to increased muscle yield and growth rate, which may be contributing to an increased rate and development of new muscle disorders in chickens. Previously reported muscle disorders and conditions are generally associated with poor meat quality traits and have a significant negative economic impact on the poultry industry. Recently, a novel myopathy phenotype has emerged which is characterized by palpably "hard" or tough breast muscle. The objective of this study is to identify the underlying biological mechanisms that contribute to this emerging muscle disorder colloquially referred to as "Wooden Breast", through the use of RNA-sequencing technology. METHODS: We constructed cDNA libraries from five affected and six unaffected breast muscle samples from a line of commercial broiler chickens. After paired-end sequencing of samples using the Illumina Hiseq platform, we used Tophat to align the resulting sequence reads to the chicken reference genome and then used Cufflinks to find significant changes in gene transcript expression between each group. By comparing our gene list to previously published histology findings on this disorder and using Ingenuity Pathways Analysis (IPA®), we aim to develop a characteristic gene expression profile for this novel disorder through analyzing genes, gene families, and predicted biological pathways. RESULTS: Over 1500 genes were differentially expressed between affected and unaffected birds. There was an average of approximately 98 million reads per sample, across all samples. Results from the IPA analysis suggested "Diseases and Disorders" such as connective tissue disorders, "Molecular and Cellular Functions" such as cellular assembly and organization, cellular function and maintenance, and cellular movement, "Physiological System Development and Function" such as tissue development, and embryonic development, and "Top Canonical Pathways" such as, coagulation system, axonal guidance signaling, and acute phase response signaling, are associated with the Wooden Breast disease. CONCLUSIONS: There is convincing evidence by RNA-seq analysis to support localized hypoxia, oxidative stress, increased intracellular calcium, as well as the possible presence of muscle fiber-type switching, as key features of Wooden Breast Disease, which are supported by reported microscopic lesions of the disease.

Single-nucleotide polymorphisms in CD8A and their associations with T lymphocyte subpopulations in pig.

Findings from previous studies suggested that the cluster of the differentiation 8 alpha (CD8A) gene plays a prominent role in human T lymphocyte subpopulations. However, the evidence from pig population is still rare. To determine whether the important role of the CD8A gene is conserved in pig, a candidate gene analysis was performed herein through genotype-phenotype associations. Five single-nucleotide polymorphisms (SNPs) locating in the regulatory region of porcine CD8A gene were detected and tested for association analysis with seven T lymphocyte subpopulations (proportion of CD4(-)CD8(-), CD4(+)CD8(+), CD4(+)CD8(-), CD4(-)CD8(+), CD4(+), CD8(+), and the ratio of CD4(+) to CD8(+) T cells in peripheral blood) in 382 Large White piglets. After Bonferroni correction for multiple testing, four SNPs were significantly associated with some or all of the seven T lymphocyte subpopulations. Analyses of pairwise D' measures of linkage disequilibrium between all SNPs were also explored. Two haplotype blocks was inferred and the association study on haplotype level revealed similar effects on T lymphocyte subpopulations. In addition, the tissue-specific RNA expression pattern and electrophoretic mobility shift assay offered further explanation of the link between the CD8A gene with porcine T lymphocyte subpopulations. The findings presented here provide strong evidence for associations of CD8A variants with T lymphocyte subpopulations and may be applied in porcine breeding programs.

Linkage disequilibrium in crossbred and pure line chickens.

BACKGROUND: Both genome-wide association (GWA) studies and genomic selection depend on the level of non-random association of alleles at different loci, i.e. linkage disequilibrium (LD), across the genome. Therefore, characterizing LD is of fundamental importance to implement both approaches. In this study, using a 60K single nucleotide polymorphism (SNP) panel, we estimated LD and haplotype structure in crossbred broiler chickens and their component pure lines (one male and two female lines) and calculated the consistency of LD between these populations. RESULTS: The average level of LD (measured by r(2)) between adjacent SNPs across the chicken autosomes studied here ranged from 0.34 to 0.40 in the pure lines but was only 0.24 in the crossbred populations, with 28.4% of adjacent SNP pairs having an r(2) higher than 0.3. Compared with the pure lines, the crossbred populations consistently showed a lower level of LD, smaller haploblock sizes and lower haplotype homozygosity on macro-, intermediate and micro-chromosomes. Furthermore, correlations of LD between markers at short distances (0 to 10 kb) were high between crossbred and pure lines (0.83 to 0.94). CONCLUSIONS: Our results suggest that using crossbred populations instead of pure lines can be advantageous for high-resolution QTL (quantitative trait loci) mapping in GWA studies and to achieve good persistence of accuracy of genomic breeding values over generations in genomic selection. These results also provide useful information for the design and implementation of GWA studies and genomic selection using crossbred populations.

Genome-wide association study for pigmentation traits in Chinese Holstein population.

With the Illumina BovineSNP50K BeadChip, we performed a genome-wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case-control design was used. Cases were the cows with PB <0.30 (n = 129) and TC <2 points (n = 140); controls were those with PB >0.90 (n = 58) and TC >4 points (n = 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome-wide significant (P < 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.

Association of the porcine cluster of differentiation 4 gene with T lymphocyte subpopulations and its expression in immune tissues.

Cluster of differentiation 4 (CD4) is mainly expressed on CD4(+) T cells, which plays an important role in immune response. The aim of this study was to detect the association between polymorphisms of the CD4 gene and T lymphocyte subpopulations in pigs, and to investigate the effects of genetic variation on the CD4 gene expression level in immune tissues. Five missense mutations in the CD4 gene were identified using DNA pooling sequencing assays, and two main haplotypes (CCTCC and AGCTG) in strong linkage disequilibrium (with frequencies of 50.26% and 46.34%, respectively) were detected in the population of Large White pigs. Our results indicated that the five SNPs and the two haplotypes were significantly associated with the proportions of CD4(-)CD8(-), CD4(+)CD8(+), CD4(+)CD8(-), CD4(+) and CD4(+)/CD8(+) in peripheral blood (p<0.05). Gene expression analysis showed the mRNA level of the CD4 gene in thymus was significantly higher than that in lymph node and spleen (p<0.05). However, no significant difference was observed between animals with CCTCC/CCTCC genotype and animals with AGCTG/AGCTG genotype in the three immune tissues (p>0.05). These results indicate that the CD4 gene may influence T lymphocyte subpopulations and can be considered as a candidate gene affecting immunity in pigs.

The Inactivated ISKNV-I Vaccine Confers Highly Effective Cross-Protection against Epidemic RSIV-I and RSIV-II from Cultured Spotted Sea Bass Lateolabrax maculatus.

The genus Megalocytivirus of the family Iridoviridae is composed of two distinct species, namely, infectious spleen and kidney necrosis virus (ISKNV) and scale drop disease virus (SDDV), and both are important causative agents in a variety of bony fish worldwide. Of them, the ISKNV species is subdivided into three genotypes, namely, red seabream iridovirus (RSIV), ISKNV, and turbot reddish body iridovirus (TRBIV), and a further six subgenotypes, RSIV-I, RSIV-II, ISKNV-I, ISKNV-II, TRBIV-I, and TRBIV-II. Commercial vaccines derived from RSIV-I , RSIV-II and ISKNV-I have been available to several fish species. However, studies regarding the cross-protection effect among different genotype or subgenotype isolates have not been fully elucidated. In this study, RSIV-I and RSIV-II were demonstrated as the causative agents in cultured spotted seabass, Lateolabrax maculatus, through serial robust evidence, including cell culture-based viral isolation, whole-genome determination and phylogeny analysis, artificial challenge, histopathology, immunohistochemistry, and immunofluorescence as well as transmission electron microscope observation. Thereafter, a formalin-killed cell (FKC) vaccine generated from an ISKNV-I isolate was prepared to evaluate the protective effects against two spotted seabass original RSIV-I and RSIV-II. The result showed that the ISKNV-I-based FKC vaccine conferred almost complete cross-protection against RSIV-I and RSIV-II as well as ISKNV-I itself. No serotype difference was observed among RSIV-I, RSIV-II, and ISKNV-I. Additionally, the mandarin fish Siniperca chuatsi is proposed as an ideal infection and vaccination fish species for the study of various megalocytiviral isolates. IMPORTANCE Red seabream iridovirus (RSIV) infects a wide mariculture bony fish and has resulted in significant annual economic loss worldwide. Previous studies showed that the phenotypic diversity of infectious RSIV isolates would lead to different virulence characteristics, viral antigenicity, and vaccine efficacy as well as host range. Importantly, it is still doubted whether a universal vaccine could confer the same highly protective effect against various genotypic isolates. Our study here presented enough experimental evidence that a water in oil (w/o) formation of inactivated ISKNV-I vaccine could confer almost complete protection against RSIV-I and RSIV-II as well as ISKNV-I itself. Our study provides valuable data for better understanding the differential infection and immunity among different genotypes of ISKNV and RSIV isolates in the genus Megalocytivirus.

A genome-wide detection of copy number variations using SNP genotyping arrays in swine.

BACKGROUND: Copy Number Variations (CNVs) have been shown important in both normal phenotypic variability and disease susceptibility, and are increasingly accepted as another important source of genetic variation complementary to single nucleotide polymorphism (SNP). Comprehensive identification and cataloging of pig CNVs would be of benefit to the functional analyses of genome variation. RESULTS: In this study, we performed a genome-wide CNV detection based on the Porcine SNP60 genotyping data of 474 pigs from three pure breed populations (Yorkshire, Landrace and Songliao Black) and one Duroc × Erhualian crossbred population. A total of 382 CNV regions (CNVRs) across genome were identified, which cover 95.76Mb of the pig genome and correspond to 4.23% of the autosomal genome sequence. The length of these CNVRs ranged from 5.03 to 2,702.7kb with an average of 250.7kb, and the frequencies of them varied from 0.42 to 20.87%. These CNVRs contains 1468 annotated genes, which possess a great variety of molecular functions, making them a promising resource for exploring the genetic basis of phenotypic variation within and among breeds. To confirmation of these findings, 18 CNVRs representing different predicted status and frequencies were chosen for validation via quantitative real time PCR (qPCR). Accordingly, 12 (66.67%) of them was successfully confirmed. CONCLUSIONS: Our results demonstrated that currently available Porcine SNP60 BeadChip can be used to capture CNVs efficiently. Our study firstly provides a comprehensive map of copy number variation in the pig genome, which would be of help for understanding the pig genome and provide preliminary foundation for investigating the association between various phenotypes and CNVs.

Genome-wide association study for T lymphocyte subpopulations in swine.

BACKGROUND: Lymphocytes act as a major component of the adaptive immune system, taking very crucial responsibility for immunity. Differences in proportions of T-cell subpopulations in peripheral blood among individuals under same conditions provide evidence of genetic control on these traits, but little is known about the genetic mechanism of them, especially in swine. Identification of the genetic control on these variants may help the genetic improvement of immune capacity through selection. RESULTS: To identify genomic regions responsible for these immune traits in swine, a genome-wide association study was conducted. A total of 675 pigs of three breeds were involved in the study. At 21 days of age, all individuals were vaccinated with modified live classical swine fever vaccine. Blood samples were collected when the piglets were 20 and 35 days of age, respectively. Seven traits, including the proportions of CD4+, CD8+, CD4+CD8+, CD4+CD8-, CD4-CD8+, CD4-CD8- and the ratio of CD4+ to CD8+ T cells were measured at the two ages. All the samples were genotyped for 62,163 single nucleotide polymorphisms (SNP) using the Illumina porcineSNP60k BeadChip. 40833 SNPs were selected after quality control for association tests between SNPs and each immune trait considered based on a single-locus regression model. To tackle the issue of multiple testing in GWAS, 10,000 permutations were performed to determine the chromosome-wise and genome-wise significance levels of association tests. In total, 61 SNPs with chromosome-wise significance level and 3 SNPs with genome-wise significance level were identified. 27 significant SNPs were located within the immune-related QTL regions reported in previous studies. Furthermore, several significant SNPs fell into the regions harboring known immunity-related genes, 14 of them fell into the regions which harbor some known T cell-related genes. CONCLUSIONS: Our study demonstrated that genome-wide association studies would be a feasible way for revealing the potential genetics variants affecting T-cell subpopulations. Results herein lay a preliminary foundation for further identifying the causal mutations underlying swine immune capacity in follow-up studies.

Improved LASSO priors for shrinkage quantitative trait loci mapping.

Recently, the Bayesian least absolute shrinkage and selection operator (LASSO) has been successfully applied to multiple quantitative trait loci (QTL) mapping, which assigns the double-exponential prior and the Student's t prior to QTL effect that lead to the shrinkage estimate of QTL effect. However, as reported by many researchers, the Bayesian LASSO usually cannot effectively shrink the effects of zero-effect QTL very close to zero. In this study, the double-exponential prior and Student's t prior are modified so that the estimate of the effect for zero-effect QTL can be effectively shrunk toward zero. It is also found that the Student's t prior is virtually the same as the Jeffreys' prior, since both the shape and scale parameters of the scaled inverse Chi-square prior involved in the Student's t prior are estimated very close to zero. Besides the two modified Bayesian Markov chain Monte Carlo (MCMC) algorithms, an expectation-maximization (EM) algorithm with use of the modified double-exponential prior is also adapted. The results shows that the three new methods perform similarly on true positive rate and false positive rate for QTL detection, and all of them outperform the Bayesian LASSO.

Tissues Expression, Polymorphisms of IFN Regulatory Factor 6 (IRF6) Gene and Their Associated with Immune Traits in Three Pig Populations.

Interferon regulatory factor 6 (IRF6) gene is a member of the IRF-family, and plays functionally diverse roles in the regulation of the immune system. In this report, the 13,720 bp porcine IRF6 genomic DNA structure was firstly identified with a putative IRF6 protein of 467 amino acids. Alignment and phylogenetic analysis of the porcine IRF6 amino acid sequences with their homologies to other species showed high identity (over 96%). Tissues expression of IRF6 mRNA was observed by RT-PCR, the results revealed IRF6 expressed widely in eight tissues. One SNP (HQ026023:1383 G>C) in exon7 and two SNPs (HQ026023:130 G>A; 232 C>T) in the 5 ' promoter region of porcine IRF6 gene were demonstrated b y DNA sequencing analysis. A further analysis of SNP genotypes associated with immune traits including IFN-γ and IL10 concentrations in serum was carried out in three pig populations including Large White, Landraces and Songliao Black pig (a Chinese indigenous breed). The results showed that the SNP (HQ026023:1383 G>C) was significantly associated with the level of IFN-γ (d 20) in serum (p = 0.038) and the ratio of IFN-γ to IL10 (d 20) in serum (p = 0.041); The other two SNPs (HQ026023:130 G>A; 232 C>T) were highly significantly associated with IL10 level in serum both at the day 20 (p = 0.005; p = 0.001) and the day 35 (p = 0.004; p = 0.006). Identification of the porcine IRF6 gene will help our further understanding of the molecular basis of the IFN regulation pathway in the porcine immune response. All these results should indicate that the IRF6 gene can be regarded as a molecular marker associated with the IL10 level in serum and used for genetic selection in the pig breeding.

Genetic Analysis of Coronary Artery Disease Using Tree-Based Automated Machine Learning Informed By Biology-Based Feature Selection.

Machine Learning (ML) approaches are increasingly being used in biomedical applications. Important challenges of ML include choosing the right algorithm and tuning the parameters for optimal performance. Automated ML (AutoML) methods, such as Tree-based Pipeline Optimization Tool (TPOT), have been developed to take some of the guesswork out of ML thus making this technology available to users from more diverse backgrounds. The goals of this study were to assess applicability of TPOT to genomics and to identify combinations of single nucleotide polymorphisms (SNPs) associated with coronary artery disease (CAD), with a focus on genes with high likelihood of being good CAD drug targets. We leveraged public functional genomic resources to group SNPs into biologically meaningful sets to be selected by TPOT. We applied this strategy to data from the U.K. Biobank, detecting a strikingly recurrent signal stemming from a group of 28 SNPs. Importance analysis of these SNPs uncovered functional relevance of the top SNPs to genes whose association with CAD is supported in the literature and other resources. Furthermore, we employed game-theory based metrics to study SNP contributions to individual-level TPOT predictions and discover distinct clusters of well-predicted CAD cases. The latter indicates a promising approach towards precision medicine.