Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Quantitative high-resolution warm season rainfall recorded in varved sediments of Lake Oeschinen, northern Swiss Alps: calibration and validation AD 1901–2008

High-resolution, well-calibrated records of lake sediments are critically important for quantitative climate reconstructions, but they remain a methodological and analytical challenge. While several comprehensive paleotemperature reconstructions have been developed across Europe, only a few quantitative high-resolution studies exist for precipitation. Here we present a calibration and verification study of lithoclastic sediment proxies from proglacial Lake Oeschinen (46°30′N, 7°44′E, 1,580 m a.s.l., north–west Swiss Alps) that are sensitive to rainfall for the period AD 1901–2008. We collected two sediment cores, one in 2007 and another in 2011. The sediments are characterized by two facies: (A) mm-laminated clastic varves and (B) turbidites. The annual character of the laminae couplets was confirmed by radiometric dating (210Pb, 137Cs) and independent flood-layer chronomarkers. Individual varves consist of a dark sand-size spring-summer layer enriched in siliciclastic minerals and a lighter clay-size calcite-rich winter layer. Three subtypes of varves are distinguished: Type I with a 1–1.5 mm fining upward sequence; Type II with a distinct fine-sand base up to 3 mm thick; and Type III containing multiple internal microlaminae caused by individual summer rainstorm deposits. Delta-fan surface samples and sediment trap data fingerprint different sediment source areas and transport processes from the watershed and confirm the instant response of sediment flux to rainfall and erosion. Based on a highly accurate, precise and reproducible chronology, we demonstrate that sediment accumulation (varve thickness) is a quantitative predictor for cumulative boreal alpine spring (May–June) and spring/summer (May–August) rainfall (rMJ = 0.71, rMJJA = 0.60, p < 0.01). Bootstrap-based verification of the calibration model reveals a root mean squared error of prediction (RMSEPMJ = 32.7 mm, RMSEPMJJA = 57.8 mm) which is on the order of 10–13 % of mean MJ and MJJA cumulative precipitation, respectively. These results highlight the potential of the Lake Oeschinen sediments for high-resolution reconstructions of past rainfall conditions in the northern Swiss Alps, central and eastern France and south-west Germany.

Co-activation of JAK/STAT and p38 MAPK pathways by type I interferon enhances apoptosis in human basophils

It has previously been published that interferon-α (type I IFN) improves clinical symptoms of asthma patients. Since human basophils are major inflammatory cells in maintaining chronic allergic asthma we investigate whether type I IFN affect human blood basophils. Furthermore, previous studies have shown that spontaneous apoptosis of human basophils is slow due to constitutive expression of anti-apoptotic BCL-2 family members. In addition, IL-3 exceptionally promotes survival of basophils by enhancing constitutive expression of BCL-2 family members and by inducing de-novo expression of Pim-1 kinase. Thus, we also assessed whether type I IFN might overcome IL-3-induced survival of human basophils. Our data show that type I IFN enhances apoptosis in purified human blood basophils compared to spontaneous apoptosis of controls or type II IFN treated cells. Furthermore, we demonstrate that both type I IFN and FasL enhance apoptosis in human basophils with similar efficiency in a rather additive than synergistic way. Analyses of signaling pathways reveal that type I IFN promote prolonged phosphorylation of STAT1/STAT2. By using a pan-JAK inhibitor the phosphorylation of STAT1/STAT2 is inhibited and most importantly the pro-apoptotic effect of type I IFN is abolished. On the other hand, type I IFN do not reduce IL-3-induced de novo expression of Pim-1 and BCL-2. This is in line with our observation that IL-3-induced survival is dominant over type I IFN-enhanced apoptosis. In addition, phosphorylation of p38 MAPK in type I IFN treated cells is comparable to non-treated cells. Particularly however, inhibition of this p-p38 activity abrogates apoptosis as well. We conclude that type I IFN-enhanced apoptosis is tightly regulated by the cooperation of JAK/STAT and p38 MAPK pathways. Our study identifies a so far unknown effect of type I IFN and may explain the improved clinical symptoms of asthma patients treated with type I IFN.

TGF-beta1 in Aplysia: role in long-term changes in the excitability of sensory neurons and distribution of TbetaR-II-like immunoreactivity.

