Investigation of candidate genes and mechanisms underlying obesity associated type 2 diabetes mellitus using bioinformatics analysis and screening of small drug molecules.

BACKGROUND: Obesity associated type 2 diabetes mellitus is a metabolic disorder ; however, the etiology of obesity associated type 2 diabetes mellitus remains largely unknown. There is an urgent need to further broaden the understanding of the molecular mechanism associated in obesity associated type 2 diabetes mellitus. METHODS: To screen the differentially expressed genes (DEGs) that might play essential roles in obesity associated type 2 diabetes mellitus, the publicly available expression profiling by high throughput sequencing data (GSE143319) was downloaded and screened for DEGs. Then, Gene Ontology (GO) and REACTOME pathway enrichment analysis were performed. The protein - protein interaction network, miRNA - target genes regulatory network and TF-target gene regulatory network were constructed and analyzed for identification of hub and target genes. The hub genes were validated by receiver operating characteristic (ROC) curve analysis and RT- PCR analysis. Finally, a molecular docking study was performed on over expressed proteins to predict the target small drug molecules. RESULTS: A total of 820 DEGs were identified between healthy obese and metabolically unhealthy obese, among 409 up regulated and 411 down regulated genes. The GO enrichment analysis results showed that these DEGs were significantly enriched in ion transmembrane transport, intrinsic component of plasma membrane, transferase activity, transferring phosphorus-containing groups, cell adhesion, integral component of plasma membrane and signaling receptor binding, whereas, the REACTOME pathway enrichment analysis results showed that these DEGs were significantly enriched in integration of energy metabolism and extracellular matrix organization. The hub genes CEBPD, TP73, ESR2, TAB1, MAP 3K5, FN1, UBD, RUNX1, PIK3R2 and TNF, which might play an essential role in obesity associated type 2 diabetes mellitus was further screened. CONCLUSIONS: The present study could deepen the understanding of the molecular mechanism of obesity associated type 2 diabetes mellitus, which could be useful in developing therapeutic targets for obesity associated type 2 diabetes mellitus.

Identification of key pathways and genes in polycystic ovary syndrome via integrated bioinformatics analysis and prediction of small therapeutic molecules.

To enhance understanding of polycystic ovary syndrome (PCOS) at the molecular level; this investigation intends to examine the genes and pathways associated with PCOS by using an integrated bioinformatics analysis. Based on the expression profiling by high throughput sequencing data GSE84958 derived from the Gene Expression Omnibus (GEO) database, the differentially expressed genes (DEGs) between PCOS samples and normal controls were identified. We performed a functional enrichment analysis. A protein-protein interaction (PPI) network, miRNA- target genes and TF - target gene networks, were constructed and visualized, with which the hub gene nodes were identified. Validation of hub genes was performed by using receiver operating characteristic (ROC) and RT-PCR. Small drug molecules were predicted by using molecular docking. A total of 739 DEGs were identified, of which 360 genes were up regulated and 379 genes were down regulated. GO enrichment analysis revealed that up regulated genes were mainly involved in peptide metabolic process, organelle envelope and RNA binding and the down regulated genes were significantly enriched in plasma membrane bounded cell projection organization, neuron projection and DNA-binding transcription factor activity, RNA polymerase II-specific. REACTOME pathway enrichment analysis revealed that the up regulated genes were mainly enriched in translation and respiratory electron transport and the down regulated genes were mainly enriched in generic transcription pathway and transmembrane transport of small molecules. The top 10 hub genes (SAA1, ADCY6, POLR2K, RPS15, RPS15A, CTNND1, ESR1, NEDD4L, KNTC1 and NGFR) were identified from PPI network, miRNA - target gene network and TF - target gene network. The modules analysis showed that genes in modules were mainly associated with the transport of respiratory electrons and signaling NGF, respectively. We find a series of crucial genes along with the pathways that were most closely related with PCOS initiation and advancement. Our investigations provide a more detailed molecular mechanism for the progression of PCOS, detail information on the potential biomarkers and therapeutic targets.

In vivo functional screening for systems-level integrative cancer genomics.

