yQTL Pipeline: A structured computational workflow for large scale quantitative trait loci discovery and downstream visualization.

Quantitative trait loci (QTL) denote regions of DNA whose variation is associated with variations in quantitative traits. QTL discovery is a powerful approach to understand how changes in molecular and clinical phenotypes may be related to DNA sequence changes. However, QTL discovery analysis encompasses multiple analytical steps and the processing of multiple input files, which can be laborious, error prone, and hard to reproduce if performed manually. To facilitate and automate large-scale QTL analysis, we developed the yQTL Pipeline, where the 'y' indicates the dependent quantitative variable being modeled. Prior to the association test, the pipeline supports the calculation or the direct input of pre-defined genome-wide principal components and genetic relationship matrix when applicable. User-specified covariates can also be provided. Depending on whether familial relatedness exists among the subjects, genome-wide association tests will be performed using either a linear mixed-effect model or a linear model. The options to run an ANOVA model or testing the interaction with a covariate are also available. Using the workflow management tool Nextflow, the pipeline parallelizes the analysis steps to optimize run-time and ensure results reproducibility. In addition, a user-friendly R Shiny App is developed to facilitate result visualization. It can generate Manhattan and Miami plots of phenotype traits, genotype-phenotype boxplots, and trait-QTL connection networks. We applied the yQTL Pipeline to analyze metabolomics profiles of blood serum from the New England Centenarians Study (NECS) participants. A total of 9.1M SNPs and 1,052 metabolites across 194 participants were analyzed. Using a p-value cutoff 5e-8, we found 14,983 mQTLs associated with 312 metabolites. The built-in parallelization of our pipeline reduced the run time from ~90 min to ~26 min. Visualization using the R Shiny App revealed multiple mQTLs shared across multiple metabolites. The yQTL Pipeline is available with documentation on GitHub at https://github.com/montilab/yQTLpipeline.

Heterogeneity analysis provides evidence for a genetically homogeneous subtype of bipolar-disorder.

Bipolar Disorder (BD) is a complex disease. It is heterogeneous, both at the phenotypic and genetic level, although the extent and impact of this heterogeneity is not fully understood. One way to assess this heterogeneity is to look for patterns in the subphenotype data, identify a more phenotypically homogeneous set of subjects, and perform a genome-wide association-study (GWAS) and subsequent secondary analyses restricted to this homogeneous subset. Because of the variability in how phenotypic data was collected by the various BD studies over the years, homogenizing the phenotypic data is a challenging task, and so is replication. As members of the Psychiatric Genomics Consortium (PGC), we have access to the raw genotypes of 18,711 BD cases and 29,738 controls. This amount of data makes it possible for us to set aside the intricacies of phenotype and allow the genetic data itself to determine which subjects define a homogeneous genetic subgroup. In this paper, we leverage recent advances in heterogeneity analysis to look for distinct homogeneous genetic BD subgroups (or biclusters) that manifest the broad phenotype we think of as Bipolar Disorder. As our data was generated by 27 studies and genotyped on a variety of platforms (OMEX, Affymetrix, Illumina), we use a biclustering algorithm capable of covariate-correction. Covariate-correction is critical if we wish to distinguish disease-related signals from those which are a byproduct of ancestry, study or genotyping platform. We rely on the raw genotyped data and do not include any data generated through imputation. We first apply this covariate-corrected biclustering algorithm to a cohort of 2524 BD cases and 4106 controls from the Bipolar Disease Research Network (BDRN: OMEX). We find evidence of genetic heterogeneity delineating a statistically significant bicluster comprising a subset of BD cases which exhibits a disease-specific pattern of differential-expression across a subset of SNPs. This pattern replicates across the remaining data-sets collected by the PGC containing 5781/8289 (OMEX), 3581/7591 (Illumina), and 6825/9752(Affymetrix) cases/controls, respectively. This bicluster includes subjects diagnosed with bipolar type-I, as well as subjects diagnosed with bipolar type-II. However, the bicluster is enriched for bipolar type-I over type-II and may represent a collection of correlated genetic risk-factors. By investigating the bicluster-informed polygenic-risk-scoring (PRS), we find that the disease-specific pattern highlighted by the bicluster can be leveraged to eliminate noise from our GWAS analyses and improve not only risk prediction, particularly when using only a relatively small subset (e.g., ~ 1%) of the available SNPs, but also SNP replication. Though our primary focus is only the analysis of disease-related signal, we also identify replicable control-related heterogeneity. Covariate-corrected biclustering of raw genetic data appears to be a promising route for untangling heterogeneity and identifying replicable homogeneous genetic subtypes of complex disease. It may also prove useful in identifying protective effects within the control group. This approach circumvents some of the difficulties presented by subphenotype data collected by meta-analyses or 23 andMe, e.g., missingness, assessment variation, and reliance on self-report.