Exogenous recombinant human transforming growth factor beta-1 (TGF-beta1) induced long-term facilitation of Aplysia sensory-motor synapses. In addition, 5-HT-induced facilitation was blocked by application of a soluble fragment of the extracellular portion of the TGF-beta1 type II receptor (TbetaR-II), which presumably acted by scavenging an endogenous TGF-beta1-like molecule. Because TbetaR-II is essential for transmembrane signaling by TGF-beta, we sought to determine whether Aplysia tissues contained TbetaR-II and specifically, whether neurons expressed the receptor. Western blot analysis of Aplysia tissue extracts demonstrated the presence of a TbetaR-II-immunoreactive protein in several tissue types. The expression and distribution of TbetaR-II-immunoreactive proteins in the central nervous system was examined by immunohistochemistry to elucidate sites that may be responsive to TGF-beta1 and thus may play a role in synaptic plasticity. Sensory neurons in the ventral-caudal cluster of the pleural ganglion were immunoreactive for TbetaR-II, as well as many neurons in the pedal, abdominal, buccal, and cerebral ganglia. Sensory neurons cultured in isolation and cocultured sensory and motor neurons were also immunoreactive. TGF-beta1 affected the biophysical properties of cultured sensory neurons, inducing an increase of excitability that persisted for at least 48 hr. Furthermore, exposure to TGF-beta1 resulted in a reduction in the firing threshold of sensory neurons. These results provide further support for the hypothesis that TGF-beta1 plays a role in long-term synaptic plasticity in Aplysia.

Molecular mechanisms of chondrocyte-specific expression on the mouse pro$\alpha$1(II) collagen gene

Type II collagen is a major chondrocyte-specific component of the cartilage extracellular matrix and it represents a typical differentiation marker of mature chondrocytes. In order to delineate cis-acting elements of the mouse pro$\alpha1$(II) collagen gene that control chondrocyte-specific expression in intact mouse embryos, we generated transgenic mice harboring chimeric constructions in which varying lengths of the promoter and intron 1 sequences were linked to a $\beta$-galactosidase reporter gene. A construction containing a 3000-bp promoter and a 3020-bp intron 1 fragment directed high levels of $\beta$-galactosidase expression specifically to chondrocytes. Successive deletions of intron 1 delineated a 48-bp fragment which targeted $\beta$-galactosidase expression to chondrocytes with the same specificity as the larger intron 1 fragment. When the Col2a1 promoter was replaced with a minimal $\beta$-globin promoter, the 48-bp intron 1 sequence was still able to target expression of the transgene to chondrocytes, specifically. Therefore a 48-bp intron 1 DNA segment of the mouse Col2a1 gene contains the necessary information to confer high-level, temporally correct, chondrocyte expression to a reporter gene in intact mouse embryos and that Col2a1 promoter sequences are dispensable for chondrocyte expression. Nuclear proteins present selectively in mouse primary chondrocytes and rat chondrosarcoma cells bind to the three putative HMG (High-Mobility-Group) domain protein binding sites in this 48-bp sequence and the chondrocyte-specific proteins likely bind the DNA through minor groove. Together, my results indicate that a 48-bp sequence in Col2a1 intron 1 controls chondrocyte-specific expression in vivo and suggest that chondrocytes contain specific nuclear proteins involved in enhancer activity. ^

The pre- and post-somatic segments of the human type I spiral ganglion neurons--structural and functional considerations related to cochlear implantation.

Human auditory nerve afferents consist of two separate systems; one is represented by the large type I cells innervating the inner hair cells and the other one by the small type II cells innervating the outer hair cells. Type I spiral ganglion neurons (SGNs) constitute 96% of the afferent nerve population and, in contrast to other mammals, their soma and pre- and post-somatic segments are unmyelinated. Type II nerve soma and fibers are unmyelinated. Histopathology and clinical experience imply that human SGNs can persist electrically excitable without dendrites, thus lacking connection to the organ of Corti. The biological background to this phenomenon remains elusive. We analyzed the pre- and post-somatic segments of the type I human SGNs using immunohistochemistry and transmission electron microscopy (TEM) in normal and pathological conditions. These segments were found surrounded by non-myelinated Schwann cells (NMSCs) showing strong intracellular expression of laminin-β2/collagen IV. These cells also bordered the perikaryal entry zone and disclosed surface rugosities outlined by a folded basement membrane (BM) expressing laminin-β2 and collagen IV. It is presumed that human large SGNs are demarcated by three cell categories: (a) myelinated Schwann cells, (b) NMSCs and (c) satellite glial cells (SGCs). Their BMs express laminin-β2/collagen IV and reaches the BM of the sensory epithelium at the habenula perforata. We speculate that the NMSCs protect SGNs from further degeneration following dendrite loss. It may give further explanation why SGNs can persist as electrically excitable monopolar cells even after long-time deafness, a blessing for the deaf treated with cochlear implantation.

