Tea consumption and incidence of type 2 diabetes in Europe: the EPIC-InterAct case-cohort study.

BACKGROUND In previous meta-analyses, tea consumption has been associated with lower incidence of type 2 diabetes. It is unclear, however, if tea is associated inversely over the entire range of intake. Therefore, we investigated the association between tea consumption and incidence of type 2 diabetes in a European population. METHODOLOGY/PRINCIPAL FINDINGS The EPIC-InterAct case-cohort study was conducted in 26 centers in 8 European countries and consists of a total of 12,403 incident type 2 diabetes cases and a stratified subcohort of 16,835 individuals from a total cohort of 340,234 participants with 3.99 million person-years of follow-up. Country-specific Hazard Ratios (HR) for incidence of type 2 diabetes were obtained after adjustment for lifestyle and dietary factors using a Cox regression adapted for a case-cohort design. Subsequently, country-specific HR were combined using a random effects meta-analysis. Tea consumption was studied as categorical variable (0, >0-<1, 1-<4, ≥ 4 cups/day). The dose-response of the association was further explored by restricted cubic spline regression. Country specific medians of tea consumption ranged from 0 cups/day in Spain to 4 cups/day in United Kingdom. Tea consumption was associated inversely with incidence of type 2 diabetes; the HR was 0.84 [95%CI 0.71, 1.00] when participants who drank ≥ 4 cups of tea per day were compared with non-drinkers (p(linear trend) = 0.04). Incidence of type 2 diabetes already tended to be lower with tea consumption of 1-<4 cups/day (HR = 0.93 [95%CI 0.81, 1.05]). Spline regression did not suggest a non-linear association (p(non-linearity) = 0.20). CONCLUSIONS/SIGNIFICANCE A linear inverse association was observed between tea consumption and incidence of type 2 diabetes. People who drink at least 4 cups of tea per day may have a 16% lower risk of developing type 2 diabetes than non-tea drinkers.

Applications of stochastic analysis in finance and statistical inference

First: A continuous-time version of Kyle's model (Kyle 1985), known as the Back's model (Back 1992), of asset pricing with asymmetric information, is studied. A larger class of price processes and of noise traders' processes are studied. The price process, as in Kyle's model, is allowed to depend on the path of the market order. The process of the noise traders' is an inhomogeneous Lévy process. Solutions are found by the Hamilton-Jacobi-Bellman equations. With the insider being risk-neutral, the price pressure is constant, and there is no equilibirium in the presence of jumps. If the insider is risk-averse, there is no equilibirium in the presence of either jumps or drifts. Also, it is analised when the release time is unknown. A general relation is established between the problem of finding an equilibrium and of enlargement of filtrations. Random announcement time is random is also considered. In such a case the market is not fully efficient and there exists equilibrium if the sensitivity of prices with respect to the global demand is time decreasing according with the distribution of the random time. Second: Power variations. it is considered, the asymptotic behavior of the power variation of processes of the form _integral_0^t u(s-)dS(s), where S_ is an alpha-stable process with index of stability 0&alpha&2 and the integral is an Itô integral. Stable convergence of corresponding fluctuations is established. These results provide statistical tools to infer the process u from discrete observations. Third: A bond market is studied where short rates r(t) evolve as an integral of g(t-s)sigma(s) with respect to W(ds), where g and sigma are deterministic and W is the stochastic Wiener measure. Processes of this type are particular cases of ambit processes. These processes are in general not of the semimartingale kind.

Lymphoblastic transformation of follicular lymphoma: a clinicopathologic and molecular analysis of 7 patients.

Approximately 30% of patients with follicular lymphoma (FL) transform to a more aggressive malignancy, most commonly diffuse large B cell lymphoma. Rarely, FL transformation results in clinical findings, histology, and immunophenotype reminiscent of B-lymphoblastic leukemia/lymphoma. We report the largest series to date with detailed analysis of 7 such patients. Lymphoblastic transformation occurred on average 2 years after initial diagnosis of FL. Five patients had prior intensive chemotherapy. Two patients developed mature high-grade lymphoma, followed by the lymphoblastic transformation. FL had BCL2 gene rearrangement in 4 of 5 cases. High-grade transformation was accompanied by MYC gene rearrangement (5 of 5). Transformation was characterized by expression of TdT, loss of Bcl6, variable loss of immunoglobulin light chain, and persistence of Pax-5, Bcl2, and CD10. Whole-exome sequencing in 1 case revealed presence of several actionable mutations (CD79B, CCND3, CDK12). FL, aggressive mature B cell lymphoma, and lymphoblastic transformation were clonally related in 6 evaluable cases. After transformation, survival ranged from 1 to 14 months. Four patients died of disease, 2 were in remission after stem cell transplant, and 1 was alive with disease.

