GENETIC PREDICTORS OF HYPERGLYCEMIA DUE TO HYDROCHLOROTHIAZIDE THERAPY

Response to pharmacological treatment is variable among individuals. Some patients respond favorably to a drug while others develop adverse reactions. Early investigations showed evidence of variation in genes that code for drug receptors, drug transporters, and drug metabolizing enzymes; and pharmacogenetics appeared as the science that studies the relationship between drug response and genetic variation. Thiazide diuretics are the recommended first-line monotherapy for hypertension (i.e. SBP>140 or DBP>90). Even so, diuretics are associated with adverse metabolic side effects, such as hyperglycemia, which increase the risk of developing type 2 diabetes. Published approaches testing variation in candidate genes (e.g. the renin-angiotensin-aldosteron system (RAAS) and salt–sensitivity genes) have met with only limited success. We conducted the first genome wide association study to identify genes influencing hyperglycemia as an adverse effect of thiazide diuretics in non-Hispanic White hypertensive patients participating in the Genetic Epidemiology of Responses to Antihypertensives (GERA) and Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) clinical trials. No SNP reached the a priori defined threshold of statistical significance (p<5x10-8). We detected 50 SNPs in 9 genomic regions with suggestive p-values (p<1x10-5). Two of them, rs6870564 (p-value=3.28 X 10-6) and rs7702121 (p-value=5.09 X 10-6), were located close to biologic candidate genes, MYO and MGAT1, and one SNP in a genomic region in chromosome 6, rs7762018 (p-value=4.59 X 10-6) has been previously related to Insulin-Dependent Diabetes Mellitus (IDDM8). I conclude that 1) there are unlikely to be common SNPs with large effects on the adverse metabolic effects to hydrochlorothiazide treatment and 2) larger sample sizes are needed for pharmacogenetic studies of inter-individual variation in response to commonly prescribed medication.

HKT2;2/1, a K+-permeable transporter identified in a salt tolerant rice cultivar through surveys of natural genetic polymorphism

We have investigated OsHKT2;1 natural variation in a collection of 49 cultivars with different levels of salt tolerance and geographical origins. The effect of identified polymorphism on OsHKT2;1 activity was analysed through heterologous expression of variants in Xenopus oocytes. OsHKT2;1 appeared to be a highly conserved protein with only five possible amino acid substitutions that have no substantial effect on functional properties. Our study, however, also identified a new HKT isoform, No-OsHKT2;2/1 in Nona Bokra, a highly salt-tolerant cultivar. No-OsHKT2;2/1 probably originated from a deletion in chromosome 6, producing a chimeric gene. Its 5¢ region corresponds to that of OsHKT2;2, whose full-length sequence is not present in Nipponbare but has been identified in Pokkali, a salt-tolerant rice cultivar. Its 3¢ region corresponds to that of OsHKT2;1. No-OsHKT2;2/1 is essentially expressed in roots and displays a significant level of expression at high Na+ concentrations, in contrast to OsHKT2;1. Expressed in Xenopus oocytes or in Saccharomyces cerevisiae, No-OsHKT2;2/1 exhibited a strong permeability to Na+ and K+, even at high external Na+ concentrations, like OsHKT2;2, and in contrast to OsHKT2;1. Our results suggest that No-OsHKT2;2/1 can contribute to Nona Bokra salt tolerance by enabling root K+ uptake under saline conditions.

Comparison of parametric and nonparametric methods to map oligogenes by linkage

A sample of 95 sib pairs affected with insulin-dependent diabetes and typed with their normal parents for 28 markers on chromosome 6 has been analyzed by several methods. When appropriate parameters are efficiently estimated, a parametric model is equivalent to the β model, which is superior to nonparametric alternatives both in single point tests (as found previously) and in multipoint tests. Theory is given for meta-analysis combined with allelic association, and problems that may be associated with errors of map location and/or marker typing are identified. Reducing by multipoint analysis the number of association tests in a dense map can give a 3-fold reduction in the critical lod, and therefore in the cost of positional cloning.

Cloning, expression, and genetic mapping of Sema W, a member of the semaphorin family

