Genome-wide association identifies nine common variants associated with fasting proinsulin levels and provides new insights into the pathophysiology of type 2 diabetes.

OBJECTIVE: Proinsulin is a precursor of mature insulin and C-peptide. Higher circulating proinsulin levels are associated with impaired β-cell function, raised glucose levels, insulin resistance, and type 2 diabetes (T2D). Studies of the insulin processing pathway could provide new insights about T2D pathophysiology. RESEARCH DESIGN AND METHODS: We have conducted a meta-analysis of genome-wide association tests of ∼2.5 million genotyped or imputed single nucleotide polymorphisms (SNPs) and fasting proinsulin levels in 10,701 nondiabetic adults of European ancestry, with follow-up of 23 loci in up to 16,378 individuals, using additive genetic models adjusted for age, sex, fasting insulin, and study-specific covariates. RESULTS: Nine SNPs at eight loci were associated with proinsulin levels (P < 5 × 10(-8)). Two loci (LARP6 and SGSM2) have not been previously related to metabolic traits, one (MADD) has been associated with fasting glucose, one (PCSK1) has been implicated in obesity, and four (TCF7L2, SLC30A8, VPS13C/C2CD4A/B, and ARAP1, formerly CENTD2) increase T2D risk. The proinsulin-raising allele of ARAP1 was associated with a lower fasting glucose (P = 1.7 × 10(-4)), improved β-cell function (P = 1.1 × 10(-5)), and lower risk of T2D (odds ratio 0.88; P = 7.8 × 10(-6)). Notably, PCSK1 encodes the protein prohormone convertase 1/3, the first enzyme in the insulin processing pathway. A genotype score composed of the nine proinsulin-raising alleles was not associated with coronary disease in two large case-control datasets. CONCLUSIONS: We have identified nine genetic variants associated with fasting proinsulin. Our findings illuminate the biology underlying glucose homeostasis and T2D development in humans and argue against a direct role of proinsulin in coronary artery disease pathogenesis.

Direct Effect of Birth Weight on Childhood Blood Pressure: a Causal Mediation Analysis

