Search

Genetic interactions (GIs), such as the synthetic lethal interaction, are promising therapeutic targets in precision medicine. However, despite extensive efforts to characterize GIs by large-scale perturbation screening, considerable false positives have been reported in multiple studies. We propose a new computational approach for improved precision in GI identification by applying constraints that consider actual biological phenomena. In this study, GIs were characterized by assessing mutation, loss of function, and expression profiles in the DEPMAP database. The expression profiles were used to exclude loss-of-function data for nonexpressed genes in GI characterization. More importantly, the characterized GIs were refined based on Kyoto Encyclopedia of Genes and Genomes (KEGG) or protein-protein interaction (PPI) networks, under the assumption that genes genetically interacting with a certain mutated gene are adjacent in the networks. As a result, the initial GIs characterized with CRISPR and RNAi screenings were refined to 65 and 23 GIs based on KEGG networks and to 183 and 142 GIs based on PPI networks. The evaluation of refined GIs showed improved precision with respect to known synthetic lethal interactions. The refining process also yielded a synthetic partner network (SPN) for each mutated gene, which provides insight into therapeutic strategies for the mutated genes; specifically, exploring the SPN of mutated BRAF revealed ELAVL1 as a potential target for treating BRAF -mutated cancer, as validated by previous research. We expect that this work will advance cancer therapeutic research.

This study investigated the potential genetic mechanisms which underlie adipose tissue mitochondrial dysfunction in Type 2 diabetes (T2D), by systematically identifying nuclear-encoded mitochondrial genes (NEMGs) among the genes regulated by T2D-associated genetic loci. The target genes of these 'disease loci' were identified by mapping genetic loci associated with both disease and gene expression levels (expression quantitative trait loci, eQTL) using high resolution genetic maps, with independent estimates co-locating to within a small genetic distance. These co-locating signals were defined as T2D-eQTL and the target genes as T2D cis- genes. In total, 763 cis- genes were associated with T2D-eQTL, of which 50 were NEMGs. Independent gene expression datasets for T2D and insulin resistant cases and controls confirmed that the cis- genes and cis- NEMGs were enriched for differential expression in cases, providing independent validation that genetic maps can identify informative functional genes. Two additional results were consistent with a potential role of T2D-eQTL in regulating the 50 identified cis- NEMGs in the context of T2D risk: (1) the 50 cis- NEMGs showed greater differential expression compared to other NEMGs and (2) other NEMGs showed a trend towards significantly decreased expression if their expression levels correlated more highly with the subset of 50 cis- NEMGs. These 50 cis -NEMGs, which are differentially expressed and associated with mapped T2D disease loci, encode proteins acting within key mitochondrial pathways, including some of current therapeutic interest such as the metabolism of branched-chain amino acids, GABA and biotin., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Maude, Lau, Maniatis and Andrew.)

The treatment of complex diseases often relies on combinatorial therapy, a strategy where drugs are used to target multiple genes simultaneously. Promising candidate genes for combinatorial perturbation often constitute epistatic genes, i.e., genes which contribute to a phenotype in a non-linear fashion. Experimental identification of the full landscape of genetic interactions by perturbing all gene combinations is prohibitive due to the exponential growth of testable hypotheses. Here we present a model for the inference of pairwise epistatic, including synthetic lethal, gene interactions from siRNA-based perturbation screens. The model exploits the combinatorial nature of siRNA-based screens resulting from the high numbers of sequence-dependent off-target effects, where each siRNA apart from its intended target knocks down hundreds of additional genes. We show that conditional and marginal epistasis can be estimated as interaction coefficients of regression models on perturbation data. We compare two methods, namely glinternet and xyz, for selecting non-zero effects in high dimensions as components of the model, and make recommendations for the appropriate use of each. For data simulated from real RNAi screening libraries, we show that glinternet successfully identifies epistatic gene pairs with high accuracy across a wide range of relevant parameters for the signal-to-noise ratio of observed phenotypes, the effect size of epistasis and the number of observations per double knockdown. xyz is also able to identify interactions from lower dimensional data sets (fewer genes), but is less accurate for many dimensions. Higher accuracy of glinternet, however, comes at the cost of longer running time compared to xyz. The general model is widely applicable and allows mining the wealth of publicly available RNAi screening data for the estimation of epistatic interactions between genes. As a proof of concept, we apply the model to search for interactions, and potential targets for treatment, among previously published sets of siRNA perturbation screens on various pathogens. The identified interactions include both known epistatic interactions as well as novel findings., Competing Interests: JJ is an employee of Novartis Institutes for BioMedical Research, which provided support in the form of the salary, but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Epistasis and cooperativity of folding both result from networks of energetic interactions in proteins. Epistasis results from energetic interactions among mutants, whereas cooperativity results from energetic interactions during folding that reduce the presence of intermediate states. The two concepts seem intuitively related, but it is unknown how they are related, particularly in terms of selection. To investigate their relationship, we simulated protein evolution under selection for cooperativity and separately under selection for epistasis. Strong selection for cooperativity created strong epistasis between contacts in the native structure but weakened epistasis between nonnative contacts. In contrast, selection for epistasis increased epistasis in both native and nonnative contacts and reduced cooperativity. Because epistasis can be used to predict protein structure only if it preferentially occurs in native contacts, this result indicates that selection for cooperativity may be key for predicting structure using epistasis. To evaluate this inference, we simulated the evolution of guanine nucleotide-binding protein (GB1) with and without cooperativity. With cooperativity, strong epistatic interactions clearly map out the native GB1 structure, while allowing the presence of intermediate states (low cooperativity) obscured the structure. This indicates that using epistasis measurements to reconstruct protein structure may be inappropriate for proteins with stable intermediates., Competing Interests: The authors declare no competing interest.

PIWI (P element induced wimpy testis) integrating RNAs (piRNAs) are small non-coding RNAs with the length of approximately 30 nucleotides that plays crucial roles in germ cells and adult stem cells. Recently, accumulating data have shown that piRNA and PIWI proteins are involved in tumorigenesis. However, the roles of PIWI proteins and piRNAs in pancreatic cancer are still elusive. Here, we showed that piR-017061 is significantly downregulated in pancreatic cancer patients' samples and pancreatic cancer cell lines. Furthermore, we studied the function of piR-017061 in pancreatic cancer and our data revealed that piR-017061 inhibits pancreatic cancer cell growth in vitro and in vivo. Moreover, we analyzed the genomic loci around piR-017061 and identified EFNA5 as a novel target of piR-017061. Importantly, our data further revealed a direct binding between piR-017061 and EFNA5 mRNA mediated by PIWIL1. Mechanically, piR-017061 cooperates with PIWIL1 to facilitate EFNA5 mRNA degradation and loss of piR-017061 results in accumulation of EFNA5 which facilitates pancreatic cancer development. Hence, our data provided novel insights into PIWI/piRNA-mediated gene regulation and their function in pancreatic cancer. Since PIWI proteins and piRNA predominately express in germline and cancer cells, our study provided novel therapeutic strategy for pancreatic cancer treatment.

Since domestication, rice has cultivated in a wide range of latitudes with different day lengths. Selection of diverse natural variations in heading date and photoperiod sensitivity is critical for adaptation of rice to different geographical environments. To unravel the genetic architecture underlying natural variation of rice flowering time, we conducted a genome wide association study (GWAS) using several association analysis strategies with a diverse worldwide collection of 529 O. sativa accessions. Heading date was investigated in three environments under long-day or short-day conditions, and photosensitivity was evaluated. By dividing the whole association panel into subpopulations and performing GWAS with both linear mixed models and multi-locus mixed-models, we revealed hundreds of significant loci harboring novel candidate genes as well as most of the known flowering time genes. In total, 127 hotspots were detected in at least two GWAS. Universal genetic heterogeneity was found across subpopulations. We further detected abundant interactions between GWAS loci, especially in indica. Functional gene families were revealed from enrichment analysis of the 127 hotspots. The results demonstrated a rich of genetic interactions in rice flowering time genes and such epistatic interactions contributed to the large portions of missing heritability in GWAS. It suggests the increased complexity of genetic heterogeneity might discount the power of increasing the sample sizes in GWAS., (Copyright © 2020 Elsevier B.V. All rights reserved.)