The semaphorins comprise a large family of membrane-bound and secreted proteins, some of which have been shown to function in axon guidance. We have cloned a transmembrane semaphorin, Sema W, that belongs to the class IV subgroup of the semaphorin family. The mouse and rat forms of Sema W show 97% amino acid sequence identity with each other, and each shows about 91% identity with the human form. The gene for Sema W is divided into 15 exons, up to 4 of which are absent in the human cDNAs that we sequenced. Unlike many other semaphorins, Sema W is expressed at low levels in the developing embryo but was found to be expressed at high levels in the adult central nervous system and lung. Functional studies with purified membrane fractions from COS7 cells transfected with a Sema W expression plasmid showed that Sema W has growth-cone collapse activity against retinal ganglion-cell axons, indicating that vertebrate transmembrane semaphorins, like secreted semaphorins, can collapse growth cones. Genetic mapping of human SEMAW with human/hamster radiation hybrids localized the gene to chromosome 2p13. Genetic mapping of mouse Semaw with mouse/hamster radiation hybrids localized the gene to chromosome 6, and physical mapping placed the gene on bacteria artificial chromosomes carrying microsatellite markers D6Mit70 and D6Mit189. This localization places Semaw within the locus for motor neuron degeneration 2, making it an attractive candidate gene for this disease.

The natural killer gene complex genetic locus Chok encodes Ly-49D, a target recognition receptor that activates natural killing

Previously, we established that natural killer (NK) cells from C57BL/6 (B6), but not BALB/c, mice lysed Chinese hamster ovary (CHO) cells, and we mapped the locus that determines this differential CHO-killing capacity to the NK gene complex on chromosome 6. The localization of Chok in the NK gene complex suggested that it may encode either an activating or an inhibitory receptor. Here, results from a lectin-facilitated lysis assay predicted that Chok is an activating B6 NK receptor. Therefore, we immunized BALB/c mice with NK cells from BALB.B6–Cmv1r congenic mice and generated a mAb, designated 4E4, that blocked B6-mediated CHO lysis. mAb 4E4 also redirected lysis of Daudi targets, indicating its reactivity with an activating NK cell receptor. Furthermore, only the 4E4+ B6 NK cell subset mediated CHO killing, and this lysis was abrogated by preincubation with mAb 4E4. Flow cytometric analysis indicated that mAb 4E4 specifically reacts with Ly-49D but not Ly-49A, B, C, E, G, H, or I transfectants. Finally, gene transfer of Ly-49DB6 into BALB/c NK cells conferred cytotoxic capacity against CHO cells, thus establishing that the Ly-49D receptor is sufficient to activate NK cells to lyse this target. Hence, Ly-49D is the Chok gene product and is a mouse NK cell receptor capable of directly triggering natural killing.

Molecular cloning of Ian4: a BCR/ABL-induced gene that encodes an outer membrane mitochondrial protein with GTP-binding activity

Using the representation difference analysis technique, we have identified a novel gene, Ian4, which is preferentially expressed in hematopoietic precursor 32D cells transfected with wild-type versus mutant forms of the Bcr/Abl oncogene. Ian4 expression was undetectable in 32D cells transfected with v-src, oncogenic Ha-ras or v-Abl. Murine Ian4 maps to chromosome 6, 25 cM from the centromere. The Ian4 mRNA contains two open reading frames (ORFs) separated by 5 nt. The first ORF has the potential to encode for a polypeptide of 67 amino acids without apparent homology to known proteins. The second ORF encodes a protein of 301 amino acids with a GTP/ATP-binding site in the N-terminus and a hydrophobic domain in the extreme C-terminus. The IAN-4 protein resides in the mitochondrial outer membrane and the last 20 amino acids are necessary for this localization. The IAN-4 protein has GTP-binding activity and shares sequence homology with a novel family of putative GTP-binding proteins: the immuno-associated nucleotide (IAN) family.

IMGT/HLA Database—a sequence database for the human major histocompatibility complex

The IMGT/HLA Database (www.ebi.ac.uk/imgt/hla/) specialises in sequences of polymorphic genes of the HLA system, the human major histocompatibility complex (MHC). The HLA complex is located within the 6p21.3 region on the short arm of human chromosome 6 and contains more than 220 genes of diverse function. Many of the genes encode proteins of the immune system and these include the 21 highly polymorphic HLA genes, which influence the outcome of clinical transplantation and confer susceptibility to a wide range of non-infectious diseases. The database contains sequences for all HLA alleles officially recognised by the WHO Nomenclature Committee for Factors of the HLA System and provides users with online tools and facilities for their retrieval and analysis. These include allele reports, alignment tools and detailed descriptions of the source cells. The online IMGT/HLA submission tool allows both new and confirmatory sequences to be submitted directly to the WHO Nomenclature Committee. The latest version (release 1.7.0 July 2000) contains 1220 HLA alleles derived from over 2700 component sequences from the EMBL/GenBank/DDBJ databases. The HLA database provides a model which will be extended to provide specialist databases for polymorphic MHC genes of other species.

Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Mérino ewes