Background: Numerous studies have shown a negative association between birth weight (BW) and blood pressure (BP) later in life. To estimate the direct effect of BW on BP, it is conventional to condition on current weight (CW). However, such conditioning can induce collider stratification bias in the estimate of the direct effect. Objective: To bound the potential bias due to U, an unmeasured common cause of CW and BP, on the estimate of the (controlled) direct effect of BW on BP. Methods: Data from a school based study in Switzerland were used (N = 4,005; 2,010 B/1,995 G; mean age: 12.3 yr [range: 10.1-14.9]). Measured common causes of BW-BP (SES, smoking, body weight, and hypertension status of the mother) and CW-BP (breastfeeding and child's physical activity and diet) were identified with DAGs. Linear regression models were fitted to estimate the association between BW and BP. Sensitivity analyses were conducted to assess the potential effect of U on the association between BW and BP. U was assumed 1) to be a binary variable that affected BP by the same magnitude in low BWand in normal BW children and 2) to have a different prevalence in low BW children and in normal BW children for a given CW. Results: A small negative association was observed between BW and BP [beta: -0.3 mmHg/kg (95% CI: -0.9 to 0.3)]. The association was strengthened upon conditioning for CW [beta: -1.5 mmHg/kg (95% CI: -2.1 to -0.9)]. Upon further conditioning on common causes of BW-BP and CW-BP, the association did not change substantially [beta: -1.4 mmHg/kg (95% CI: -2.0 to -0.8)]. The negative association could be explained by U only if U was strongly associated with BP and if there was a large difference in the prevalence of U between low BWand normal BW children. Conclusion: The observed negative association between BW and BP upon adjustment for CW was not easily explained by an unmeasured common cause of CWand BP.

Comparative genomic and phylogeographic analysis of Mycobacterium leprae.

Reductive evolution and massive pseudogene formation have shaped the 3.31-Mb genome of Mycobacterium leprae, an unculturable obligate pathogen that causes leprosy in humans. The complete genome sequence of M. leprae strain Br4923 from Brazil was obtained by conventional methods (6x coverage), and Illumina resequencing technology was used to obtain the sequences of strains Thai53 (38x coverage) and NHDP63 (46x coverage) from Thailand and the United States, respectively. Whole-genome comparisons with the previously sequenced TN strain from India revealed that the four strains share 99.995% sequence identity and differ only in 215 polymorphic sites, mainly SNPs, and by 5 pseudogenes. Sixteen interrelated SNP subtypes were defined by genotyping both extant and extinct strains of M. leprae from around the world. The 16 SNP subtypes showed a strong geographical association that reflects the migration patterns of early humans and trade routes, with the Silk Road linking Europe to China having contributed to the spread of leprosy.

Analysis and spatial modeling of the indoor radon Swiss data

The present research deals with an important public health threat, which is the pollution created by radon gas accumulation inside dwellings. The spatial modeling of indoor radon in Switzerland is particularly complex and challenging because of many influencing factors that should be taken into account. Indoor radon data analysis must be addressed from both a statistical and a spatial point of view. As a multivariate process, it was important at first to define the influence of each factor. In particular, it was important to define the influence of geology as being closely associated to indoor radon. This association was indeed observed for the Swiss data but not probed to be the sole determinant for the spatial modeling. The statistical analysis of data, both at univariate and multivariate level, was followed by an exploratory spatial analysis. Many tools proposed in the literature were tested and adapted, including fractality, declustering and moving windows methods. The use of Quan-tité Morisita Index (QMI) as a procedure to evaluate data clustering in function of the radon level was proposed. The existing methods of declustering were revised and applied in an attempt to approach the global histogram parameters. The exploratory phase comes along with the definition of multiple scales of interest for indoor radon mapping in Switzerland. The analysis was done with a top-to-down resolution approach, from regional to local lev¬els in order to find the appropriate scales for modeling. In this sense, data partition was optimized in order to cope with stationary conditions of geostatistical models. Common methods of spatial modeling such as Κ Nearest Neighbors (KNN), variography and General Regression Neural Networks (GRNN) were proposed as exploratory tools. In the following section, different spatial interpolation methods were applied for a par-ticular dataset. A bottom to top method complexity approach was adopted and the results were analyzed together in order to find common definitions of continuity and neighborhood parameters. Additionally, a data filter based on cross-validation was tested with the purpose of reducing noise at local scale (the CVMF). At the end of the chapter, a series of test for data consistency and methods robustness were performed. This lead to conclude about the importance of data splitting and the limitation of generalization methods for reproducing statistical distributions. The last section was dedicated to modeling methods with probabilistic interpretations. Data transformation and simulations thus allowed the use of multigaussian models and helped take the indoor radon pollution data uncertainty into consideration. The catego-rization transform was presented as a solution for extreme values modeling through clas-sification. Simulation scenarios were proposed, including an alternative proposal for the reproduction of the global histogram based on the sampling domain. The sequential Gaussian simulation (SGS) was presented as the method giving the most complete information, while classification performed in a more robust way. An error measure was defined in relation to the decision function for data classification hardening. Within the classification methods, probabilistic neural networks (PNN) show to be better adapted for modeling of high threshold categorization and for automation. Support vector machines (SVM) on the contrary performed well under balanced category conditions. In general, it was concluded that a particular prediction or estimation method is not better under all conditions of scale and neighborhood definitions. Simulations should be the basis, while other methods can provide complementary information to accomplish an efficient indoor radon decision making.

Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways.

Through genome-wide association meta-analyses of up to 133,010 individuals of European ancestry without diabetes, including individuals newly genotyped using the Metabochip, we have increased the number of confirmed loci influencing glycemic traits to 53, of which 33 also increase type 2 diabetes risk (q < 0.05). Loci influencing fasting insulin concentration showed association with lipid levels and fat distribution, suggesting impact on insulin resistance. Gene-based analyses identified further biologically plausible loci, suggesting that additional loci beyond those reaching genome-wide significance are likely to represent real associations. This conclusion is supported by an excess of directionally consistent and nominally significant signals between discovery and follow-up studies. Functional analysis of these newly discovered loci will further improve our understanding of glycemic control.

Clear detection of ADIPOQ locus as the major gene for plasma adiponectin: results of genome-wide association analyses including 4659 European individuals.

OBJECTIVE: Plasma adiponectin is strongly associated with various components of metabolic syndrome, type 2 diabetes and cardiovascular outcomes. Concentrations are highly heritable and differ between men and women. We therefore aimed to investigate the genetics of plasma adiponectin in men and women. METHODS: We combined genome-wide association scans of three population-based studies including 4659 persons. For the replication stage in 13795 subjects, we selected the 20 top signals of the combined analysis, as well as the 10 top signals with p-values less than 1.0 x 10(-4) for each the men- and the women-specific analyses. We further selected 73 SNPs that were consistently associated with metabolic syndrome parameters in previous genome-wide association studies to check for their association with plasma adiponectin. RESULTS: The ADIPOQ locus showed genome-wide significant p-values in the combined (p=4.3 x 10(-24)) as well as in both women- and men-specific analyses (p=8.7 x 10(-17) and p=2.5 x 10(-11), respectively). None of the other 39 top signal SNPs showed evidence for association in the replication analysis. None of 73 SNPs from metabolic syndrome loci exhibited association with plasma adiponectin (p>0.01). CONCLUSIONS: We demonstrated the ADIPOQ gene as the only major gene for plasma adiponectin, which explains 6.7% of the phenotypic variance. We further found that neither this gene nor any of the metabolic syndrome loci explained the sex differences observed for plasma adiponectin. Larger studies are needed to identify more moderate genetic determinants of plasma adiponectin.

Common variants in UMOD associate with urinary uromodulin levels: a meta-analysis.

Uromodulin is expressed exclusively in the thick ascending limb and is the most abundant protein excreted in normal urine. Variants in UMOD, which encodes uromodulin, are associated with renal function, and urinary uromodulin levels may be a biomarker for kidney disease. However, the genetic factors regulating uromodulin excretion are unknown. We conducted a meta-analysis of urinary uromodulin levels to identify associated common genetic variants in the general population. We included 10,884 individuals of European descent from three genetic isolates and three urban cohorts. Each study measured uromodulin indexed to creatinine and conducted linear regression analysis of approximately 2.5 million single nucleotide polymorphisms using an additive model. We also tested whether variants in genes expressed in the thick ascending limb associate with uromodulin levels. rs12917707, located near UMOD and previously associated with renal function and CKD, had the strongest association with urinary uromodulin levels (P<0.001). In all cohorts, carriers of a G allele of this variant had higher uromodulin levels than noncarriers did (geometric means 10.24, 14.05, and 17.67 μg/g creatinine for zero, one, or two copies of the G allele). rs12446492 in the adjacent gene PDILT (protein disulfide isomerase-like, testis expressed) also reached genome-wide significance (P<0.001). Regarding genes expressed in the thick ascending limb, variants in KCNJ1, SORL1, and CAB39 associated with urinary uromodulin levels. These data indicate that common variants in the UMOD promoter region may influence urinary uromodulin levels. They also provide insights into uromodulin biology and the association of UMOD variants with renal function.

Role of hepatitis C virus genotype 3 in liver fibrosis progression : a systematic review and meta-analysis ; et, Impact of nurse vaccination program on hepatitis B immunity in a Swiss HIV clinic

Rôle du génotype 3 du virus de l'hépatite C dans la progression de la fibrose hépatique, une revue systématique avec méta-analyse. On estime à 170 millions le nombre de personnes atteintes d'hépatite C chronique dans le monde. La principale conséquence de cette maladie est la fibrose du foie, qui évolue plus ou moins rapidement, pour aboutir au développement d'une cirrhose et/ou d'un hépatocarcinome. Certains des facteurs accélérateurs de la fibrose, comme l'âge avancé au moment de l'infection, le sexe masculin, la consommation d'alcool, sont bien connus. On a longtemps considéré que les six différents génotypes viraux n'influençaient pas la progression de la fibrose. Des études récentes ont cependant suggéré que certains génotypes, en particulier ie génotype 3, pouvaient entraîner une fibrose plus rapide. Le but de ce travail de thèse était de déterminer à l'aide d'une méta-analyse le rôle du génotype viral dans la progression de la fibrose dans l'infection chronique au virus de l'hépatite C. Les études ont été sélectionnées dans la littérature médicale à partir d'une série de mots-clés. Le degré de fibrose a été estimé par biopsie, en utilisant le score Metavir. Deux types d'études ont décrits de manière différente la durée d'infection. Les premières ont calculé la progression de la fibrose depuis le moment estimée de l'infection (« études avec une biopsie »), les secondes ont exprimés cette durée comme étant l'intervalle entre deux biopsies (« études avec deux biopsies »). L'analyse a permis d'identifier 8 études avec une biopsie pour un collectif total de 3182 patients ainsi que 8 études avec deux biopsies pour un collectif de 896 patients. Dans une méta-analyse de type « random effect », le rapport de cote pour l'association du génotype 3 avec une fibrose accélérée est de 1.52 (95% IC 1.12-2.07, p=0.007) pour les études à une biopsie. Pour les études à deux biopsies, le rapport de cote pour cette association est de 1.37 (95% IC 0.87-2.17, P=0.17). Cette étude montre que les patients avec une hépatite C chronique due au génotype 3 ont une progression de fibrose plus rapide que ceux qui sont infectés par les autres génotypes. Alors que la méta-analyse des études avec une biopsie est clairement significative, celle des études avec deux biopsies est au-dessous du seuil de significativité. Les études à deux biopsies peuvent être limitées par plusieurs facteurs, comprenant un « biais d'indication » (seuls les patients évoluant rapidement vers la cirrhose ont plus de risque d'avoir une deuxième biopsie), une durée d'observation très courte (5 années comparée à 13 années pour les études à 2 biopsies), et un nombre de patient limité (896 pour le études à 2 biopsies comparé à 3182 pour les études à 1 biopsie). Impact d'un programme de vaccination sur l'immunité contre l'hépatite Β dans une clinique suisse du VIH Le virus de l'hépatite Β cause une infection aigûe dont la symptomatologie varie d'une présentation subclinique à une progression fulminante. Dans une minorité de cas, l'infection aigiie est suivie d'une infection chronique pouvant évoluer vers une cirrhose hépatique et/ou un hépatocarcinome. La prévalence de l'hépatite Β aiguë et chronique chez les personnes vivant avec le virus d'immunodéficience humaine (VIH) est supérieure à celle de la population générale. Par ailleurs la co-infection avec le virus du VIH entraine une progression plus rapide de l'hépatite B. Dès lors, l'immunité pour le virus de l'hépatite Β représente un facteur primordial de prévention dans la population infectée par le virus de l'HIV. Bien que l'administration d'un vaccin contre l'hépatite Β soit particulièrement recommandée chez tous les individus infectés par le VIH, la couverture vaccinale dans cette population est souvent insuffisante. Le but de cette étude était de déterminer l'état d'immunisation contre le virus de l'hépatite Β dans la population infectée par le VIH de la cohorte Suisse HIV et d'analyser l'efficacité d'un programme de vaccination administré par le personnel soignant. L'immunité avant et après intervention dans notre centre a été comparée aux autres centres de la cohorte HIV en Suisse. L'immunité pour le centre d'intervention a passé de 32% avant intervention à 76% après intervention alors que pour les autres centres, l'immunité n'a progressé que de 33% à 39% dans le même laps de temps (n=2712, P=0.001). Cette étude montre qu'un contrôle systématique de l'immunité par du personnel soignant augmente de manière significative l'immunité pour le vaccin de l'hépatite Β dans la population HIV.