While opioids constitute the major component of perioperative analgesic regimens for surgery in general, a variety of evidence points to an association between perioperative opioid exposure and longer term oncologic outcomes. The mechanistic details underlying these effects are not well understood. In this study, we focused on clear cell renal cell carcinoma (ccRCC) and utilized RNA sequencing and outcome data from both The Cancer Genome Atlas, as well as a local patient cohort to identify survival-associated gene coexpression networks. We then projected drug-induced transcriptional profiles from in vitro cancer cells to predict drug effects on these networks and recurrence-free, cancer-specific, and overall survival. The opioid receptor agonist, leu-enkephalin, was predicted to have antisurvival effects in ccRCC, primarily through Th2 immune- and NRF2 -dependent macrophage networks. Conversely, the antagonist, naloxone, was predicted to have prosurvival effects, primarily through angiogenesis, fatty acid metabolism, and hemopoesis pathways. Eight coexpression networks associated with survival endpoints in ccRCC were identified, and master regulators of the transition from the normal to disease state were inferred, a number of which are linked to opioid pathways. These results are the first to suggest a mechanism for opioid effects on cancer outcomes through modulation of survival-associated coexpression networks. While we focus on ccRCC, this methodology may be employed to predict opioid effects on other cancer types and to personalize analgesic regimens in patients with cancer for optimal outcomes. SIGNIFICANCE: This study suggests a possible molecular mechanism for opioid effects on cancer outcomes generally, with implications for personalization of analgesic regimens., (©2020 American Association for Cancer Research.)

The neurodevelopmental hypothesis of schizophrenia has been widely accepted. In light of our previous microarray data, two neurodevelopment-related genes were focused on inclduing the N-acylsphingosine amidohydrolase 1 gene (ASAH1) and the nerve growth factor gene (NGF). The evidence that ASAH1 and NGF are associated with schizophrenia is far from conclusive. Furthermore, their interactions in schizophrenia have not been investigated. Total 413 patients and 578 controls were included. Eleven single-nucleotide polymorphisms (SNPs) in ASAH1 and NGF were selected. A multifactor dimensionality reduction (MDR) was applied to investigate gene-gene interactions in schizophrenia, and the traditional odds ratio methods was applied to validate it. The effects of ASAH1, NGF and their interaction on the severity of the disease were analyzed by 3 × 3 covariance analysis of (ANCOVA). The biological interaction between ASAH1 and NGF was examined. KEGG was used to identify the related signaling pathways. After correction by Bonferroni, there were no differences in the genotypic, allelic, or haplotypic frequencies of 11 SNPs between patients and controls. However, the interaction of certain SNPs had effect on susceptibility to schizophrenia, including two high-risk and one low-risk genotypic combinations (OR = 1.49 [1.11-2.00]; OR = 1.45 [1.09-1.92], and OR = 0.64 [0.41-0.98]). ASAH1-rs7830490 and its interaction with NGF-rs4332358 were associated with the general psychopathological subscale score (F adjusted  = 3.94, p adjusted  = 0.01; F adjusted  = 2.36, p adjusted  = 0.03). We also found that ASAH1 and NGF interacted with CaMK2B involving in the neurotrophin signaling pathway. Our results suggest that the interaction of ASAH1 and NGF with CaMK2B involved in neurotrophin signaling pathway may contribute to schizophrenia susceptibility and psychopathology., (Copyright © 2020 Elsevier Inc. All rights reserved.)

The coordination between mitochondrial and nuclear genes is crucial to eukaryotic organisms. Predicting the nature of these epistatic interactions can be difficult because of the transmission asymmetry of the genes involved. While autosomes and X-linked genes are transmitted through both sexes, genes on the Y chromosome and in the mitochondrial genome are uniparentally transmitted through males and females, respectively. Here, we generate 36 otherwise isogenic Drosophila melanogaster strains differing only in the geographical origin of their mitochondrial genome and Y chromosome, to experimentally examine the effects of the uniparentally inherited parts of the genome, as well as their interaction, in males. We assay longevity and gene expression through RNA-sequencing. We detect an important role for both mitochondrial and Y-linked genes, as well as extensive mitochondrial-Y chromosome epistasis. In particular, genes involved in male reproduction appear to be especially sensitive to such interactions, and variation on the Y chromosome is associated with differences in longevity. Despite these interactions, we find no evidence that the mitochondrial genome and Y chromosome are co-adapted within a geographical region. Overall, our study demonstrates a key role for the uniparentally inherited parts of the genome for male biology, but also that mito-nuclear interactions are complex and not easily predicted from simple transmission asymmetries.

Background: This study aimed to investigate the roles of lncRNA SNHG10 (SNHG10) and miR-218 in osteosarcoma (OS)., Methods: Paired OS and non-tumor tissues were collected from 58 OS patients. The expression of SNHG10 and miR-218 in tissue samples were determined by RT-qPCR. The interaction between SNHG10 and miR-218 was evaluated by overexpression experiment. Methylation-specific PCR was performed to assess the methylation status of miR-218. Glucose uptake in OS cells was analyzed by glucose uptake assay. Cell proliferation was detected by cell proliferation assay., Results: SNHG10 was upregulated in OS, while miR-218 was downregulated in OS. The expression of SNHG10 and miR-218 were inversely correlated. In OS cells, high glucose induced the upregulation of SNHG10 and downregulation of miR-218. In OS cells, SNHG10 positively, and miR-218 negatively regulated glucose uptake. Overexpression of SNHG10 increased miR-218 gene methylation and decreased the expression of miR-218. In addition, overexpression of SNHG10 also suppressed the inhibitory effects of overexpression of miR-218 on cell proliferation., Conclusions: SNHG10 increases the methylation of miR-218 gene to promote glucose uptake and cell proliferation in OS.

The role of histone modifications in the pathogenesis of schizophrenia has been proposed previously. H3F3B is a member of the histone 3. NSD2 is a histone methyltransferase that mediates dimethylation of Histone 3 lysine 36 (H3K36me2). The aim of the current study was to explore the associations between SNPs within H3F3B gene (rs60700976, rs3214028) and NSD2 gene (rs13148597, rs75820801) and the susceptibility to schizophrenia in a Chinese population. A total of 810 patients and 490 healthy controls were recruited and genetic association analyses were performed. The H3F3B gene polymorphisms rs3214028 and rs60700976 were significantly associated with schizophrenia. Rs60700976 was also associated with psychotic symptoms in schizophrenia patients. Furthermore, we found the interaction between NSD2 gene and H3F3B gene was related to the susceptibility to schizophrenia. The corresponding best three-locus model was H3F3B (rs60700976) - NSD2 (rs75820801, rs13148597), and the high-risk genotype combination was rs13148597(CC)- rs60700976(GG)-rs75820801(TT) (OR = 1.388[1.091-1.766], P = .007). The low-risk genotype combination was rs13148597(CC)-rs60700976(GG)-rs75820801(CT) (OR = 0.57 [0.330-0.985], P = .042). Our findings provided the preliminary evidence that the histone modification related H3F3B and NSD2 genes may confer the susceptibility to schizophrenia in a Chinese population., Competing Interests: Declaration of Competing Interest The authors declare that they have no conflict of interest., (Copyright © 2020. Published by Elsevier Inc.)

