Nonadditive Effects of Genes in Human Metabolomics

Genome-wide association studies (GWAS) are widely applied to analyze the genetic effects on phenotypes. With the availability of high-throughput technologies for metabolite measurements, GWAS successfully identified loci that affect metabolite concentrations and underlying pathways. In most GWAS, the effect of each SNP on the phenotype is assumed to be additive. Other genetic models such as recessive, dominant, or overdominant were considered only by very few studies. In contrast to this, there are theories that emphasize the relevance of nonadditive effects as a consequence of physiologic mechanisms. This might be especially important for metabolites because these intermediate phenotypes are closer to the underlying pathways than other traits or diseases. In this study we analyzed systematically nonadditive effects on a large panel of serum metabolites and all possible ratios (22,801 total) in a population-based study [Cooperative Health Research in the Region of Augsburg (KORA) F4, N = 1,785]. We applied four different 1-degree-of-freedom (1-df) tests corresponding to an additive, dominant, recessive, and overdominant trait model as well as a genotypic model with two degree-of-freedom (2-df) that allows a more general consideration of genetic effects. Twenty-three loci were found to be genome-wide significantly associated (Bonferroni corrected P ≤ 2.19 × 10⁻¹²) with at least one metabolite or ratio. For five of them, we show the evidence of nonadditive effects. We replicated 17 loci, including 3 loci with nonadditive effects, in an independent study (TwinsUK, N = 846). In conclusion, we found that most genetic effects on metabolite concentrations and ratios were indeed additive, which verifies the practice of using the additive model for analyzing SNP effects on metabolites.     

Detecting low frequent loss-of-function alleles in genome wide association studies with red hair color as example

Multiple loss-of-function (LOF) alleles at the same gene may influence a phenotype not only in the homozygote state when alleles are considered individually, but also in the compound heterozygote (CH) state. Such LOF alleles typically have low frequencies and moderate to large effects. Detecting such variants is of interest to the genetics community, and relevant statistical methods for detecting and quantifying their effects are sorely needed. We present a collapsed double heterozygosity (CDH) test to detect the presence of multiple LOF alleles at a gene. When causal SNPs are available, which may be the case in next generation genome sequencing studies, this CDH test has overwhelmingly higher power than single SNP analysis. When causal SNPs are not directly available such as in current GWA settings, we show the CDH test has higher power than standard single SNP analysis if tagging SNPs are in linkage disequilibrium with the underlying causal SNPs to at least a moderate degree (r²>0.1). The test is implemented for genome-wide analysis in the publically available software package GenABEL which is based on a sliding window approach. We provide the proof of principle by conducting a genome-wide CDH analysis of red hair color, a trait known to be influenced by multiple loss-of-function alleles, in a total of 7,732 Dutch individuals with hair color ascertained. The association signals at the MC1R gene locus from CDH were uniformly more significant than traditional GWA analyses (the most significant P for CDH = 3.11×10⁻¹⁴² vs. P for rs258322 = 1.33×10⁻⁶⁶). The CDH test will contribute towards finding rare LOF variants in GWAS and sequencing studies.     

Increase of histidine content in Brassica rapa subsp. oleifera by over-expression of histidine-rich fusion proteins

An approach that enables the increase of the quantity of a specific amino acid in crop plants is reported. Oleosin gene from Arabidopsis thaliana or 30K movement protein gene of Tobacco mosaic virus (TMV; genus Tobamovirus) were cloned under the control of napin or hybrid promoters, and in fusion to synthetic poly-histidine (poly-His) sequences for transformation into spring turnip rape (Brassica rapa subsp. oleifera; synonym to B. campestris). The most stable expression cassettes for the poly-His production prior to the plant transformation were selected by analyzing the protein expression in in vitro translation and in transient plant expression systems using GFP as marker. Expression of the poly-His-constructs in transgenic Brassica rapa plants was analyzed using dot and western blotting and PCR. The constructs were stably expressed in the third generation of the transgenic plant lines. Histidine content was measured from the seeds of the transgenic plants, and some plant lines had more than 20% increase in histidine content compared to wild type. The methodology may be widely applicable to increase the content of any amino acid in crop plants including those encoded by rare codons.     

Effective size of subpopulations in early-run sockeye salmon Oncorhynchus nerka from Azabach'e Lake (Kamchatka): the effect of relative reproductive success of different-year cohorts

The effect of variation in reproductive success of cohorts of different year of birth (within generation) on the effective subpopulation (breeding group) size in early-run sockeye salmon Oncorhynchus nerka from Azabach'e Lake (Kamchatka). The annual variation in census size and overlapping of year classes reduced the ratio of the effective subpopulation size to the census size by 7 to 88% in different subpopulations. The total effect of the variance of reproductive success in individual years and the variance of reproductive success of different cohorts reduced the effective size/census size ratio by 68-96%.     

Influence of human B-casomorphin-7 on specific binding of 3H-spiperone to the 5-HT2-receptors of rat brain frontal cortex

