Interpreting metabolic complexity via isotope-assisted metabolic flux analysis.

Isotope-assisted metabolic flux analysis (iMFA) is a powerful method to mathematically determine the metabolic fluxome from experimental isotope labeling data and a metabolic network model. While iMFA was originally developed for industrial biotechnological applications, it is increasingly used to analyze eukaryotic cell metabolism in physiological and pathological states. In this review, we explain how iMFA estimates the intracellular fluxome, including data and network model (inputs), the optimization-based data fitting (process), and the flux map (output). We then describe how iMFA enables analysis of metabolic complexities and discovery of metabolic pathways. Our goal is to expand the use of iMFA in metabolism research, which is essential to maximizing the impact of metabolic experiments and continuing to advance iMFA and biocomputational techniques.

Induced pluripotent stem cell-derived cells model brain microvascular endothelial cell glucose metabolism.

Glucose transport from the blood into the brain is tightly regulated by brain microvascular endothelial cells (BMEC), which also use glucose as their primary energy source. To study how BMEC glucose transport contributes to cerebral glucose hypometabolism in diseases such as Alzheimer's disease, it is essential to understand how these cells metabolize glucose. Human primary BMEC (hpBMEC) can be used for BMEC metabolism studies; however, they have poor barrier function and may not recapitulate in vivo BMEC function. iPSC-derived BMEC-like cells (hiBMEC) are readily available and have good barrier function but may have an underlying epithelial signature. In this study, we examined differences between hpBMEC and hiBMEC glucose metabolism using a combination of dynamic metabolic measurements, metabolic mass spectrometry, RNA sequencing, and Western blots. hiBMEC had decreased glycolytic flux relative to hpBMEC, and the overall metabolomes and metabolic enzyme levels were different between the two cell types. However, hpBMEC and hiBMEC had similar glucose metabolism, including nearly identical glucose labeled fractions of glycolytic and TCA cycle metabolites. Treatment with astrocyte conditioned media and high glucose increased glycolysis in both hpBMEC and hiBMEC, though hpBMEC decreased glycolysis in response to fluvastatin while hiBMEC did not. Together, these results suggest that hiBMEC can be used to model cerebral vascular glucose metabolism, which expands their use beyond barrier models.

Isotope-Assisted Metabolic Flux Analysis: A Powerful Technique to Gain New Insights into the Human Metabolome in Health and Disease.

Cell metabolism represents the coordinated changes in genes, proteins, and metabolites that occur in health and disease. The metabolic fluxome, which includes both intracellular and extracellular metabolic reaction rates (fluxes), therefore provides a powerful, integrated description of cellular phenotype. However, intracellular fluxes cannot be directly measured. Instead, flux quantification requires sophisticated mathematical and computational analysis of data from isotope labeling experiments. In this review, we describe isotope-assisted metabolic flux analysis (iMFA), a rigorous computational approach to fluxome quantification that integrates metabolic network models and experimental data to generate quantitative metabolic flux maps. We highlight practical considerations for implementing iMFA in mammalian models, as well as iMFA applications in in vitro and in vivo studies of physiology and disease. Finally, we identify promising new frontiers in iMFA which may enable us to fully unlock the potential of iMFA in biomedical research.

Genetic Variation and Recurrent Haplotypes on Chromosome 6q23-25 Risk Locus in Familial Lung Cancer.

Although lung cancer is known to be caused by environmental factors, it has also been shown to have genetic components, and the genetic etiology of lung cancer remains understudied. We previously identified a lung cancer risk locus on 6q23-25 using microsatellite data in families with a history of lung cancer. To further elucidate that signal, we performed targeted sequencing on nine of our most strongly linked families. Two-point linkage analysis of the sequencing data revealed that the signal was heterogeneous and that different families likely had different risk variants. Three specific haplotypes were shared by some of the families: 6q25.3-26 in families 42 and 44, 6q25.2-25.3 in families 47 and 59, and 6q24.2-25.1 in families 30, 33, and 35. Region-based logarithm of the odds scores and expression data identified the likely candidate genes for each haplotype overlap: ARID1B at 6q25.3, MAP3K4 at 6q26, and UTRN (6q24.1) and PHACTR2 (6q24.2). Further annotation was used to zero in on potential risk variants in those genes. All four genes are good candidate genes for lung cancer risk, having been linked to either lung cancer specifically or other cancers. However, this is the first time any of these genes has been implicated in germline risk. Functional analysis of these four genes is planned for future work. SIGNIFICANCE: This study identifies four genes associated with lung cancer risk, which could help guide future lung cancer prevention and treatment approaches.