Growth hormone (GH) is a peptide hormone predominantly produced by the anterior pituitary and is essential for normal growth and metabolism. The GH locus contains five evolutionarily related genes under the control of an upstream locus control region that coordinates tissue-specific expression of these genes. Compromised GH signalling and genetic variation in these genes has been implicated in various disorders including cancer. We hypothesised that regulatory regions within the GH locus coordinate expression of a gene network that extends the impact of the GH locus control region. We used the CoDeS3D algorithm to analyse 529 common single nucleotide polymorphisms (SNPs) across the GH locus. This algorithm identifies colocalised Hi-C and eQTL associations to determine which SNPs are associated with a change in gene expression at loci that physically interact within the nucleus. One hundred and eighty-one common SNPs were identified that interacted with 292 eGenes across 48 different tissues. One hundred and forty-five eGenes were regulated in trans. eGenes were found to be enriched in GH/GHR-related cellular signalling pathways including MAPK, PI3K-AKT-mTOR, ERBB and insulin signalling, suggesting that these pathways may be co-regulated with GH signalling. Enrichment was also observed in the Wnt and Hippo signalling pathways and in pathways associated with hepatocellular, colorectal, breast and non-small cell lung carcinoma. Thirty-three eQTL SNPs identified in our study were found to be of regulatory importance in a genome-wide Survey of Regulatory Elements reporter screen. Our data suggest that the GH locus functions as a complex regulatory region that coordinates expression of numerous genes in cis and trans, many of which may be involved in modulating GH function in normal and disease states.

Background: Recent studies of synapse form and function highlight the importance of the actin cytoskeleton in regulating multiple aspects of morphogenesis, neurotransmission, and neural plasticity. The conserved actin-associated protein Enabled (Ena) is known to regulate development of the Drosophila larval neuromuscular junction through a postsynaptic mechanism. However, the functions and regulation of Ena within the presynaptic terminal has not been determined., Methods: Here, we use a conditional genetic approach to address a presynaptic role for Ena on presynaptic morphology and ultrastructure, and also examine the pathway in which Ena functions through epistasis experiments., Results: We find that Ena is required to promote the morphogenesis of presynaptic boutons and branches, in contrast to its inhibitory role in muscle. Moreover, while postsynaptic Ena is regulated by microRNA-mediated mechanisms, presynaptic Ena relays the output of the highly conserved receptor protein tyrosine phosphatase Dlar and associated proteins including the heparan sulfate proteoglycan Syndecan, and the non-receptor Abelson tyrosine kinase to regulate addition of presynaptic varicosities. Interestingly, Ena also influences active zones, where it restricts active zone size, regulates the recruitment of synaptic vesicles, and controls the amplitude and frequency of spontaneous glutamate release., Conclusion: We thus show that Ena, under control of the Dlar pathway, is required for presynaptic terminal morphogenesis and bouton addition and that Ena has active zone and neurotransmission phenotypes. Notably, in contrast to Dlar, Ena appears to integrate multiple pathways that regulate synapse form and function.

Epithelial-mesenchymal transition (EMT) is a complex biological program between physiology and pathology. Here, amniotic epithelial cells (AEC) were used as in vitro model of transiently inducible EMT in order to evaluate the transcriptional insights underlying this process. Therefore, RNA-seq was used to identify the differentially expressed genes and enrichment analyses were carried out to assess the intracellular pathways involved. As a result, molecules exclusively expressed in AEC that experienced EMT (GSTA1-1 and GSTM3) or when this process is inhibited (KLHL14 and KCNE3) were identified. Lastly, the network theory was used to obtain a computational model able to recognize putative controller genes involved in the induction and in the prevention of EMT. The results suggested an opposite role of lysophosphatidic acid (LPA) synthesis and degradation enzymes in the regulation of EMT process. In conclusion, these molecules may represent novel EMT regulators and also targets for developing new therapeutic strategies.

Aims/hypothesis: Peptide hormones are first synthesised as larger, inactive precursors that are converted to their active forms by endopeptidase cleavage and post-translational modifications, such as amidation. Recent, large-scale genome-wide studies have suggested that two coding variants of the amidating enzyme, peptidylglycine α-amidating monooxygenase (PAM), are associated with impaired insulin secretion and increased type 2 diabetes risk. We aimed to elucidate the role of PAM in modulating beta cell peptide amidation, beta cell function and the development of diabetes., Methods: PAM transcript and protein levels were analysed in mouse islets following induction of endoplasmic reticulum (ER) or cytokine stress, and PAM expression patterns were examined in human islets. To study whether haploinsufficiency of PAM accelerates the development of diabetes, Pam +/- and Pam +/+ mice were fed a low-fat diet (LFD) or high-fat diet (HFD) and glucose homeostasis was assessed. Since aggregates of the PAM substrate human islet amyloid polypeptide (hIAPP) lead to islet inflammation and beta cell failure, we also investigated whether PAM haploinsufficiency accelerated hIAPP-induced diabetes and islet amyloid formation in Pam +/- and Pam +/+ mice with beta cell expression of hIAPP., Results: Immunostaining revealed high expression of PAM in alpha, beta and delta cells in human pancreatic islets. Pam mRNA and PAM protein expression were reduced in mouse islets following administration of an HFD, and in isolated islets following induction of ER stress with thapsigargin, or cytokine stress with IL-1β, IFN-γ and TFN-α. Despite Pam +/- only having 50% PAM expression and enzyme activity as compared with Pam +/+ mice, glucose tolerance and body mass composition were comparable in the two models. After 24 weeks of HFD, both Pam +/- and Pam +/+ mice had insulin resistance and impaired glucose tolerance, but no differences in glucose tolerance, insulin sensitivity or plasma insulin levels were observed in PAM haploinsufficient mice. Islet amyloid formation and beta cell function were also similar in Pam +/- and Pam +/+ mice with beta cell expression of hIAPP., Conclusions/interpretation: Haploinsufficiency of PAM in mice does not accelerate the development of diet-induced obesity or hIAPP transgene-induced diabetes.

Although cystic fibrosis is a monogenic disease, a considerable clinical phenotypic variability is observed in patients with the same CFTR mutations. Thanks to the development of new and powerful tools for carrying out genetic studies, several genes called "modifier genes" have been identified as being associated with the severity of the lung function disorder in cystic fibrosis patients. Among these genes, SLC6A14 may modulate the anti-infective response and epithelial integrity of the airways, thus providing a potential therapeutic target to improve the patient's lung function., (Copyright © 2020 SPLF. Published by Elsevier Masson SAS. All rights reserved.)

Background: Genetic risk factors play an important role in the pathogenesis of Alzheimer's disease (AD). However, the gene-gene interaction (epistasis) between specific allelic variants is only partially understood., Objective: In our study, we examined the presence of the ɛ4 allele of apolipoprotein E (APOE) and the presence of C677T and A1298C (rs1801133 and rs1801131) polymorphisms in the methylenetetrahydrofolate reductase (MTHFR) gene in patients with AD and controls. We also evaluated the epistatic interaction between MTHFR and the APOE variants., Methods: A total of 564 patients with AD and 534 cognitively unimpaired age-matched controls were involved in the study., Results: The presence of the ɛ4 allele of APOE increases the risk of developing AD in a dose-dependent manner (OR 32.7: homozygotes, 15.6: homozygotes + heterozygotes, 14.3: heterozygotes). The combination of genotypes also increases the risk of developing AD in a dose-dependent manner: OR 18.3 (APOE 4/X and 4/4 + CT rs1801133), OR 19.4 (APOE 4/X and 4/4 + CT rs1801133 + AC rs1801131), OR 22.4 (APOE 4/X and 4/4 + TT rs1801133), and OR 21.2 (APOE 4/X and 4/4 + CC rs1801131). Homozygotes for variant alleles of MTHFR as well as patients with AD had significantly higher levels of homocysteine than homozygotes for standard alleles or controls., Conclusion: Homozygotes for APOE4 and carriers of APOE4 with TT genotype of rs1801133 were found to be at the highest risk of developing AD. These findings suggest that the epistatic interaction of specific gene variants can have a significant effect on the development of AD.