Interactions of dietary whole-grain intake with fasting glucose- and insulin-related genetic loci in individuals of European descent: a meta-analysis of 14 cohort studies.

Whole-grain foods are touted for multiple health benefits, including enhancing insulin sensitivity and reducing type 2 diabetes risk. Recent genome-wide association studies (GWAS) have identified several single nucleotide polymorphisms (SNPs) associated with fasting glucose and insulin concentrations in individuals free of diabetes. We tested the hypothesis that whole-grain food intake and genetic variation interact to influence concentrations of fasting glucose and insulin. Via meta-analysis of data from 14 cohorts comprising ∼ 48,000 participants of European descent, we studied interactions of whole-grain intake with loci previously associated in GWAS with fasting glucose (16 loci) and/or insulin (2 loci) concentrations. For tests of interaction, we considered a P value <0.0028 (0.05 of 18 tests) as statistically significant. Greater whole-grain food intake was associated with lower fasting glucose and insulin concentrations independent of demographics, other dietary and lifestyle factors, and BMI (β [95% CI] per 1-serving-greater whole-grain intake: -0.009 mmol/l glucose [-0.013 to -0.005], P < 0.0001 and -0.011 pmol/l [ln] insulin [-0.015 to -0.007], P = 0.0003). No interactions met our multiple testing-adjusted statistical significance threshold. The strongest SNP interaction with whole-grain intake was rs780094 (GCKR) for fasting insulin (P = 0.006), where greater whole-grain intake was associated with a smaller reduction in fasting insulin concentrations in those with the insulin-raising allele. Our results support the favorable association of whole-grain intake with fasting glucose and insulin and suggest a potential interaction between variation in GCKR and whole-grain intake in influencing fasting insulin concentrations.