Single-Cell Patch-Clamp/Proteomics of Human Alzheimer's Disease iPSC-Derived Excitatory Neurons Versus Isogenic Wild-Type Controls Suggests Novel Causation and Therapeutic Targets.

Standard single-cell (sc) proteomics of disease states inferred from multicellular organs or organoids cannot currently be related to single-cell physiology. Here, a scPatch-Clamp/Proteomics platform is developed on single neurons generated from hiPSCs bearing an Alzheimer's disease (AD) genetic mutation and compares them to isogenic wild-type controls. This approach provides both current and voltage electrophysiological data plus detailed proteomics information on single-cells. With this new method, the authors are able to observe hyperelectrical activity in the AD hiPSC-neurons, similar to that observed in the human AD brain, and correlate it to ≈1400 proteins detected at the single neuron level. Using linear regression and mediation analyses to explore the relationship between the abundance of individual proteins and the neuron's mutational and electrophysiological status, this approach yields new information on therapeutic targets in excitatory neurons not attainable by traditional methods. This combined patch-proteomics technique creates a new proteogenetic-therapeutic strategy to correlate genotypic alterations to physiology with protein expression in single-cells.

The relationship between 11 different polygenic longevity scores, parental lifespan, and disease diagnosis in the UK Biobank.

Large-scale genome-wide association studies (GWAS) strongly suggest that most traits and diseases have a polygenic component. This observation has motivated the development of disease-specific "polygenic scores (PGS)" that are weighted sums of the effects of disease-associated variants identified from GWAS that correlate with an individual's likelihood of expressing a specific phenotype. Although most GWAS have been pursued on disease traits, leading to the creation of refined "Polygenic Risk Scores" (PRS) that quantify risk to diseases, many GWAS have also been pursued on extreme human longevity, general fitness, health span, and other health-positive traits. These GWAS have discovered many genetic variants seemingly protective from disease and are often different from disease-associated variants (i.e., they are not just alternative alleles at disease-associated loci) and suggest that many health-positive traits also have a polygenic basis. This observation has led to an interest in "polygenic longevity scores (PLS)" that quantify the "risk" or genetic predisposition of an individual towards health. We derived 11 different PLS from 4 different available GWAS on lifespan and then investigated the properties of these PLS using data from the UK Biobank (UKB). Tests of association between the PLS and population structure, parental lifespan, and several cancerous and non-cancerous diseases, including death from COVID-19, were performed. Based on the results of our analyses, we argue that PLS are made up of variants not only robustly associated with parental lifespan, but that also contribute to the genetic architecture of disease susceptibility, morbidity, and mortality.

Emerging therapeutic drug monitoring technologies: considerations and opportunities in precision medicine.

In recent years, the development of sensor and wearable technologies have led to their increased adoption in clinical and health monitoring settings. One area that is in early, but promising, stages of development is the use of biosensors for therapeutic drug monitoring (TDM). Traditionally, TDM could only be performed in certified laboratories and was used in specific scenarios to optimize drug dosage based on measurement of plasma/blood drug concentrations. Although TDM has been typically pursued in settings involving medications that are challenging to manage, the basic approach is useful for characterizing drug activity. TDM is based on the idea that there is likely a clear relationship between plasma/blood drug concentration (or concentration in other matrices) and clinical efficacy. However, these relationships may vary across individuals and may be affected by genetic factors, comorbidities, lifestyle, and diet. TDM technologies will be valuable for enabling precision medicine strategies to determine the clinical efficacy of drugs in individuals, as well as optimizing personalized dosing, especially since therapeutic windows may vary inter-individually. In this mini-review, we discuss emerging TDM technologies and their applications, and factors that influence TDM including drug interactions, polypharmacy, and supplement use. We also discuss how using TDM within single subject (N-of-1) and aggregated N-of-1 clinical trial designs provides opportunities to better capture drug response and activity at the individual level. Individualized TDM solutions have the potential to help optimize treatment selection and dosing regimens so that the right drug and right dose may be matched to the right person and in the right context.