: A numerous factor can cause infertility, but around one of four reproductive failure cases remain unexplained and diagnosed as idiopathic infertility. In the past few decades, analysis of gene polymorphisms takes a significant place in pathogenesis of infertility. The aim of this study was to evaluate the possible role of hemostasis-related gene polymorphisms in unexplained infertility. The study includes 117 female patients with idiopathic infertility and 130 fertile women with at least one born child. Eight polymorphisms important for hemostasis (ITGB3 1565T>C, FV 1691G>A, FII 20210G>A, MTHFR 677C>T and 1298A>C, ATIII 786G>A, PAI-14G/5G and ACE I/D) were genotyped by real-time PCR system. The frequencies of alleles and genotypes of examined polymorphisms were analyzed in SPSS statistical program, whereas gene interactions were identified using the GMDR software. Examination of etiological factors has shown that family history is a significant factor in assessing individual risk for infertility. The alleles and genotypes frequency of FV 1691G>A and FII 20210G>A polymorphisms were statistically different between control and patient group leading to a greater risk for infertility. The analysis of epistatic relationship between examined hemostasis-related gene polymorphisms identified more complex high-risk genotypes associated with infertility. Our results suggest that positive family history could be important predictive factor for fertility problems, pointing to the potential hereditary basis of this condition. Polymorphisms FVL and FII prothrombin are independent risk factors for idiopathic infertility, whereas multilocus interactions approach should be taken into consideration for the future research.

We previously reported that the Bcl-2 family member human Bcl-rambo, also known as BCL2L13, induces apoptosis in human embryonic kidney 293T cells. Mouse Bcl-rambo has recently been reported to mediate mitochondrial fragmentation and mitophagy. In the present study, we showed that the transfection of human Bcl-rambo and its microtubule-associated protein light chain 3-interacting region motif mutant (W276A/I279A) caused mitochondrial fragmentation and the perinuclear accumulation of fragmented mitochondria in human lung adenocarcinoma A549 cells. In comprehensive screening using the Drosophila model in which human Bcl-rambo was ectopically expressed in eye imaginal discs, voltage-dependent anion channels (VDAC), also known as mitochondrial porin, were found to manifest a genetic interaction with human Bcl-rambo. In addition to human adenine nucleotide translocase (ANT) 1 and ANT2, the human Bcl-rambo protein bound to human VDAC1, albeit to a lesser extent than ANT2. Moreover, human VDAC1 and human VDAC2 in particular promoted the activation of effector caspases only when they were co-expressed with human Bcl-rambo in 293T cells. Bcl-rambo induced the perinuclear accumulation of fragmented mitochondria by the knockdown of VDAC1, VDAC2, and VDAC3 in A549 cells. Thus, the present study revealed that human Bcl-rambo and VDAC cooperatively promote the activation of effector caspases in human cultured cells., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Background: Recent genome-wide association studies (GWAS) have identified several genetic variants that influence the risk of dyslipidemia and coronary artery disease (CAD). In this study, we have examined the potential association of five SNPs variants related to lipid pathway, previously identified in GWAS studies (ZNF259 C>G, CETP I405VA/G, LPA C>T, LPLS447X and PSRC1 A>G) with CAD., Methods: Two hundred and ninety subjects including 194 patients with coronary artery disease and 96 controls were enrolled, followed by the analyses of anthropometric/biochemical parameters. Genotyping was carried out using Taq-Man real-time PCR based method. The association of the genetic polymorphisms with CAD was determined using univariate and multivariate analyses., Results: CAD patients had a higher (p < 0.05) fasting blood glucose (FBG), total cholesterol (TC), high sensitivity C-reactive protein (hs-CRP), low-density lipoprotein cholesterol (LDL-C) and waist circumference. Results showed that subjects with CETP rs5882 genetic variant, AA&AG genotypes, had a higher risk of developing Coronary artery disease [OR: 2.1, 95% CI (1.2-4.1), p value = 0.015]. Also subjects who carried the G allele of the ZNF259 polymorphism were at an increased the risk of developing CAD [OR 1.86, 95% CI: 1.06-3.25, p value = 0.029] and had an increased TC, LDL and TG levels (p < 0.05). Furthermore, no statistically significant association was found between genetic polymorphisms of PSRC1 A>G, LPL S447X and LPA C>T and CAD., Conclusion: We identified a relationship between a genetic variant in CETP and ZNF259 gene with CAD and CAD and lipid profile, respectively. Further investigation in a larger population may help to investigate the value of emerging marker as a risk stratification marker in CAD and its risk factors., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Recent studies of the tumor genome seek to identify cancer pathways as groups of genes in which mutations are epistatic with one another or, specifically, "mutually exclusive." Here, we show that most mutations are mutually exclusive not due to pathway structure but to interactions with disease subtype and tumor mutation load. In particular, many cancer driver genes are mutated preferentially in tumors with few mutations overall, causing mutations in these cancer genes to appear mutually exclusive with numerous others. Researchers should view current epistasis maps with caution until we better understand the multiple cause-and-effect relationships among factors such as tumor subtype, positive selection for mutations, and gross tumor characteristics including mutational signatures and load., (Copyright © 2019. Published by Elsevier Inc.)

The understanding that differences in biological epistasis may impact disease risk, diagnosis, or disease management stands in wide contrast to the unavailability of widely accepted large-scale epistasis analysis protocols. Several choices in the analysis workflow will impact false-positive and false-negative rates. One of these choices relates to the exploitation of particular modelling or testing strategies. The strengths and limitations of these need to be well understood, as well as the contexts in which these hold. This will contribute to determining the potentially complementary value of epistasis detection workflows and is expected to increase replication success with biological relevance. In this contribution, we take a recently introduced regression-based epistasis detection tool as a leading example to review the key elements that need to be considered to fully appreciate the value of analytical epistasis detection performance assessments. We point out unresolved hurdles and give our perspectives towards overcoming these.

The serotonin transporter (5-HTTLPR) and brain-derived neurotrophic factor (BDNF) gene polymorphisms have been associated with risk for affective disorders and functional variability of the amygdala. We examined whether the two genotypes interactively influence intrinsic functional connectivity (FC) of the amygdala and whether FC mediates the genetic association with anxiety. Eighty genotyped healthy adults underwent resting state fMRI and completed the self-reported State-Trait Anxiety Inventory. Interactive genetic association with anxiety was observed such that effects of 5-HTTLPR depended on the BDNF Val66Met polymorphism (rs6265 variant), with higher anxiety scores in short and Met carriers compared to the other allelic groups. Voxel-wise FC with left and right amygdala seeds identified regions that were sensitive to variability in anxiety scores. A significant moderated mediation model demonstrated that the effect of 5-HTTLPR genotype on anxiety, moderated by BDNF Val66Met genotype, was fully mediated by FC between the left amygdala and the right dorsolateral prefrontal cortex, a cognitive control-related region, during a task-free state. FC was highest in carriers of the 5-HTTLPR short allele and BDNF Met allele. These findings establish intrinsic amygdala-prefrontal functional connectivity as a potential intermediate phenotype for anxiety, an important step toward identification of causal pathways for vulnerability to affective disorders., (Copyright © 2019. Published by Elsevier B.V.)

