Accurate Phenotype Prediction Research Articles

Residual networks without pooling layers improve the accuracy of genomic predictions.

Residual neural network genomic selection is the first GS algorithm to reach 35 layers, and its prediction accuracy surpasses previous algorithms. With the decrease in DNA sequencing costs and the development of deep learning, phenotype prediction accuracy by genomic selection (GS) continues to improve. Residual networks, a widely validated deep learning technique, are introduced to deep learning for GS. Since each locus has a different weighted impact on the phenotype, strided convolutions are more suitable for GS problems than pooling layers. Through the above technological innovations, we propose a GS deep learning algorithm, residual neural network for genomic selection (ResGS). ResGS is the first neural network to reach 35 layers in GS. In 15 cases from four public data, the prediction accuracy of ResGS is higher than that of ridge-regression best linear unbiased prediction, support vector regression, random forest, gradient boosting regressor, and deep neural network genomic prediction in most cases. ResGS performs well in dealing with gene-environment interaction. Phenotypes from other environments are imported into ResGS along with genetic data. The prediction results are much better than just providing genetic data as input, which demonstrates the effectiveness of GS multi-modal learning. Standard deviation is recommended as an auxiliary GS evaluation metric, which could improve the distribution of predicted results. Deep learning for GS, such as ResGS, is becoming more accurate in phenotype prediction.

Multi-omics data integration and drug screening of AML cancer using Generative Adversarial Network

In the era of precision medicine, accurate disease phenotype prediction for heterogeneous diseases, such as cancer, is emerging due to advanced technologies that link genotypes and phenotypes. However, it is difficult to integrate different types of biological data because they are so varied. In this study, we focused on predicting the traits of a blood cancer called Acute Myeloid Leukemia (AML) by combining different kinds of biological data. We used a recently developed method called Omics Generative Adversarial Network (GAN) to better classify cancer outcomes. The primary advantages of a GAN include its ability to create synthetic data that is nearly indistinguishable from real data, its high flexibility, and its wide range of applications, including multi-omics data analysis. In addition, the GAN was effective at combining two types of biological data. We created synthetic datasets for gene activity and DNA methylation. Our method was more accurate in predicting disease traits than using the original data alone. The experimental results provided evidence that the creation of synthetic data through interacting multi-omics data analysis using GANs improves the overall prediction quality. Furthermore, we identified the top-ranked significant genes through statistical methods and pinpointed potential candidate drug agents through in-silico studies. The proposed drugs, also supported by other independent studies, might play a crucial role in the treatment of AML cancer. The code is available on GitHub; https://github.com/SabrinAfroz/omicsGAN_codes?fbclid=IwAR1-/stuffmlE0hyWgSu2wlXo6dYlKUei3faLdlvpxTOOUPVlmYCloXf4Uk9ejK4I

Accessory genes define species-specific routes to antibiotic resistance.

A deeper understanding of the relationship between the antimicrobial resistance (AMR) gene carriage and phenotype is necessary to develop effective response strategies against this global burden. AMR phenotype is often a result of multi-gene interactions; therefore, we need approaches that go beyond current simple AMR gene identification tools. Machine-learning (ML) methods may meet this challenge and allow the development of rapid computational approaches for AMR phenotype classification. To examine this, we applied multiple ML techniques to 16,950 bacterial genomes across 28 genera, with corresponding MICs for 23 antibiotics with the aim of training models to accurately determine the AMR phenotype from sequenced genomes. This resulted in a >1.5-fold increase in AMR phenotype prediction accuracy over AMR gene identification alone. Furthermore, we revealed 528 unique (often species-specific) genomic routes to antibiotic resistance, including genes not previously linked to the AMR phenotype. Our study demonstrates the utility of ML in predicting AMR phenotypes across diverse clinically relevant organisms and antibiotics. This research proposes a rapid computational method to support laboratory-based identification of the AMR phenotype in pathogens.

