Interpreting Human Genetic Variation through Genomics and In Vivo Models

Improvements in genomic technology, both in the increased speed and reduced cost of sequencing, have expanded the appreciation of the abundance of human genetic variation. However the sheer amount of variation, as well as the varying type and genomic content of variation, poses a challenge in understanding the clinical consequence of a single mutation. This work uses several methodologies to interpret the observed variation in the human genome, and presents novel strategies for the prediction of allele pathogenicity.

Using the zebrafish model system as an in vivo assay of allele function, we identified a novel driver of Bardet-Biedl Syndrome (BBS) in CEP76. A combination of targeted sequencing of 785 cilia-associated genes in a cohort of BBS patients and subsequent in vivo functional assays recapitulating the human phenotype gave strong evidence for the role of CEP76 mutations in the pathology of an affected family. This portion of the work demonstrated the necessity of functional testing in validating disease-associated mutations, and added to the catalogue of known BBS disease genes.

Further study into the role of copy-number variations (CNVs) in a cohort of BBS patients showed the significant contribution of CNVs to disease pathology. Using high-density array comparative genomic hybridization (aCGH) we were able to identify pathogenic CNVs as small as several hundred bp. Dissection of constituent gene and in vivo experiments investigating epistatic interactions between affected genes allowed for an appreciation of several paradigms by which CNVs can contribute to disease. This study revealed that the contribution of CNVs to disease in BBS patients is much higher than previously expected, and demonstrated the necessity of consideration of CNV contribution in future (and retrospective) investigations of human genetic disease.

Finally, we used a combination of comparative genomics and in vivo complementation assays to identify second-site compensatory modification of pathogenic alleles. These pathogenic alleles, which are found compensated in other species (termed compensated pathogenic deviations [CPDs]), represent a significant fraction (from 3 – 10%) of human disease-associated alleles. In silico pathogenicity prediction algorithms, a valuable method of allele prioritization, often misrepresent these alleles as benign, leading to omission of possibly informative variants in studies of human genetic disease. We created a mathematical model that was able to predict CPDs and putative compensatory sites, and functionally showed in vivo that second-site mutation can mitigate the pathogenicity of disease alleles. Additionally, we made publically available an in silico module for the prediction of CPDs and modifier sites.

These studies have advanced the ability to interpret the pathogenicity of multiple types of human variation, as well as made available tools for others to do so as well.

Genetic interactions of the Ashergillus nidulans atmA(ATM) homolog with different components of the DNA damage response pathway

Ataxia telangiectasia mutated (ATM) is a phosphatidyl-3-kinase-related protein kinase that functions as a central regulator of the DNA damage response in eukaryotic cells. In humans, mutations in ATM cause the devastating neurodegenerative disease ataxia telangiectasia. Previously, we characterized the homolog of ATM (AtmA) in the filamentous fungus Aspergillus nidulans. In addition to its expected role in the DNA damage response, we found that AtmA is also required for polarized hyphal growth. Here, we extended these studies by investigating which components of the DNA damage response pathway are interacting with AtmA. The AtmA(ATM) loss of function caused synthetic lethality when combined with mutation in UvsB(ATR). Our results suggest that AtmA and UvsB are interacting and they are probably partially redundant in terms of DNA damage sensing and/or repairing and polar growth. We identified and inactivated A. nidulans chkA(CHK1) and chkB(CHK2) genes. These genes are also redundantly involved in A. nidulans DNA damage response. We constructed several combinations of double mutants for Delta atmA, Delta uvsB, Delta chkA, and Delta chkB. We observed a complex genetic relationship with these mutations during the DNA replication checkpoint and DNA damage response. Finally, we observed epistatic and synergistic interactions between AtmA, and bimE(APCI), ankA(WEE1) and the cdc2-related kinase npkA, at S-phase checkpoint and in response to DNA-damaging agents.

A bi-dimensional genome scan for prolificacy traits in pigs shows the existence of multiple epistatic QTL