Protein evolution occurs via restricted evolutionary paths that are influenced by both previous and subsequent mutations. This effect, termed epistasis, is critical in population genetics, drug resistance, and immune escape; however, the effect of epistasis on the level of protein fitness is less well characterized. We generated and characterized a 6615-member library of all two-amino acid combinations in a highly mutable loop of a virus-like particle. This particle is a model of protein self-assembly and a promising vehicle for drug delivery and imaging. In addition to characterizing the effect of all double mutants on assembly, thermostability, and acid stability, we observed many instances of epistasis, in which combinations of mutations are either more deleterious or more beneficial than expected. These results were used to generate rules governing the effects of multiple mutations on the self-assembly of the virus-like particle.

Experimental evolution has become an increasingly common approach for studying evolutionary phenomena, as well as uncovering physiological connections in a manner complementary to traditional genetics. Here I describe the development of Methylobacterium as a model system for using experimental evolution to study questions at the intersection of metabolism and evolution. Each experiment was initiated to address a particular question inspired by patterns in natural methylotrophs, such as tradeoffs between single-carbon and multi-carbon growth, or the challenges involved in incorporating novel metabolic pathways or genes with poor codon usage that are acquired via horizontal gene transfer. What I could not have appreciated initially, however, was just how many fortuitous surprise findings would emerge. These have ranged from the repeatability of evolution, complex dynamics within populations, epistasis between beneficial mutations, and even the ability to use simple mathematical models to generate testable, quantitative hypotheses about the fitness landscape.

Identification of genetic variants with effects on trait variability can provide insights into the biological mechanisms that control variation and can identify potential interactions. We propose a two-degree-of-freedom test for jointly testing mean and variance effects to identify such variants. We implement the test in a linear mixed model, for which we provide an efficient algorithm and software. To focus on biologically interesting settings, we develop a test for dispersion effects, that is, variance effects not driven solely by mean effects when the trait distribution is non-normal. We apply our approach to body mass index in the subsample of the UK Biobank population with British ancestry (n ~408,000) and show that our approach can increase the power to detect associated loci. We identify and replicate novel associations with significant variance effects that cannot be explained by the non-normality of body mass index, and we provide suggestive evidence for a connection between leptin levels and body mass index variability.

Understanding the relationship between protein sequence, function, and stability is a fundamental problem in biology. The essential function of many proteins that fold into a specific structure is their ability to bind to a ligand, which can be assayed for thousands of mutated variants. However, binding assays do not distinguish whether mutations affect the stability of the binding interface or the overall fold. Here, we introduce a statistical method to infer a detailed energy landscape of how a protein folds and binds to a ligand by combining information from many mutated variants. We fit a thermodynamic model describing the bound, unbound, and unfolded states to high quality data of protein G domain B1 binding to IgG-Fc. We infer distinct folding and binding energies for each mutation providing a detailed view of how mutations affect binding and stability across the protein. We accurately infer the folding energy of each variant in physical units, validated by independent data, whereas previous high-throughput methods could only measure indirect changes in stability. While we assume an additive sequence-energy relationship, the binding fraction is epistatic due its nonlinear relation to energy. Despite having no epistasis in energy, our model explains much of the observed epistasis in binding fraction, with the remaining epistasis identifying conformationally dynamic regions.

Low temperature stress is one of the major abiotic stresses limiting the productivity of Geng (japonica) rice grown the temperate regions as well as in tropical high lands worldwide. To develop rice varieties with improved cold tolerance (CT) at the reproductive stage, 84 BC2 CT introgression lines (ILs) were developed from five populations through backcross breeding. These CT ILs plus 310 random ILs from the same BC populations were used for dissecting genetic networks underlying CT in rice by detecting QTLs and functional genetic units (FGUs) contributing to CT. Seventeen major QTLs for CT were identified using five selective introgression populations and the method of segregation distortion. Of them, three QTLs were confirmed using the random populations and seven others locate in the regions with previously reported CT QTLs/genes. Using multi-locus probability tests and linkage disequilibrium (LD) analyses, 46 functional genetic units (FGUs) (37 single loci and 9 association groups or AGs) distributed in 37 bins (~20%) across the rice genome for CT were detected. Together, each of the CT loci (bins) was detected in 1.7 populations, including 18 loci detected in two or more populations. Putative genetic networks (multi-locus structures) underlying CT were constructed based on strong non-random associations between or among donor alleles at the unlinked CT loci/FGUs identified in the CT ILs, suggesting the presence of strong epistasis among the detected CT loci. Our results demonstrated the power and usefulness of using selective introgression for simultaneous improvement and genetic dissection of complex traits such as CT in rice., Competing Interests: The authors have declared that no competing interests exist.

Elucidating the relationships between antimicrobial resistance and virulence is key to understanding the evolution and population dynamics of resistant pathogens. Here, we show that the susceptibility of the gram-positive bacterium Listeria monocytogenes to the antibiotic fosfomycin is a complex trait involving interactions between resistance and virulence genes and the environment. We found that a FosX enzyme encoded in the listerial core genome confers intrinsic fosfomycin resistance to both pathogenic and non-pathogenic Listeria spp. However, in the genomic context of the pathogenic L. monocytogenes, FosX-mediated resistance is epistatically suppressed by two members of the PrfA virulence regulon, hpt and prfA, which upon activation by host signals induce increased fosfomycin influx into the bacterial cell. Consequently, in infection conditions, most L. monocytogenes isolates become susceptible to fosfomycin despite possessing a gene that confers high-level resistance to the drug. Our study establishes the molecular basis of an epistatic interaction between virulence and resistance genes controlling bacterial susceptibility to an antibiotic. The reported findings provide the rationale for the introduction of fosfomycin in the treatment of Listeria infections even though these bacteria are intrinsically resistant to the antibiotic in vitro., Competing Interests: The authors have declared that no competing interests exist.

The Red Queen hypothesis predicts that host-parasite coevolutionary dynamics can select for host resistance through increased genetic diversity, recombination and evolutionary rates. However, in haplodiploid organisms such as the honeybee (Apis mellifera), models suggest the selective pressure is weaker than in diploids. Haplodiploid sex determination, found in A. mellifera, can allow deleterious recessive alleles to persist in the population through the diploid sex with negative effects predominantly expressed in the haploid sex. To overcome these negative effects in haploid genomes, epistatic interactions have been hypothesized to play an important role. Here, we use the interaction between A. mellifera and the parasitic mite Varroa destructor to test epistasis in the expression of resistance, through the inhibition of parasite reproduction, in haploid drones. We find novel loci on three chromosomes which explain over 45% of the resistance phenotype. Two of these loci interact only additively, suggesting their expression is independent of each other, but both loci interact epistatically with the third locus. With drone offspring inheriting only one copy of the queen's chromosomes, the drones will only possess one of two queen alleles throughout the years-long lifetime of the honeybee colony. Varroa, in comparison, completes its highly inbred reproductive cycle in a matter of weeks, allowing it to rapidly evolve resistance. Faced with the rapidly evolving Varroa, a diversity of pathways and epistatic interactions for the inhibition of Varroa reproduction could therefore provide a selective advantage to the high levels of recombination seen in A. mellifera. This allows for the remixing of phenotypes despite a fixed queen genotype., (© 2018 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2018 European Society For Evolutionary Biology.)

Transcription factors (TFs) direct developmental transitions by binding to target DNA sequences, influencing gene expression and establishing complex gene-regultory networks. To systematically determine the molecular components that enable or constrain TF activity, we investigated the genomic occupancy of FOXA2, GATA4 and OCT4 in several cell types. Despite their classification as pioneer factors, all three TFs exhibit cell-type-specific binding, even when supraphysiologically and ectopically expressed. However, FOXA2 and GATA4 can be distinguished by low enrichment at loci that are highly occupied by these factors in alternative cell types. We find that expression of additional cofactors increases enrichment at a subset of these sites. Finally, FOXA2 occupancy and changes to DNA accessibility can occur in G 1 -arrested cells, but subsequent loss of DNA methylation requires DNA replication.

