Assigning function to SNPs: Considerations when interpreting genetic variation.

Assigning function to single nucleotide polymorphisms (SNPs) to understand the mechanisms that link genetic and phenotypic variation and disease is an area of intensive research that is necessary to contribute to the continuing development of precision medicine. However, despite the apparent simplicity that is captured in the name SNP - 'single nucleotide' changes are not easy to functionally characterize. This complexity arises from multiple features of the genome including the fact that function is development and environment specific. As such, we are often fooled by our terminology and underlying assumptions that there is a single function for a SNP. Here we discuss some of what is known about SNPs, their functions and how we can go about characterizing them.

Spatially constrained gene regulation identifies key genetic contributions of preeclampsia, hypertension, and proteinuria†.

Preeclampsia (PE) is a relatively common but severe pregnancy disorder (with very limited effective treatments) characterized by hypertension (HTN) and usually proteinuria (PRO) or other organ damage. Genome-wide association studies (GWAS) of PE, HTN, and PRO have mostly identified risk loci single nucleotide polymorphisms (SNPs) located in noncoding genomic regions, likely impacting the regulation of distal gene expression. The latest GWAS associated (P < 1 × 10-6) SNPs to PE (n = 25), HTN (n = 1926), and PRO (n = 170). Our algorithmic analysis (CoDeS3D) used chromatin connection data (Hi-C) derived from 70 cell lines followed by analysis of two expression quantitative trail loci (eQTL) cohorts: GTEx (838 donors, 54 tissues, totaling 15 253 samples) and DICE (91 donors, 13 blood tissue types). We identified spatially constrained eQTLs which implicate gene targets in PE (n = 16), HTN (n = 3561), and PRO (n = 335). By overlapping these target genes and their molecular pathways (protein-protein interaction networks), we identified shared functional impacts between PE and HTN, which are significantly enriched for regulatory interactions which target genes intolerant to loss-of-function mutations. While the disease-associated SNP loci mostly do not overlap, the regulatory signals (target genes and pathways) overlap, informing on PE risk mechanisms. This demonstrates a model in which genetic predisposition to HTN and PRO lays a molecular groundwork toward risk for PE pathogenesis. This overlap at the gene regulatory network level identifies possible shared therapeutic targets for future study.

Establishing gene regulatory networks from Parkinson's disease risk loci.

The latest meta-analysis of genome-wide association studies identified 90 independent variants across 78 genomic regions associated with Parkinson's disease, yet the mechanisms by which these variants influence the development of the disease remains largely elusive. To establish the functional gene regulatory networks associated with Parkinson's disease risk variants, we utilized an approach combining spatial (chromosomal conformation capture) and functional (expression quantitative trait loci) data. We identified 518 genes subject to regulation by 76 Parkinson's variants across 49 tissues, whicih encompass 36 peripheral and 13 CNS tissues. Notably, one-third of these genes were regulated via trans-acting mechanisms (distal; risk locus-gene separated by >1 Mb, or on different chromosomes). Of particular interest is the identification of a novel trans-expression quantitative trait loci-gene connection between rs10847864 and SYNJ1 in the adult brain cortex, highlighting a convergence between familial studies and Parkinson's disease genome-wide association studies loci for SYNJ1 (PARK20) for the first time. Furthermore, we identified 16 neurodevelopment-specific expression quantitative trait loci-gene regulatory connections within the foetal cortex, consistent with hypotheses suggesting a neurodevelopmental involvement in the pathogenesis of Parkinson's disease. Through utilizing Louvain clustering we extracted nine significant and highly intraconnected clusters within the entire gene regulatory network. The nine clusters are enriched for specific biological processes and pathways, some of which have not previously been associated with Parkinson's disease. Together, our results not only contribute to an overall understanding of the mechanisms and impact of specific combinations of Parkinson's disease variants, but also highlight the potential impact gene regulatory networks may have when elucidating aetiological subtypes of Parkinson's disease.

Do Cutting Kinematics Change as Boys Mature? A Longitudinal Cohort Study of High-School Athletes.