Crystallization and preliminary X-ray analysis of monalysin, a novel β-pore-forming toxin from the entomopathogen Pseudomonas entomophila.

Monalysin was recently described as a novel pore-forming toxin (PFT) secreted by the Drosophila pathogen Pseudomonas entomophila. Recombinant monalysin is multimeric in solution, whereas PFTs are supposed to be monomeric until target membrane association. Monalysin crystals were obtained by the hanging-drop vapour-diffusion method using PEG 8000 as precipitant. Preliminary X-ray diffraction analysis revealed that monalysin crystals belonged to the monoclinic space group C2, with unit-cell parameters a = 162.4, b = 146.2, c = 144.4 Å, β = 122.8°, and diffracted to 2.85 Å resolution using synchrotron radiation. Patterson self-rotation analysis and Matthews coefficient calculation indicate that the asymmetric unit contains nine copies of monalysin. Heavy-atom derivative data were collected and a Ta6Br14 cluster derivative data set confirmed the presence of ninefold noncrystallographic symmetry.

Persistent association of nailfold capillaroscopy changes and skin involvement over thirty-six months with duration of untreated disease in patients with juvenile dermatomyositis.

OBJECTIVE: To determine the association of changes on nailfold capillaroscopy with clinical findings and genotype in children with juvenile dermatomyositis (DM), in order to identify potential differences in disease course over 36 months. METHODS: At diagnosis of juvenile DM in 61 children prior to the initiation of treatment, tumor necrosis factor alpha (TNFalpha) -308 allele and DQA1*0501 status was determined, juvenile DM Disease Activity Scores (DAS) were obtained, and nailfold capillaroscopy was performed. The disease course was monitored for 36 months. Variations within and between patients were assessed by regression analysis. RESULTS: At diagnosis, shorter duration of untreated disease (P = 0.05) and a lower juvenile DM skin DAS (P = 0.035) were associated with a unicyclic disease course. Over 36 months, end-row loop (ERL) regeneration was associated with lower skin DAS (P &lt; 0.001) but not muscle DAS (P = 0.98); ERL regeneration and decreased bushy loops were associated with a shorter duration of untreated disease (P = 0.04 for both). At 36 months, increased ERL regeneration (P = 0.007) and improvement of skin DAS (P &lt; 0.001) and muscle DAS (P = 0.025) were associated with a unicyclic disease course. CONCLUSION: Early treatment of juvenile DM may lead to a unicyclic disease course. The non-unicyclic disease course usually involves continuing skin manifestations with persistent nailfold capillaroscopy changes. The correlation of nailfold capillaroscopy results with cutaneous but not with musculoskeletal signs of juvenile DM over a 36-month period suggests that the cutaneous and muscle vasculopathies have different pathophysiologic mechanisms. These findings indicate that efforts to identify the optimal treatment of cutaneous features in juvenile DM require greater attention.