With the genetic portraits of all major human malignancies now available, we next face the challenge of characterizing the function of mutated genes, their downstream targets, interactions and molecular networks. Moreover, poorly understood at the functional level are also non-mutated but dysregulated genomes, epigenomes or transcriptomes. Breakthroughs in manipulative mouse genetics offer new opportunities to probe the interplay of molecules, cells and systemic signals underlying disease pathogenesis in higher organisms. Herein, we review functional screening strategies in mice using genetic perturbation and chemical mutagenesis. We outline the spectrum of genetic tools that exist, such as transposons, CRISPR and RNAi and describe discoveries emerging from their use. Genome-wide or targeted screens are being used to uncover genomic and regulatory landscapes in oncogenesis, metastasis or drug resistance. Versatile screening systems support experimentation in diverse genetic and spatio-temporal settings to integrate molecular, cellular or environmental context-dependencies. We also review the combination of in vivo screening and barcoding strategies to study genetic interactions and quantitative cancer dynamics during tumour evolution. These scalable functional genomics approaches are transforming our ability to interrogate complex biological systems.

Identification and analysis of 35 genes associated with vitamin D deficiency: A systematic review to identify genetic variants.

Vitamin D deficiency is a public health concern associated with, but not limited to, skeletal anomalies, chronic diseases, immune conditions, and cancer, among others. Hypovitaminosis D is mainly associated with environmental and lifestyle factors that affect sunlight exposure. However, genetic factors also influence 25-hydroxyvitamin D (25[OH]D) serum concentration. Although there is available information of genes with clear biological relevance or markers identified by Genome-Wide Association Studies, an overall view and screening tool to identify known genetic causes of altered serum levels of 25(OH)D is lacking. Moreover, there are no studies including the total genetic evidence associated with abnormal serum concentration of 25(OH)D. Therefore, we conducted a de-novo systematic literature review to propose a set of genes comprehensive of all genetic variants reported to be associated with deficiency of vitamin D. Abstracts retrieved from PubMed search were organized by gene and curated one-by-one using the PubTerm web tool. The genes identified were classified according to the type of genetic evidence associated with serum 25(OH)D levels and were also compared with the few commonly screened genes related to vitamin D status. This strategy allowed the identification of 35 genes associated with serum 25(OH)D concentrations, 27 (75%) of which are not commercially available and are not, therefore, analyzed in clinical practice for genetic counseling, nor are they sufficiently studied for research purposes. Functional analysis of the genes identified confirmed their role in vitamin D pathways and diseases. Thus, the list of genes is an important source to understand the genetic determinants of 25(OH)D levels. To further support our findings, we provide a map of the reported functional variants and SNPs not included in ClinVar, minor allelic frequencies, SNP effect sizes, associated diseases, and an integrated overview of the biological role of the genes. In conclusion, we identified a comprehensive candidate list of genes associated with serum 25(OH)D concentrations, most of which are not commercially available, but would prove of importance in clinical practice in screening for patients that should respond to supplementation because of alterations in absorption, patients that would have little benefit because alterations in the downstream metabolism of vitamin D, and to study non-responsiveness to supplementation with vitamin D.

Molecular Dissection of Quantitative Variation in Bermudagrass Hybrids (Cynodon dactylon x transvaalensis): Morphological Traits.

Bermudagrass (Cynodon (L.)) is the most important warm-season grass grown for forage or turf. It shows extensive variation in morphological characteristics and growth attributes, but the genetic basis of this variation is little understood. Detection and tagging of quantitative trait loci (QTL) affecting above-ground morphology with diagnostic DNA markers would provide a foundation for genetic and molecular breeding applications in bermudagrass. Here, we report early findings regarding genetic architecture of foliage (canopy height, HT), stolon (stolon internode length, ILEN and length of the longest stolon LLS), and leaf traits (leaf blade length, LLEN and leaf blade width, LW) in 110 F1 individuals derived from a cross between Cynodon dactylon (T89) and C. transvaalensis (T574). Separate and joint environment analyses were performed on trait data collected across two to five environments (locations, and/or years, or time), finding significant differences (P < 0.001) among the hybrid progeny for all traits. Analysis of marker-trait associations detected 74 QTL and 135 epistatic interactions. Composite interval mapping (CIM) and mixed-model CIM (MCIM) identified 32 main effect QTL (M-QTL) and 13 interacting QTL (int-QTL). Colocalization of QTL for plant morphology partially explained significant correlations among traits. M-QTL qILEN-3-2 (for ILEN; R2 = 11-19%), qLLS-7-1 (for LLS; R2 = 13-27%), qLEN-1-1 (for LLEN; R2 = 10-11%), and qLW-3-2 (for LW; R2 = 10-12%) were 'stable' across multiple environments, representing candidates for fine mapping and applied breeding applications. QTL correspondence between bermudagrass and divergent grass lineages suggests opportunities to accelerate progress by predictive breeding of bermudagrass.