OBJECTIVE: Examine longitudinal changes in trunk, hip, and knee kinematics in maturing boys during an unanticipated cutting task. DESIGN: Prospective cohort study. SETTING: Biomechanical laboratory. PARTICIPANTS: Forty-two high-school male basketball, volleyball, and soccer athletes. ASSESSMENT OF RISK FACTORS: Trunk, hip, and knee range-of-motion (RoM), peak angles, and angles at initial contact during an unanticipated 45 degrees sidestep cutting task were estimated using laboratory-based three-dimensional optoelectronic motion capture. Maturation was classified using a modified Pubertal Maturational Observational Scale (PMOS) into prepubertal, midpubertal, or postpubertal stages. MAIN OUTCOME MEASURES: Trunk total RoM in frontal, sagittal, and transverse planes; peak trunk flexion, right lateral flexion and right rotation angles; hip total RoM in frontal, sagittal, and transverse planes; hip flexion angle at initial contact; peak hip flexion and adduction angles; knee total RoM in frontal, sagittal, and transverse planes; knee flexion angle at initial contact; peak knee flexion and abduction angles. RESULTS: As boys matured, there was a decrease in hip sagittal-plane RoM (49.02 degrees to 43.45 degrees, Benjamini-Hochberg adjusted P = 0.027), hip flexion at initial contact (29.33 degrees to 23.08 degrees, P = 0.018), and peak hip flexion (38.66 degrees to 32.71 degrees, P = 0.046), and an increase in trunk contralateral rotation (17.47 degrees to 25.05 degrees, P = 0.027). CONCLUSIONS: Maturing male athletes adopted a more erect cutting strategy that is associated with greater knee joint loading. Knee kinematic changes that increase knee joint loading were not observed in this cohort.

Copper-mediated oxidation of imidazopyrazinones inhibits marine luciferase activity.

The development of bioluminescence-based tools has seen steady growth in the field of chemical biology over the past few decades ranging in uses from reporter genes to assay development and targeted imaging. More recently, coelenterazine-utilizing luciferases such as Gaussia, Renilla, and the engineered nano-luciferases have been utilized due to their intense luminescence relative to firefly luciferin/luciferase. The emerging importance of these systems warrants investigations into the components that affect their light production. Previous work has reported that one marine luciferase, Gaussia, is potently inhibited by copper salt. The mechanism for inhibition was not elucidated but was hypothesized to occur via binding to the enzyme. In this study, we provide the first report of a group of nonhomologous marine luciferases also exhibiting marked decreases in light emission in the presence of copper (II). We investigate the mechanism of action behind this inhibition and demonstrate that the observed copper inhibition does not stem from a luciferase interaction but rather the chemical oxidation of imidazopyrazinone luciferins generating inert, dehydrated luciferins.

A new method for determining ribosomal DNA copy number shows differences between Saccharomyces cerevisiae populations.

Ribosomal DNA genes (rDNA) encode the major ribosomal RNAs and in eukaryotes typically form tandem repeat arrays. Species have characteristic rDNA copy numbers, but there is substantial intra-species variation in copy number that results from frequent rDNA recombination. Copy number differences can have phenotypic consequences, however difficulties in quantifying copy number mean we lack a comprehensive understanding of how copy number evolves and the consequences. Here we present a genomic sequence read approach to estimate rDNA copy number based on modal coverage to help overcome limitations with existing mean coverage-based approaches. We validated our method using Saccharomyces cerevisiae strains with known rDNA copy numbers. Application of our pipeline to a global sample of S. cerevisiae isolates showed that different populations have different rDNA copy numbers. Our results demonstrate the utility of the modal coverage method, and highlight the high level of rDNA copy number variation within and between populations.

Genomic dissection of 43 serum urate-associated loci provides multiple insights into molecular mechanisms of urate control.

High serum urate is a prerequisite for gout and associated with metabolic disease. Genome-wide association studies (GWAS) have reported dozens of loci associated with serum urate control; however, there has been little progress in understanding the molecular basis of the associated loci. Here, we employed trans-ancestral meta-analysis using data from European and East Asian populations to identify 10 new loci for serum urate levels. Genome-wide colocalization with cis-expression quantitative trait loci (eQTL) identified a further five new candidate loci. By cis- and trans-eQTL colocalization analysis, we identified 34 and 20 genes, respectively, where the causal eQTL variant has a high likelihood that it is shared with the serum urate-associated locus. One new locus identified was SLC22A9 that encodes organic anion transporter 7 (OAT7). We demonstrate that OAT7 is a very weak urate-butyrate exchanger. Newly implicated genes identified in the eQTL analysis include those encoding proteins that make up the dystrophin complex, a scaffold for signaling proteins and transporters at the cell membrane; MLXIP that, with the previously identified MLXIPL, is a transcription factor that may regulate serum urate via the pentose-phosphate pathway and MRPS7 and IDH2 that encode proteins necessary for mitochondrial function. Functional fine mapping identified six loci (RREB1, INHBC, HLF, UBE2Q2, SFMBT1 and HNF4G) with colocalized eQTL containing putative causal SNPs. This systematic analysis of serum urate GWAS loci identified candidate causal genes at 24 loci and a network of previously unidentified genes likely involved in control of serum urate levels, further illuminating the molecular mechanisms of urate control.