13C Metabolic Flux Analysis Indicates Endothelial Cells Attenuate Metabolic Perturbations by Modulating TCA Activity.

Disrupted endothelial metabolism is linked to endothelial dysfunction and cardiovascular disease. Targeted metabolic inhibitors are potential therapeutics; however, their systemic impact on endothelial metabolism remains unknown. In this study, we combined stable isotope labeling with 13C metabolic flux analysis (13C MFA) to determine how targeted inhibition of the polyol (fidarestat), pentose phosphate (DHEA), and hexosamine biosynthetic (azaserine) pathways alters endothelial metabolism. Glucose, glutamine, and a four-carbon input to the malate shuttle were important carbon sources in the baseline human umbilical vein endothelial cell (HUVEC) 13C MFA model. We observed two to three times higher glutamine uptake in fidarestat and azaserine-treated cells. Fidarestat and DHEA-treated HUVEC showed decreased 13C enrichment of glycolytic and TCA metabolites and amino acids. Azaserine-treated HUVEC primarily showed 13C enrichment differences in UDP-GlcNAc. 13C MFA estimated decreased pentose phosphate pathway flux and increased TCA activity with reversed malate shuttle direction in fidarestat and DHEA-treated HUVEC. In contrast, 13C MFA estimated increases in both pentose phosphate pathway and TCA activity in azaserine-treated cells. These data show the potential importance of endothelial malate shuttle activity and suggest that inhibiting glycolytic side branch pathways can change the metabolic network, highlighting the need to study systemic metabolic therapeutic effects.

Optimizing Stem Cell Therapy after Ischemic Brain Injury.

Stem cells have been used for regenerative and therapeutic purposes in a variety of diseases. In ischemic brain injury, preclinical studies have been promising, but have failed to translate results to clinical trials. We aimed to explore the application of stem cells after ischemic brain injury by focusing on topics such as delivery routes, regeneration efficacy, adverse effects, and in vivo potential optimization. PUBMED and Web of Science were searched for the latest studies examining stem cell therapy applications in ischemic brain injury, particularly after stroke or cardiac arrest, with a focus on studies addressing delivery optimization, stem cell type comparison, or translational aspects. Other studies providing further understanding or potential contributions to ischemic brain injury treatment were also included. Multiple stem cell types have been investigated in ischemic brain injury treatment, with a strong literature base in the treatment of stroke. Studies have suggested that stem cell administration after ischemic brain injury exerts paracrine effects via growth factor release, blood-brain barrier integrity protection, and allows for exosome release for ischemic injury mitigation. To date, limited studies have investigated these therapeutic mechanisms in the setting of cardiac arrest or therapeutic hypothermia. Several delivery modalities are available, each with limitations regarding invasiveness and safety outcomes. Intranasal delivery presents a potentially improved mechanism, and hypoxic conditioning offers a potential stem cell therapy optimization strategy for ischemic brain injury. The use of stem cells to treat ischemic brain injury in clinical trials is in its early phase; however, increasing preclinical evidence suggests that stem cells can contribute to the down-regulation of inflammatory phenotypes and regeneration following injury. The safety and the tolerability profile of stem cells have been confirmed, and their potent therapeutic effects make them powerful therapeutic agents for ischemic brain injury patients.

Whole Exome Sequencing of Highly Aggregated Lung Cancer Families Reveals Linked Loci for Increased Cancer Risk on Chromosomes 12q, 7p, and 4q.