Background: Congenital hyperinsulinism (CHI) is frequently caused by mutations in one of the KATP channel subunits encoded by the genes ABCC8 and KCNJ11. The effect of simultaneous mutations in both of these genes on the pancreatic β-cell function is not known and patients with CHI carrying both ABCC8 and KCNJ11 mutations have not yet been reported. We questioned if a combination of heterozygous mutations in the ABCC8 and KCNJ11 genes could also lead to β-cell dysfunction presenting as CHI., Methods: As a model, we used a patient with transient CHI that paternally inherited novel heterozygous mutations in ABCC8 (p.Tyr1293Asp) and KCNJ11 (p.Arg50Trp) genes. The pathogenic effects on the pancreatic β-cells function were examined in an in vitro functional study using radioactive rubidium efflux assay., Results: We showed that the activation of the mutated KATP channels by diazoxide was decreased by 60.9% in the channels with the heterozygous combination of both mutations compared to the wild type channels. This could indicate the pathogenic effect on the pancreatic β-cell function leading to CHI although conclusive evidence is needed to be added., Conclusions: Our findings may widen the spectrum of genetic causes of CHI and suggest a novel pathogenic mechanism of CHI that must however, be further investigated.

Wnt signaling through the low-density lipoprotein-related receptor 5 (Lrp5) coreceptor regulates osteoblast maturation, matrix mineralization, and intermediary metabolism. In the mature osteoblast, signals downstream of Lrp5 are required for normal long-chain fatty acid β-oxidation. Mice rendered deficient for this coreceptor in osteoblasts and osteocytes accumulate body fat with elevated serum lipid levels but retain normal insulin sensitivity. In the present study, we challenged Lrp5-mutant mice with a high-fat diet (HFD) to determine whether they were more susceptible to diet-induced disturbances in glucose homeostasis. The HFD-fed Lrp5 mutant mice maintained a low bone mass phenotype with an increase in adipose tissue mass and hypertriglyceridemia and hypercholesterolemia. Examination of glucose metabolism revealed that Lrp5 deficiency in the osteoblast also resulted in hyperglycemia and hyperinsulinemia, with reductions in glucose tolerance, insulin sensitivity, and serum undercarboxylated osteocalcin. The results from in vivo genetic epistasis and in vitro studies suggest that this phenotype proceeds via the accumulation of diacylglycerol species and impaired insulin signaling in Lrp5-deficient osteoblasts. In turn, glucose uptake and osteocalcin production are diminished in mutant osteoblasts. Taken together, these data identify a link between Wnt-Lrp5 signaling and insulin signaling in the osteoblast that has the potential to influence energy balance and compound the detrimental effects of a HFD on whole-body metabolism., (Copyright © 2017 Endocrine Society.)

To further our understanding of the complexity and genetic heterogeneity of rare diseases, it has become essential to shed light on how combinations of variants in different genes are responsible for a disease phenotype. With the appearance of a resource on digenic diseases, it has become possible to evaluate how digenic combinations differ in terms of the phenotypes they produce. All instances in this resource were assigned to two classes of digenic effects, annotated as true digenic and composite classes. Whereas in the true digenic class variants in both genes are required for developing the disease, in the composite class, a variant in one gene is sufficient to produce the phenotype, but an additional variant in a second gene impacts the disease phenotype or alters the age of onset. We show that a combination of variant, gene and higher-level features can differentiate between these two classes with high accuracy. Moreover, we show via the analysis of three digenic disorders that a digenic effect decision profile, extracted from the predictive model, motivates why an instance was assigned to either of the two classes. Together, our results show that digenic disease data generates novel insights, providing a glimpse into the oligogenic realm., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Previously, we showed that cortical mineralization is coordinately adjusted to mechanically offset external bone size differences between A/J (narrow) and C57BL/6J (wide) mouse femora to achieve whole bone strength equivalence at adulthood. The identity of the genes and their interactions that are responsible for establishing this homeostatic state (ie, canalization) remain unknown. We hypothesize that these inbred strains, whose interindividual differences in bone structure and material properties mimic that observed among humans, achieve functional homeostasis by differentially adjusting key molecular pathways regulating external bone size and mineralization throughout growth. The cortices of A/J and C57BL/6J male mouse femora were phenotyped and gene expression levels were assessed across growth (ie, ages 2, 4, 6, 8, 12, 16 weeks). A difference in total cross-sectional area (p < 0.01) and cortical tissue mineral density were apparent between mouse strains by age 2 weeks and maintained at adulthood (p < 0.01). These phenotypic dissimilarities corresponded to gene expression level differences among key regulatory pathways throughout growth. A/J mice had a 1.55- to 7.65-fold greater expression among genes inhibitory to Wnt pathway induction, whereas genes involved in cortical mineralization were largely upregulated 1.50- to 3.77-fold to compensate for their narrow diaphysis. Additionally, both mouse strains showed an upregulation among Wnt pathway antagonists corresponding to the onset of adult ambulation (ie, increased physiological loads). This contrasts with other studies showing an increase in Wnt pathway activation after functionally isolated, experimental in vivo loading regimens. A/J and C57BL/6J long bones provide a model to develop a systems-based approach to identify individual genes and the gene-gene interactions that contribute to trait differences between the strains while being involved in the process by which these traits are coordinately adjusted to establish similar levels of mechanical function, thus providing insight into the process of canalization. © 2017 American Society for Bone and Mineral Research., (© 2017 American Society for Bone and Mineral Research.)

The replication protein A (RPA) complex binds single-stranded DNA generated at stalled replication forks and recruits other DNA repair proteins to promote recovery of these forks. Here, we identify Ewing tumor-associated antigen 1 (ETAA1), which has been linked to susceptibility to pancreatic cancer, as a new repair protein that is recruited to stalled forks by RPA. We demonstrate that ETAA1 interacts with RPA through two regions, each of which resembles two previously identified RPA-binding domains, RPA70N-binding motif and RPA32C-binding motif, respectively. In response to replication stress, ETAA1 is recruited to stalled forks where it colocalizes with RPA, and this recruitment is diminished when RPA is depleted. Notably, inactivation of the ETAA1 gene increases the collapse level of the stalled replication forks and decreases the recovery efficiency of these forks. Moreover, epistasis analysis shows that ETAA1 stabilizes stalled replication forks in an ataxia telangiectasia and Rad3-related protein (ATR)-independent manner. Thus, our results reveal that ETAA1 is a novel RPA-interacting protein that promotes restart of stalled replication forks., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

Soybean [Glycine max (L.) Merr.] is an economically important crop that is grown worldwide. Sudden death syndrome (SDS), caused by Fusarium virguliforme, is one of the top yield-limiting diseases in soybean. However, the genetic basis of SDS resistance, especially with respect to epistatic interactions, is still unclear. To better understand the genetic architecture of soybean SDS resistance, genome-wide association and epistasis studies were performed using a population of 214 germplasm accessions and 31,914 SNPs from the SoySNP50K Illumina Infinium BeadChip. Twelve loci and 12 SNP-SNP interactions associated with SDS resistance were identified at various time points after inoculation. These additive and epistatic loci together explained 24-52% of the phenotypic variance. Disease-resistant, pathogenesis-related and chitin- and wound-responsive genes were identified in the proximity of peak SNPs, including stress-induced receptor-like kinase gene 1 (SIK1), which is pinpointed by a trait-associated SNP and encodes a leucine-rich repeat-containing protein. We report that the proportion of phenotypic variance explained by identified loci may be considerably improved by taking epistatic effects into account. This study shows the necessity of considering epistatic effects in soybean SDS resistance breeding using marker-assisted and genomic selection approaches. Based on our findings, we propose a model for soybean root defense against the SDS pathogen. Our results facilitate identification of the molecular mechanism underlying SDS resistance in soybean, and provide a genetic basis for improvement of soybean SDS resistance through breeding strategies based on additive and epistatic effects., (© 2015 The Authors The Plant Journal © 2015 John Wiley & Sons Ltd.)