Ewes from the Booroola strain of Australian Mérino sheep are characterized by high ovulation rate and litter size. This phenotype is due to the action of the FecBB allele of a major gene named FecB, as determined by statistical analysis of phenotypic data. By genetic analysis of 31 informative half-sib families from heterozygous sires, we showed that the FecB locus is situated in the region of ovine chromosome 6 corresponding to the human chromosome 4q22–23 that contains the bone morphogenetic protein receptor IB (BMPR-IB) gene encoding a member of the transforming growth factor-β (TGF-β) receptor family. A nonconservative substitution (Q249R) in the BMPR-IB coding sequence was found to be associated fully with the hyperprolificacy phenotype of Booroola ewes. In vitro, ovarian granulosa cells from FecBB/FecBB ewes were less responsive than granulosa cells from FecB+/FecB+ ewes to the inhibitory effect on steroidogenesis of GDF-5 and BMP-4, natural ligands of BMPR-IB. It is suggested that in FecBB/FecBB ewes, BMPR-IB would be inactivated partially, leading to an advanced differentiation of granulosa cells and an advanced maturation of ovulatory follicles.

Identification and characterization of a second melanin-concentrating hormone receptor, MCH-2R

Melanin-concentrating hormone (MCH) is a 19-aa cyclic neuropeptide originally isolated from chum salmon pituitaries. Besides its effects on the aggregation of melanophores in fish several lines of evidence suggest that in mammals MCH functions as a regulator of energy homeostasis. Recently, several groups reported the identification of an orphan G protein-coupled receptor as a receptor for MCH (MCH-1R). We hereby report the identification of a second human MCH receptor termed MCH-2R, which shares about 38% amino acid identity with MCH-1R. MCH-2R displayed high-affinity MCH binding, resulting in inositol phosphate turnover and release of intracellular calcium in mammalian cells. In contrast to MCH-1R, MCH-2R signaling is not sensitive to pertussis toxin and MCH-2R cannot reduce forskolin-stimulated cAMP production, suggesting an exclusive Gαq coupling of the MCH-2R in cell-based systems. Northern blot and in situ hybridization analysis of human and monkey tissue shows that expression of MCH-2R mRNA is restricted to several regions of the brain, including the arcuate nucleus and the ventral medial hypothalamus, areas implicated in regulation of body weight. In addition, the human MCH-2R gene was mapped to the long arm of chromosome 6 at band 6q16.2–16.3, a region reported to be associated with cytogenetic abnormalities of obese patients. The characterization of a second mammalian G protein-coupled receptor for MCH potentially indicates that the control of energy homeostasis in mammals by the MCH neuropeptide system may be more complex than initially anticipated.

Cloning and comparative analysis of the human pre-T-cell receptor alpha-chain gene.

In immature T cells the T-cell receptor (TCR) beta-chain gene is rearranged and expressed before the TCR alpha-chain gene. At this stage TCR beta chain can form disulfide-linked heterodimers with the pre-T-cell receptor alpha chain (pTalpha). Using the recently isolated murine pTalpha cDNA as a probe, we have isolated the human pTalpha cDNA. The complete nucleotide sequence predicts a mature protein of 282 aa consisting of an extracellular immunoglobulin-like domain, a connecting peptide, a transmembrane region, and a long cytoplasmic tail. Amino acid sequence comparison of human pTalpha with the mouse pTalpha molecule reveals high sequence homology in the extracellular as well as the transmembrane region. In contrast, the cytoplasmic region differs in amino acid composition and in length from the murine homologue. The human pTalpha gene is expressed in immature but not mature T cells and is located at the p21.2-p12 region of the short arm of chromosome 6.

Genetic analysis of type 1 diabetes using whole genome approaches.

Whole genome linkage analysis of type 1 diabetes using affected sib pair families and semi-automated genotyping and data capture procedures has shown how type 1 diabetes is inherited. A major proportion of clustering of the disease in families can be accounted for by sharing of alleles at susceptibility loci in the major histocompatibility complex on chromosome 6 (IDDM1) and at a minimum of 11 other loci on nine chromosomes. Primary etiological components of IDDM1, the HLA-DQB1 and -DRB1 class II immune response genes, and of IDDM2, the minisatellite repeat sequence in the 5' regulatory region of the insulin gene on chromosome 11p15, have been identified. Identification of the other loci will involve linkage disequilibrium mapping and sequencing of candidate genes in regions of linkage.

Structural organization of mouse peroxisome proliferator-activated receptor gamma (mPPAR gamma) gene: alternative promoter use and different splicing yield two mPPAR gamma isoforms.