yQTL Pipeline: a structured computational workflow for large scale quantitative trait loci discovery and downstream visualization.

Quantitative trait loci (QTL) denote regions of DNA whose variation is associated with variations in quantitative traits. QTL discovery is a powerful approach to understand how changes in molecular and clinical phenotypes may be related to DNA sequence changes. However, QTL discovery analysis encompasses multiple analytical steps and the processing of multiple input files, which can be laborious, error prone, and hard to reproduce if performed manually. In order to facilitate and automate large-scale QTL analysis, we developed the yQTL Pipeline, where the 'y' indicates the dependent quantitative variable being modeled. Prior to genome-wide association test, the pipeline supports the calculation or the direct input of pre-defined genome-wide principal components and genetic relationship matrix when applicable. User-specified covariates can also be provided. Depending on whether familial relatedness exists among the subjects, genome-wide association tests will be performed using either a linear mixed-effect model or a linear model. Using the workflow management tool Nextflow, the pipeline parallelizes the analysis steps to optimize run-time and ensure results reproducibility. In addition, a user-friendly R Shiny App is developed to facilitate result visualization. Upon uploading the result file, it can generate Manhattan plots of user-selected phenotype traits and trait-QTL connection networks based on user-specified p-value thresholds. We applied the yQTL Pipeline to analyze metabolomics profiles of blood serum from the New England Centenarians Study (NECS) participants. A total of 9.1M SNPs and 1,052 metabolites across 194 participants were analyzed. Using a p-value cutoff 5e-8, we found 14,983 mQTLs cumulatively associated with 312 metabolites. The built-in parallelization of our pipeline reduced the run time from ~90 min to ~26 min. Visualization using the R Shiny App revealed multiple mQTLs shared across multiple metabolites. The yQTL Pipeline is available with documentation on GitHub at https://github.com/montilab/yQTL-Pipeline.

Pentadecanoic Acid (C15:0), an Essential Fatty Acid, Shares Clinically Relevant Cell-Based Activities with Leading Longevity-Enhancing Compounds.

Pentadecanoic acid (C15:0) is an essential odd-chain saturated fatty acid with broad activities relevant to protecting cardiometabolic, immune, and liver health. C15:0 activates AMPK and inhibits mTOR, both of which are core components of the human longevity pathway. To assess the potential for C15:0 to enhance processes associated with longevity and healthspan, we used human cell-based molecular phenotyping assays to compare C15:0 with three longevity-enhancing candidates: acarbose, metformin, and rapamycin. C15:0 (n = 36 activities in 10 of 12 cell systems) and rapamycin (n = 32 activities in 12 of 12 systems) had the most clinically relevant, dose-dependent activities. At their optimal doses, C15:0 (17 µM) and rapamycin (9 µM) shared 24 activities across 10 cell systems, including anti-inflammatory (e.g., lowered MCP-1, TNFα, IL-10, IL-17A/F), antifibrotic, and anticancer activities, which are further supported by previously published in vitro and in vivo studies. Paired with prior demonstrated abilities for C15:0 to target longevity pathways, hallmarks of aging, aging rate biomarkers, and core components of type 2 diabetes, heart disease, cancer, and nonalcoholic fatty liver disease, our results support C15:0 as an essential nutrient with activities equivalent to, or surpassing, leading longevity-enhancing candidate compounds.

Pathway-Specific Polygenic Risk Scores Correlate with Clinical Status and Alzheimer's Disease-Related Biomarkers.