Analysis of BMP signalling through the receptor type IA in embryogenesis using conditional gene inactivation in mice

Although bone morphogenetic proteins (BMPs) were initially identified for their potent bone-inducing activity, their precise roles in processes of endochondral and intramembranous bone formation are far from being clear. Tissue-specific loss-of-function experiments using the BMP receptor type IA (BMPR-IA) are particularly attractive since this receptor is thought to be essential for signaling by the closely related BMPs -2, 4, and 7. To ablate signaling through this receptor during chondrogenesis, we have generated transgenic mice expressing Cre recombinase under the control of the collagen type II (Col2a1) gene regulatory sequences. Mice lacking BMPR-IA function in chondrocytes display a number of skeletal abnormalities, including defects in bones of the chondrocranium, abnormal dorsal vertebral processes, scapulae with severe hypoplasia of dorsal elements, and shortening of the long bones. Alterations in the growth plate of long bones in mutants suggest that BMPR-IA is not required for early steps of the chondrocyte specification, but is rather important in regulation of terminal differentiation. Molecular analysis revealed noticeable downregulation of the Ihh/Ptch signalling pathway, decreased chondrocyte proliferation rate and deregulation of hypertrophy. ^ In order to elucidate the role of BMP signalling in development of the limb and intramembranous ossification, we have used mice expressing Cre recombinase under control of the Prx1 (MHox) regulatory elements (M. Logan, pers comm.). Cre activity was found in those mice in the developing limb bud mesenchyme, as well as in a subset of cranial neural crest cells. Prx1-Cre-induced conditional mutants display prominent defects in distal limb outgrowth, as well as ossification defects in a number of neural crest-derived calvarial bones. Intriguingly, mutant limbs displayed alterations in patterning along all three axes. Molecular analysis revealed ectopic anterior Shh/Ptch signalling pathway activation and expression of some Hox genes. Observed loss of Msx1 and Msx2 expression in the progress zone correlates with downregulation of Cyclin D1 and decreased distal outgrowth. Abnormal ventral localization of Lmx1b-expressing cells along with observed later morphological abnormalities suggest a novel role for BMP signalling in establishment or maintaining of the dorso-ventral polarity in the limb mesoderm. ^

The magnesium isotope composition of diagenetic dolomites from ODP Hole 112-685A and 201-1230A of the Peru Margin

The magnesium isotope composition of diagenetic dolomites and their adjacent pore fluids were studied in a 250 m thick sedimentary section drilled into the Peru Margin during Ocean Drilling Program (ODP) Leg 201 (Site 1230) and Leg 112 (Site 685). Previous studies revealed the presence of two types of dolomite: type I dolomite forms at ~ 6 m below seafloor (mbsf) due to an increase in alkalinity associated with anaerobic methane oxidation, and type II dolomite forms at focused sites below ~ 230 mbsf due to episodic inflow of deep-sourced fluids into an intense methanogenesis zone. The pore fluid delta 26Mg composition becomes progressively enriched in 26Mg with depth from values similar to seawater (i.e. -0.8 per mil, relative to DSM3 Mg reference material) in the top few meters below seafloor (mbsf) to 0.8 ± 0.2 per mil within the sediments located below 100 mbsf. Type I dolomites have a delta 26Mg of -3.5 per mil, and exhibit apparent dolomite-pore fluid fractionation factors of about -2.6 per mil consistent with previous studies of dolomite precipitation from seawater. In contrast, type II dolomites have delta 26Mg values ranging from -2.5 to -3.0 per mil and are up to -3.6 per mil lighter than the modern pore fluid Mg isotope composition. The enrichment of pore fluids in 26Mg and depletion in total Mg concentration below ~ 200 mbsf is likely the result of Mg isotope fractionation during dolomite formation, The 26Mg enrichment of pore fluids in the upper ~ 200 mbsf of the sediment sequence can be attributed to desorption of Mg from clay mineral surfaces. The obtained results indicate that Mg isotopes recorded in the diagenetic carbonate record can distinguish near surface versus deep formed dolomite demonstrating their usefulness as a paleo-diagenetic proxy.