Association analysis of human leukocyte antigen class II (DRB1) alleles with leprosy in individuals from São Luís, state of Maranhão, Brazil

Epidemiological studies have demonstrated that the variability of the clinical response to infection caused by Mycobacterium leprae is associated with host genetic factors. The present study investigated the frequency of human leukocyte antigen (HLA) class II (DRB1) alleles in patients with leprosy from São Luís, Maranhão, Brazil. A case-control study was performed in 85 individuals with leprosy and 85 healthy subjects. All samples were analysed via polymerase chain reaction-sequence specific oligonucleotide probes. The HLA-DRB1*16 allele showed a higher frequency in the group with leprosy [(9.41% vs. 4.12%) odds ratio (OR) = 2.41 95% confidence interval (CI) (0.96-6.08) p = 0.05], whereas the HLA-DRB1*11 allele was less frequent in the group with leprosy [(6.47% vs. 11.76%) OR = 0.51 95% CI (0.23-1.12) p = 0.09]. The frequency of HLA-DRB1* alleles between the control group and leprosy patient subgroups presenting different forms of the disease showed that the HLA-DRB1*16 (16.13% vs. 8.24%, OR = 4.10, CI = 1.27-13.27, p = 0.010) and HLA-DRB1*14 (5% vs. 3.53%, OR = 4.63, CI = 1.00-21.08, p = 0.032) alleles were significantly more frequent in patients with different clinical subtypes of leprosy. The sample size was a limitation in this study. Nevertheless, the results demonstrated the existence of a genetic susceptibility associated with the clinical forms of leprosy. The low frequency of the HLA-DRB1*11 allele should be further studied to investigate the possible protective effect of this allele.

Vertebral fractures in patients with inflammatory bowel disease compared with a healthy population: a prospective case-control study.