BACKGROUND: APOE is the largest genetic risk factor for Alzheimer's disease (AD), but there is a substantial polygenic component. Polygenic risk scores (PRS) can summarize small effects across the genome but may obscure differential risk across molecular processes and pathways that contribute to heterogeneity of disease presentation. OBJECTIVE: We examined polygenic risk impacting specific AD-associated pathways and its relationship with clinical status and biomarkers of amyloid, tau, and neurodegeneration (A/T/N). METHODS: We analyzed data from 1,411 participants from the Alzheimer's Disease Neuroimaging Initiative (ADNI). We applied pathway analysis and clustering to identify AD-associated "pathway clusters" and construct pathway-specific PRSs (excluding the APOE region). We tested associations with diagnostic status, abnormal levels of amyloid and ptau, and hippocampal volume. RESULTS: Thirteen pathway clusters were identified, and eight pathway-specific PRSs were significantly associated with AD diagnosis. Amyloid-positivity was associated with endocytosis and fibril formation, response misfolded protein, and regulation protein tyrosine PRSs. Ptau positivity and hippocampal volume were both related to protein localization and mitophagy PRS, and ptau-positivity was also associated with an immune signaling PRS. A global AD PRS showed stronger associations with diagnosis and all biomarkers compared to pathway PRSs. CONCLUSIONS: Pathway PRS may contribute to understanding separable disease processes, but do not add significant power for predictive purposes. These findings demonstrate that AD-phenotypes may be preferentially associated with risk in specific pathways, and defining genetic risk along multiple dimensions may clarify etiological heterogeneity in AD. This approach to delineate pathway-specific PRS can be used to study other complex diseases.

Borrelia PeptideAtlas: A proteome resource of common Borrelia burgdorferi isolates for Lyme research.

Lyme disease, caused by an infection with the spirochete Borrelia burgdorferi, is the most common vector-borne disease in North America. B. burgdorferi strains harbor extensive genomic and proteomic variability and further comparison is key to understanding the spirochetes infectivity and biological impacts of identified sequence variants. To achieve this goal, both transcript and mass spectrometry (MS)-based proteomics was applied to assemble peptide datasets of laboratory strains B31, MM1, B31-ML23, infective isolates B31-5A4, B31-A3, and 297, and other public datasets, to provide a publicly available Borrelia PeptideAtlas http://www.peptideatlas.org/builds/borrelia/. Included is information on total proteome, secretome, and membrane proteome of these B. burgdorferi strains. Proteomic data collected from 35 different experiment datasets, with a total of 855 mass spectrometry runs, identified 76,936 distinct peptides at a 0.1% peptide false-discovery-rate, which map to 1,221 canonical proteins (924 core canonical and 297 noncore canonical) and covers 86% of the total base B31 proteome. The diverse proteomic information from multiple isolates with credible data presented by the Borrelia PeptideAtlas can be useful to pinpoint potential protein targets which are common to infective isolates and may be key in the infection process.

Lifespan-extending interventions induce consistent patterns of fatty acid oxidation in mouse livers.

Aging manifests as progressive deteriorations in homeostasis, requiring systems-level perspectives to investigate the gradual molecular dysregulation of underlying biological processes. Here, we report systemic changes in the molecular regulation of biological processes under multiple lifespan-extending interventions. Differential Rank Conservation (DIRAC) analyses of mouse liver proteomics and transcriptomics data show that mechanistically distinct lifespan-extending interventions (acarbose, 17α-estradiol, rapamycin, and calorie restriction) generally tighten the regulation of biological modules. These tightening patterns are similar across the interventions, particularly in processes such as fatty acid oxidation, immune response, and stress response. Differences in DIRAC patterns between proteins and transcripts highlight specific modules which may be tightened via augmented cap-independent translation. Moreover, the systemic shifts in fatty acid metabolism are supported through integrated analysis of liver transcriptomics data with a mouse genome-scale metabolic model. Our findings highlight the power of systems-level approaches for identifying and characterizing the biological processes involved in aging and longevity.

Genetic and Protein Network Underlying the Convergence of Rett-Syndrome-like (RTT-L) Phenotype in Neurodevelopmental Disorders.