DCT Type-III for Multicarrier Modulation

In this paper we propose the use of Discrete Cosine Transform Type-III (DCT3) for multicarrier modulation. There are two DCT3 (even and odd) and, for each of them, we derive the expressions for both prefix and suffix to be appended into each data symbol to be transmitted. Moreover, DCT3 are closely related to the corresponding inverse DCT Type-II even and odd. Furthermore, we give explicit expressions for the 1-tap per subcarrier equalizers that must be implemented at the receiver to perform the channel equalization in the frequency-domain. As a result, the proposed DCT3-based multicarrier modulator can be used as an alternative to DFT-based systems to perform Orthogonal Frequency-Division Multiplexing or Discrete Multitone Modulation

Accumulation of Major Histocompatibility Complex Class II Molecules in Mast Cell Secretory Granules and Their Release upon Degranulation

To investigate the relationship between major histocompatibility complex (MHC) class II compartments, secretory granules, and secretory lysosomes, we analyzed the localization and fate of MHC class II molecules in mast cells. In bone marrow-derived mast cells, the bulk of MHC class II molecules is contained in two distinct compartments, with features of both lysosomal compartments and secretory granules defined by their protein content and their accessibility to endocytic tracers. Type I granules display internal membrane vesicles and are accessed by exogenous molecules after a time lag of 20 min; type II granules are reached by the endocytic tracer later and possess a serotonin-rich electron-dense core surrounded by a multivesicular domain. In these type I and type II granules, MHC class II molecules, mannose-6-phosphate receptors and lysosomal membrane proteins (lamp1 and lamp2) localize to small intralumenal vesicles. These 60–80-nm vesicles are released along with inflammatory mediators during mast cell degranulation triggered by IgE-antigen complexes. These observations emphasize the intimate connection between the endocytic and secretory pathways in cells of the hematopoietic lineage which allows regulated secretion of the contents of secretory lysosomes, including membrane proteins associated with small vesicles.

Blockade of type β transforming growth factor signaling prevents liver fibrosis and dysfunction in the rat

We eliminated type β transforming growth factor (TGF-β) signaling by adenovirus-mediated local expression of a dominant-negative type II TGF-β receptor (AdCATβ-TR) in the liver of rats treated with dimethylnitrosamine, a model of persistent liver fibrosis. In rats that received a single application of AdCATβ-TR via the portal vein, liver fibrosis as assessed by histology and hydroxyproline content was markedly attenuated. All AdCATβ-TR-treated rats remained alive, and their serum levels of hyaluronic acid and transaminases remained at low levels, whereas all the AdCATβ-TR-untreated rats died of liver dysfunction. The results demonstrate that TGF-β does play a central role in liver fibrogenesis and indicate clearly in a persistent fibrosis model that prevention of fibrosis by anti-TGF-β intervention could be therapeutically useful.

A missense mutation of the Na+ channel αII subunit gene Nav1.2 in a patient with febrile and afebrile seizures causes channel dysfunction

Generalized epilepsy with febrile seizures plus (GEFS+), a clinical subset of febrile seizures (FS), is characterized by frequent episodes beyond 6 years of age (FS+) and various types of subsequent epilepsy. Mutations in β1 and αI-subunit genes of voltage-gated Na+ channels have been associated with GEFS+1 and 2, respectively. Here, we report a mutation resulting in an amino acid exchange (R187W) in the gene encoding the α-subunit of neuronal voltage-gated Na+ channel type II (Nav1.2) in a patient with FS associated with afebrile seizures. The mutation R187W occurring on Arg187, a highly conserved residue among voltage-gated Na+ channels, was not found in 224 alleles of unaffected individuals. Whole-cell patch clamp recordings on human embryonic kidney (HEK) cells expressing a rat wild-type (rNav1.2) and the corresponding mutant channels showed that the mutant channel inactivated more slowly than wild-type whereas the Na+ channel conductance was not affected. Prolonged residence in the open state of the R187W mutant channel may augment Na+ influx and thereby underlie the neuronal hyperexcitability that induces seizure activity. Even though a small pedigree could not show clear cosegregation with the disease phenotype, these findings strongly suggest the involvement of Nav1.2 in a human disease and propose the R187W mutation as the genetic defect responsible for febrile seizures associated with afebrile seizures.