Multi-omics Deep-learning Prediction of Homologous Recombination Deficiency-like Phenotype Improved Risk Stratification and Guided Therapeutic Decisions in Gynecological Cancers.

Homologous recombination deficiency (HRD) is a well-recognized important biomarker in determining the clinical benefits of platinum-based chemotherapy and PARP inhibitor therapy for patients diagnosed with gynecologic cancers. Accurate prediction of HRD phenotype remains challenging. Here, we proposed a novel Multi-Omics integrative Deep-learning framework named MODeepHRD for detecting HRD-positive phenotype. MODeepHRD utilizes a convolutional attention autoencoder that effectively leverages omics-specific and cross-omics complementary knowledge learning. We trained MODeepHRD on 351 ovarian cancer (OV) patients using transcriptomic, DNA methylation and mutation data, and validated it in 2133 OV samples of 22 datasets. The predicted HRD-positive tumors were significantly associated with improved survival (HR = 0.68; 95% CI, 0.60-0.77; log-rank p < 0.001 for meta-cohort; HR = 0.5; 95% CI, 0.29-0.86; log-rank p = 0.01 for ICGC-OV cohort) and higher responsetoplatinum-basedchemotherapy compared to predicted HRD-negative tumors. The translational potential of MODeepHRDs was further validated in multicenter breast and endometrial cancer cohorts. Furthermore, MODeepHRD outperforms conventional machine-learning methods and other similar task approaches. In conclusion, our study demonstrates the promising value of deep learning as a solution for HRD testing in the clinical setting. MODeepHRD holds potential clinical applicability in guiding patient risk stratification and therapeutic decisions, providing valuable insights for precision oncology and personalized treatment strategies.

Morphological Signal Processing for Phenotype Recognition of Human Pluripotent Stem Cells Using Machine Learning Methods.

Human pluripotent stem cells have the potential for unlimited proliferation and controlled differentiation into various somatic cells, making them a unique tool for regenerative and personalized medicine. Determining the best clone selection is a challenging problem in this field and requires new sensing instruments and methods able to automatically assess the state of a growing colony ('phenotype') and make decisions about its destiny. One possible solution for such label-free, non-invasive assessment is to make phase-contrast images and/or videos of growing stem cell colonies, process the morphological parameters ('morphological portrait', or signal), link this information to the colony phenotype, and initiate an automated protocol for the colony selection. As a step in implementing this strategy, we used machine learning methods to find an effective model for classifying the human pluripotent stem cell colonies of three lines according to their morphological phenotype ('good' or 'bad'), using morphological parameters from the previously published data as predictors. We found that the model using cellular morphological parameters as predictors and artificial neural networks as the classification method produced the best average accuracy of phenotype prediction (67%). When morphological parameters of colonies were used as predictors, logistic regression was the most effective classification method (75% average accuracy). Combining the morphological parameters of cells and colonies resulted in the most effective model, with a 99% average accuracy of phenotype prediction. Random forest was the most efficient classification method for the combined data. We applied feature selection methods and showed that different morphological parameters were important for phenotype recognition via either cellular or colonial parameters. Our results indicate a necessity for retaining both cellular and colonial morphological information for predicting the phenotype and provide an optimal choice for the machine learning method. The classification models reported in this study could be used as a basis for developing and/or improving automated solutions to control the quality of human pluripotent stem cells for medical purposes.

Current Uses and Future Perspectives of Genomic Technologies in Clinical Microbiology.