Mutations of the X-linked gene encoding methyl-CpG-binding protein 2 (MECP2) cause classical forms of Rett syndrome (RTT) in girls. A subset of patients who are recognized to have an overlapping neurological phenotype with RTT but are lacking a mutation in a gene that causes classical or atypical RTT can be described as having a 'Rett-syndrome-like phenotype (RTT-L). Here, we report eight patients from our cohort diagnosed as having RTT-L who carry mutations in genes unrelated to RTT. We annotated the list of genes associated with RTT-L from our patient cohort, considered them in the light of peer-reviewed articles on the genetics of RTT-L, and constructed an integrated protein-protein interaction network (PPIN) consisting of 2871 interactions connecting 2192 neighboring proteins among RTT- and RTT-L-associated genes. Functional enrichment analysis of RTT and RTT-L genes identified a number of intuitive biological processes. We also identified transcription factors (TFs) whose binding sites are common across the set of RTT and RTT-L genes and appear as important regulatory motifs for them. Investigation of the most significant over-represented pathway analysis suggests that HDAC1 and CHD4 likely play a central role in the interactome between RTT and RTT-L genes.

The human pangenome reference anticipates equitable and fundamental genomic insights.

For the past few years, researchers in the Human Pangenome Reference Consortium (HPRC) have been working to catalog almost all human genomic diversity. Frazer and Schork preview an article recently published in Nature, "A draft human pangenome reference,"1 which represents the initial release of 47 fully phased diploid assemblies of genomes of individuals with diverse ancestries.

A CRISPR-enhanced metagenomic NGS test to improve pandemic preparedness.

The lack of preparedness for detecting and responding to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pathogen (i.e., COVID-19) has caused enormous harm to public health and the economy. Testing strategies deployed on a population scale at day zero, i.e., the time of the first reported case, would be of significant value. Next-generation sequencing (NGS) has such capabilities; however, it has limited detection sensitivity for low-copy-number pathogens. Here, we leverage the CRISPR-Cas9 system to effectively remove abundant sequences not contributing to pathogen detection and show that NGS detection sensitivity of SARS-CoV-2 approaches that of RT-qPCR. The resulting sequence data can also be used for variant strain typing, co-infection detection, and individual human host response assessment, all in a single molecular and analysis workflow. This NGS work flow is pathogen agnostic and, therefore, has the potential to transform how large-scale pandemic response and focused clinical infectious disease testing are pursued in the future.

Distinct longevity mechanisms across and within species and their association with aging.

Lifespan varies within and across species, but the general principles of its control remain unclear. Here, we conducted multi-tissue RNA-seq analyses across 41 mammalian species, identifying longevity signatures and examining their relationship with transcriptomic biomarkers of aging and established lifespan-extending interventions. An integrative analysis uncovered shared longevity mechanisms within and across species, including downregulated Igf1 and upregulated mitochondrial translation genes, and unique features, such as distinct regulation of the innate immune response and cellular respiration. Signatures of long-lived species were positively correlated with age-related changes and enriched for evolutionarily ancient essential genes, involved in proteolysis and PI3K-Akt signaling. Conversely, lifespan-extending interventions counteracted aging patterns and affected younger, mutable genes enriched for energy metabolism. The identified biomarkers revealed longevity interventions, including KU0063794, which extended mouse lifespan and healthspan. Overall, this study uncovers universal and distinct strategies of lifespan regulation within and across species and provides tools for discovering longevity interventions.

Exploring the genetic heterogeneity of Alzheimer's disease: Evidence for genetic subtypes.