Association naevus et papillome de la caroncule. A propos d'un cas observation anatomo-clinique [Papillomatous and nevus of the caruncle].

A 63-year-old male patient was operated of a single conjunctival papillomatous lesion located on the left caruncle. The curative tumoral excision was followed by its histological analysis. Histologically, the tumor associated an epithelial squamous papilloma with a subepithelial nevus, confirming a combination of two different lesions.

Predictive analysis and mapping of indoor radon concentrations in Switzerland

It is estimated that around 230 people die each year due to radon (222Rn) exposure in Switzerland. 222Rn occurs mainly in closed environments like buildings and originates primarily from the subjacent ground. Therefore it depends strongly on geology and shows substantial regional variations. Correct identification of these regional variations would lead to substantial reduction of 222Rn exposure of the population based on appropriate construction of new and mitigation of already existing buildings. Prediction of indoor 222Rn concentrations (IRC) and identification of 222Rn prone areas is however difficult since IRC depend on a variety of different variables like building characteristics, meteorology, geology and anthropogenic factors. The present work aims at the development of predictive models and the understanding of IRC in Switzerland, taking into account a maximum of information in order to minimize the prediction uncertainty. The predictive maps will be used as a decision-support tool for 222Rn risk management. The construction of these models is based on different data-driven statistical methods, in combination with geographical information systems (GIS). In a first phase we performed univariate analysis of IRC for different variables, namely the detector type, building category, foundation, year of construction, the average outdoor temperature during measurement, altitude and lithology. All variables showed significant associations to IRC. Buildings constructed after 1900 showed significantly lower IRC compared to earlier constructions. We observed a further drop of IRC after 1970. In addition to that, we found an association of IRC with altitude. With regard to lithology, we observed the lowest IRC in sedimentary rocks (excluding carbonates) and sediments and the highest IRC in the Jura carbonates and igneous rock. The IRC data was systematically analyzed for potential bias due to spatially unbalanced sampling of measurements. In order to facilitate the modeling and the interpretation of the influence of geology on IRC, we developed an algorithm based on k-medoids clustering which permits to define coherent geological classes in terms of IRC. We performed a soil gas 222Rn concentration (SRC) measurement campaign in order to determine the predictive power of SRC with respect to IRC. We found that the use of SRC is limited for IRC prediction. The second part of the project was dedicated to predictive mapping of IRC using models which take into account the multidimensionality of the process of 222Rn entry into buildings. We used kernel regression and ensemble regression tree for this purpose. We could explain up to 33% of the variance of the log transformed IRC all over Switzerland. This is a good performance compared to former attempts of IRC modeling in Switzerland. As predictor variables we considered geographical coordinates, altitude, outdoor temperature, building type, foundation, year of construction and detector type. Ensemble regression trees like random forests allow to determine the role of each IRC predictor in a multidimensional setting. We found spatial information like geology, altitude and coordinates to have stronger influences on IRC than building related variables like foundation type, building type and year of construction. Based on kernel estimation we developed an approach to determine the local probability of IRC to exceed 300 Bq/m3. In addition to that we developed a confidence index in order to provide an estimate of uncertainty of the map. All methods allow an easy creation of tailor-made maps for different building characteristics. Our work is an essential step towards a 222Rn risk assessment which accounts at the same time for different architectural situations as well as geological and geographical conditions. For the communication of 222Rn hazard to the population we recommend to make use of the probability map based on kernel estimation. The communication of 222Rn hazard could for example be implemented via a web interface where the users specify the characteristics and coordinates of their home in order to obtain the probability to be above a given IRC with a corresponding index of confidence. Taking into account the health effects of 222Rn, our results have the potential to substantially improve the estimation of the effective dose from 222Rn delivered to the Swiss population.