Recent advancements in sequencing technology and data analytics have led to a transformative era in pathogen detection and typing. These developments not only expedite the process, but also render it more cost-effective. Genomic analyses of infectious diseases are swiftly becoming the standard for pathogen analysis and control. Additionally, national surveillance systems can derive substantial benefits from genomic data, as they offer profound insights into pathogen epidemiology and the emergence of antimicrobial-resistant strains. Antimicrobial resistance (AMR) is a pressing global public health issue. While clinical laboratories have traditionally relied on culture-based antimicrobial susceptibility testing, the integration of genomic data into AMR analysis holds immense promise. Genomic-based AMR data can furnish swift, consistent, and highly accurate predictions of resistance phenotypes for specific strains or populations, all while contributing invaluable insights for surveillance. Moreover, genome sequencing assumes a pivotal role in the investigation of hospital outbreaks. It aids in the identification of infection sources, unveils genetic connections among isolates, and informs strategies for infection control. The One Health initiative, with its focus on the intricate interconnectedness of humans, animals, and the environment, seeks to develop comprehensive approaches for disease surveillance, control, and prevention. When integrated with epidemiological data from surveillance systems, genomic data can forecast the expansion of bacterial populations and species transmissions. Consequently, this provides profound insights into the evolution and genetic relationships of AMR in pathogens, hosts, and the environment.

GSCtool: A Novel Descriptor that Characterizes the Genome for Applying Machine Learning in Genomics

Machine learning (ML) is one of the core driving forces for the next breeding stage, and Breeding 4.0. Genotype matrix based on single‐nucleotide polymorphisms (SNPs) is often used in ML for genome‐to‐phenotype prediction. Genotype matrix has an inherent defect, as the feature spaces it generates across different individuals or groups are inconsistent, and this hinders the application of ML. To overcome the challenge, a genome descriptor, Genic SNPs Composition Tool (GSCtool) is developed, which counts the number of SNPs in each gene of the genome so the dimension of the feature vectors equals the number of annotated genes in a species. Compared to using the genotype matrix, using GSCtool significantly decreases the model training time and has a higher accuracy of phenotype prediction. GSCtool also achieves good performance in variety identification, which is useful in crop variety protection. In general, GSCtool will help facilitate the application and study of genomic ML. The source code and test data of GSCtool are freely available at https://github.com/SZJhacker/GSCtool and https://gitee.com/shenzijie/GSCtool.

Identification and validation of obesity related genes signature based on microenvironment phenotypes in prostate adenocarcinoma.

The role of obesity related genes (ORGs) in the immune checkpoint inhibitors (ICIs) treatment of prostate adenocarcinoma (PRAD) has not yet been proved by research. We comprehensively evaluated the ORGs patterns in PRAD based on tumor microenvironment (TME) phenotypes and immunotherapy efficacies. Then we constructed a ORGs risk score for prognosis and a ORGs signature for accurate prediction of TME phenotype and immunotherapy efficacy in order to evaluate individual patients. Two distinct ORGs patterns were generated. The two ORGs patterns were consistent with inflammatory and non-inflammatory TME phenotypes. ORGs patterns had an important role for predicting immunotherapy efficacies. Next, we constructed a ORGs risk score for predicting each patient's prognosis with high performance in TCGA-PRAD. The ORGs risk score could be well verified in the external cohorts including GSE70769 and GSE21034. Then, we developed a ORGs signature and found it was significantly positively correlated with tumor-infiltrating lymphocytes in TCGA-PRAD. We found that each patient in the high-risk ORGs signature group represented a non-inflamed TME phenotype on the single cell level. The patients with high ORGs signature had more sensitivity to immunotherapy. And those ORGs were verified. ORGs pattern depicts different TME phenotypes in PRAD. The ORGs risk score and ORGs signature have an important role for predicting prognosis and immunotherapy efficacies.

Local flux coordination and global gene expression regulation in metabolic modeling

Genome-scale metabolic networks (GSMs) are fundamental systems biology representations of a cell’s entire set of stoichiometrically balanced reactions. However, such static GSMs do not incorporate the functional organization of metabolic genes and their dynamic regulation (e.g., operons and regulons). Specifically, there are numerous topologically coupled local reactions through which fluxes are coordinated; the global growth state often dynamically regulates many gene expression of metabolic reactions via global transcription factor regulators. Here, we develop a GSM reconstruction method, Decrem, by integrating locally coupled reactions and global transcriptional regulation of metabolism by cell state. Decrem produces predictions of flux and growth rates, which are highly correlated with those experimentally measured in both wild-type and mutants of three model microorganisms Escherichia coli, Saccharomyces cerevisiae, and Bacillus subtilis under various conditions. More importantly, Decrem can also explain the observed growth rates by capturing the experimentally measured flux changes between wild-types and mutants. Overall, by identifying and incorporating locally organized and regulated functional modules into GSMs, Decrem achieves accurate predictions of phenotypes and has broad applications in bioengineering, synthetic biology, and microbial pathology.