Background: Alzheimer's disease (AD) exhibits heterogeneity in cognitive impairment, atrophy, and pathological accumulation, suggesting the potential existence of subtypes. AD is under substantial genetic influence, thus identifying systematic variation in genetic risk may provide insights into disease origins. Objective: We investigated genetic heterogeneity in AD risk through a multi-step analysis. Methods: We performed principal component analysis (PCA) on AD-associated variants in the UK Biobank (AD cases=2,739, controls=5,478) to assess the presence of structured genetic heterogeneity. Subsequently, a biclustering algorithm searched for distinct disease-specific genetic signatures among subsets of cases. Replication tests were conducted using the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset (AD cases=500, controls=470). We categorized a separate set of ADNI individuals with mild cognitive impairment (MCI; n=399) into genetic subtypes and examined cognitive, amyloid, and tau trajectories. Results: PCA revealed three distinct clusters ("constellations") within AD-associated variants containing a mixture of cases and controls, reflecting disease-relevant structure. We found two disease-specific biclusters among AD cases. Pathway analysis linked bicluster-associated variants to neuron morphogenesis and outgrowth, including genes related to cellular components and development-modulating factors. Both disease-relevant and disease-specific structure replicated in ADNI. Individuals with genetic signatures resembling bicluster 2 exhibited increased CSF p-tau and cognitive decline over time. Conclusions: This study unveils a hierarchical structure of AD genetic risk. Disease-relevant constellations may represent differential biological vulnerability that is itself not sufficient to increase risk. Biclusters may represent distinct AD genetic subtypes. This structure replicates in an independent dataset and relates to differential pathological accumulation and cognitive decline over time.

Pathway-specific polygenic risk scores correlate with clinical status and Alzheimer's-related biomarkers.

Background: APOE is the largest genetic risk factor for sporadic Alzheimer's disease (AD), but there is a substantial polygenic component as well. Polygenic risk scores (PRS) can summarize small effects across the genome but may obscure differential risk associated with different molecular processes and pathways. Variability at the genetic level may contribute to the extensive phenotypic heterogeneity of Alzheimer's disease (AD). Here, we examine polygenic risk impacting specific pathways associated with AD and examined its relationship with clinical status and AD biomarkers of amyloid, tau, and neurodegeneration (A/T/N). Methods: A total of 1,411 participants from the Alzheimer's Disease Neuroimaging Initiative (ADNI) with genotyping data were included. Sets of variants identified from a pathway analysis of AD GWAS summary statistics were combined into clusters based on their assigned pathway. We constructed pathway-specific PRSs for each participant and tested their associations with diagnostic status (AD vs cognitively normal), abnormal levels of amyloid and ptau (positive vs negative), and hippocampal volume. The APOE region was excluded from all PRSs, and analyses controlled for APOE -ε4 carrier status. Results: Thirteen pathway clusters were identified relating to categories such as immune response, amyloid precursor processing, protein localization, lipid transport and binding, tyrosine kinase, and endocytosis. Eight pathway-specific PRSs were significantly associated with AD dementia diagnosis. Amyloid-positivity was associated with endocytosis and fibril formation, response misfolded protein, and regulation protein tyrosine PRSs. Ptau positivity and hippocampal volume were both related to protein localization and mitophagy PRS, and ptau positivity was additionally associated with an immune signaling PRS. A global AD PRS showed stronger associations with diagnosis and all biomarkers compared to pathway PRSs, suggesting a strong synergistic effect of all loci contributing to the global AD PRS. Conclusions: Pathway PRS may contribute to understanding separable disease processes, but do not appear to add significant power for predictive purposes. These findings demonstrate that, although genetic risk for AD is widely distributed, AD-phenotypes may be preferentially associated with risk in specific pathways. Defining genetic risk along multiple dimensions at the individual level may help clarify the etiological heterogeneity in AD.

A time-series analysis of blood-based biomarkers within a 25-year longitudinal dolphin cohort.

Causal interactions and correlations between clinically-relevant biomarkers are important to understand, both for informing potential medical interventions as well as predicting the likely health trajectory of any individual as they age. These interactions and correlations can be hard to establish in humans, due to the difficulties of routine sampling and controlling for individual differences (e.g., diet, socio-economic status, medication). Because bottlenose dolphins are long-lived mammals that exhibit several age-related phenomena similar to humans, we analyzed data from a well controlled 25-year longitudinal cohort of 144 dolphins. The data from this study has been reported on earlier, and consists of 44 clinically relevant biomarkers. This time-series data exhibits three starkly different influences: (A) directed interactions between biomarkers, (B) sources of biological variation that can either correlate or decorrelate different biomarkers, and (C) random observation-noise which combines measurement error and very rapid fluctuations in the dolphin's biomarkers. Importantly, the sources of biological variation (type-B) are large in magnitude, often comparable to the observation errors (type-C) and larger than the effect of the directed interactions (type-A). Attempting to recover the type-A interactions without accounting for the type-B and type-C variation can result in an abundance of false-positives and false-negatives. Using a generalized regression which fits the longitudinal data with a linear model accounting for all three influences, we demonstrate that the dolphins exhibit many significant directed interactions (type-A), as well as strong correlated variation (type-B), between several pairs of biomarkers. Moreover, many of these interactions are associated with advanced age, suggesting that these interactions can be monitored and/or targeted to predict and potentially affect aging.