BACKGROUND A prospective study was performed to compare the prevalence of morphometric vertebral fractures (MVF) between patients with inflammatory bowel disease (IBD) and healthy subjects and to identify predictive factors of fracture. METHODS A total of 107 patients with IBD (53 with Crohn's disease and 54 with ulcerative colitis) and 51 healthy subjects participated in the study. Information about anthropometric parameters, toxins, previous fractures, and parameters related to this disease were evaluated. The index of vertebral deformity, bone mass density (BMD), and biochemical parameters were calculated. RESULTS A total of 72 fractures were detected in 38.32% of patients with IBD, and 10 fractures were detected in 13.73% of healthy subjects; the risk of fracture in patients with IBD was higher than that in control subjects (OR, 4.03; 95% CI, 1.652-9.847; p < 0.002). We found no correlation between fracture and BMD in patients with IBD (lumbar spine, r = -0.103, p = 0.17 and femoral neck, r = -0.138, p = 0.07). Corticosteroid treatment was not associated with prevalent vertebral fractures nor with taking corticosteroids (r = 0.135, p = 0.14) or the duration for which they were taken (r = 0.08, p = 0.38), whereas this relationship was present in the controls (r = -0.365, p = 0.01). In the multivariate analysis, none of the measured parameters were significantly predictive of fracture, only to manifested IBD. Hypovitaminosis D was observed in 55.14% of patients with IBD. CONCLUSIONS The prevalence of morphometric vertebral fractures is higher in patients with IBD than in the healthy population, without association with BMD or corticoid treatment. Simply having IBD was proven to be a predictive factor of fracture. We observed a high incidence of hypovitaminosis D in patients with IBD.

Constraints on the upper Permian to upper Triassic marine carbon isotope curve : case studies from the Tethys

Secular variations of the seawater carbon isotopic composition provide evidence for paleoceanographic and paleoclimatic changes and may serve for chemiostratigraphic correlations. The present study aimed to improve the current knowledge on the Upper Permian and Triassic segment of the Phanerozoic marine carbon isotope curve, whose Triassic part was poorly constrained by previous studies. Profiles of inorganic carbon isotopes are provided for sections from Himalaya (Salt Range, Kashmir, Spiti and Nepal), Oman and North Dobrogea (Romania) on the basis of whole-rock carbonate analysis. The data acquired, together with a literature compilation confirmed that most of the Upper Permian is characterized by high δ13C values (averaging +40/00) but failed to detect a positive excursion as suggested by recent compilations. In the light of these observations, the large drop in δ13C values associated with the end-Permian mass extinction appears to be driven by a sudden transfer of previously stocked 13C depleted carbon, rather than by the overturn of a Late Permian stratified ocean. The Triassic data-set outlines significant secular variations. The best documented is a carbon isotope positive excursion just across the Lower-Middle Triassic boundary, globally developed since it was detected in various paleogeographic settings. It is interpreted to reflect variations in surface ocean chemistry, possibly related to increased primary productivity, at times when the biotic recovery after the end-Permian mass-extinction began to accelerate significantly and when a sharp rise in seawater δ34S values occurred globally. Strontium isotope data obtained from well preserved biogenic phosphates allow a refinement of the Middle Triassic segment of the seawater strontium isotope curve and show a major inflexion point of the seawater strontium isotope curve also near the Lower Triassic - Middle Triassic boundary. These facts suggest that the transition from the Early to the Middle Triassic was a time of revolutionary global change which represented an important step in the evolution of Mesozoic marine environments. A tentative carbon isotope curve for the Upper Permian to Upper Triassic time interval is proposed. Its major features are: ? high but constant δ13C values during the Late Permian ? a sharp drop in δ13C values in the latest Permian ? subsequent recovery of δ13C values ? a short-lived positive excursion across the Early-Middle Triassic boundary ? a gradual rise in δ13C values starting in the Late Ladinian or in the Early Carnian It is foreseen that these fluctuations of the carbon isotope curve may serve as chronostratigraphic markers and further assist in the correlation of Permian and Triassic carbonate deposits.

Family history of cancer: pooled analysis in the International Head and Neck Cancer Epidemiology Consortium.

Alcohol and tobacco consumption are well-recognized risk factors for head and neck cancer (HNC). Evidence suggests that genetic predisposition may also play a role. Only a few epidemiologic studies, however, have considered the relation between HNC risk and family history of HNC and other cancers. We pooled individual-level data across 12 case-control studies including 8,967 HNC cases and 13,627 controls. We obtained pooled odds ratios (OR) using fixed and random effect models and adjusting for potential confounding factors. All statistical tests were two-sided. A family history of HNC in first-degree relatives increased the risk of HNC (OR=1.7, 95% confidence interval, CI, 1.2-2.3). The risk was higher when the affected relative was a sibling (OR=2.2, 95% CI 1.6-3.1) rather than a parent (OR=1.5, 95% CI 1.1-1.8) and for more distal HNC anatomic sites (hypopharynx and larynx). The risk was also higher, or limited to, in subjects exposed to tobacco. The OR rose to 7.2 (95% CI 5.5-9.5) among subjects with family history, who were alcohol and tobacco users. A weak but significant association (OR=1.1, 95% CI 1.0-1.2) emerged for family history of other tobacco-related neoplasms, particularly with laryngeal cancer (OR=1.3, 95% CI 1.1-1.5). No association was observed for family history of nontobacco-related neoplasms and the risk of HNC (OR=1.0, 95% CI 0.9-1.1). Familial factors play a role in the etiology of HNC. In both subjects with and without family history of HNC, avoidance of tobacco and alcohol exposure may be the best way to avoid HNC.

Binary proposal for assessing quality of Open Access Institutional Repositories : the case of Spanish repositories

The reason for this study is to propose a new quantitative approach on how to assess the quality of Open Access University Institutional Repositories. The results of this new approach are tested in the Spanish University Repositories. The assessment method is based in a binary codification of a proposal of features that objectively describes the repositories. The purposes of this method are assessing the quality and an almost automatically system for updating the data of the characteristics. First of all a database was created with the 38 Spanish institutional repositories. The variables of analysis are presented and explained either if they are coming from bibliography or are a set of new variables. Among the characteristics analyzed are the features of the software, the services of the repository, the features of the information system, the Internet visibility and the licenses of use. Results from Spanish universities ARE provided as a practical example of the assessment and for having a picture of the state of the development of the open access movement in Spain.

Preinfarction angina prior to first myocardial infarction does not influence long-term prognosis: a retrospective study with subgroup analysis in elderly and diabetic patients.

BACKGROUND AND HYPOTHESIS Although prodromal angina occurring shortly before an acute myocardial infarction (MI) has protective effects against in-hospital complications, this effect has not been well documented after initial hospitalization, especially in older or diabetic patients. We examined whether angina 1 week before a first MI provides protection in these patients. METHODS A total of 290 consecutive patients, 143 elderly (>64 years of age) and 147 adults (<65 years of age), 68 of whom were diabetic (23.4%) and 222 nondiabetic (76.6%), were examined to assess the effect of preceding angina on long-term prognosis (56 months) after initial hospitalization for a first MI. RESULTS No significant differences were found in long-term complications after initial hospitalization in these adult and elderly patients according to whether or not they had prodromal angina (44.4% with angina vs 45.4% without in adults; 45.5% vs 58% in elderly, P < 0.2). Nor were differences found according to their diabetic status (61.5% with angina vs 72.7% without in diabetics; 37.3% vs 38.3% in nondiabetics; P = 0.4). CONCLUSION The occurrence of angina 1 week before a first MI does not confer long-term protection against cardiovascular complications after initial hospitalization in adult or elderly patients, whether or not they have diabetes.

Gene-lifestyle interaction and type 2 diabetes: the EPIC interact case-cohort study.

BACKGROUND Understanding of the genetic basis of type 2 diabetes (T2D) has progressed rapidly, but the interactions between common genetic variants and lifestyle risk factors have not been systematically investigated in studies with adequate statistical power. Therefore, we aimed to quantify the combined effects of genetic and lifestyle factors on risk of T2D in order to inform strategies for prevention. METHODS AND FINDINGS The InterAct study includes 12,403 incident T2D cases and a representative sub-cohort of 16,154 individuals from a cohort of 340,234 European participants with 3.99 million person-years of follow-up. We studied the combined effects of an additive genetic T2D risk score and modifiable and non-modifiable risk factors using Prentice-weighted Cox regression and random effects meta-analysis methods. The effect of the genetic score was significantly greater in younger individuals (p for interaction  = 1.20×10-4). Relative genetic risk (per standard deviation [4.4 risk alleles]) was also larger in participants who were leaner, both in terms of body mass index (p for interaction  = 1.50×10-3) and waist circumference (p for interaction  = 7.49×10-9). Examination of absolute risks by strata showed the importance of obesity for T2D risk. The 10-y cumulative incidence of T2D rose from 0.25% to 0.89% across extreme quartiles of the genetic score in normal weight individuals, compared to 4.22% to 7.99% in obese individuals. We detected no significant interactions between the genetic score and sex, diabetes family history, physical activity, or dietary habits assessed by a Mediterranean diet score. CONCLUSIONS The relative effect of a T2D genetic risk score is greater in younger and leaner participants. However, this sub-group is at low absolute risk and would not be a logical target for preventive interventions. The high absolute risk associated with obesity at any level of genetic risk highlights the importance of universal rather than targeted approaches to lifestyle intervention.

Determinants of school efficiency : The case of primary schools in the State of Geneva, Switzerland

The public primary school system in the State of Geneva, Switzerland, is characterized by centrally evaluated pupil performance measured with the use of standardized tests. As a result, consistent data are collected among the system. The 2010-2011 dataset is used to develop a two-stage data envelopment analysis (DEA) of school efficiency. In the first stage, DEA is employed to calculate an individual efficiency score for each school. It shows that, on average, each school could reduce its inputs by 7% whilst maintaining the same quality of pupil performance. The cause of inefficiency lies in perfectible management. In the second stage, efficiency is regressed on school characteristics and environmental variables;external factors outside of the control of headteachers. The model is tested for multicollinearity, heteroskedasticity and endogeneity. Four variables are identified as statistically significant. School efficiency is negatively influenced by (1) the provision of special education, (2) the proportion of disadvantaged pupils enrolled at the school and (3) operations being held on multiple sites, but positively influenced by school size (captured by the number of pupils). The proportion of allophone pupils; schools located in urban areas and the provision of reception classes for immigrant pupils are not significant. Although the significant variables influencing school efficiency are outside of the control of headteachers, it is still possible to either boost the positive impact or curb the negative impact. Dans le canton de Genève (Suisse), les écoles publiques primaires sont caractérisées par un financement assuré par les collectivités publiques (canton et communes) et par une évaluation des élèves à l'aide d'épreuves standardisées à trois moments distincts de leur scolarité. Cela permet de réunir des informations statistiques consistantes. La base de données de l'année 2010-2011 est utilisée dans une analyse en deux étapes de l'efficience des écoles. Dans une première étape, la méthode d'analyse des données par enveloppement (DEA) est utilisée pour calculer un score d'efficience pour chaque école. Cette analyse démontre que l'efficience moyenne des écoles s'élève à 93%. Chaque école pourrait, en moyenne, réduire ses ressources de 7% tout en conservant constants les résultats des élèves aux épreuves standardisées. La source de l'inefficience réside dans un management des écoles perfectible. Dans une seconde étape, les scores d'efficience sont régressés sur les caractéristiques des écoles et sur des variables environnementales. Ces variables ne sont pas sous le contrôle (ou l'influence) des directeurs d'école. Le modèle est testé pour la multicolinéartié, l'hétéroscédasticité et l'endogénéité. Quatre variables sont statistiquement significatives. L'efficience des écoles est influencée négativement par (1) le fait d'offrir un enseignement spécialisé en classe séparée, (2) la proporition d'élèves défavorisés et (3) le fait d'opérer sur plusieurs sites différents. L'efficience des écoles est influencée positivement par la taille de l'école, mesurée par le nombre d'élèves. La proporition d'élèves allophones, le fait d'être situé dans une zone urbaine et d'offrir des classes d'accueil pour les élèves immigrants constituent autant de variables non significatives. Le fait que les variables qui influencent l'efficience des écoles ne soient pas sous le contrôle des directeurs ne signifie pas qu'il faille céder au fatalisme. Différentes pistes sont proposées pour permettre soit de réduire l'impact négatif soit de tirer parti de l'impact positif des variables significatives.

Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

A combined analysis of the short-term effects of photochemical air pollutants on mortality within the EMECAM project.

In recent years, some epidemiologic studies have attributed adverse effects of air pollutants on health not only to particles and sulfur dioxide but also to photochemical air pollutants (nitrogen dioxide and ozone). The effects are usually small, leading to some inconsistencies in the results of the studies. Furthermore, the different methodologic approaches of the studies used has made it difficult to derive generic conclusions. We provide here a quantitative summary of the short-term effects of photochemical air pollutants on mortality in seven Spanish cities involved in the EMECAM project, using generalized additive models from analyses of single and multiple pollutants. Nitrogen dioxide and ozone data were provided by seven EMECAM cities (Barcelona, Gijón, Huelva, Madrid, Oviedo, Seville, and Valencia). Mortality indicators included daily total mortality from all causes excluding external causes, daily cardiovascular mortality, and daily respiratory mortality. Individual estimates, obtained from city-specific generalized additive Poisson autoregressive models, were combined by means of fixed effects models and, if significant heterogeneity among local estimates was found, also by random effects models. Significant positive associations were found between daily mortality (all causes and cardiovascular) and NO(2), once the rest of air pollutants were taken into account. A 10 microg/m(3) increase in the 24-hr average 1-day NO(2)level was associated with an increase in the daily number of deaths of 0.43% [95% confidence interval (CI), -0.003-0.86%] for all causes excluding external. In the case of significant relationships, relative risks for cause-specific mortality were nearly twice as much as that for total mortality for all the photochemical pollutants. Ozone was independently related only to cardiovascular daily mortality. No independent statistically significant relationship between photochemical air pollutants and respiratory mortality was found. The results in this study suggest that, given the present levels of photochemical pollutants, people living in Spanish cities are exposed to health risks derived from air pollution.

Spatial analysis of the relationship between mortality from cardiovascular and cerebrovascular disease and drinking water hardness.

Previously published scientific papers have reported a negative correlation between drinking water hardness and cardiovascular mortality. Some ecologic and case-control studies suggest the protective effect of calcium and magnesium concentration in drinking water. In this article we present an analysis of this protective relationship in 538 municipalities of Comunidad Valenciana (Spain) from 1991-1998. We used the Spanish version of the Rapid Inquiry Facility (RIF) developed under the European Environment and Health Information System (EUROHEIS) research project. The strategy of analysis used in our study conforms to the exploratory nature of the RIF that is used as a tool to obtain quick and flexible insight into epidemiologic surveillance problems. This article describes the use of the RIF to explore possible associations between disease indicators and environmental factors. We used exposure analysis to assess the effect of both protective factors--calcium and magnesium--on mortality from cerebrovascular (ICD-9 430-438) and ischemic heart (ICD-9 410-414) diseases. This study provides statistical evidence of the relationship between mortality from cardiovascular diseases and hardness of drinking water. This relationship is stronger in cerebrovascular disease than in ischemic heart disease, is more pronounced for women than for men, and is more apparent with magnesium than with calcium concentration levels. Nevertheless, the protective nature of these two factors is not clearly established. Our results suggest the possibility of protectiveness but cannot be claimed as conclusive. The weak effects of these covariates make it difficult to separate them from the influence of socioeconomic and environmental factors. We have also performed disease mapping of standardized mortality ratios to detect clusters of municipalities with high risk. Further standardization by levels of calcium and magnesium in drinking water shows changes in the maps when we remove the effect of these covariates.

Evaluation of the efficacy and safety of lanreotide in combination with targeted therapies in patients with neuroendocrine tumours in clinical practice: a retrospective cross-sectional analysis.

BACKGROUND Based on the mechanism of action, combining somatostatin analogues (SSAs) with mTOR inhibitors or antiangiogenic agents may provide synergistic effects for the treatment of patients with neuroendocrine tumours (NETs). Herein, we investigate the use of these treatment combinations in clinical practice. METHODS This retrospective cross-sectional analysis of patients with NETs treated with the SSA lanreotide and targeted therapies at 35 Spanish hospitals evaluated the efficacy and safety of lanreotide treatment combinations in clinical practice. The data of 159 treatment combinations with lanreotide in 133 patients was retrospectively collected. RESULTS Of the 133 patients, with a median age of 59.4 (16-83) years, 70 (52.6 %) patients were male, 64 (48.1 %) had pancreatic NET, 23 (17.3 %) had ECOG PS ≥2, 41 (30.8 %) had functioning tumours, 63 (47.7 %) underwent surgery of the primary tumour, 45 (33.8 %) had received prior chemotherapy, and 115 (86.5 %) had received prior SSAs. 115 patients received 1 lanreotide treatment combination and 18 patients received between 2 and 5 combinations. Lanreotide was mainly administered in combination with everolimus (73 combinations) or sunitinib (61 combinations). The probability of being progression-free was 78.5 % (6 months), 68.6 % (12 months) and 57.0 % (18 months) for patients who only received everolimus plus lanreotide (n = 57) and 89.3 % (6 months), 73.0 % (12 months), and 67.4 % (18 months) for patients who only received sunitinib and lanreotide (n = 50). In patients who only received everolimus plus lanreotide the median time-to-progression from the initiation of lanreotide combination treatment was 25.8 months (95 % CI, 11.3, 40.3) and it had not yet been reached among the subgroup of patients only receiving sunitinib plus lanreotide. The safety profile of the combination treatment was comparable to that of the targeted agent alone. CONCLUSIONS The combination of lanreotide and targeted therapies, mainly everolimus and sunitinib, is widely used in clinical practice without unexpected toxicities and suggests efficacy that should be explored in randomized prospective clinical trials.