Can adult polygenic scores improve prediction of body mass index in childhood?

BACKGROUND/OBJECTIVES: Modelling genetic pre-disposition may identify children at risk of obesity. However, most polygenic scores (PGSs) have been derived in adults, and lack validation during childhood. This study compared the utility of existing large-scale adult-derived PGSs to predict common anthropometric traits (body mass index (BMI), waist circumference, and body fat) in children and adults, and examined whether childhood BMI prediction could be improved by combining PGSs and non-genetic factors (maternal and earlier child BMI). SUBJECTS/METHODS: Participants (n = 1365 children, and n = 2094 adults made up of their parents) were drawn from the Longitudinal Study of Australian Children. Children were weighed and measured every two years from 0-1 to 12-13 years, and adults were measured or self-reported measurements were obtained concurrently (average analysed). Participants were genotyped from blood or oral samples, and PGSs were derived based on published genome-wide association studies. We used linear regression to compare the relative utility of these PGSs to predict their respective traits at different ages. RESULTS: BMI PGSs explained up to 12% of child BMI z-score variance in 10-13 year olds, compared with up to 15% in adults. PGSs for waist circumference and body fat explained less variance (up to 8%). An interaction between BMI PGSs and puberty (p = 0.001-0.002) suggests the effect of some variants may differ across the life course. Individual BMI measures across childhood predicted 10-60% of the variance in BMI at 12-13 years, and maternal BMI and BMI PGS each added 1-9% above this. CONCLUSION: Adult-derived PGSs for BMI, particularly those derived by modelling between-variant interactions, may be useful for predicting BMI during adolescence with similar accuracy to that obtained in adulthood. The level of precision presented here to predict BMI during childhood may be relevant to public health, but is likely to be less useful for individual clinical purposes.

A de novo ACTB gene pathogenic variant in identical twins with phenotypic variation for hydrops and jejunal atresia.