Background: Prolificacy is the most important trait influencing the reproductive efficiency of pig production systems. The low heritability and sex-limited expression of prolificacy have hindered to some extent the improvement of this trait through artificial selection. Moreover, the relative contributions of additive, dominant and epistatic QTL to the genetic variance of pig prolificacy remain to be defined. In this work, we have undertaken this issue by performing one-dimensional and bi-dimensional genome scans for number of piglets born alive (NBA) and total number of piglets born (TNB) in a three generation Iberian by Meishan F2 intercross. Results: The one-dimensional genome scan for NBA and TNB revealed the existence of two genome-wide highly significant QTL located on SSC13 (P < 0.001) and SSC17 (P < 0.01) with effects on both traits. This relative paucity of significant results contrasted very strongly with the wide array of highly significant epistatic QTL that emerged in the bi-dimensional genome-wide scan analysis. As much as 18 epistatic QTL were found for NBA (four at P < 0.01 and five at P < 0.05) and TNB (three at P < 0.01 and six at P < 0.05), respectively. These epistatic QTL were distributed in multiple genomic regions, which covered 13 of the 18 pig autosomes, and they had small individual effects that ranged between 3 to 4% of the phenotypic variance. Different patterns of interactions (a × a, a × d, d × a and d × d) were found amongst the epistatic QTL pairs identified in the current work.Conclusions: The complex inheritance of prolificacy traits in pigs has been evidenced by identifying multiple additive (SSC13 and SSC17), dominant and epistatic QTL in an Iberian × Meishan F2 intercross. Our results demonstrate that a significant fraction of the phenotypic variance of swine prolificacy traits can be attributed to first-order gene-by-gene interactions emphasizing that the phenotypic effects of alleles might be strongly modulated by the genetic background where they segregate.

Maintenance of pre-mRNA secondary structure by epistatic selection.

Linkage disequilibrium between polymorphisms in a natural population may result from various evolutionary forces, including random genetic drift due to sampling of gametes during reproduction, restricted migration between subpopulations in a subdivided population, or epistatic selection. In this report, we present evidence that the majority of significant linkage disequilibria observed in introns of the alcohol dehydrogenase locus (Adh) of Drosophila pseudoobscura are due to epistatic selection maintaining secondary structure of precursor mRNA (pre-mRNA). Based on phylogenetic-comparative analysis and a likelihood approach, we propose secondary structure models of Adh pre-mRNA for the regions of the adult intron and intron 2 where clustering of linkage disequilibria has been observed. Furthermore, we applied the likelihood ratio test to the phylogenetically predicted secondary structure in intron 1. In contrast to the other two structures, polymorphisms associated with the more conserved stem-loop structure of intron 1 are in low frequency, and linkage disequilibria have not been observed. These findings are qualitatively consistent with a model of compensatory fitness interactions. This model assumes that mutations disrupting pairing in a secondary structural element are individually deleterious if they destabilize a functionally important structure; a second "compensatory" mutation, however, may restabilize the structure and restore fitness.

Interactions between quantitative trait loci in soybean in which trait variation at one locus is conditional upon a specific allele at another.

A large recombinant inbred population of soybean has been characterized for 220 restriction fragment-length polymorphism (RFLP) markers. Values for agronomic traits also have been measured. Quantitative trait loci (QTL) for height, yield, and maturity were located by their linkage to RFLP markers. QTL controlling large amounts of trait variation were analyzed for the dependence of trait variation on particular alleles at a second locus by comparing cumulative distributions of the trait for each genotype (four genotypes per pair of loci). Interesting pairs of loci were analyzed statistically with maximum likelihood and Monte Carlo comparison of additive and epistatic models. For each locus affecting height, variation was conditional upon the presence of a particular allele at a second unlinked locus that itself explained little or no trait variation. The results show that interactions between QTL are frequent and control large effects. Interactions distinguished between different QTL in a single linkage group and between QTL that affect different traits closely linked to one RFLP marker--i.e., distinguished between pleiotropy and closely linked genes. The implications for the evolution of inbreeding plants and for the construction of agronomic breeding strategies are discussed.

Interactions between endothelin receptor B and transcription factors Sox10 and Pax3 in the melanocyte lineage