BACKGROUND: Lung cancer kills more people than any other cancer in the United States. In addition to environmental factors, lung cancer has genetic risk factors as well, though the genetic etiology is still not well understood. We have performed whole exome sequencing on 262 individuals from 28 extended families with a family history of lung cancer. METHODS: Parametric genetic linkage analysis was performed on these samples using two distinct analyses-the lung cancer only (LCO) analysis, where only patients with lung cancer were coded as affected, and the all aggregated cancers (AAC) analysis, where other cancers seen in the pedigree were coded as affected. RESULTS: The AAC analysis yielded a genome-wide significant result at rs61943670 in POLR3B at 12q23.3. POLR3B has been implicated somatically in lung cancer, but this germline finding is novel and is a significant expression quantitative trait locus in lung tissue. Interesting genome-wide suggestive haplotypes were also found within individual families, particularly near SSPO at 7p36.1 in one family and a large linked haplotype spanning 4q21.3-28.3 in a different family. The 4q haplotype contains potential causal rare variants in DSPP at 4q22.1 and PTPN13 at 4q21.3. CONCLUSIONS: Regions on 12q, 7p, and 4q are linked to increased cancer risk in highly aggregated lung cancer families, 12q across families and 7p and 4q within a single family. POLR3B, SSPO, DSPP, and PTPN13 are currently the best candidate genes. IMPACT: Functional work on these genes is planned for future studies and if confirmed would lead to potential biomarkers for risk in cancer.

Small posterior fossa in Chiari I malformation affected families is significantly linked to 1q43-44 and 12q23-24.11 using whole exome sequencing.

The posterior fossa of the cranium contains the cerebellum and brainstem. Processes that reduce the volume of the posterior fossa squeeze the cerebellum and brainstem caudally, resulting in Chiari I malformation (CM1). CM1 causes neck pain, balance issues, decreased motor skills and headaches in those affected. We have posterior fossa measurements and whole exome sequence data on individuals from 7 extended families from Russia that have a family history of CM1. We performed parametric linkage analyses using an autosomal dominant inheritance model with a disease allele frequency of 0.01 and a penetrance of 0.8 for carriers and 0.0 for non-carriers. Variant-based two-point linkage analysis and gene-based linkage analysis was performed. Our results found a genome-wide significant signal on chromosome 1q43-44 (max HLOD = 3.3) in the variant-based analysis and 12q23 (max HLOD = 4.2) in the gene-based analysis. In both cases, the signal was driven by a single (different) family that contained a long, linked haplotype across the region in question. Using functional annotation, we were able to identify several rare nonsynonymous variants that were enriched in each family. The best candidate genes were rs765865412:G>A in MYBPC1 for the 12q haplotype and rs61749963:A>G in COX20 for the 1q haplotype. Good candidate variants in the 1q haplotype were also identified in CEP170 and AKT. Further laboratory work is planned to verify the causality of these genes.

Exome genotyping and linkage analysis identifies two novel linked regions and replicates two others for myopia in Ashkenazi Jewish families.

BACKGROUND: Myopia is one of most common eye diseases in the world and affects 1 in 4 Americans. It is a complex disease caused by both environmental and genetics effects; the genetics effects are still not well understood. In this study, we performed genetic linkage analyses on Ashkenazi Jewish families with a strong familial history of myopia to elucidate any potential causal genes. METHODS: Sixty-four extended Ashkenazi Jewish families were previously collected from New Jersey. Genotypes from the Illumina ExomePlus array were merged with prior microsatellite linkage data from these families. Additional custom markers were added for candidate regions reported in literature for myopia or refractive error. Myopia was defined as mean spherical equivalent (MSE) of -1D or worse and parametric two-point linkage analyses (using TwoPointLods) and multi-point linkage analyses (using SimWalk2) were performed as well as collapsed haplotype pattern (CHP) analysis in SEQLinkage and association analyses performed with FBAT and rv-TDT. RESULTS: Strongest evidence of linkage was on 1p36(two-point LOD = 4.47) a region previously linked to refractive error (MYP14) but not myopia. Another genome-wide significant locus was found on 8q24.22 with a maximum two-point LOD score of 3.75. CHP analysis also detected the signal on 1p36, localized to the LINC00339 gene with a maximum HLOD of 3.47, as well as genome-wide significant signals on 7q36.1 and 11p15, which overlaps with the MYP7 locus. CONCLUSIONS: We identified 2 novel linkage peaks for myopia on chromosomes 7 and 8 in these Ashkenazi Jewish families and replicated 2 more loci on chromosomes 1 and 11, one previously reported in refractive error but not myopia in these families and the other locus previously reported in the literature. Strong candidate genes have been identified within these linkage peaks in our families. Targeted sequencing in these regions will be necessary to definitively identify causal variants under these linkage peaks.