Towards frailty biomarkers: Candidates from genes and pathways regulated in aging and age-related diseases.

OBJECTIVE: Use of the frailty index to measure an accumulation of deficits has been proven a valuable method for identifying elderly people at risk for increased vulnerability, disease, injury, and mortality. However, complementary molecular frailty biomarkers or ideally biomarker panels have not yet been identified. We conducted a systematic search to identify biomarker candidates for a frailty biomarker panel. METHODS: Gene expression databases were searched (http://genomics.senescence.info/genes including GenAge, AnAge, LongevityMap, CellAge, DrugAge, Digital Aging Atlas) to identify genes regulated in aging, longevity, and age-related diseases with a focus on secreted factors or molecules detectable in body fluids as potential frailty biomarkers. Factors broadly expressed, related to several "hallmark of aging" pathways as well as used or predicted as biomarkers in other disease settings, particularly age-related pathologies, were identified. This set of biomarkers was further expanded according to the expertise and experience of the authors. In the next step, biomarkers were assigned to six "hallmark of aging" pathways, namely (1) inflammation, (2) mitochondria and apoptosis, (3) calcium homeostasis, (4) fibrosis, (5) NMJ (neuromuscular junction) and neurons, (6) cytoskeleton and hormones, or (7) other principles and an extensive literature search was performed for each candidate to explore their potential and priority as frailty biomarkers. RESULTS: A total of 44 markers were evaluated in the seven categories listed above, and 19 were awarded a high priority score, 22 identified as medium priority and three were low priority. In each category high and medium priority markers were identified. CONCLUSION: Biomarker panels for frailty would be of high value and better than single markers. Based on our search we would propose a core panel of frailty biomarkers consisting of (1) CXCL10 (C-X-C motif chemokine ligand 10), IL-6 (interleukin 6), CX3CL1 (C-X3-C motif chemokine ligand 1), (2) GDF15 (growth differentiation factor 15), FNDC5 (fibronectin type III domain containing 5), vimentin (VIM), (3) regucalcin (RGN/SMP30), calreticulin, (4) PLAU (plasminogen activator, urokinase), AGT (angiotensinogen), (5) BDNF (brain derived neurotrophic factor), progranulin (PGRN), (6) α-klotho (KL), FGF23 (fibroblast growth factor 23), FGF21, leptin (LEP), (7) miRNA (micro Ribonucleic acid) panel (to be further defined), AHCY (adenosylhomocysteinase) and KRT18 (keratin 18). An expanded panel would also include (1) pentraxin (PTX3), sVCAM/ICAM (soluble vascular cell adhesion molecule 1/Intercellular adhesion molecule 1), defensin α, (2) APP (amyloid beta precursor protein), LDH (lactate dehydrogenase), (3) S100B (S100 calcium binding protein B), (4) TGFß (transforming growth factor beta), PAI-1 (plasminogen activator inhibitor 1), TGM2 (transglutaminase 2), (5) sRAGE (soluble receptor for advanced glycosylation end products), HMGB1 (high mobility group box 1), C3/C1Q (complement factor 3/1Q), ST2 (Interleukin 1 receptor like 1), agrin (AGRN), (6) IGF-1 (insulin-like growth factor 1), resistin (RETN), adiponectin (ADIPOQ), ghrelin (GHRL), growth hormone (GH), (7) microparticle panel (to be further defined), GpnmB (glycoprotein nonmetastatic melanoma protein B) and lactoferrin (LTF). We believe that these predicted panels need to be experimentally explored in animal models and frail cohorts in order to ascertain their diagnostic, prognostic and therapeutic potential.

SLC39A8 deficiency: biochemical correction and major clinical improvement by manganese therapy.

PurposeSLC39A8 deficiency is a severe inborn error of metabolism that is caused by impaired function of manganese metabolism in humans. Mutations in SLC39A8 lead to impaired function of the manganese transporter ZIP8 and thus manganese deficiency. Due to the important role of Mn2+ as a cofactor for a variety of enzymes, the resulting phenotype is complex and severe. The manganese-dependence of ß-1,4-galactosyltransferases leads to secondary hypoglycosylation, making SLC39A8 deficiency both a disorder of trace element metabolism and a congenital disorder of glycosylation. Some hypoglycosylation disorders have previously been treated with galactose administration. The development of an effective treatment of the disorder by high-dose manganese substitution aims at correcting biochemical, and hopefully, clinical abnormalities.MethodsTwo SCL39A8 deficient patients were treated with 15 and 20 mg MnSO4/kg bodyweight per day. Glycosylation and blood manganese were monitored closely. In addition, magnetic resonance imaging was performed to detect potential toxic effects of manganese.ResultsAll measured enzyme dysfunctions resolved completely and considerable clinical improvement regarding motor abilities, hearing, and other neurological manifestations was observed.ConclusionHigh-dose manganese substitution was effective in two patients with SLC39A8 deficiency. Close therapy monitoring by glycosylation assays and blood manganese measurements is necessary to prevent manganese toxicity.

Discovery of targetable genetic alterations in advanced non-small cell lung cancer using a next-generation sequencing-based circulating tumor DNA assay.

Next-generation sequencing (NGS)-based circulating tumor DNA (ctDNA) assays have provided a new method of identifying tumor-driving genes in patients with advanced non-small cell lung carcinoma (NSCLC), especially in those whose cancer tissues are unavailable or in those that have acquired treatment resistance. Here, we describe a total of 119 patients with advanced EGFR-TKI-naive NSCLC and 15 EGFR-TKI-resistant patients to identify somatic SNVs, small indels, CNVs and gene fusions in 508 tumor-related genes. Somatic ctDNA mutations were detected in 82.8% (111/134) of patients in the total cohort. Of the 119 patients with advanced NSCLC, 27.7% (33/119) were suitable for treatment with National Comprehensive Cancer Network (NCCN) guideline-approved targeted drugs. Actionable genetic alterations included 25 EGFR mutations, 5 BRAF mutations, and 1 MET mutation, as well as 1 EML4-ALK gene fusion and 1 KIF5B-RET gene fusion. In 19.3% (23/119) of the patients, we also identified genomic alterations with that could be targeted by agents that are in clinical trials, such as mTOR inhibitors, PARP inhibitors, and CDK4/6 inhibitors. Additionally, the EGFR T790M mutation was found in 46.7% (7/15) of the patients with EGFR-TKI-resistant NSCLC, suggesting that the NGS-based ctDNA assay might be an optional method to monitor EGFR-TKI resistance and to discover mechanisms of drug resistance.

Adapting simultaneous analysis phylogenomic techniques to study complex disease gene relationships.

The characterization of complex diseases remains a great challenge for biomedical researchers due to the myriad interactions of genetic and environmental factors. Network medicine approaches strive to accommodate these factors holistically. Phylogenomic techniques that can leverage available genomic data may provide an evolutionary perspective that may elucidate knowledge for gene networks of complex diseases and provide another source of information for network medicine approaches. Here, an automated method is presented that leverages publicly available genomic data and phylogenomic techniques, resulting in a gene network. The potential of approach is demonstrated based on a case study of nine genes associated with Alzheimer Disease, a complex neurodegenerative syndrome. The developed technique, which is incorporated into an update to a previously described Perl script called "ASAP," was implemented through a suite of Ruby scripts entitled "ASAP2," first compiles a list of sequence-similarity based orthologues using PSI-BLAST and a recursive NCBI BLAST+ search strategy, then constructs maximum parsimony phylogenetic trees for each set of nucleotide and protein sequences, and calculates phylogenetic metrics (Incongruence Length Difference between orthologue sets, partitioned Bremer support values, combined branch scores, and Robinson-Foulds distance) to provide an empirical assessment of evolutionary conservation within a given genetic network. In addition to the individual phylogenetic metrics, ASAP2 provides results in a way that can be used to generate a gene network that represents evolutionary similarity based on topological similarity (the Robinson-Foulds distance). The results of this study demonstrate the potential for using phylogenomic approaches that enable the study of multiple genes simultaneously to provide insights about potential gene relationships that can be studied within a network medicine framework that may not have been apparent using traditional, single-gene methods. Furthermore, the results provide an initial integrated evolutionary history of an Alzheimer Disease gene network and identify potentially important co-evolutionary clustering that may warrant further investigation.

Genetic variants in renalase and blood pressure responses to dietary salt and potassium interventions: a family-based association study.

BACKGROUND/AIMS: Renalase (gene name RNLS), a recently discovered enzyme with monoamine oxidase activity, is implicated in the degradation of catecholamines. Recent studies indicate that common variations in the gene with RNLS are associated with hypertension. The aim of this study was to examine the association between genetic variants in RNLS and blood pressure (BP) responses to strict dietary interventions of salt and potassium intake. METHODS: A total of 334 subjects from 124 families were selected and sequentially maintained on a low-salt diet for 7 days (3.0 g/day, NaCl), then a high-salt diet for 7 days (18.0 g/day, NaCl), high-salt diet with potassium supplementation for another 7 days (4.5 g/day, KCl). RESULTS: SNPs rs919115 and rs792205 of the RNLS gene were significantly associated with diastolic BP (DBP) and mean arterial pressure (MAP) responses to high-salt intervention. In addition, rs12356177 was significantly associated with systolic BP (SBP) and DBP responses to low-salt diet, and SBP, DBP or MAP during the high-salt intervention. Unfortunately, no associations for the 7 RNLS SNPs with BP response to high-salt diet with potassium supplementation reached nominal statistical significance. CONCLUSIONS: This family-based study indicates that genetic variants in the RNLS gene are significantly associated with BP responses to dietary salt intake.

Association study of gene polymorphisms and depression with abnormal humor in traditional Uighur medicine.

BACKGROUND: According to the humor theory of Traditional Uighur Medicine (TUM), a same disease is classified into different abnormal humor types and corresponding methods are applied to treat the diseases according to the type of abnormal humor characteristics. To date the biological foundation of classification of diseases by humor theory has been little studied and the mechanism of action is still unclear. In the present study, we aimed to investigate the association between some related gene polymorphisms and depression with abnormal humor in TUM. METHODS: 201 cases of depression patients in a Uighur population were divided into two groups as: 107 cases of depression patients with abnormal black bile (ABB), 94 cases of depression patients with none abnormal black bile (nABB), and 50 healthy people were served as control group. Venous blood was used to isolate DNA samples, and the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) technique was used for genotyping of single nucleotide polymorphisms (SNPs). Polymorphisms in the serotonin 2A (5-HT2A) receptor gene, brain derived neurotrophic factor (BDNF), serotonin 1A (5-HT1A) receptor gene were investigated in each groups, respectively. RESULTS: The 5-HT2A A-1438G, 5-HT2A T102C, BDNF Val66Met, and 5-HT1A C-1019G gene polymorphisms showed significant association with ABB. However, no difference between nABB and controls was found for those genotype distribution and allele frequency. Moreover, the T102C and A1438G SNPs in the 5-HT2A receptor gene polymorphisms were in linkage disequilibrium. In addition, the OR associated with the combination of Val66Met-Val/Val genotype plus the presence of -1019C allele was 8.393 for ABB compared with controls (OR 8.393; 95% CI 1.807 ~ 38.991; P = 0.003). Moreover, the OR associated with the presence of -Met plus -1019C alleles was 12.194 for ABB compared with controls (OR 12.194; 95% CI 1.433 ~ 103.776; P = 0.005). The OR associated with the presence of -1438C/C plus Val/Val genotypes was 7.738 for ABB compared with controls (OR 7.738; 95% CI 1.566 ~ 38.241; P = 0.005). CONCLUSION: It was concluded that there were significant relationship between the gene polymorphisms and classification of depression with abnormal humor in TUM. The 5-HT2A A-1438G, 5-HT2A T102C, BDNF Val66Met, and 5-HT1A C-1019G gene polymorphisms might predict the incidence of depression with ABB.