Genetic interactions that underlie developmental processes such as cell differentiation and pattern formation are complex and difficult to elucidate. Neural Crest (NC) cells and their derivatives offer an optimal system in which to probe for these complex interactions as they acquire different cell fates and constitute a variety of structures. The transcription factors Sox10 and Pax3 as well as the transmembrane receptor Endothelin receptor b (Ednrb) are temporally and spatially co-expressed early in NC cells and mutations in these genes lead to similar hypopigmentation phenotypes due to a reduced number of NC-derived melanocyte precursors, the melanoblasts. The goal of this study was to establish whether Sox10 and Ednrb or Pax3 and Ednrb interact to promote normal murine melanocyte development. Crosses of Sox10 or Pax3 with Ednrb heterozygous mutants showed that the double heterozygous hypopigmentation phenotype was significantly more pronounced than phenotypes of single heterozygotes, implying that a synergistic interaction exists between Sox10 and Ednrb and Pax3 and Ednrb. This interaction was further explored by the attempt to rescue the Sox10 and Pax3 hypopigmentation phenotypes by the transgenic addition of Ednrb to melanoblasts. Pigmentation was completely restored in the Sox10 and partially restored in the Pax3 mutant mice. The comparison of the number of melanoblasts in transgenic and non-transgenic Sox10 mutant embryos showed that the transgenic rescue occurred as early as E11.5, a critical time for melanoblast population expansion. Cell survival assays indicated that the rescue was not due to an effect of the transgene on melanoblast survival. A novel phenotype arose when studying the interaction between Ednrb and Pax3. Newborns appeared normal but by 3.5 weeks of age, the affected pups were smaller than normal littermates and developed a dome-shaped head; some also developed thoracic kyphosis. Affected pups were dead by 4 weeks of age: 80% were Pax3Sp/+ and 75% were female. When compared to normal littermates, affected mice had brains with enlarged 4th ventricles and more glia while skeletal staining showed kyphosis, wider rib cages and pelvic differences. An epistatic interaction resulting from the mixing of genetic backgrounds that is exacerbated in the presence of Pax3 heterozygosity is suspected.

Interactions between Endothelin Receptor B and Transcription Factors Sox10 and Pax3 in the Melanocyte Lineage

Genetic interactions that underlie developmental processes such as cell differentiation and pattern formation are complex and difficult to elucidate. Neural Crest (NC) cells and their derivatives offer an optimal system in which to probe for these complex interactions as they acquire different cell fates and constitute a variety of structures. The transcription factors Sox10 and Pax3 as well as the transmembrane receptor Endothelin receptor b (Ednrb) are temporally and spatially co-expressed early in NC cells and mutations in these genes lead to similar hypopigmentation phenotypes due to a reduced number of NC-derived melanocyte precursors, the melanoblasts. The goal of this study was to establish whether Sox10 and Ednrb or Pax3 and Ednrb interact to promote normal murine melanocyte development. Crosses of Sox10 or Pax3 with Ednrb heterozygous mutants showed that the double heterozygous hypopigmentation phenotype was significantly more pronounced than phenotypes of single heterozygotes, implying that a synergistic interaction exists between Sox10 and Ednrb and Pax3 and Ednrb. This interaction was further explored by the attempt to rescue the Sox10 and Pax3 hypopigmentation phenotypes by the transgenic addition of Ednrb to melanoblasts. Pigmentation was completely restored in the Sox10 and partially restored in the Pax3 mutant mice. The comparison of the number of melanoblasts in transgenic and non-transgenic Sox10 mutant embryos showed that the transgenic rescue occurred as early as E11.5, a critical time for melanoblast population expansion. Cell survival assays indicated that the rescue was not due to an effect of the transgene on melanoblast survival. A novel phenotype arose when studying the interaction between Ednrb and Pax3. Newborns appeared normal but by 3.5 weeks of age, the affected pups were smaller than normal littermates and developed a dome-shaped head; some also developed thoracic kyphosis. Affected pups were dead by 4 weeks of age: 80% were Pax3Sp/+ and 75% were female. When compared to normal littermates, affected mice had brains with enlarged 4th ventricles and more glia while skeletal staining showed kyphosis, wider rib cages and pelvic differences. An epistatic interaction resulting from the mixing of genetic backgrounds that is exacerbated in the presence of Pax3 heterozygosity is suspected.