The beta-actin gene (ACTB) encodes a ubiquitous cytoskeletal protein, essential for embryonic development in humans. De novo heterozygous missense variants in the ACTB are implicated in causing Baraitser-Winter cerebrofrontofacial syndrome (BWCFFS; MIM#243310). ACTB pathogenic variants are rarely associated with intestinal malformations. We report on a rare case of monozygotic twins presenting with proximal small bowel atresia and hydrops in one, and apple-peel bowel atresia and laryngeal dysgenesis in the other. The twin with hydrops could not be resuscitated. Intensive and surgical care was provided to the surviving twin. Rapid trio genome sequencing identified a de novo missense variant in ACTB (NM_00101.3:c.1043C>T; p.(Ser348Leu)) that guided the care plan. The identical variant subsequently was identified in the demised twin. To characterize the functional effect, the variant was recreated as a pseudoheterozygote in a haploid wild-type S. cerevisiae strain. There was an obvious growth defect of the yACT1S348L/WT pseudoheterozygote compared to a yACT1WT/WT strain when grown at 22°C but not when grown at 30°C, consistent with the yACT1 S348L variant having a functional defect that is dominant over the wild-type allele. The functional results provide supporting evidence that the Ser348Leu variant is likely to be a pathogenic variant, including being associated with intestinal malformations in BWCFFS, and can demonstrate variable expressivity within monozygotic twins.

Development of an ATP-independent bioluminescent probe for detection of extracellular hydrogen peroxide.

This work reports a new ATP-independent bioluminescent probe (bor-DTZ) for detecting hydrogen peroxide that is compatible with the Nanoluciferase enzyme. The probe is designed with an arylboronate ester protecting group appended to a diphenylterazine core via a self-immolative phenolate linker. Reaction with hydrogen peroxide reveals diphenylterazine, which can then react with Nanoluciferase to produce a detectable bioluminescent signal. Bor-DTZ shows a dose-dependent response to hydrogen peroxide and selectivity over other biologically relevant reactive oxygen species and can be applied to detect either intra- or extracellular species. We further demonstrate the ability of this platform to monitor fluxes in extracellular hydrogen peroxide in a breast cancer cell line in response to the anticancer treatment, cisplatin.

Low tolerance for transcriptional variation at cohesin genes is accompanied by functional links to disease-relevant pathways.

BACKGROUND: The cohesin complex plays an essential role in genome organisation and cell division. A full complement of the cohesin complex and its regulators is important for normal development, since heterozygous mutations in genes encoding these components can be sufficient to produce a disease phenotype. The implication that genes encoding the cohesin subunits or cohesin regulators must be tightly controlled and resistant to variability in expression has not yet been formally tested. METHODS: Here, we identify spatial-regulatory connections with potential to regulate expression of cohesin loci (Mitotic: SMC1A, SMC3, STAG1, STAG2, RAD21/RAD21-AS; Meiotic: SMC1B, STAG3, REC8, RAD21L1), cohesin-ring support genes (NIPBL, MAU2, WAPL, PDS5A, PDS5B) and CTCF, including linking their expression to that of other genes. We searched the genome-wide association studies (GWAS) catalogue for SNPs mapped or attributed to cohesin genes by GWAS (GWAS-attributed) and the GTEx catalogue for SNPs mapped to cohesin genes by cis-regulatory variants in one or more of 44 tissues across the human body (expression quantitative trail locus-attributed). RESULTS: Connections that centre on the cohesin ring subunits provide evidence of coordinated regulation that has little tolerance for perturbation. We used the CoDeS3D SNP-gene attribution methodology to identify transcriptional changes across a set of genes coregulated with the cohesin loci that include biological pathways such as extracellular matrix production and proteasome-mediated protein degradation. Remarkably, many of the genes that are coregulated with cohesin loci are themselves intolerant to loss-of-function. CONCLUSIONS: The results highlight the importance of robust regulation of cohesin genes and implicate novel pathways that may be important in the human cohesinopathy disorders.

Needs assessment for health service design for people with back pain in a hospital setting: A qualitative study.

BACKGROUND: There is a need for effective health service solutions to provide greater structure and support for implementing evidence-based practice in back pain care. Patient involvement in developing these solutions is crucial to increase relevance, acceptability and uptake. OBJECTIVES: To determine patients' perceived needs and barriers to best-practice back pain care, and potential solutions to better address care needs. The study is the third in a series of needs assessment studies feeding into the 'idea generation' for service design in a large teaching hospital in a culturally and linguistically diverse community in metropolitan Sydney, Australia. DESIGN: We conducted a combination of focus groups and in-depth interviews using an interpretive description approach. We used inductive thematic analysis to identify the main themes. SETTING AND PARTICIPANTS: We purposively sampled patients with diverse characteristics from the neurosurgery and physiotherapy outpatient clinics, in particular those whose primary language was English, Arabic, Persian or Mandarin. Non-English audio recordings were translated and transcribed by bilingual researchers. RESULTS: There were 24 participants (focus groups = 9; individual interviews = 15) when data saturation was reached. The analysis identified three key themes with several subthemes around what service designers needed to understand in helping people with back pain in this setting: (1) This is who I am; (2) It's not working for me; and (3) What I think I need. DISCUSSION AND CONCLUSION: This study highlights that perceived unmet needs of patients are underpinned by unhelpful beliefs about the causes of and solutions for back pain, misaligned care expectations, unclear expectations of the hospital role and fragmentations in the health system. To design and implement a service that can deliver better back pain care, several solutions need to be integrated around: developing new resources that challenge unhelpful beliefs and set realistic expectations; improving access to education and self-management resources; focusing on individualized care; using a collaborative multidisciplinary approach within the hospital; and better connecting with and directing primary health care services. PATIENT OR PUBLIC CONTRIBUTION: A consumer representative of the Western Sydney Local Health District provided input during study conceptualisation and is duly recognized in the Acknowledgements section.

Cohesin facilitates zygotic genome activation in zebrafish.

At zygotic genome activation (ZGA), changes in chromatin structure are associated with new transcription immediately following the maternal-to-zygotic transition (MZT). The nuclear architectural proteins cohesin and CCCTC-binding factor (CTCF) contribute to chromatin structure and gene regulation. We show here that normal cohesin function is important for ZGA in zebrafish. Depletion of the cohesin subunit Rad21 delays ZGA without affecting cell cycle progression. In contrast, CTCF depletion has little effect on ZGA, whereas complete abrogation is lethal. Genome-wide analysis of Rad21 binding reveals a change in distribution from pericentromeric satellite DNA and other locations, including the miR-430 locus (the products of which are responsible for maternal transcript degradation), to genes, as embryos progress through the MZT. After MZT, a subset of Rad21 binding overlaps the pioneer factor Pou5f3, which activates early expressed genes. Rad21 depletion disrupts the formation of nucleoli and RNA polymerase II foci, suggestive of global defects in chromosome architecture. We propose that Rad21/cohesin redistribution to active areas of the genome is key to the establishment of chromosome organization and the embryonic developmental program.

Lower insulin sensitivity remains a feature of children born very preterm.

BACKGROUND: The first report of children born very preterm (<32 weeks of gestation) having insulin resistance was made 16 years ago. However, neonatal care has improved since. Thus, we aimed to assess whether children born very preterm still have lower insulin sensitivity than term controls. METHODS: Participants were prepubertal children aged 5 to 11 years born very preterm (<32 weeks of gestation; n = 51; 61% boys) or at term (37-41 weeks; n = 50; 62% boys). Frequently sampled intravenous glucose tolerance tests were performed, and insulin sensitivity was calculated using Bergman's minimal model. Additional clinical assessments included anthropometry, body composition using whole-body dual-energy X-ray absorptiometry scans, clinic blood pressure, and 24-hour ambulatory blood pressure monitoring. RESULTS: Children born very preterm were 0.69 standard deviation score (SDS) lighter (P < .001), 0.53 SDS shorter (P = .003), and had body mass index 0.57 SDS lower (P = .003) than children born at term. Notably, children born very preterm had insulin sensitivity that was 25% lower than term controls (9.4 vs 12.6 × 10-4 minutes-1 ·[mU/L]; P = .001). Other parameters of glucose metabolism, including fasting insulin levels, were similar in the two groups. The awake systolic blood pressure (from 24-hour monitoring) tended to be 3.1 mm Hg higher on average in children born very preterm (P = .054), while the clinic systolic blood pressure was 5.4 mm Hg higher (P = .002). CONCLUSIONS: Lower insulin sensitivity remains a feature of children born very preterm, despite improvements in neonatal intensive care. As reported in our original study, our findings suggest the defect in insulin action in prepubertal children born very pretermis primarily peripheral and not hepatic.

Hi-C detects novel structural variants in HL-60 and HL-60/S4 cell lines.

Cancer cell lines often have large structural variants (SVs) that evolve over time. There are many reported differences in large scale SVs between HL-60 and HL-60/S4, two cell lines derived from the same acute myeloid leukemia sample. However, the stability and variability of inter- and intra-chromosomal structural variants between different sources of the same cell line is unknown. Here, we used Hi-C and RNA-seq to identify and compare large SVs in HL-60 and HL-60/S4 cell lines. Comparisons with previously published karyotypes identified novel SVs in both cell lines. Hi-C was used to characterize the known expansion centered on the MYC locus. The MYC expansion was integrated into known locations in HL-60/S4, and a novel location (chr4) in HL-60. The HL-60 cell line has more within-line structural variation than the HL-60/S4 derivative cell line. Collectively we demonstrate the usefulness of Hi-C and with RNA-seq data for the identification and characterization of SVs.

A non-coding genetic variant maximally associated with serum urate levels is functionally linked to HNF4A-dependent PDZK1 expression.

The precise molecular mechanisms by which urate-associated genetic variants affect urate levels are unknown. Here, we tested for functional linkage of the maximally associated genetic variant rs1967017 at the PDZK1 locus to elevated PDZK1 expression. We performed expression quantitative trait loci (eQTL) and likelihood analyses and gene expression assays. Zebrafish were used to evaluate tissue-specific gene expression. Luciferase assays in HEK293 and HepG2 cells measured the effect of rs1967017 on transcription amplitude. Probabilistic Annotation Integrator analysis revealed rs1967017 as most likely to be causal and rs1967017 was an eQTL for PDZK1 in the intestine. The region harboring rs1967017 was capable of directly driving green fluorescent protein expression in the kidney, liver and intestine of zebrafish embryos, consistent with a conserved ability to confer tissue-specific expression. Small interfering RNA depletion of HNF4A reduced endogenous PDZK1 expression in HepG2 cells. Luciferase assays showed that the T allele of rs1967017 gains enhancer activity relative to the urate-decreasing C allele, with T allele enhancer activity abrogated by HNF4A depletion. HNF4A physically binds the rs1967017 region, suggesting direct transcriptional regulation of PDZK1 by HNF4A. Computational prediction of increased motif strength, together with our functional assays, suggests that the urate-increasing T allele of rs1967017 strengthens a binding site for the transcription factor HNF4A. Our and other data predict that the urate-raising T allele of rs1967017 enhances HNF4A binding to the PDZK1 promoter, thereby increasing PDZK1 expression. As PDZK1 is a scaffold protein for many ion channel transporters, increased expression can be predicted to increase activity of urate transporters and alter excretion of urate.

Gut microbiome transfer-Finding the perfect fit.

Gut microbiome transfer (GMT; also referred to as faecal microbiota transplantation or FMT) has been propelled from fringe therapy to mainstream science as a highly effective treatment for recurrent Clostridioides difficile infection. As a result, there has been great interest in the potential efficacy and safety of GMT in treating other medical conditions, for example inflammatory bowel disease, and more recently as a novel therapy for obesity and metabolic diseases. For these chronic conditions, the results from clinical trials have been mixed. Further, specifically in obesity and metabolic diseases, there are limited available data, with only a few published studies with a small number of participants and short duration of follow-up. Therefore, this review aims to explore the human, microbial and formulation factors that may affect the success of GMT. This includes various aspects in the preparation and administration of GMT, such as stool processing, modes of delivery, pretreatment with antibiotics and/or bowel lavage, frequency of GMT and possible use of precision bacteriotherapy. In addition, we examine the potential use of GMT in obesity, type 2 diabetes and metabolic diseases based on current available literature, highlighting some recent advances in GMT research in this area, as well as potential adverse effects after GMT therapy.

Common Variants Coregulate Expression of GBA and Modifier Genes to Delay Parkinson's Disease Onset.

BACKGROUND: GBA mutations are numerically the most significant genetic risk factor for Parkinson's disease (PD), yet these mutations have low penetrance, suggesting additional mechanisms. OBJECTIVES: The objective of this study was to determine if the penetrance of GBA in PD can be explained by regulatory effects on GBA and modifier genes. METHODS: Genetic variants associated with the regulation of GBA were identified by screening 128 common single nucleotide polymorphisms (SNPs) in the GBA locus for spatial cis-expression quantitative trail locus (supported by chromatin interactions). RESULTS: We identified common noncoding SNPs within GBA that (1) regulate GBA expression in peripheral tissues, some of which display α-synuclein pathology and (2) coregulate potential modifier genes in the central nervous system and/or peripheral tissues. Haplotypes based on 3 of these SNPs delay disease onset by 5 years. In addition, SNPs on 6 separate chromosomes coregulate GBA expression specifically in either the substantia nigra or cortex, and their combined effect potentially modulates motor and cognitive symptoms, respectively. CONCLUSIONS: This work provides a new perspective on the haplotype-specific effects of GBA and the genetic etiology of PD, expanding the role of GBA from the gene encoding the ß-glucocerebrosidase (GCase) to that of a central regulator and modifier of PD onset, with GBA expression itself subject to distant regulation. Some idiopathic patients might possess insufficient GBA-encoded GCase activity in the substantia nigra as the result of distant regulatory variants and therefore might benefit from GBA-targeting therapeutics. The SNPs' regulatory impacts provide a plausible explanation for the variable phenotypes also observed in GBA-centric Gaucher's disease and dementia with Lewy bodies. © 2020 The Authors. Movement Disorders published by Wiley Periodicals, LLC on behalf of International Parkinson and Movement Disorder Society.

Robotic automation of a UHPLC/MS-MS method profiling one-carbon metabolites, amino acids, and precursors in plasma.

Amino acids (AAs) and one-carbon (1-C) metabolism compounds are involved in a range of key metabolic pathways, and mediate numerous health and disease processes in the human body. Previous assays have quantified a limited selection of these compounds and typically require extensive manual handling. Here, we describe the robotic automation of an analytical method for the simultaneous quantification of 37 1-C metabolites, amino acids, and precursors using reversed-phase ultra-high-pressure liquid chromatography coupled with tandem mass spectrometry (UHPLC/MS-MS). Compound extraction from human plasma was tested manually before being robotically automated. The final automated analytical panel was validated on human plasma samples. Our automated and multiplexed method holds promise for application to large cohort studies.