High affinity type I interleukin 1 receptor antagonists discovered by screening recombinant peptide libraries.

Two families of peptides that specifically bind the extracellular domain of the human type I interleukin I (IL-1) receptor were identified from recombinant peptide display libraries. Peptides from one of these families blocked binding of IL-lalpha to the type I IL-1 receptor with IC50 values of 45-140 microM. Affinity-selective screening of variants of these peptides produced ligands of much higher affinity (IC50 approximately 2 nM). These peptides block IL-1-driven responses in human and monkey cells; they do not bind the human type II IL-1 receptor or the murine type I IL-1 receptor. This is the first example (that we know of) of a high affinity peptide that binds to a cytokine receptor and acts as a cytokine antagonist.

Structure and conformational changes of DNA topoisomerase II visualized by electron microscopy.

Type II DNA topoisomerases, which create a transient gate in duplex DNA and transfer a second duplex DNA through this gate, are essential for topological transformations of DNA in prokaryotic and eukaryotic cells and are of interest not only from a mechanistic perspective but also because they are targets of agents for anticancer and antimicrobial chemotherapy. Here we describe the structure of the molecule of human topoisomerase II [DNA topoisomerase (ATP-hydrolyzing), EC 5.99.1.3] as seen by scanning transmission electron microscopy. A globular approximately 90-angstrom diameter core is connected by linkers to two approximately 50-angstrom domains, which were shown by comparison with genetically truncated Saccharomyces cerevisiae topoisomerase II to contain the N-terminal region of the approximately 170-kDa subunits and that are seen in different orientations. When the ATP-binding site is occupied by a nonhydrolyzable ATP analog, a quite different structure is seen that results from a major conformational change and consists of two domains approximately 90 angstrom and approximately 60 angstrom in diameter connected by a linker, and in which the N-terminal domains have interacted. About two-thirds of the molecules show an approximately 25 A tunnel in the apical part of the large domain, and the remainder contain an internal cavity approximately 30 A wide in the large domain close to the linker region. We propose that structural rearrangements lead to this displacement of an internal tunnel. The tunnel is likely to represent the channel through which one DNA duplex, after capture in the clamp formed by the N-terminal domains, is transferred across the interface between the enzyme's subunits. These images are consistent with biochemical observations and provide a structural basis for understanding the reaction of topoisomerase II.

Expression of a renal type I sodium/phosphate transporter (NaPi-1) induces a conductance in Xenopus oocytes permeable for organic and inorganic anions.

Two distinct molecular types (I and II) of renal proximal tubular brush border Na+/Pi cotransporters have been identified by expression cloning on the basis of their capacity to induce Na+-dependent Pi influx in tracer experiments. Whereas the type II transporters (e.g., NaPi-2 and NaPi-3) resemble well known characteristics of brush border Na+/Pi cotransport, little is known about the properties of the type I transporter (NaPi-1). In contrast to type II, type I transporters produced electrogenic transport only at high extracellular Pi concentrations (> or =3 mM). On the other hand, expression of NaPi-1 induced a Cl- conductance in Xenopus laevis oocytes, which was inhibited by Cl- channel blockers [5-nitro-2-(3-phenylpropylamino)benzoic acid (NPPB) > niflumic acid >> 4,4'-diisothiocyanatostilbene-2,2'-disulfonic acid]. Further, the Cl- conductance was inhibited by the organic anions phenol red, benzylpenicillin (penicillin G), and probenecid. These organic anions induced outwardly directed currents in the absence of Cl-. In tracer studies, we observed uptake of benzylpenicillin with a Km of 0.22 mM; benzylpenicillin uptake was inhibited by NPPB and niflumic acid. These findings suggest that the type I Na+/Pi cotransporter functions also as a novel type of anion channel permeable not only for Cl- but also for organic anions. Such an apical anion channel could serve an important role in the transport of Cl- and the excretion of anionic xenobiotics.