To gain insight into the regulation of expression of peroxisome proliferator-activated receptor (PPAR) isoforms, we have determined the structural organization of the mouse PPAR gamma (mPPAR gamma) gene. This gene extends > 105 kb and gives rise to two mRNAs (mPPAR gamma 1 and mPPAR gamma 2) that differ at their 5' ends. The mPPAR gamma 2 cDNA encodes an additional 30 amino acids N-terminal to the first ATG codon of mPPAR gamma 1 and reveals a different 5' untranslated sequence. We show that mPPAR gamma 1 mRNA is encoded by eight exons, whereas the mPPAR gamma 2 mRNA is encoded by seven exons. Most of the 5' untranslated sequence of mPPAR gamma 1 mRNA is encoded by two exons, whereas the 5' untranslated sequence and the extra 30 N-terminal amino acids of mPPAR gamma 2 are encoded by one exon, which is located between the second and third exons coding for mPPAR gamma 1. The last six exons of mPPAR gamma gene code for identical sequences in mPPAR gamma 1 and mPPAR gamma 2 isoforms. The mPPAR gamma 1 and mPPAR gamma 2 isoforms are transcribed from different promoters. The mPPAR gamma gene has been mapped to chromosome 6 E3-F1 by in situ hybridization using a biotin-labeled probe. These results establish that at least one of the PPAR genes yields more than one protein product, similar to that encountered with retinoid X receptor and retinoic acid receptor genes. The existence of multiple PPAR isoforms transcribed from different promoters could increase the diversity of ligand and tissue-specific transcriptional responses.

Investigação de mecanismos genéticos e epigenéticos em distúrbios do crescimento humano