From integrative genomics to therapeutic targets.

Combinatorial approaches that integrate conventional pathology with genomic profiling and functional genomics have begun to enhance our understanding of the genetic basis of breast cancer. These methods have identified key genotypic-phenotypic correlations in different breast cancer subtypes that have led to the discovery of genetic dependencies that drive their behavior. Moreover, this knowledge has been applied to define novel tailored therapies for these groups of patients with cancer. With the current emphasis on characterizing the mutational repertoire of breast cancers by next-generation sequencing, the question remains as to what constitutes a driver event. By focusing efforts on homogenous subgroups of breast cancer and integrating orthogonal data-types combined with functional approaches, we can begin to unravel the heterogeneity and identify aberrations that can be therapeutically targeted.

Investigating highly replicated asthma genes as candidate genes for allergic rhinitis.

BACKGROUND: Asthma genetics has been extensively studied and many genes have been associated with the development or severity of this disease. In contrast, the genetic basis of allergic rhinitis (AR) has not been evaluated as extensively. It is well known that asthma is closely related with AR since a large proportion of individuals with asthma also present symptoms of AR, and patients with AR have a 5-6 fold increased risk of developing asthma. Thus, the relevance of asthma candidate genes as predisposing factors for AR is worth investigating. The present study was designed to investigate if SNPs in highly replicated asthma genes are associated with the occurrence of AR. METHODS: A total of 192 SNPs from 21 asthma candidate genes reported to be associated with asthma in 6 or more unrelated studies were genotyped in a Swedish population with 246 AR patients and 431 controls. Genotypes for 429 SNPs from the same set of genes were also extracted from a Singapore Chinese genome-wide dataset which consisted of 456 AR cases and 486 controls. All SNPs were subsequently analyzed for association with AR and their influence on allergic sensitization to common allergens. RESULTS: A limited number of potential associations were observed and the overall pattern of P-values corresponds well to the expectations in the absence of an effect. However, in the tests of allele effects in the Chinese population the number of significant P-values exceeds the expectations. The strongest signals were found for SNPs in NPSR1 and CTLA4. In these genes, a total of nine SNPs showed P-values <0.001 with corresponding Q-values <0.05. In the NPSR1 gene some P-values were lower than the Bonferroni correction level. Reanalysis after elimination of all patients with asthmatic symptoms excluded asthma as a confounding factor in our results. Weaker indications were found for IL13 and GSTP1 with respect to sensitization to birch pollen in the Swedish population. CONCLUSIONS: Genetic variation in the majority of the highly replicated asthma genes were not associated to AR in our populations which suggest that asthma and AR could have less in common than previously anticipated. However, NPSR1 and CTLA4 can be genetic links between AR and asthma and associations of polymorphisms in NPSR1 with AR have not been reported previously.