Noninvasive screening of fetal RHD genotype in Chinese pregnant women with serologic RhD-negative phenotype.

Noninvasive fetal RHD genotyping has been provided to nonimmunized RhD-negative pregnant women to guide anti-D prophylaxis. Among the Chinese, more than 30% of the RhD-negative phenotype is associated with variant RHD alleles, which would limit the accuracy of fetal RHD status prediction; thus, more targeting and proper programs need to be developed. Fluorescence quantitative polymerase chain reaction PCR (qPCR) or Sanger sequencing on all RHD exons was used to detect maternal RHD genotypes. For pregnant women with RHD*01N.01 or RHD*01N.03 alleles, the presence of RHD exons 5 and 10 in cell-free DNA was determined by qPCR. For pregnant women with the RHD(1227G>A) allele, high-throughput sequencing on exon 9 of the RHD gene and RHCE gene was used to predict fetal RhD phenotype. Among 65 cases of Chinese pregnant women with the serologic RhD-negative phenotype, three major genotypes were identified: RHD*01N.01/RHD*01N.01 (61.5%), RHD*01N.01/RHD(1227G>A) or RHD*01N.03/RHD(1227G>A) (20%), and RHD*01N.01/RHD*01N.03 (13.8%), along with three cases of minor genotypes (4.6%). For 43 pregnant women with the RHD*01N.01 or RHD*01N.03 alleles, qPCR on maternal cell-free DNA yielded a 98.5% (42/43) accuracy rate and 100% successful prediction rate. High-throughput sequencing was successfully used to predict fetal RhD phenotypes for 13 pregnant women with RHD(1227G>A). On the basis of maternal RHD genotyping, fetal genotyping through qPCR or high-throughput sequencing can improve the accuracy and success rate of prenatal fetal RhD phenotype prediction among Chinese pregnant women. It plays a potential role in guiding anti-D prophylaxis and pregnancy management in Chinese pregnant women.

Interdisciplinary approaches to predicting evolutionary biology.

Interdisciplinary approaches to predicting evolutionary biology.

The CellPhe toolkit for cell phenotyping using time-lapse imaging and pattern recognition

With phenotypic heterogeneity in whole cell populations widely recognised, the demand for quantitative and temporal analysis approaches to characterise single cell morphology and dynamics has increased. We present CellPhe, a pattern recognition toolkit for the unbiased characterisation of cellular phenotypes within time-lapse videos. CellPhe imports tracking information from multiple segmentation and tracking algorithms to provide automated cell phenotyping from different imaging modalities, including fluorescence. To maximise data quality for downstream analysis, our toolkit includes automated recognition and removal of erroneous cell boundaries induced by inaccurate tracking and segmentation. We provide an extensive list of features extracted from individual cell time series, with custom feature selection to identify variables that provide greatest discrimination for the analysis in question. Using ensemble classification for accurate prediction of cellular phenotype and clustering algorithms for the characterisation of heterogeneous subsets, we validate and prove adaptability using different cell types and experimental conditions.

A tractable physical model for the yeast polarity predicts epistasis and fitness.