Gene-environment interactions in Autism Spectrum Disorders (ASD)

Spectrum Disorder (ASD), is a heterogeneous neurodevelopmental disorder with na estimated global prevalence rate of 17:10000, and a male to female ratio of 4:1. Patients with ASD presente language and communication difficulties and stereotyped behaviours. Comorbidity with other disorders, such as Intelectual Disability, Fragile-X syndrome (FXS) epilepsy and tuberous sclerosis frequently occurs. ASD presents amultifactorial etiopathology, and genetic factos alone are not suficiente to explain how the syndrome arises, with recente studies establishing ASD heritability at approximately 50%. Pre-, peri- and post-natal exposure to toxic environmental factos has been implicated in the development of ASD. Involvement of epigenetic regulatory mechanisms has been suggested, supported by the occurrence of autistic symptoms in patients with disorders aris ing from epigenetic mutations, such as FXS. A polygenic and epistatic model is a strong hypothesis to explain ASD. The main goal of this project is to identify specific exposure patterns to environmental toxicants in children diagnosed with ASD and integrate the results with genetic and epigenetic data.

Iis--integrated Interactome System: A Web-based Platform For The Annotation, Analysis And Visualization Of Protein-metabolite-gene-drug Interactions By Integrating A Variety Of Data Sources And Tools.

High-throughput screening of physical, genetic and chemical-genetic interactions brings important perspectives in the Systems Biology field, as the analysis of these interactions provides new insights into protein/gene function, cellular metabolic variations and the validation of therapeutic targets and drug design. However, such analysis depends on a pipeline connecting different tools that can automatically integrate data from diverse sources and result in a more comprehensive dataset that can be properly interpreted. We describe here the Integrated Interactome System (IIS), an integrative platform with a web-based interface for the annotation, analysis and visualization of the interaction profiles of proteins/genes, metabolites and drugs of interest. IIS works in four connected modules: (i) Submission module, which receives raw data derived from Sanger sequencing (e.g. two-hybrid system); (ii) Search module, which enables the user to search for the processed reads to be assembled into contigs/singlets, or for lists of proteins/genes, metabolites and drugs of interest, and add them to the project; (iii) Annotation module, which assigns annotations from several databases for the contigs/singlets or lists of proteins/genes, generating tables with automatic annotation that can be manually curated; and (iv) Interactome module, which maps the contigs/singlets or the uploaded lists to entries in our integrated database, building networks that gather novel identified interactions, protein and metabolite expression/concentration levels, subcellular localization and computed topological metrics, GO biological processes and KEGG pathways enrichment. This module generates a XGMML file that can be imported into Cytoscape or be visualized directly on the web. We have developed IIS by the integration of diverse databases following the need of appropriate tools for a systematic analysis of physical, genetic and chemical-genetic interactions. IIS was validated with yeast two-hybrid, proteomics and metabolomics datasets, but it is also extendable to other datasets. IIS is freely available online at: http://www.lge.ibi.unicamp.br/lnbio/IIS/.

Ant-caterpillar Antagonism At The Community Level: Interhabitat Variation Of Tritrophic Interactions In A Neotropical Savanna.

