Architecture of population-differentiated polymorphisms in the human genome

Reliable RNA-seq analysis from FFPE specimens as a means to accelerate cancer-related health disparities research

Whole transcriptome sequencing (WTS/ RNA-Seq) is a ubiquitous tool for investigating cancer biology. RNA isolated from frozen sources limits possible studies for analysis of associations with phenotypes or clinical variables requiring long-term follow-up. Although good correlations are reported in RNA-Seq data from paired frozen and formalin fixed paraffin embedded (FFPE) samples, uncertainties regarding RNA quality, methods of extraction, and data reliability are hurdles to utilization of archival samples. We compared three different platforms for performing RNA-seq using archival FFPE oropharyngeal squamous carcinoma (OPSCC) specimens stored up to 20 years, as part of an investigation of transcriptional profiles related to health disparities. We developed guidelines to purify DNA and RNA from FFPE tissue and perform downstream RNA-seq and DNA SNP arrays. RNA was extracted from 150 specimens, with an average yield of 401.8 ng/cm2 of tissue. Most samples yielded sufficient RNA reads >13,000 protein coding genes which could be used to differentiate HPV-associated from HPV-independent OPSCCs. Co-isolated DNA was used to identify reliably define patient ancestry which correlated well with patient-reported race. Utilizing the methods described in this study provides a robust, reliable, and standardized means of DNA & RNA extraction from FFPE as well as a means by which to assure the quality of the data generated. Optimized RNA extraction techniques, combined with robust bioinformatic approaches designed to optimize data homogenization, analysis and biological validation can revolutionize our ability to transcriptomically profile large solid tumor sets derived from ancestrally varied patient populations.

Reliable RNA-seq analysis from FFPE specimens as a means to accelerate cancer-related health disparities research

Whole transcriptome sequencing (WTS/ RNA-Seq) is a ubiquitous tool for investigating cancer biology. RNA isolated from frozen sources limits possible studies for analysis of associations with phenotypes or clinical variables requiring long-term follow-up. Although good correlations are reported in RNA-Seq data from paired frozen and formalin fixed paraffin embedded (FFPE) samples, uncertainties regarding RNA quality, methods of extraction, and data reliability are hurdles to utilization of archival samples. We compared three different platforms for performing RNA-seq using archival FFPE oropharyngeal squamous carcinoma (OPSCC) specimens stored up to 20 years, as part of an investigation of transcriptional profiles related to health disparities. We developed guidelines to purify DNA and RNA from FFPE tissue and perform downstream RNA-seq and DNA SNP arrays. RNA was extracted from 150 specimens, with an average yield of 401.8 ng/cm 2 of tissue. Most samples yielded sufficient RNA reads >13,000 protein coding genes which could be used to differentiate HPV-associated from HPV-independent OPSCCs. Co-isolated DNA was used to identify patient ancestry. Utilizing the methods described in this study provides a robust, reliable, and standardized means of DNA & RNA extraction from FFPE as well as a means by which to assure the quality of the data generated.

Artificial-Intelligence-Assisted Discovery of Genetic Factors for Precision Medicine of Antiplatelet Therapy in Diabetic Peripheral Artery Disease

An increased risk of cardiovascular events was identified in patients with peripheral artery disease (PAD). Clopidogrel is one of the most widely used antiplatelet medications. However, there are heterogeneous outcomes when clopidogrel is used to prevent cardiovascular events in PAD patients. Here, we use an artificial intelligence (AI)-assisted methodology to identify genetic factors potentially involved in the clopidogrel-resistant mechanism, which is currently unclear. Several discoveries can be pinpointed. Firstly, a high proportion (>50%) of clopidogrel resistance was found among diabetic PAD patients in Taiwan. Interestingly, our result suggests that platelet function test-guided antiplatelet therapy appears to reduce the post-interventional occurrence of major adverse cerebrovascular and cardiac events in diabetic PAD patients. Secondly, AI-assisted genome-wide association study of a single-nucleotide polymorphism (SNP) database identified a SNP signature composed of 20 SNPs, which are mapped into 9 protein-coding genes (SLC37A2, IQSEC1, WASHC3, PSD3, BTBD7, GLIS3, PRDM11, LRBA1, and CNR1). Finally, analysis of the protein connectivity map revealed that LRBA, GLIS3, BTBD7, IQSEC1, and PSD3 appear to form a protein interaction network. Intriguingly, the genetic factors seem to pinpoint a pathway related to endocytosis and recycling of P2Y12 receptor, which is the drug target of clopidogrel. Our findings reveal that a combination of AI-assisted discovery of SNP signatures and clinical parameters has the potential to develop an ethnic-specific precision medicine for antiplatelet therapy in diabetic PAD patients.