Accurate phenotype prediction based on genetic information has numerous societal applications, such as crop design or cellular factories. Epistasis, when biological components interact, complicates modelling phenotypes from genotypes. Here we show an approach to mitigate this complication for polarity establishment in budding yeast, where mechanistic information is abundant. We coarse-grain molecular interactions into a so-called mesotype, which we combine with gene expression noise into a physical cell cycle model. First, we show with computer simulations that the mesotype allows validation of the most current biochemical polarity models by quantitatively matching doubling times. Second, the mesotype elucidates epistasis emergence as exemplified by evaluating the predicted mutational effect of key polarity protein Bem1p when combined with known interactors or under different growth conditions. This example also illustrates how unlikely evolutionary trajectories can become more accessible. The tractability of our biophysically justifiable approach inspires a road-map towards bottom-up modelling complementary to statistical inferences. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Predicting phenotypes of beef eating quality traits.

Introduction: Phenotype predictions of beef eating quality for individual animals could be used to allocate animals to longer and more expensive feeding regimes as they enter the feedlot if they are predicted to have higher eating quality, and to sort carcasses into consumer or market value categories. Phenotype predictions can include genetic effects (breed effects, heterosis and breeding value), predicted from genetic markers, as well as fixed effects such as days aged and carcass weight, hump height, ossification, and hormone growth promotant (HGP) status. Methods: Here we assessed accuracy of phenotype predictions for five eating quality traits (tenderness, juiciness, flavour, overall liking and MQ4) in striploins from 1701 animals from a wide variety of backgrounds, including Bos indicus and Bos taurus breeds, using genotypes and simple fixed effects including days aged and carcass weight. The genetic components were predicted based on 709k single nucleotide polymorphism (SNP) using BayesR model, which assumes some markers may have a moderate to large effect. Fixed effects in the prediction included principal components of the genomic relationship matrix, to account for breed effects, heterosis, days aged and carcass weight. Results and Discussion: A model which allowed breed effects to be captured in the SNP effects (e.g., not explicitly fitting these effects) tended to have slightly higher accuracies (0.43-0.50) compared to when these effects were explicitly fitted as fixed effects (0.42-0.49), perhaps because breed effects when explicitly fitted were estimated with more error than when incorporated into the (random) SNP effects. Adding estimates of effects of days aged and carcass weight did not increase the accuracy of phenotype predictions in this particular analysis. The accuracy of phenotype prediction for beef eating quality traits was sufficiently high that such predictions could be useful in predicting eating quality from DNA samples taken from an animal/carcass as it enters the processing plant, to enable optimal supply chain value extraction by sorting product into markets with different quality. The BayesR predictions identified several novel genes potentially associated with beef eating quality.

Transferability of the PRS estimates for height and BMI obtained from the European ethnic groups to the Western Russian populations.

Genetic data plays an increasingly important role in modern medicine. Decrease in the cost of sequencing with subsequent increase in imputation accuracy, and the accumulation of large amounts of high-quality genetic data enable the creation of polygenic risk scores (PRSs) to perform genotype-phenotype associations. The accuracy of phenotype prediction primarily depends on the overall trait heritability, Genome-wide association studies cohort size, and the similarity of genetic background between the base and the target cohort. Here we utilized 8,664 high coverage genomic samples collected across Russia by "Evogen", a Russian biomedical company, to evaluate the predictive power of PRSs based on summary statistics established on cohorts of European ancestry for basic phenotypic traits, namely height and BMI. We have demonstrated that the PRSs calculated for selected traits in three distinct Russian populations, recapitulate the predictive power from the original studies. This is evidence that GWAS summary statistics calculated on cohorts of European ancestry are transferable onto at least some ethnic groups in Russia.

Estimating the In Vivo Function of CYP2D6 Alleles through Population Pharmacokinetic Modeling of Brexpiprazole.