Ant foraging on foliage can substantially affect how phytophagous insects use host plants and represents a high predation risk for caterpillars, which are important folivores. Ant-plant-herbivore interactions are especially pervasive in cerrado savanna due to continuous ant visitation to liquid food sources on foliage (extrafloral nectaries, insect honeydew). While searching for liquid rewards on plants, aggressive ants frequently attack or kill insect herbivores, decreasing their numbers. Because ants vary in diet and aggressiveness, their effect on herbivores also varies. Additionally, the differential occurrence of ant attractants (plant and insect exudates) on foliage produces variable levels of ant foraging within local floras and among localities. Here, we investigate how variation of ant communities and of traits among host plant species (presence or absence of ant attractants) can change the effect of carnivores (predatory ants) on herbivore communities (caterpillars) in a cerrado savanna landscape. We sampled caterpillars and foliage-foraging ants in four cerrado localities (70-460 km apart). We found that: (i) caterpillar infestation was negatively related with ant visitation to plants; (ii) this relationship depended on local ant abundance and species composition, and on local preference by ants for plants with liquid attractants; (iii) this was not related to local plant richness or plant size; (iv) the relationship between the presence of ant attractants and caterpillar abundance varied among sites from negative to neutral; and (v) caterpillars feeding on plants with ant attractants are more resistant to ant predation than those feeding on plants lacking attractants. Liquid food on foliage mediates host plant quality for lepidopterans by promoting generalized ant-caterpillar antagonism. Our study in cerrado shows that the negative effects of generalist predatory ants on herbivores are detectable at a community level, affecting patterns of abundance and host plant use by lepidopterans. The magnitude of ant-induced effects on caterpillar occurrence across the cerrado landscape may depend on how ants use plants locally and how they respond to liquid food on plants at different habitats. This study enhances the relevance of plant-ant and ant-herbivore interactions in cerrado and highlights the importance of a tritrophic perspective in this ant-rich environment.

Contrasting Effects Of Land Use Intensity And Exotic Host Plants On The Specialization Of Interactions In Plant-herbivore Networks.

Human land use tends to decrease the diversity of native plant species and facilitate the invasion and establishment of exotic ones. Such changes in land use and plant community composition usually have negative impacts on the assemblages of native herbivorous insects. Highly specialized herbivores are expected to be especially sensitive to land use intensification and the presence of exotic plant species because they are neither capable of consuming alternative plant species of the native flora nor exotic plant species. Therefore, higher levels of land use intensity might reduce the proportion of highly specialized herbivores, which ultimately would lead to changes in the specialization of interactions in plant-herbivore networks. This study investigates the community-wide effects of land use intensity on the degree of specialization of 72 plant-herbivore networks, including effects mediated by the increase in the proportion of exotic plant species. Contrary to our expectation, the net effect of land use intensity on network specialization was positive. However, this positive effect of land use intensity was partially canceled by an opposite effect of the proportion of exotic plant species on network specialization. When we analyzed networks composed exclusively of endophagous herbivores separately from those composed exclusively of exophagous herbivores, we found that only endophages showed a consistent change in network specialization at higher land use levels. Altogether, these results indicate that land use intensity is an important ecological driver of network specialization, by way of reducing the local host range of herbivore guilds with highly specialized feeding habits. However, because the effect of land use intensity is offset by an opposite effect owing to the proportion of exotic host species, the net effect of land use in a given herbivore assemblage will likely depend on the extent of the replacement of native host species with exotic ones.

Graphene Oxide: A Carrier For Pharmaceuticals And A Scaffold For Cell Interactions.

During the last ten years, graphene oxide has been explored in many applications due to its remarkable electroconductivity, thermal properties and mobility of charge carriers, among other properties. As discussed in this review, the literature suggests that a total characterization of graphene oxide must be conducted because oxidation debris (synthesis impurities) present in the graphene oxides could act as a graphene oxide surfactant, stabilizing aqueous dispersions. It is also important to note that the structure models of graphene oxide need to be revisited because of significant implications for its chemical composition and its direct covalent functionalization. Another aspect that is discussed is the need to consider graphene oxide surface chemistry. The hemolysis assay is recommended as a reliable test for the preliminary assessment of graphene oxide toxicity, biocompatibility and cell membrane interaction. More recently, graphene oxide has been extensively explored for drug delivery applications. An important increase in research efforts in this emerging field is clearly represented by the hundreds of related publications per year, including some reviews. Many studies have been performed to explore the graphene oxide properties that enable it to deliver more than one activity simultaneously and to combine multidrug systems with photothermal therapy, indicating that graphene oxide is an attractive tool to overcome hurdles in cancer therapies. Some strategic aspects of the application of these materials in cancer treatment are also discussed. In vitro studies have indicated that graphene oxide can also promote stem cell adhesion, growth and differentiation, and this review discusses the recent and pertinent findings regarding graphene oxide as a valuable nanomaterial for stem cell research in medicine. The protein corona is a key concept in nanomedicine and nanotoxicology because it provides a biomolecular identity for nanomaterials in a biological environment. Understanding protein corona-nanomaterial interactions and their influence on cellular responses is a challenging task at the nanobiointerface. New aspects and developments in this area are discussed.