Variability in porcine microRNA genes and its association with mRNA expression and lipid phenotypes

BACKGROUND

Mature microRNAs (miRNAs) play an important role in repressing the expression of a wide range of mRNAs. The presence of polymorphic sites in miRNA genes and their corresponding 3'UTR binding sites can disrupt canonical conserved miRNA-mRNA pairings, and thus modify gene expression patterns. However, to date such polymorphic sites in miRNA genes and their association with gene expression phenotypes and complex traits are poorly characterized in pigs.

RESULTS

By analyzing whole-genome sequences from 120 pigs and wild boars from Europe and Asia, we identified 285 single nucleotide polymorphisms (SNPs) that map to miRNA loci, and 109,724 SNPs that are located in predicted 7mer-m8 miRNA binding sites within porcine 3'UTR. In porcine miRNA genes, SNP density is reduced compared with their flanking non-miRNA regions. By sequencing the genomes of five Duroc boars, we identified 12 miRNA SNPs that were subsequently genotyped in their offspring (N = 345, Lipgen population). Association analyses of miRNA SNPs with 38 lipid-related traits and hepatic and muscle microarray expression phenotypes recorded in the Lipgen population were performed. The most relevant detected association was between the genotype of the rs319154814 (G/A) SNP located in the apical loop of the ssc-miR-326 hairpin precursor and PPP1CC mRNA levels in the liver (q-value = 0.058). This result was subsequently confirmed by qPCR (P-value = 0.027). The rs319154814 (G/A) genotype was also associated with several fatty acid composition traits.

CONCLUSIONS

Our findings show a reduced variability of porcine miRNA genes, which is consistent with strong purifying selection, particularly in the seed region that plays a critical role in miRNA binding. Although it is generally assumed that SNPs mapping to the seed region are those with the most pronounced consequences on mRNA expression, we show that a SNP mapping to the apical region of ssc-miR-326 is significantly associated with hepatic mRNA levels of the PPP1CC gene, one of its predicted targets. Although experimental confirmation of such an interaction is reported in humans but not in pigs, this result highlights the need to further investigate the functional effects of miRNA polymorphisms that are located outside the seed region on gene expression in pigs.

Robust Performance of Potentially Functional SNPs in Machine Learning Models for the Prediction of Atorvastatin-Induced Myalgia

Statins can cause muscle symptoms resulting in poor adherence to therapy and increased cardiovascular risk. We hypothesize that combinations of potentially functional SNPs (pfSNPs), rather than individual SNPs, better predict myalgia in patients on atorvastatin. This study assesses the value of potentially functional single nucleotide polymorphisms (pfSNPs) and employs six machine learning algorithms to identify the combination of SNPs that best predict myalgia.

Methods: Whole genome sequencing of 183 Chinese, Malay and Indian patients from Singapore was conducted to identify genetic variants associated with atorvastatin induced myalgia. To adjust for confounding factors, demographic and clinical characteristics were also examined for their association with myalgia. The top factor, sex, was then used as a covariate in the whole genome association analyses. Variants that were highly associated with myalgia from this and previous studies were extracted, assessed for potential functionality (pfSNPs) and incorporated into six machine learning models. Predictive performance of a combination of different models and inputs were compared using the average cross validation area under ROC curve (AUC). The minimum combination of SNPs to achieve maximum sensitivity and specificity as determined by AUC, that predict atorvastatin-induced myalgia in most, if not all the six machine learning models was determined.

Results: Through whole genome association analyses using sex as a covariate, a larger proportion of pfSNPs compared to non-pf SNPs were found to be highly associated with myalgia. Although none of the individual SNPs achieved genome wide significance in univariate analyses, machine learning models identified a combination of 15 SNPs that predict myalgia with good predictive performance (AUC >0.9). SNPs within genes identified in this study significantly outperformed SNPs within genes previously reported to be associated with myalgia. pfSNPs were found to be more robust in predicting myalgia, outperforming non-pf SNPs in the majority of machine learning models tested.

Conclusion: Combinations of pfSNPs that were consistently identified by different machine learning models to have high predictive performance have good potential to be clinically useful for predicting atorvastatin-induced myalgia once validated against an independent cohort of patients.