Accurate prediction of CYP2D6 phenotype from genotype information is important to support safe and efficacious pharmacotherapy with CYP2D6 substrates. To facilitate accurate CYP2D6 genotype-phenotype translation, there remains a need to investigate the enzyme activity associated with individual CYP2D6 alleles using large clinical data sets. This study aimed to quantify and compare the in vivo function of different CYP2D6 alleles through population pharmacokinetic (PopPK) modeling of brexpiprazole using data from 13 clinical studies. A PopPK model of brexpiprazole and its two metabolites, DM-3411 and DM-3412, was developed based on plasma concentration samples from 826 individuals. As the minor metabolite, DM-3412, is formed via CYP2D6, the metabolic ratio of DM-3412:brexpiprazole calculated from the PopPK parameter estimates was used as a surrogate measure of CYP2D6 activity. A CYP2D6 genotype-phenotype analysis based on 496 subjects showed that the CYP2D6*2 allele (n=183) was associated with only 10% enzyme activity relative to the wild-type allele (CYP2D6*1) and a low enzyme activity was consistently observed across genotypes containing CYP2D6*2. Among the decreased function alleles, the following enzyme activities relative to CYP2D6*1 were estimated: 23% for CYP2D6*9 (n=20), 32% for CYP2D6*10 (n=62), 64% for CYP2D6*14 (n=1), 4% for CYP2D6*17 (n=37), 4% for CYP2D6*29 (n=13), and 9% for CYP2D6*41 (n=64). These findings imply that a lower functional value would more accurately reflect the in vivo function of many reduced function CYP2D6 alleles in the metabolism of brexpiprazole. The low enzyme activity observed for CYP2D6*2, which has also been reported by others, suggests that the allele exhibits substrate-specific enzyme activity.

Genome-wide identification of genes enabling accurate prediction of hybrid performance from parents across environments and populations for gene-based breeding in maize

Accurate prediction of hybrid offspring complex trait phenotype from parents is paramount to enhanced plant breeding, animal breeding, and human medicine. Here we report genome-wide identification of genes enabling accurate prediction of hybrid offspring complex traits from parents using maize grain yield as the target trait. We identified 181 ZmF1GY genes enabling prediction of maize (Zea mays L.) F1 hybrid grain yield from parents and tested their utility and efficiency for predicting F1 hybrid grain yields from parents using their expressions, genic SNPs, and number of favorable alleles (NFAs), respectively. The ZmF1GY genes predicted hybrid grain yields from parents at an accuracy of 0.86, presented by correlation coefficient between predicted and observed phenotypes, within an environment, 0.74 across environments, and 0.64 across populations, outperforming genomic prediction by 27–406%, 23%, and 40%, respectively. Furthermore, we identified nine of the ZmF1GY genes containing SNPs or InDels in parents that increased or decreased hybrid grain yields by 14–46%. When the NFAs of these nine ZmF1GY genes were used for hybrid grain yield prediction from parents, they predicted hybrid grain yields at an accuracy of 0.79, outperforming genomic prediction by 21% that was based on up to tens of thousands of genome-wide SNPs. These results demonstrate the feasibility of developing a gene toolkit for a species enabling gene-based breeding across environments and populations that is much more powerful and efficient than current breeding, thereby helping secure the world’s food production. The methodology is applicable to all crops, livestock, and humans.

Accurate predictions of barley phenotypes using genomewide markers and environmental covariates

AbstractPredicting the performance of new plant genotypes under new environmental conditions could accelerate the development of locally adapted and climate resilient cultivars. Enabling these predictions may rely on extending the genomewide prediction framework to include environmental covariates (EC), but such models have generally been tested under limited, less breeding‐realistic contexts. Using a barley (Hordeum vulgare L.) multi‐environment dataset, our objectives were to compare multi‐environment prediction models and scenarios that target genotypes from different breeding generations, use different levels and timescales of ECs, and are applied to different agronomic and quality traits. When predicting the phenotypes of previously tested genotypes in untested environments, models that included the interaction of genomewide markers and pre‐selected in‐season ECs resulted in more accurate predictions (rMG or rMP) within (rMG = 0.56–0.94) and across (rMP = 0.63–0.87) environments; similar accuracy was achieved within (rMP = 0.46–0.89) and across (rMP = 0.87–0.95) locations when using only ECs from realistically available historical climate data. Shifting the prediction target to a distinct, untested offspring population slightly reduced model performance within environments or locations, but rMP across environments (rMP = 0.60–0.86) or locations (rMP = 0.87–0.94) remained very high. Though we achieved moderately high rMP for most traits in the challenging scenario of predicting the offspring population in holdout environments, the similarity between training and target environments, like that between populations, will be a limiting factor for enabling accurate predictions of new genotypes under new growing conditions.