Interplay Of Weak Interactions In The Atom-by-atom Condensation Of Xenon Within Quantum Boxes.

Condensation processes are of key importance in nature and play a fundamental role in chemistry and physics. Owing to size effects at the nanoscale, it is conceptually desired to experimentally probe the dependence of condensate structure on the number of constituents one by one. Here we present an approach to study a condensation process atom-by-atom with the scanning tunnelling microscope, which provides a direct real-space access with atomic precision to the aggregates formed in atomically defined 'quantum boxes'. Our analysis reveals the subtle interplay of competing directional and nondirectional interactions in the emergence of structure and provides unprecedented input for the structural comparison with quantum mechanical models. This approach focuses on-but is not limited to-the model case of xenon condensation and goes significantly beyond the well-established statistical size analysis of clusters in atomic or molecular beams by mass spectrometry.

Alcohol dehydrogenase activities and ethanol tolerance in Anastrepha (Diptera, Tephritidae) fruit-fly species and their hybrids

The ADH (alcohol dehydrogenase) system is one of the earliest known models of molecular evolution, and is still the most studied in Drosophila. Herein, we studied this model in the genus Anastrepha (Diptera, Tephritidae). Due to the remarkable advantages it presents, it is possible to cross species with different Adh genotypes and with different phenotype traits related to ethanol tolerance. The two species studied here each have a different number of Adh gene copies, whereby crosses generate polymorphisms in gene number and in composition of the genetic background. We measured certain traits related to ethanol metabolism and tolerance. ADH specific enzyme activity presented gene by environment interactions, and the larval protein content showed an additive pattern of inheritance, whilst ADH enzyme activity per larva presented a complex behavior that may be explained by epistatic effects. Regression models suggest that there are heritable factors acting on ethanol tolerance, which may be related to enzymatic activity of the ADHs and to larval mass, although a pronounced environmental effect on ethanol tolerance was also observed. By using these data, we speculated on the mechanisms of ethanol tolerance and its inheritance as well as of associated traits.

Effect of zinc insertion and hydrophobicity on the membrane interactions and PDT activity of porphyrin photosensitizers

A series of photosensitizers (PS), which are meso-substituted tetra-cationic porphyrins, was synthesized in order to study the role of amphiphilicity and zinc insertion in photodynamic therapy (PDT) efficacy. Several properties of the PS were evaluated and compared within the series including photophysical properties (absorption spectra, fluorescence quantum yield Phi(f), and singlet oxygen quantum yield Phi(Delta)), uptake by vesicles, mitochondria and HeLa cells, dark and phototoxicity in HeLa cells. The photophysical properties of all compounds are quite similar (Phi(f) <= 0.02; Phi(Delta) similar to 0.8). An increase in lipophilicity and the presence of zinc in the porphyrin ring result in higher vesicle and cell uptake. Binding in mitochondria is dependent on the PS lipophilicity and on the electrochemical membrane potential, i.e., in uncoupled mitochondria PS binding decreases by up to 53%. The porphyrin substituted with octyl groups (TC8PyP) is the compound that is most enriched in mitochondria, and its zinc derivative (ZnTC8PyP) has the highest global uptake. The stronger membrane interaction of the zinc-substituted porphyrins is attributed to a complexing effect with phosphate groups of the phospholipids. Zinc insertion was also shown to decrease the interaction with isolated mitochondria and with the mitochondria of HeLa cells, an effect that has been explained by the particular characteristics of the mitochondrial internal membrane. Phototoxicity was shown to increase proportionally with membrane binding efficiency, which is attributed to favorable membrane interactions which allow more efficient membrane photooxidation. For this series of compounds, photodynamic efficiency is directly proportional to the membrane binding and cell uptake, but it is not totally related to mitochondrial targeting.