Precise prediction for genetic architecture of complex traits is impeded by the limited understanding on genetic effects of complex traits, especially on gene-by-gene (GxG) and gene-by-environment (GxE) interaction. In the past decades, an explosion of high throughput technologies enables omics studies at multiple levels (such as genomics, transcriptomics, proteomics, and metabolomics). The analyses of large omics data, especially two-loci interaction analysis, are very time intensive. Integrating the diverse omics data and environmental effects in the analyses also remain challenges. We proposed mixed linear model approaches using GPU (Graphic Processing Unit) computation to simultaneously dissect various genetic effects. Analyses can be performed for estimating genetic main effects, GxG epistasis effects, and GxE environment interaction effects on large-scale omics data for complex traits, and for estimating heritability of specific genetic effects. Both mouse data analyses and Monte Carlo simulations demonstrated that genetic effects and environment interaction effects could be unbiasedly estimated with high statistical power by using the proposed approaches.

Context: Genetic factors account for 60-80% of the areal bone mineral density (aBMD) variance, whereas the heritability of bone microstructure is not clearly established. aBMD and microstructure are under the control of osteocytes, which regulate bone formation through the expression of molecules such as sclerostin (SOST) and periostin (POSTN)., Objective: We hypothesized that additive genetic effects contribute to serum levels of SOST and POSTN and thereby to the individual variance of bone microstructure., Subjects and Methods: In a retrospective analysis of 432 subjects from the Geneva Retiree Cohort age 64.9 ± 1.4 years and 96 of their offspring age 37.9 ± 5.7 years, we measured serum SOST (sSOST) and serum POSTN (sPOSTN), distal radius and tibia microstructure, hip and lumbar spine aBMD, and bone turnover markers, Heritability (h(2), %) was calculated as twice the slope of the regression (β) between parents and offspring., Results: cPOSTN levels were significantly higher in men than women and in offspring than parents. h(2) values for bone microstructural traits ranged from 22-64% depending on the envelope (trabecular [Tb] or cortical [Ct]) and skeletal site (radius or tibia), whereas h(2) for sPOSTN and sSOST was 50% and 40%, respectively. sPOSTN was positively associated with Tb bone volume on total volume and Ct thickness, and negatively with Ct porosity. The associations for Ct parameters remain significant after adjustment for propetide of type-I procollagen, cross-linked telopeptide of type I collagen, femoral neck aBMD, sex or age. After adjustment of bone traits for sPOSTN, h(2) values decreased for several Tb and Ct bone parameters, but not for aBMD. In contrast, adjusting for sSOST did not alter h(2) values for bone traits., Conclusions: Additive genetic effects account for a substantial proportion of the individual variance of bone microstructure, sPOSTN, and sSOST. sPOSTN is largely inherited as a sex-related trait and carries an important contribution to the heritability of bone microstructure, indicating that these traits are at least partly determined by common genetic effects.

Directional cell migration is fundamental for neural development, and extracellular factors are pivotal for this process. Heparan sulfate proteoglycans (HSPGs) that carry long chains of differentially modified sugar residues contribute to extracellular matrix; however, the functions of HSPG in guiding cell migration remain elusive. Here, we used the Caenorhabditis elegans mutant pool from the Million Mutation Project and isolated a mutant allele of the heparan sulfate-modifying enzyme glucuronyl C5-epimerase HSE-5. Loss-of-function of this enzyme resulted in defective Q neuroblast migration. We showed that hse-5 controlled Q cell migration in a cell non-autonomous manner. By performing live cell imaging in hse-5 mutant animals, we found that hse-5 controlled initial polarization during Q neuroblast migration. Furthermore, our genetic epistasis analysis demonstrated that lon-2 might act downstream of hse-5. Finally, rescue of the hse-5 mutant phenotype by expression of human and mouse hse-5 homologs suggested a conserved function for this gene in neural development. Taken together, our results indicated that proper HSPG modification in the extracellular matrix by HSE-5 is essential for neuroblast polarity during migration., (Copyright © 2015 Elsevier Inc. All rights reserved.)

The study of biological pathways is key to a large number of systems analyses. However, many relevant tools consider a limited number of pathway sources, missing out on many genes and gene-to-gene connections. Simply pooling several pathways sources would result in redundancy and the lack of systematic pathway interrelations. To address this, we exercised a combination of hierarchical clustering and nearest neighbor graph representation, with judiciously selected cutoff values, thereby consolidating 3215 human pathways from 12 sources into a set of 1073 SuperPaths. Our unification algorithm finds a balance between reducing redundancy and optimizing the level of pathway-related informativeness for individual genes. We show a substantial enhancement of the SuperPaths' capacity to infer gene-to-gene relationships when compared with individual pathway sources, separately or taken together. Further, we demonstrate that the chosen 12 sources entail nearly exhaustive gene coverage. The computed SuperPaths are presented in a new online database, PathCards, showing each SuperPath, its constituent network of pathways, and its contained genes. This provides researchers with a rich, searchable systems analysis resource. Database URL: http://pathcards.genecards.org/, (© The Author(s) 2015. Published by Oxford University Press.)

Economically important traits in many species generally show polygenic, quantitative inheritance. The components of genetic variation (additive, dominant and epistatic effects) of these traits conferred by multiple genes in shared biological pathways remain to be defined. Here, we investigated 11 full-length genes in cellulose biosynthesis, on 10 growth and wood-property traits, within a population of 460 unrelated Populus tomentosa individuals, via multi-gene association. To validate positive associations, we conducted single-marker analysis in a linkage population of 1,200 individuals. We identified 118, 121, and 43 associations (P< 0.01) corresponding to additive, dominant, and epistatic effects, respectively, with low to moderate proportions of phenotypic variance (R(2)). Epistatic interaction models uncovered a combination of three non-synonymous sites from three unique genes, representing a significant epistasis for diameter at breast height and stem volume. Single-marker analysis validated 61 associations (false discovery rate, Q ≤ 0.10), representing 38 SNPs from nine genes, and its average effect (R(2) = 3.8%) nearly 2-fold higher than that identified with multi-gene association, suggesting that multi-gene association can capture smaller individual variants. Moreover, a structural gene-gene network based on tissue-specific transcript abundances provides a better understanding of the multi-gene pathway affecting tree growth and lignocellulose biosynthesis. Our study highlights the importance of pathway-based multiple gene associations to uncover the nature of genetic variance for quantitative traits and may drive novel progress in molecular breeding., (© The Author 2014. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.)

