Recovering Power in Association Mapping Panels with Variable Levels of Linkage Disequilibrium

Association mapping has permitted the discovery of major QTL in many species. It can be applied to existing populations and, as a consequence, it is generally necessary to take into account structure and relatedness among individuals in the statistical model to control false positives. We analytically studied power in association studies by computing noncentrality parameter of the tests and its relationship with parameters characterizing diversity (genetic differentiation between groups and allele frequencies) and kinship between individuals. Investigation of three different maize diversity panels genotyped with the 50k SNPs array highlighted contrasted average power among panels and revealed gaps of power of classical mixed models in regions with high linkage disequilibrium (LD). These gaps could be related to the fact that markers are used for both testing association and estimating relatedness. We thus considered two alternative approaches to estimating the kinship matrix to recover power in regions of high LD. In the first one, we estimated the kinship with all the markers that are not located on the same chromosome than the tested SNP. In the second one, correlation between markers was taken into account to weight the contribution of each marker to the kinship. Simulations revealed that these two approaches were efficient to control false positives and were more powerful than classical models.     

Resequencing the susceptibility gene,ITGAM, identifies two functionally deleterious rare variants in systemic lupus erythematosus cases

Introduction

The majority of the genetic variance of systemic lupus erythematosus (SLE) remains unexplained by the common disease-common variant hypothesis. Rare variants, which are not detectable by genome-wide association studies because of their low frequencies, are predicted to explain part of this ”missing heritability.” However, recent studies identifying rare variants within known disease-susceptibility loci have failed to show genetic associations because of their extremely low frequencies, leading to the questioning of the contribution of rare variants to disease susceptibility. A common (minor allele frequency = 17.4% in cases) nonsynonymous coding variant rs1143679 (R77H) in ITGAM (CD11b), which forms half of the heterodimeric integrin receptor, complement receptor 3 (CR3), is robustly associated with SLE and has been shown to impair CR3-mediated phagocytosis.

Methods

We resequenced ITGAM in 73 SLE cases and identified two previously unidentified, case-specific nonsynonymous variants, F941V and G1145S. Both variants were genotyped in 2,107 and 949 additional SLE cases, respectively, to estimate their frequencies in a disease population. An in vitro model was used to assess the impact of F941V and G1145S, together with two nonsynonymous ITGAM polymorphisms, A858V (rs1143683) and M441T (rs11861251), on CR3-mediated phagocytosis. A paired two-tailed t test was used to compare the phagocytic capabilities of each variant with that of wild-type CR3.

Results

Both rare variants, F941V and G1145S, significantly impair CR3-mediated phagocytosis in an in vitro model (61% reduction, P = 0.006; 26% reduction, P = 0.0232). However, neither of the common variants, M441T and A858V, had an effect on phagocytosis. Neither rare variant was observed again in the genotyping of additional SLE cases, suggesting that there frequencies are extremely low.

Conclusions

Our results add further evidence to the functional importance of ITGAM in SLE pathogenesis through impaired phagocytosis. Additionally, this study provides a new example of the identification of rare variants in common-allele-associated loci, which, because of their extremely low frequencies, are not statistically associated. However, the demonstration of their functional effects adds support to their contribution to disease risk, and questions the current notion of dismissing the contribution of very rare variants on purely statistical analyses.     

Analysis of a genome-wide association study-linked locus (CCR6) in Asian rheumatoid arthritis

A genome-wide association study in Japan identified the C-C chemokine receptor type 6 gene (CCR6) as associated with rheumatoid arthritis (RA). This finding has not been validated in other Asian populations. A case-control study involving 996 subjects, comprising 440 controls and 556 RA patients, was done to determine their anticyclic citrullinated peptide (anti-CCP) antibody status and CCR6 polymorphism (rs3093024) genotype. Three hundred eighty-seven patients were anti-CCP positive and 153 anti-CCP negative. Logistic regression showed that allele A was likely to increase the risk of developing RA among females via a recessive model (odds ratio [OR]=1.55, 95% confidence interval [CI]=1.01, 2.39), whereas the risk effect appeared to be reduced among males via an additive model (OR=0.60, 95% CI=0.42, 0.85). Considering only subjects who are anti-CCP positive, allele A increased RA risk among females via a recessive model (OR=1.68, 95% CI=1.07, 2.64) but decreased the risk among males via an additive model (OR=0.59, 95% CI=0.39, 0.89). We showed that CCR6 polymorphism was a risk factor among females but a protective factor among males. Functional studies are warranted to unravel the pathophysiological relevance of the gene variant and other linked variants with RA.     