Parcela considerável de pacientes com distúrbios de crescimento não têm a causa de seus quadros clínicos estabelecida, incluindo aproximadamente 50% dos pacientes com diagnóstico clínico de síndrome de Silver−Russell (SRS) e 10-20% dos pacientes com síndrome de Beckwith-Wiedemann (BWS). O objetivo deste estudo foi investigar as causas genéticas e epigenéticas de distúrbios de crescimento, de etiologia desconhecida, numa contribuição para o entendimento de mecanismos que regulam o crescimento. O estudo compreendeu: (1) a investigação de microdesequilíbrios cromossômicos, por aCGH; (2) a análise do perfil de expressão alelo-específica de genes sujeitos a imprinting (IG), por pirossequenciamento (PSQ) ou sequenciamento de Sanger; (3) a investigação do padrão de metilação global em pacientes com restrição de crescimento, utilizando microarray de metilação. A casuística constituiu-se de 41 pacientes não aparentados, com distúrbios de crescimento, de etiologia desconhecida: (1) 25, com hipótese diagnóstica de SRS; (2) seis, com restrição de crescimento intrauterino e peso ao nascimento abaixo do 10º percentil, associados a outros sinais clínicos; (3) sete, com hipótese diagnóstica de BWS; e (4) três, com macrossomia pré-natal ou pós-natal, associada a outros sinais. A investigação de microdesequilíbrios cromossômicos foi realizada em 40 pacientes. Foram detectadas 58 variantes raras em 30/40 pacientes (75%): 40 foram consideradas provavelmente benignas (18 pacientes, 45%), 12, com efeito patogênico desconhecido (11 pacientes, 27,5%), duas, provavelmente patogênicas (um paciente, 2,5%) e quatro, patogênicas (três pacientes, 7,5%). Essas frequências são comparáveis àquelas descritas em estudos que investigaram CNV em grupos de pacientes com distúrbios de crescimento e outras alterações congênitas, incluindo SRS, e mostram a importância da investigação de microdesequilíbrios cromossômicos nesses pacientes. A diversidade dos microdesequilíbrios cromossômicos identificados é reflexo da heterogeneidade clínica das casuísticas. Neste estudo, muitos dos pacientes com hipótese diagnóstica de SRS e BWS apresentavam sinais clínicos atípicos, explicando a ausência neles das alterações (epi)genéticas que causam essas síndromes. A identificação de CNV características de outras síndromes reflete a sobreposição de sinais clínicos com BWS e SRS. A análise do perfil de expressão alelo-específica de IG foi realizada em um subgrupo de 18 pacientes com restrição de crescimento. Trinta IG com função em proliferação celular, crescimento fetal ou neurodesenvolvimento foram inicialmente selecionados. Após seleção de SNP transcritos com alta frequência na população, genotipagem de pacientes, genitores e indivíduos controle, determinação da expressão dos IG em sangue periférico e seu padrão de expressão (mono ou bialélico), 13 IG, expressos no sangue, tiveram a expressão alelo-específica avaliada, sete deles por PSQ e seis por sequenciamento de Sanger. Alterações no perfil de expressão de dois genes, de expressão normalmente paterna, foram detectadas em 4/18 pacientes (22%). Este estudo é o primeiro a utilizar pirossequenciamento e sequenciamento de Sanger na avaliação do perfil de expressão alelo-específica de IG, em pacientes com restrição de crescimento. Apesar de terem limitações, ambas as técnicas mostraram-se robustas e revelaram alterações de expressão alélica interessantes; entretanto, a relação dessas alterações com o quadro clínico dos pacientes permanece por esclarecer. A investigação da metilação global do DNA foi realizada em subgrupo de 21 pacientes com restrição de crescimento e em 24 indivíduos controle. Dois tipos de análise foram realizados: (1) análise diferencial de grupo e (2) análise diferencial individual. Na primeira análise, em que foi comparado o padrão de metilação do grupo de pacientes com quadro clínico sugestivo de SRS (n=16) com o do grupo controle (n=24), não houve indicação de hipo ou hipermetilação global no grupo SRS. Na segunda análise, foi comparado o padrão de metilação de cada um dos 21 pacientes com restrição de crescimento e dos 24 indivíduos controle, com o padrão de metilação do grupo controle. O número médio de CpG hipermetilados e de segmentos diferencialmente metilados (SDM) foi significativamente maior nos pacientes. Foram identificados 82 SDM hipermetilados, estando 57 associados a gene(s) (69,5%), em 16 pacientes, e 51 SDM hipometilados, 41 deles associados a gene(s) (80,4%), em 10 pacientes. A análise de ontologia genética dos 61 genes associados aos SDM hipo ou hipermetilados nos pacientes destacou genes que atuam no desenvolvimento e na morfogênese do sistema esquelético e de órgãos fetais, e na regulação da transcrição gênica e de processos metabólicos. Alterações de metilação em genes que atuam em processos de proliferação e diferenciação celulares e crescimento foram identificadas em 9/20 dos pacientes (45%), sugerindo implicação clínica. Não foi detectada alteração epigenética comum aos pacientes com diagnóstico clínico de SRS, explicável provavelmente pela heterogeneidade clínica. A investigação de metilação global, utilizando microarray, produziu novos dados que podem contribuir para a compreensão de mecanismos moleculares que influenciam o crescimento pré- e pós-natal. Na translocação aparentemente equilibrada - t(5;6)(q35.2;p22.3)dn, detectada em paciente com suspeita clínica de SRS, a interrupção de um gene, pela quebra no cromossomo 6, pode ser a causa do quadro clínico; alternativamente, a translocação pode ter impactado a regulação de genes de desenvolvimento localizados próximos aos pontos de quebra. A análise de expressão em sangue periférico mostrou que os níveis de cDNA do gene, interrompido pelo ponto de quebra da translocação, estavam reduzidos à metade. Além de sinais típicos da SRS, a paciente apresentava algumas características clínicas sugestivas de displasia cleidocraniana. Assim, a translocação t(5;6) pode ter alterado a interação de genes de desenvolvimento e seus elementos reguladores, levando à desregulação de sua expressão espaço-temporal, e resultando num fenótipo atípico, com características sobrepostas de mais de uma síndrome genética

Methylomic and Transcriptomic Insights Linking Air Pollution and Atherosclerosis

Thesis (Ph.D.)--University of Washington, 2016-06

HLA and Genomewide Allele Sharing in Dizygotic Twins

Gametic selection during fertilization or the effects of specific genotypes on the viability of embryos may cause a skewed transmission of chromosomes to surviving offspring. A recent analysis of transmission distortion in humans reported significant excess sharing among full siblings. Dizygotic (DZ) twin pairs are a special case of the simultaneous survival of two genotypes, and there have been reports of DZ pairs with excess allele sharing around the HLA locus, a candidate locus for embryo survival. We performed an allele-sharing study of 1,592 DZ twin pairs from two independent Australian cohorts, of which 1,561 pairs were informative for linkage on chromosome 6. We also analyzed allele sharing in 336 DZ twin pairs from The Netherlands. We found no evidence of excess allele sharing, either at the HLA locus or in the rest of the genome. In contrast, we found evidence of a small but significant (P = .003 for the Australian sample) genomewide deficit in the proportion of two alleles shared identical by descent among DZ twin pairs. We reconciled conflicting evidence in the literature for excess genomewide allele sharing by performing a simulation study that shows how undetected genotyping errors can lead to an apparent deficit or excess of allele sharing among sibling pairs, dependent on whether parental genotypes are known. Our results imply that gene-mapping studies based on affected sibling pairs that include DZ pairs will not suffer from false-positive results due to loci involved in embryo survival.