The adaptive evolution of large asexual populations is generally characterized by competition between clones carrying different beneficial mutations. Interference slows down the adaptation speed and makes the theoretical description of the dynamics more complex with respect to the successional occurrence and fixation of beneficial mutations typical of small populations. A simplified modeling framework considering multiple beneficial mutations with equal and constant fitness advantage is known to capture some of the essential features of laboratory evolution experiments. However, in these experiments the relative advantage of a beneficial mutation is generally dependent on the genetic background. In particular, the general pattern is that, as mutations in different loci accumulate, the relative advantage of new mutations decreases, a trend often referred to as "diminishing return" epistasis. Here, we propose a phenomenological model that generalizes the fixed-advantage framework to include this negative epistasis in a simple way. We evaluate analytically as well as with direct simulations the quantitative consequences of diminishing returns on the evolutionary dynamics. The speed of adaptation decreases in time and reaches a limit value corresponding to neutral evolution in the long time limit. This corresponds to an increase of the diversity in terms of "classes of mutation" in the population. Finally, we show how the model can be compared with dynamic data on fitness and number of beneficial mutations from laboratory evolution experiments., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Background: Respiratory system development is regulated by a complex series of endoderm-mesoderm interactions that are not fully understood. Recently Xenopus has emerged as an alternative model to investigate early respiratory system development, but the extent to which the morphogenesis and molecular pathways involved are conserved between Xenopus and mammals has not been systematically documented., Results: In this study, we provide a histological and molecular atlas of Xenopus respiratory system development, focusing on Nkx2.1+ respiratory cell fate specification in the developing foregut. We document the expression patterns of Wnt/β-catenin, fibroblast growth factor (FGF), and bone morphogenetic protein (BMP) signaling components in the foregut and show that the molecular mechanisms of respiratory lineage induction are remarkably conserved between Xenopus and mice. Finally, using several functional experiments we refine the epistatic relationships among FGF, Wnt, and BMP signaling in early Xenopus respiratory system development., Conclusions: We demonstrate that Xenopus trachea and lung development, before metamorphosis, is comparable at the cellular and molecular levels to embryonic stages of mouse respiratory system development between embryonic days 8.5 and 10.5. This molecular atlas provides a fundamental starting point for further studies using Xenopus as a model to define the conserved genetic programs controlling early respiratory system development., (© 2014 Wiley Periodicals, Inc.)

Hypomorphic mutation of apoptosis-inducing factor (AIF) in the whole body or organ-specific knockout of AIF compromises the activity of respiratory chain complexes I and IV, as it confers resistance to obesity and diabetes induced by high-fat diet. The mitochondrial defect induced by AIF deficiency can be explained by reduced AIF-dependent mitochondrial import of CHCHD4, which in turn is required for optimal import and assembly of respiratory chain complexes. Here we show that, as compared to wild type control littermates, mice with a heterozygous knockout of CHCHD4 exhibit reduced weight gain when fed with a Western style high-fat diet. This finding suggests widespread metabolic epistasis among AIF and CHCHD4. Targeting either of these proteins or their functional interaction might constitute a novel strategy to combat obesity.

UNC-6/Netrin is a conserved axon guidance cue that can mediate both attraction and repulsion. We previously discovered that attractive UNC-40/DCC receptor signaling stimulates growth cone filopodial protrusion and that repulsive UNC-40-UNC-5 heterodimers inhibit filopodial protrusion in C. elegans. Here, we identify cytoplasmic signaling molecules required for UNC-6-mediated inhibition of filopodial protrusion involved in axon repulsion. We show that the Rac-like GTPases CED-10 and MIG-2, the Rac GTP exchange factor UNC-73/Trio, UNC-44/Ankyrin and UNC-33/CRMP act in inhibitory UNC-6 signaling. These molecules were required for the normal limitation of filopodial protrusion in developing growth cones and for inhibition of growth cone filopodial protrusion caused by activated MYR::UNC-40 and MYR::UNC-5 receptor signaling. Epistasis studies using activated CED-10 and MIG-2 indicated that UNC-44 and UNC-33 act downstream of the Rac-like GTPases in filopodial inhibition. UNC-73, UNC-33 and UNC-44 did not affect the accumulation of full-length UNC-5::GFP and UNC-40::GFP in growth cones, consistent with a model in which UNC-73, UNC-33 and UNC-44 influence cytoskeletal function during growth cone filopodial inhibition., (© 2014. Published by The Company of Biologists Ltd.)

The 5'-3' resection of DNA ends is a prerequisite for the repair of DNA double strand breaks by homologous recombination, microhomology-mediated end joining, and single strand annealing. Recent studies in yeast have shown that, following initial DNA end processing by the Mre11-Rad50-Xrs2 complex and Sae2, the extension of resection tracts is mediated either by exonuclease 1 or by combined activities of the RecQ family DNA helicase Sgs1 and the helicase/endonuclease Dna2. Although human DNA2 has been shown to cooperate with the BLM helicase to catalyze the resection of DNA ends, it remains a matter of debate whether another human RecQ helicase, WRN, can substitute for BLM in DNA2-catalyzed resection. Here we present evidence that WRN and BLM act epistatically with DNA2 to promote the long-range resection of double strand break ends in human cells. Our biochemical experiments show that WRN and DNA2 interact physically and coordinate their enzymatic activities to mediate 5'-3' DNA end resection in a reaction dependent on RPA. In addition, we present in vitro and in vivo data suggesting that BLM promotes DNA end resection as part of the BLM-TOPOIIIα-RMI1-RMI2 complex. Our study provides new mechanistic insights into the process of DNA end resection in mammalian cells., (© 2014 by The American Society for Biochemistry and Molecular Biology, Inc.)

Higher cognitive performance, maintenance of mental health and psychological well-being require adequate prefrontal cortex (PFC) function. "Inverted U-shaped" dopamine model indicates optimal PFC dopamine level is important to attain its function while high or low levels have adverse effects. Catechol-O-methyltransferase (COMT) and methylenetetrahydrofolate reductase (MTHFR) may be involved in this complex non-linear PFC dopamine regulation. We addressed whether genetic variation reflecting COMT and MTHFR activities can explain the inter-individual mental health differences in healthy Japanese men (n=188). The mental health was measured by Mental Health Inventory (MHI)-5 score. The rs4633-rs4818-rs4680 haplotypes were used to represent the multilevel COMT activities, while for MTHFR, the functional single polymorphism, rs1801133 (C677T), was used. We examined the effectiveness of haplotype-based association analysis of COMT on mental health together with studying its interaction with MTHFR-C677T. As a result, the relation between activity-ranked COMT genotype and MHI-5 score showed a tendency to fit into an "inverted U-shaped" quadratic curve (P=0.054). This curvilinear correlation was significant in the subjects with MTHFR-CC (P<0.001), but not with MTHFR T-allele carriers (P=0.793). Our pilot study implies a potential influence of COMT and MTHFR genotypic combination on normal variation of mental health., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Because increased oxidative stress may mediate the detrimental actions of enhanced sympathetic nervous activity on renal function and vice versa, we investigated the effect of the polymorphic Arg16Gly in the β2 -adrenoceptor (ADRB2) gene, Trp64Arg in the β3 -adrenoceptor (ADRB3) gene and C242T in the NADPH oxidase p22phox (CYBA) gene on estimated glomerular filtration rate (eGFR) in a Chinese population. Initially recruited from different outpatient services of HeBei General Hospital in northern China, 668 individuals were finally included in the study, with complete demographic information. Laboratory tests were performed and estimated glomerular filtration rate (eGFR) was derived from the Modification of Diet in Renal Disease (MDRD) equation for the Chinese population. Plasma noradrenaline levels and genotype were determined by HPLC and the TaqMan method, respectively. Only across the Arg16Gly polymorphism did eGFR show significant difference: it was lower in individuals with the Gly16Gly variation, who also had the highest plasma noradrenaline levels. This polymorphism remained a significant determinant of eGFR after multivariate analysis. Of importance, the multifactor dimensionality reduction method further detected a significant synergism between the Arg16Gly and C242T polymorphisms in reducing eGFR. These observations clarify the effects of the studied polymorphisms on eGFR and exemplify gene-gene interactions influencing renal function., (© 2014 Wiley Publishing Asia Pty Ltd.)