Transgene-mediated and elicitor-induced perturbation of metabolic channeling at the entry point into the phenylpropanoid pathway

3H-l-Phenylalanine is incorporated into a range of phenylpropanoid compounds when fed to tobacco cell cultures. A significant proportion of (3)H-trans-cinnamic acid formed from (3)H-l-phenylalanine did not equilibrate with exogenous trans-cinnamic acid and therefore may be rapidly channeled through the cinnamate 4-hydroxylase (C4H) reaction to 4-coumaric acid. Such compartmentalization of trans-cinnamic acid was not observed after elicitation or in cell cultures constitutively expressing a bean phenylalanine ammonia-lyase (PAL) transgene. Channeling between PAL and C4H was confirmed in vitro in isolated microsomes from tobacco stems or cell suspension cultures. This channeling was strongly reduced in microsomes from stems or cell cultures of transgenic PAL-overexpressing plants or after elicitation of wild-type cell cultures. Protein gel blot analysis showed that tobacco PAL1 and bean PAL were localized in both soluble and microsomal fractions, whereas tobacco PAL2 was found only in the soluble fraction. We propose that metabolic channeling of trans-cinnamic acid requires the close association of specific forms of PAL with C4H on microsomal membranes.     

Effect of variety and location on growth and yield components of Roselle,HIbiscus sabdariffa L. and its infestation with the spiny bollworm Earias insulana (BOISD.)

Three roselle, Hibiscus sabdariffa L. varieties (Sudani, Masri and White) were cultivated at three different locations to recognize the transportation ability of roselle cultivation from the narrow old valley land to broad new land in Egypt. Qena as origin in situ old land, El-Kanater as ex situ old land and Nubaria as ex situ new land were the considered locations. Six growth quantitative characters and bolls infestation by spiny bollworm, Earias insulana were evaluated. Growth characters of roselle plants were affected significantly by either variety or location. Qena region was more suitable for roselle plant growth as judged with plant height, number of branches, number of fruits and sepals dry weight, followed by Nubaria followed by El-Kanater. Whereas, plants grown at Nubaria produced more fresh sepals weight than Qena or El-Kanater grown plants. As for Sudani, Nubaria exhibited the tallest plants, with the highest number of fruits and the heaviest fresh sepals as compared with the corresponding plants in Qena or El-Kanater. Values of broad sense heritability were highest for all characters in Qena. While the number of fruits per plant had the highest heritability in all locations. Dry sepals yield had highly significant correlation with all studied characters except percentage of water loss in Qena and Nubaria. Path coefficient analysis confirmed that fresh sepals yield had the highest direct and indirect effects on dried sepals yield. Chemical constituents responsible to sepal quality tended to produce significant variations due to the changes in varieties or locations. The highest levels of anthocyanins and sugars were achieved by Sudani variety, but the highest levels of free amino acids and total soluble solids were recorded for Masri variety. Moreover, Nubaria region was the most favourable for the accumulation of more anthocyanins in the sepals of all varieties followed by Qena. Plants grown at Qena produced sepals with the highest levels of sugars, free amino acids, organic acids and total soluble solids, followed by Nubaria followed by El-Kanater plants. Infestation with spiny bollworm Earias insulana was increased from Sudani up to Masri up to White varieties. Plants grown at Nubaria had the lowest number of attacks by bolls in all varieties, followed by those at El-Kanater followed by Qena plants. Spiny bollworm infestation was positively correlated with the number of branches and dry sepals weight, but negatively correlated with sepal moisture loss and anthocyanin contents. These findings clearly indicated that the Nubaria region was considered as a promising reclaimed area suitable for roselle cultivation, especially for Sudani, the most economic variety.     