Using radio-receptor analysis, it has been demonstrated that human beta-casomorphin-7 (Tyr-Pro-Phe-Val-Glu-Pro-Ile) displaces 3H-spiperone from 5-HT2-receptors of rat brain frontal cortex. IC50 of human beta-casomorphin-7 was 8 microM. These data suggest that one of the mechanisms of neurotropic action of beta-casomorphin-7 is might be associated with its influence on the serotoninergic system.     

The rule of "one migrant per generation" and genetic differentiation in subdivided populations

The genetic differentiation in population with migration according island and two-dimensional stepping stone models was studied with simulation methods. It is shown that migration of one or few individuals per generation is insufficient for leveling differences between subpopulations in allele frequencies. Even if migration is estimated as m = 0.5 (exchange of 50 and 500 individuals per generation) statistically significant differences remain at least in the half of populations with insular structure. Spatial heterogeneity disappears completely only if m = 0.7-0.8. In case of two-dimensional step model the level of genetic differentiation is higher and statistically significant heterogeneity remains at all levels of genetic exchange including that which was estimated as m = 1.     

Time-resolved microspectroscopy on a single crystal of bacteriorhodopsin reveals lattice-induced differences in the photocycle kinetics

The determination of the intermediate state structures of the bacteriorhodopsin photocycle has lead to an unprecedented level of understanding of the catalytic process exerted by a membrane protein. However, the crystallographic structures of the intermediate states are only relevant if the working cycle is not impaired by the crystal lattice. Therefore, we applied visible and Fourier transform infrared spectroscopy (FTIR) microspectroscopy with microsecond time resolution to compare the photoreaction of a single bacteriorhodopsin crystal to that of bacteriorhodopsin residing in the native purple membrane. The analysis of the FTIR difference spectra of the resolved intermediate states reveals great similarity in structural changes taking place in the crystal and in PM. However, the kinetics of the photocycle are significantly altered in the three-dimensional crystal as compared to PM. Strikingly, the L state decay is accelerated in the crystal, whereas the M decay is delayed. The physical origin of this deviation and the implications for trapping of intermediate states are discussed. As a methodological advance, time-resolved step-scan FTIR spectroscopy on a single protein crystal is demonstrated for the first time which may be used in the future to gauge the functionality of other crystallized proteins with the molecular resolution of vibrational spectroscopy.     

Caldesmon restricts the movement of both C- and N-termini of tropomyosin on F-actin in ghost fibers during the actomyosin ATPase cycle

New data on the movements of tropomyosin singly labeled at alpha- or beta-chain during the ATP hydrolysis cycle in reconstituted ghost fibers have been obtained by using the polarized fluorescence technique which allowed us following the azimuthal movements of tropomyosin on actin filaments. Pronounced structural changes in tropomyosin evoked by myosin heads suggested the "rolling" of the tropomyosin molecule on F-actin surface during the ATP hydrolysis cycle. The movements of actin-bound tropomyosin correlated to the strength of S1 to actin binding. Weak binding of myosin to actin led to an increase in the affinity of the tropomyosin N-terminus to actin with simultaneous decrease in the affinity of the C-terminus. On the contrary, strong binding of myosin to actin resulted in the opposite changes of the affinity to actin of both ends of the tropomyosin molecule. Caldesmon inhibited the "rolling" of tropomyosin on the surface of the thin filament during the ATP hydrolysis cycle, drastically decreased the affinity of the whole tropomyosin molecule to actin, and "freezed" tropomyosin in the position characteristic of the weak binding of myosin to actin.     

Proteomic analysis of infiltrating ductal carcinoma tissues by coupled 2-D DIGE/MS/MS analysis

There is a growing interest in protein expression profiling aiming to identify novel diagnostic markers in breast cancer. Proteomic approaches such as two-dimensional differential gel electrophoresis coupled with tandem mass spectrometry analysis (2-D DIGE/MS/MS) have been used successfully for the identification of candidate biomarkers for screening, diagnosis, prognosis and monitoring of treatment response in various types of cancer. Identifying previously unknown proteins of potential clinical relevance will ultimately help in reaching effective ways to manage the disease. We analyzed breast cancer tissues from five tumor and five normal tissue samples from ten breast cancer subjects with infiltrating ductal carcinoma (IDC) by 2-D DIGE using two types of immobilized pH gradient (IPG) strips: pH 3-10 and pH 4-7. From all the spots detected, differentially expressed (p < 0.05 and ratio > 2) were 50 spots. Of these, 39 proteins were successfully identified by MS, representing 29 different proteins. Ten proteins were overexpressed in the tumor samples. The 2-D DIGE/MS/MS analysis revealed an increase in the expression levels in tumor samples of several proteins not previously associated with breast cancer, such as: macrophage-capping protein (CAPG), phosphomannomutase 2 (PMM2), ATPase ASN1, methylthioribose-1-phosphate isomerase (MRI1), peptidyl-prolyl cis-trans isomerase FKBP4, cellular retinoic acid-binding protein 2 (CRABP2), lamin B1 and keratin, type II cytoskeletal 8 (KRT8). Ingenuity Pathway Analysis (IPA) revealed highly significant (p = 10(-26)) interactions between the identified proteins and their association with cancer. These proteins are involved in many diverse pathways and have established roles in cellular metabolism. It remains the goal of future work to test the suitability of the identified proteins in samples of larger and independent patient groups.