Inbreeding Calculated with Runs of Homozygosity Suggests Chromosome-Specific Inbreeding Depression Regions in Line 1 Hereford.

Simple SummaryInbreeding depression is the decreased fitness of offspring of closely related individuals. It is a common problem in the livestock industry as a consequence of intense selection favoring few animals with higher genetic merit. The literature shows that inbred individuals may also exhibit decreased performance for economically relevant traits, such as growth in beef cattle. However, this relationship is not simple. The use of genetic markers can improve the detection of homozygous regions and improve the calculation of inbreeding coefficients. The working hypothesis of this study is that genomic homozygous regions can be used to predict animals’ performance. Moreover, homozygous regions across the genome differ in their effects. Our results indicate that the relationship between homozygosity and decreased performance varies by chromosome. Additionally, the inbreeding coefficient can be used to predict animals’ phenotype when calculated within a chromosome. This research warrants the development of a new approach in genomic selection to increase prediction accuracy by including information on homozygous regions. Additionally, management tools can be developed with the goal of breeding individuals with a purpose to decrease homozygosity only in the relevant regions of the genome.The goal of this study was to evaluate inbreeding in a closed beef cattle population and assess phenotype prediction accuracy using inbreeding information. Effects of inbreeding on average daily gain phenotype in the Line 1 Hereford cattle population were assessed in this study. Genomic data were used to calculate inbreeding based on runs of homozygosity (ROH), and pedigree information was used to calculate the probability of an allele being identical by descent. Prediction ability of phenotypes using inbreeding coefficients calculated based on pedigree information and runs of homozygosity over the whole genome was close to 0, even in the case of significant inbreeding coefficient effects. On the other hand, inbreeding calculated per individual chromosomes’ ROH yielded higher accuracies of prediction. Additionally, including only ROH from chromosomes with higher predicting ability further increased prediction accuracy. Phenotype prediction accuracy, inbreeding depression, and the effects of chromosome-specific ROHs varied widely across the genome. The results of this study suggest that inbreeding should be evaluated per individual regions of the genome. Moreover, mating schemes to avoid inbreeding depression should focus more on specific ROH with negative effects. Finally, using ROH as added information may increase prediction of the genetic merit of animals in a genomic selection program.

Multi-omics data integration by generative adversarial network.

Accurate disease phenotype prediction plays an important role in the treatment of heterogeneous diseases like cancer in the era of precision medicine. With the advent of high throughput technologies, more comprehensive multi-omics data is now available that can effectively link the genotype to phenotype. However, the interactive relation of multi-omics datasets makes it particularly challenging to incorporate different biological layers to discover the coherent biological signatures and predict phenotypic outcomes. In this study, we introduce omicsGAN, a generative adversarial network model to integrate two omics data and their interaction network. The model captures information from the interaction network as well as the two omics datasets and fuse them to generate synthetic data with better predictive signals. Large-scale experiments on The Cancer Genome Atlas breast cancer, lung cancer and ovarian cancer datasets validate that (i) the model can effectively integrate two omics data (e.g. mRNA and microRNA expression data) and their interaction network (e.g. microRNA-mRNA interaction network). The synthetic omics data generated by the proposed model has a better performance on cancer outcome classification and patients survival prediction compared to original omics datasets. (ii) The integrity of the interaction network plays a vital role in the generation of synthetic data with higher predictive quality. Using a random interaction network does not allow the framework to learn meaningful information from the omics datasets; therefore, results in synthetic data with weaker predictive signals. Source code is available at: https://github.com/CompbioLabUCF/omicsGAN. Supplementary data are available at Bioinformatics online.