Nonalcoholic fatty liver disease risk and histologic severity are associated with genetic polymorphisms in children.

BACKGROUND AND AIMS: NAFLD is the most common chronic liver disease in children. Large pediatric studies identifying single nucleotide polymorphisms (SNPs) associated with risk and histologic severity of NAFLD are limited. Study aims included investigating SNPs associated with risk for NAFLD using family trios and association of candidate alleles with histologic severity. APPROACH AND RESULTS: Children with biopsy-confirmed NAFLD were enrolled from the NASH Clinical Research Network. The Expert Pathology Committee reviewed liver histology. Genotyping was conducted with allele-specific primers for 60 candidate SNPs. Parents were enrolled for trio analysis. To assess risk for NAFLD, the transmission disequilibrium test was conducted in trios. Among cases, regression analysis assessed associations with histologic severity. A total of 822 children with NAFLD had mean age 13.2 years (SD 2.7) and mean ALT 101 U/L (SD 90). PNPLA3 (rs738409) demonstrated the strongest risk ( p = 2.24 × 10 -14 ) for NAFLD. Among children with NAFLD, stratifying by PNPLA3 s738409 genotype, the variant genotype associated with steatosis ( p = 0.005), lobular ( p = 0.03) and portal inflammation ( p = 0.002). Steatosis grade associated with TM6SF2 ( p = 0.0009), GCKR ( p = 0.0032), PNPLA3 rs738409 ( p = 0.0053), and MTTP ( p = 0.0051). Fibrosis stage associated with PARVB rs6006473 ( p = 0.0001), NR1I2 ( p = 0.0021), ADIPOR2 ( p = 0.0038), and OXTR ( p = 0.0065). PNPLA3 rs738409 ( p = 0.0002) associated with borderline zone 1 NASH. CONCLUSIONS: This study demonstrated disease-associated SNPs in children with NAFLD. In particular, rs6006473 was highly associated with severity of fibrosis. These hypothesis-generating results support future mechanistic studies of development of adverse outcomes such as fibrosis and generation of therapeutic targets for NAFLD in children.

Interrogating the Human Diplome: Computational Methods, Emerging Applications, and Challenges.

Human DNA sequencing protocols have revolutionized human biology, biomedical science, and clinical practice, but still have very important limitations. One limitation is that most protocols do not separate or assemble (i.e., "phase") the nucleotide content of each of the maternally and paternally derived chromosomal homologs making up the 22 autosomal pairs and the chromosomal pair making up the pseudo-autosomal region of the sex chromosomes. This has led to a dearth of studies and a consequent underappreciation of many phenomena of fundamental importance to basic and clinical genomic science. We discuss a few protocols for obtaining phase information as well as their limitations, including those that could be used in tumor phasing settings. We then describe a number of biological and clinical phenomena that require phase information. These include phenomena that require precise knowledge of the nucleotide sequence in a chromosomal segment from germline or somatic cells, such as DNA binding events, and insight into unique cis vs. trans-acting functionally impactful variant combinations-for example, variants implicated in a phenotype governed by compound heterozygosity. In addition, we also comment on the need for reliable and consensus-based diploid-context computational workflows for variant identification as well as the need for laboratory-based functional verification strategies for validating cis vs. trans effects of variant combinations. We also briefly describe available resources, example studies, as well as areas of further research, and ultimately argue that the science behind the study of human diploidy, referred to as "diplomics," which will be enabled by nucleotide-level resolution of phased genomes, is a logical next step in the analysis of human genome biology.