Myopia in Chinese families shows linkage to 10q26.13.

Purpose: To determine genetic linkage between myopia and Han Chinese patients with a family history of the disease. Methods: One hundred seventy-six Han Chinese patients from 34 extended families were given eye examinations, and mean spherical equivalent (MSE) in diopters (D) was calculated by adding the spherical component of the refraction to one-half the cylindrical component and taking the average of both eyes. The MSE was converted to a binary phenotype, where all patients with an MSE of -1.00 D or less were coded as affected. Unaffected individuals had an MSE greater than 0.00 D (ages 21 years and up), +1.50 (ages 11-20), or +2.00 D (ages 6-10 years). Individuals between the given upper threshold and -1.00 were coded as unknown. Patients were genotyped on an exome chip. Three types of linkage analyses were performed: single-variant two-point, multipoint, and collapsed haplotype pattern (CHP) variant two-point. Results: The CHP variant two-point results identified a significant peak (heterogeneity logarithm of the odds [HLOD] = 3.73) at 10q26.13 in TACC2. The single-variant two-point and multipoint analyses showed highly suggestive linkage to the same region. The single-variant two-point results identified 25 suggestive variants at HTRA1, also at 10q26.13. Conclusions: We report a significant genetic linkage between myopia and Han Chinese patients at 10q26.13. 10q26.13 contains several good candidate genes, such as TACC2 and the known age-related macular degeneration gene HTRA1. Targeted sequencing of the region is planned to identify the causal variant(s).

Caucasian Families Exhibit Significant Linkage of Myopia to Chromosome 11p.

Purpose: Myopia is a common visual disorder caused by eye overgrowth, resulting in blurry vision. It affects one in four Americans, and its prevalence is increasing. The genetic mechanisms that underpin myopia are not completely understood. Here, we use genotype data and linkage analyses to identify high-risk genetic loci that are significantly linked to myopia. Methods: Individuals from 56 Caucasian families with a history of myopia were genotyped on an exome-based array, and the single nucleotide polymorphism (SNP) data were merged with microsatellite genotype data. Refractive error measures on the samples were converted into binary phenotypes consisting of affected, unaffected, or unknown myopia status. Parametric linkage analyses assuming an autosomal dominant model with 90% penetrance and 10% phenocopy rate were performed. Results: Single variant two-point analyses yielded three significantly linked SNPs at 11p14.1 and 11p11.2; a further 45 SNPs at 11p were found to be suggestive. No other chromosome had any significant SNPs or more than seven suggestive linkages. Two of the significant SNPs were located in BBOX1-AS1 and one in the intergenic region between ORA47 and TRIM49B. Collapsed haplotype pattern two-point analysis and multipoint analyses also yielded multiple suggestively linked genes at 11p. Multipoint analysis also identified suggestive evidence of linkage on 20q13. Conclusions: We identified three genome-wide significant linked variants on 11p for myopia in Caucasians. Although the novel specific signals still need to be replicated, 11p is a promising region that has been identified by other linkage studies with a number of potentially interesting candidate genes. We hope that the identification of these regions on 11p as potential causal regions for myopia will lead to more focus on these regions and maybe possible replication of our specific linkage peaks in other studies. We further plan targeted sequencing on 11p for our most highly linked families to more clearly understand the source of the linkage in this region.