Ambulatory monitoring in the genetics of psychosomatic medicine.

Psychosomatic disorders are composed of an array of psychological, biologic, and environmental features. The existing evidence points to a role for genetic factors in explaining individual differences in the development and maintenance of a variety of disorders, but studies to date have not shown consistent and replicable effects. As such, the attempt to uncover individual differences in the expression of psychosomatic disorders as a function of genetic architecture requires careful attention to their phenotypic architecture or the various intermediate phenotypes that make up a heterogeneous disorder. Ambulatory monitoring offers a novel approach to measuring time-variant and situation-dependent intermediate phenotypes. Recent examples of the use of ambulatory monitoring in genetic studies of stress reactivity, chronic pain, alcohol use disorders, and psychosocial resilience are reviewed in an effort to highlight the benefits of ambulatory monitoring for genetic study designs.

The genetics of anorexia nervosa.

Anorexia nervosa (AN) is a disease defined by inappropriate weight loss and maintenance of body weight <85% of that expected for weight and height; it is most common in adolescent women aged 15-19 years. Numerous studies have highlighted the familial aggregation of the disease, suggesting a significant genetic component to its etiology. The purpose of this review is to discuss the different fields of genetic research--both in humans and animals--that have contributed to the understanding of this complex disorder. Candidate gene studies focusing on genes involved in the hypothalamic control of appetite and energy regulation have found genetic risk variants that increase risk for AN. A recent genome-wide association study has highlighted novel loci for further investigation in AN. Animal models and epigenetic studies are also considered; the most recent advances in each field and their contributions to the understanding of AN are emphasized.

Genome-wide association mapping reveals epistasis and genetic interaction networks in sugar beet.

Epistasis is defined as interactions between alleles of two or more genetic loci. Detection of epistatic interactions is the key to understand the genetic architecture and gene networks underlying complex traits. Here, we examined the extent of epistasis for seven quantitative traits with an association mapping approach in a large population of elite sugar beet lines. We found that correction for population stratification is required and that in terms of reducing the false-positive rate the mixed model approach including the kinship matrix performed best. In genome-wide scans, we detected both main effects and epistatic QTL. For physiological traits, the detected digenic and higher-order epistasis explained a considerable proportion of the genotypic variance. We illustrate that the identified epistatic interactions define comprehensive genetic networks, which may serve as starting points towards a systems-oriented approach to understand the regulation of complex traits.