Quantitative analysis of sugar constituents of glycoproteins by capillary electrophoresis

A method for quantitative analysis of monosaccharides including N- acetylneuraminic acid derived from sialic acid-containing oligosaccharides and glycoproteins is presented. The analysis is based on the combination of chemical and enzymatic methods coupled with capillary electrophoretic (CE) separation and laser-induced fluorescence (LIF) detection. The present method utilizes a simplified acid hydrolysis procedure consisting of mild hydrolysis (0.1 M TFA) to release sialic acid and strong acid hydrolysis (2.0 N TFA) to produce amino and neutral sugars. Amino sugars released from strong acid hydrolysis of oligosaccharides and glycoproteins were reacetylated and derivatized with 8-aminopyrene-1,3,6-trisulfonate (APTS) along with neutral sugars in the presence of sodium cyanoborohydride to yield quantitatively the highly stable fluorescent APTS adducts. N- acetylneuraminic acid (Neu5Ac), a major component of most mammalian glycoproteins, was converted in a fast specific reaction by the action of neuraminic acid aldolase (N-acylneuraminate pyruvate-lyase EC 4.1.3.3) to N-acetylmannosamine (ManNAc) and pyruvate. ManNAc was then derivatized with APTS in the same manner as the other monosaccharides. This method was demonstrated for the quantitation of pure Neu5Ac and the species derived from mild acid hydrolysis of 6'-sialyl-N- acetyllactosamine and bovine fetuin glycan. Quantitative recovery of the N-acetylmannosamine was obtained from a known amount of Neu5Ac in a mixture of seven other monosaccharides or from the sialylated oligosaccharides occurring in glycoproteins. The sequence of procedures consists of acid hydrolysis, enzymatic conversion and APTS derivatization which produced quantitative recovery of APTS- monosaccharide adducts. The detection limits for sugars derivatized with APTS and detected by CE-LIF are 100 pmol for Neu5Ac and 50 pmol for the other sugars.      

Role of the reticulum in the stability and shape of the isolated human erythrocyte membrane

In order to examine the widely held hypothesis that the reticulum of proteins which covers the cytoplamsic surface of the human erythrocyte membrane controls cell stability and shape, we have assessed some of its properties. The reticulum, freed of the bilayer by extraction with Triton X-100, was found to be mechanically stable at physiological ionic strength but physically unstable at low ionic strength. The reticulum broke down after a characteristic lag period which decreased 500-fold between 0 degrees and 37 degrees C. The release of polypeptide band 4.1 from the reticulum preceded that of spectrin and actin, suggesting that band 4.1 might stabilize the ensemble but is not essential to its integrity. The time-course of breakdown was similar for ghosts, the reticulum inside of ghosts, and the isolated reticulum. However, at very low ionic strength, the reticulum was less stable within the ghost than when free; at higher ionic strength, the reverse was true. Over a wide range of conditions the membrane broke down to vesicles just as the reticulum disintegrated, presumably because the bilayer was mechanically stabilized by this network. The volume of both ghosts and naked reticula varied inversely and reversibly with ionic strength. The volume of the naked reticulum varied far more widely than the ghost, suggesting that its deformation was normally limited by the less extensible bilayer. The contour of the isolated reticulum was discoid and often dimpled or indented, as visualized in the fluorescence microscope after labeling of the ghosts with fluoroscein isothiocyanate. Reticula derived from ghosts which had lost the ability to crenate in isotonic saline were shriveled, even though the bilayer was smooth and expanded. Conversly, ghosts crenated by dinitrophenol yielded smooth, expanded reticula. We conclude that the reticulum is a durable, flexible, and elastic network which assumes and stabilizes the contour of the membrane but is not responsible for its crenation.