Antiretroviral therapy during pregnancy and premature birth: analysis of Swiss data.

BACKGROUND: There is an ongoing debate as to whether combined antiretroviral treatment (cART) during pregnancy is an independent risk factor for prematurity in HIV-1-infected women. OBJECTIVE: The aim of the study was to examine (1) crude effects of different ART regimens on prematurity, (2) the association between duration of cART and duration of pregnancy, and (3) the role of possibly confounding risk factors for prematurity. METHOD: We analysed data from 1180 pregnancies prospectively collected by the Swiss Mother and Child HIV Cohort Study (MoCHiV) and the Swiss HIV Cohort Study (SHCS). RESULTS: Odds ratios for prematurity in women receiving mono/dual therapy and cART were 1.8 [95% confidence interval (CI) 0.85-3.6] and 2.5 (95% CI 1.4-4.3) compared with women not receiving ART during pregnancy (P=0.004). In a subgroup of 365 pregnancies with comprehensive information on maternal clinical, demographic and lifestyle characteristics, there was no indication that maternal viral load, age, ethnicity or history of injecting drug use affected prematurity rates associated with the use of cART. Duration of cART before delivery was also not associated with duration of pregnancy. CONCLUSION: Our study indicates that confounding by maternal risk factors or duration of cART exposure is not a likely explanation for the effects of ART on prematurity in HIV-1-infected women.