An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.

CYP2C19 and ABCB1 gene polymorphisms are differently distributed according to ethnicity in the Brazilian general population

Background: Recent studies have reported the clinical importance of CYP2C19 and ABCB1 polymorphisms in an individualized approach to clopidogrel treatment. The aims of this study were to evaluate the frequencies of CYP2C19 and ABCB1 polymorphisms and to identify the clopidogrel-predicted metabolic phenotypes according to ethnic groups in a sample of individuals representative of a highly admixtured population. Methods: One hundred and eighty-three Amerindians and 1,029 subjects of the general population of 4 regions of the country were included. Genotypes for the ABCB1c.C3435T (rs1045642), CYP2C19*2 (rs4244285), CYP2C19*3 (rs4986893), CYP2C19*4 (rs28399504), CYP2C19*5 (rs56337013), and CYP2C19*17 (rs12248560) polymorphisms were detected by polymerase chain reaction followed by high resolution melting analysis. The CYP2C19*3, CYP2C19*4 and CYP2C19*5 variants were genotyped in a subsample of subjects (300 samples randomly selected). Results: The CYP2C19*3 and CYP2C19*5 variant alleles were not detected and the CYP2C19*4 variant allele presented a frequency of 0.3%. The allelic frequencies for the ABCB1c.C3435T, CYP2C19*2 and CYP2C19*17 polymorphisms were differently distributed according to ethnicity: Amerindian (51.4%, 10.4%, 15.8%); Caucasian descent (43.2%, 16.9%, 18.0%); Mulatto (35.9%, 16.5%, 21.3%); and African descent (32.8%, 20.2%, 26.3%) individuals, respectively. As a result, self-referred ethnicity was able to predict significantly different clopidogrel-predicted metabolic phenotypes prevalence even for a highly admixtured population. Conclusion: Our findings indicate the existence of inter-ethnic differences in the ABCB1 and CYP2C19 variant allele frequencies in the Brazilian general population plus Amerindians. This information could help in stratifying individuals from this population regarding clopidogrel-predicted metabolic phenotypes and design more cost-effective programs towards individualization of clopidogrel therapy.

Disequilibrium Coefficient: A Bayesian Perspective

Hardy-Weinberg Equilibrium (HWE) is an important genetic property that populations should have whenever they are not observing adverse situations as complete lack of panmixia, excess of mutations, excess of selection pressure, etc. HWE for decades has been evaluated; both frequentist and Bayesian methods are in use today. While historically the HWE formula was developed to examine the transmission of alleles in a population from one generation to the next, use of HWE concepts has expanded in human diseases studies to detect genotyping error and disease susceptibility (association); Ryckman and Williams (2008). Most analyses focus on trying to answer the question of whether a population is in HWE. They do not try to quantify how far from the equilibrium the population is. In this paper, we propose the use of a simple disequilibrium coefficient to a locus with two alleles. Based on the posterior density of this disequilibrium coefficient, we show how one can conduct a Bayesian analysis to verify how far from HWE a population is. There are other coefficients introduced in the literature and the advantage of the one introduced in this paper is the fact that, just like the standard correlation coefficients, its range is bounded and it is symmetric around zero (equilibrium) when comparing the positive and the negative values. To test the hypothesis of equilibrium, we use a simple Bayesian significance test, the Full Bayesian Significance Test (FBST); see Pereira, Stern andWechsler (2008) for a complete review. The disequilibrium coefficient proposed provides an easy and efficient way to make the analyses, especially if one uses Bayesian statistics. A routine in R programs (R Development Core Team, 2009) that implements the calculations is provided for the readers.

Gastroenterology training in Latin America

Latin America is characterized by ethnic, geographical, cultural, and economic diversity; therefore, training in gastroenterology in the region must be considered in this context. The continent's medical education is characterized by a lack of standards and the volume of research continues to be relatively small. There is a multiplicity of events in general gastroenterology and in sub-disciplines, both at regional and local levels, which ensure that many colleagues have access to information. Medical education programs must be based on a clinical vision and be considered in close contact with the patients. The programs should be properly supervised, appropriately defined, and evaluated on a regular basis. The disparity between the patients' needs, the scarce resources available, and the pressures exerted by the health systems on doctors are frequent cited by those complaining of poor professionalism. Teaching development can play a critical role in ensuring the quality of teaching and learning in universities. Continuing professional development programs activities must be planned on the basis of the doctors' needs, with clearly defined objectives and using proper learning methodologies designed for adults. They must be evaluated and accredited by a competent body, so that they may become the basis of a professional regulatory system. The specialty has made progress in the last decades, offering doctors various possibilities for professional development. The world gastroenterology organization has contributed to the speciality through three distinctive, but closely inter-related, programs: Training Centers, Train-the-Trainers, and Global Guidelines, in which Latin America is deeply involved. (C) 2011 Baishideng. All rights reserved.

The role of the observed tropical convection in the generation of frost events in the southern cone of South America

Based on previous results obtained from observations and linear wave theory analysis, the hypothesis that large-scale patterns can generate extreme cold events in southeast South America through the propagation of remotely excited Rossby waves was already suggested. This work will confirm these findings and extend their analysis through a series of numerical experiments using a primitive equation model where waves are excited by a thermal forcing situated in positions chosen according to observed convection anomalies over the equatorial region. The basic state used for these experiments is a composite of austral winters with maximum and minimum frequency of occurrence of generalized frosts that can affect a large area known as the Wet Pampas located in the central and eastern part of Argentina. The results suggest that stationary Rossby waves may be one important mechanism linking anomalous tropical convection with the extreme cold events in the Wet Pampas. The combination of tropical convection and a specific basic state can generate the right environment to guide the Rossby waves trigged by the tropical forcing towards South America. Depending on the phase of the waves entering the South American continent, they can favour the advection of anomalous wind at low levels from the south carrying cold and dry air over the whole southern extreme of the continent, producing a generalized frost in the Wet Pampa region. On the other hand, when a basic state based on the composites of minimum frosts is used, an anomalous anticyclone over the southern part of the continent generates a circulation with a south-southeast wind which brings maritime air and therefore humidity over the Wet Pampas region, creating negative temperature anomalies only over the northeastern part of the region. Under these conditions even if frosts occur they would not be generalized, as observed for the other basic state with maximum frequency of occurrence of generalized frosts.

Chemical potential and the nature of dark energy: The case of a phantom field

The influence of a possible nonzero chemical potential mu on the nature of dark energy is investigated by assuming that the dark energy is a relativistic perfect simple fluid obeying the equation of state, p=omega rho (omega < 0, constant). The entropy condition, S >= 0, implies that the possible values of omega are heavily dependent on the magnitude, as well as on the sign of the chemical potential. For mu > 0, the omega parameter must be greater than -1 (vacuum is forbidden) while for mu < 0 not only the vacuum but even a phantomlike behavior (omega <-1) is allowed. In any case, the ratio between the chemical potential and temperature remains constant, that is, mu/T=mu(0)/T(0). Assuming that the dark energy constituents have either a bosonic or fermionic nature, the general form of the spectrum is also proposed. For bosons mu is always negative and the extended Wien's law allows only a dark component with omega <-1/2, which includes vacuum and the phantomlike cases. The same happens in the fermionic branch for mu < 0. However, fermionic particles with mu > 0 are permitted only if -1 equation-of-state parameter to be omega <-1/2, a result surprisingly close to the maximal value required to accelerate a Friedmann-Robertson-Walker-type universe dominated by matter and dark energy (omega less than or similar to-10/21).

Dynamics and constraints of the massive graviton dark matter flat cosmologies

We discuss the dynamics of the Universe within the framework of the massive graviton cold dark matter scenario (MGCDM) in which gravitons are geometrically treated as massive particles. In this modified gravity theory, the main effect of the gravitons is to alter the density evolution of the cold dark matter component in such a way that the Universe evolves to an accelerating expanding regime, as presently observed. Tight constraints on the main cosmological parameters of the MGCDM model are derived by performing a joint likelihood analysis involving the recent supernovae type Ia data, the cosmic microwave background shift parameter, and the baryonic acoustic oscillations as traced by the Sloan Digital Sky Survey red luminous galaxies. The linear evolution of small density fluctuations is also analyzed in detail. It is found that the growth factor of the MGCDM model is slightly different (similar to 1-4%) from the one provided by the conventional flat Lambda CDM cosmology. The growth rate of clustering predicted by MGCDM and Lambda CDM models are confronted to the observations and the corresponding best fit values of the growth index (gamma) are also determined. By using the expectations of realistic future x-ray and Sunyaev-Zeldovich cluster surveys we derive the dark matter halo mass function and the corresponding redshift distribution of cluster-size halos for the MGCDM model. Finally, we also show that the Hubble flow differences between the MGCDM and the Lambda CDM models provide a halo redshift distribution departing significantly from the those predicted by other dark energy models. These results suggest that the MGCDM model can observationally be distinguished from Lambda CDM and also from a large number of dark energy models recently proposed in the literature.

Two formation channels of ultra-compact dwarf galaxies in Hickson compact groups

Context. The formation of ultra-compact dwarf galaxies (UCDs) is believed to be driven by interaction, and UCDs are abundant in the cores of galaxy clusters, environments that mark the end-point of galaxy evolution. Nothing is known about the properties of UCDs in compact groups of galaxies, environments where most of galaxy evolution and interaction is believed to occur and where UCDs in an intermediate stage in their evolution may be expected. Aims. The main goal of this study is to detect and characterize, for the first time, the UCD population of compact groups of galaxies. For that, two nearby groups in different evolutionary stages, HCG22 and HCG90, were targeted. Methods. We selected about 40 UCD candidates from pre-existing photometry of both groups, and obtained spectra of these candidates using the VLT FORS2 instrument in MXU mode. Archival HST/ACS imaging was used to measure their structural parameters. Results. We detect 16 and 5 objects belonging to HCG22 and HCG90, respectively, covering the magnitude range -10.0 > M(R) > -11.5 mag. Their integrated colours are consistent with old ages covering a broad range in metallicities (metallicities confirmed by the spectroscopic measurements). Photometric mass estimates put 4 objects in HCG90 and 9 in HCG22 in the mass range of UCDs (> 2 x 10(6) M(circle dot)) for an assumed age of 12Gyr. These UCDs are on average 2-3 times larger than the typical size of Galactic GCs, covering a range of 2 less than or similar to r(h) less than or similar to 21 pc. The UCDs in HCG22 are more concentrated around the central galaxy than in HCG90, at the 99% confidence level. They cover a broad range in [alpha/Fe] abundances from sub-to super-solar. The spectra of 3 UCDs (2 in HCG22, 1 in HCG90) show tentative evidence of intermediate age stellar populations. The clearest example is the largest and most massive UCD (similar to 10(7) M(circle dot)) in our sample, which is detected in HCG22. Its properties are most consistent with a stripped dwarf galaxy nucleus. We calculate the specific frequency (S(N)) of UCDs for both groups, finding that HCG22 has about three times higher S(N) than HCG90. Conclusions. The ensemble properties of the detected UCDs supports two co-existing formation channels: a star cluster origin (low-luminosity, compact sizes, old ages, super-solar alpha/Fe), and an origin as tidally stripped dwarf nuclei (more extended and younger stellar populations). Our results imply that the UCDs detected in both groups do not, in their majority, originate from relatively recent galaxy interactions. Most of the detected UCDs have likely been brought into the group along with their host galaxies.

TESTING THE DISTANCE-DUALITY RELATION WITH GALAXY CLUSTERS AND TYPE Ia SUPERNOVAE

In this Letter, we propose a new and model-independent cosmological test for the distance-duality (DD) relation, eta = D(L)(z)(1 + z)(-2)/D(A)(z) = 1, where D(L) and D(A) are, respectively, the luminosity and angular diameter distances. For D(L) we consider two sub-samples of Type Ia supernovae (SNe Ia) taken from Constitution data whereas D(A) distances are provided by two samples of galaxy clusters compiled by De Filippis et al. and Bonamente et al. by combining Sunyaev-Zeldovich effect and X-ray surface brightness. The SNe Ia redshifts of each sub-sample were carefully chosen to coincide with the ones of the associated galaxy cluster sample (Delta z < 0.005), thereby allowing a direct test of the DD relation. Since for very low redshifts, D(A)(z) approximate to D(L)(z), we have tested the DD relation by assuming that. is a function of the redshift parameterized by two different expressions: eta(z) = 1 + eta(0)z and eta(z) = 1 +eta(0)z/(1 + z), where eta(0) is a constant parameter quantifying a possible departure from the strict validity of the reciprocity relation (eta(0) = 0). In the best scenario (linear parameterization), we obtain eta(0) = -0.28(-0.44)(+0.44) (2 sigma, statistical + systematic errors) for the De Filippis et al. sample (elliptical geometry), a result only marginally compatible with the DD relation. However, for the Bonamente et al. sample (spherical geometry) the constraint is eta(0) = -0.42(-0.34)(+0.34) (3 sigma, statistical + systematic errors), which is clearly incompatible with the duality-distance relation.

Transiting exoplanets from the CoRoT space mission XI. CoRoT-8b: a hot and dense sub-Saturn around a K1 dwarf

Aims. We report the discovery of CoRoT-8b, a dense small Saturn-class exoplanet that orbits a K1 dwarf in 6.2 days, and we derive its orbital parameters, mass, and radius. Methods. We analyzed two complementary data sets: the photometric transit curve of CoRoT-8b as measured by CoRoT and the radial velocity curve of CoRoT-8 as measured by the HARPS spectrometer**. Results. We find that CoRoT-8b is on a circular orbit with a semi-major axis of 0.063 +/- 0.001 AU. It has a radius of 0.57 +/- 0.02 R(J), a mass of 0.22 +/- 0.03 M(J), and therefore a mean density of 1.6 +/- 0.1 g cm(-3). Conclusions. With 67% of the size of Saturn and 72% of its mass, CoRoT-8b has a density comparable to that of Neptune (1.76 g cm(-3)). We estimate its content in heavy elements to be 47-63 M(circle plus), and the mass of its hydrogen-helium envelope to be 7-23 M(circle plus). At 0.063 AU, the thermal loss of hydrogen of CoRoT-8b should be no more than similar to 0.1% over an assumed integrated lifetime of 3 Ga.

Age distribution of the central stars of galactic disk planetary nebulae

Context. Determination of the ages of central stars of planetary nebulae (CSPN) is a complex problem, and there is presently no single method that can be generally applied. We have developed several methods of estimating the ages of CSPN, based on both the observed nebular properties and some properties of the stars themselves. Aims. Our aim is to estimate the ages and the age distribution of CSPN and to compare the derived results with mass and age determinations of CSPN and white dwarfs based on empirical determinations of these quantities. Methods. We considered a sample of planetary nebulae in the galactic disk, most of which (similar to 69%) are located in the solar neighbourhood, within 3 kpc from the Sun. We discuss several methods of deriving the age distribution of CSPN, namely; (i) the use of an age-metallicity relation that also depends on the galactocentric distance; (ii) the use of an age-metallicity relation obtained for the galactic disk; and (iii) the determination of ages from the central star masses obtained from the observed nitrogen abundances. Results. We estimated the age distribution of CSPN with average uncertainties of 1-2 Gyr, and compared our results with the expected distribution based both on the observed mass distribution of white dwarfs and on the age distribution derived from available mass distributions of CSPN. Based on our derived age distributions, we conclude that most CSPN in the galactic disk have ages under 6 Gyr, and that the age distribution is peaked around 2-4 Gyr.

Kinematic analysis and membership status of TWA22 AB

Context. TWA22 was initially regarded as a member of the TW Hydrae association (TWA). In addition to being one of the youngest (approximate to 8 Myr) and nearest (approximate to 20 pc) stars to Earth, TWA22 has proven to be very interesting after being resolved as a tight, very low-mass binary. This binary can serve as a very useful dynamical calibrator for pre-main sequence evolutionary models. However, its membership in the TWA has been recently questioned despite due to the lack of accurate kinematic measurements. Aims. Based on proper motion, radial velocity, and trigonometric parallax measurements, we aim here to re-analyze the membership of TWA22 to young, nearby associations. Methods. Using the ESO NTT/SUSI2 telescope, we observed TWA22 AB during 5 different observing runs over 1.2 years to measure its trigonometric parallax and proper motion. This is a part of a larger project measuring trigonometric parallaxes and proper motions of most known TWA members at a sub-milliarcsec level. HARPS at the ESO 3.6 m telescope was also used to measure the system's radial velocity over 2 years. Results. We report an absolute trigonometric parallax of TWA22 AB, pi = 57.0 +/- 0.7 mas, corresponding to a distance 17.5 +/- 0.2 pc from Earth. Measured proper motions of TWA 22AB are mu(alpha) cos(delta) = -175.8 +/- 0.8 mas/yr and mu delta = -21.3 +/- 0.8 mas/yr. Finally, from HARPS measurements, we obtain a radial velocity V(rad) = 14.8 +/- 2.1 km s(-1). Conclusions. A kinematic analysis of TWA22 AB space motion and position implies that a membership of TWA22 AB to known young, nearby associations can be excluded except for the beta Pictoris and TW Hydrae associations. Membership probabilities based on the system's Galactic space motion and/or the trace-back technique support a higher chance of being a member to the beta Pictoris association. Membership of TWA22 in the TWA cannot be fully excluded because of large uncertainties in parallax measurements and radial velocities and to the uncertain internal velocity dispersion of its members.

DERIVING METALLICITIES FROM THE INTEGRATED SPECTRA OF EXTRAGALACTIC GLOBULAR CLUSTERS USING THE NEAR-INFRARED CALCIUM TRIPLET

The Ca II triplet (CaT) feature in the near-infrared has been employed as a metallicity indicator for individual stars as well as integrated light of Galactic globular clusters (GCs) and galaxies with varying degrees of success, and sometimes puzzling results. Using the DEIMOS multi-object spectrograph on Keck we obtain a sample of 144 integrated light spectra of GCs around the brightest group galaxy NGC 1407 to test whether the CaT index can be used as ametallicity indicator for extragalactic GCs. Different sets of single stellar population models make different predictions for the behavior of the CaT as a function of metallicity. In this work, the metallicities of the GCs around NGC 1407 are obtained from CaT index values using an empirical conversion. The measured CaT/metallicity distributions show unexpected features, the most remarkable being that the brightest red and blue GCs have similar CaT values despite their large difference in mean color. Suggested explanations for this behavior in the NGC 1407 GC system are (1) the CaT may be affected by a population of hot blue stars, (2) the CaT may saturate earlier than predicted by the models, and/or (3) color may not trace metallicity linearly. Until these possibilities are understood, the use of the CaT as a metallicity indicator for the integrated spectra of extragalactic GCs will remain problematic.

Kinematic constraints to the transition redshift from supernovae type Ia union data

The kinematic approach to cosmological tests provides direct evidence to the present accelerating stage of the Universe that does not depend on the validity of general relativity, as well as on the matter-energy content of the Universe. In this context, we consider here a linear two-parameter expansion for the decelerating parameter, q(z)=q(0)+q(1)z, where q(0) and q(1) are arbitrary constants to be constrained by the union supernovae data. By assuming a flat Universe we find that the best fit to the pair of free parameters is (q(0),q(1))=(-0.73,1.5) whereas the transition redshift is z(t)=0.49(-0.07)(+0.14)(1 sigma) +0.54-0.12(2 sigma). This kinematic result is in agreement with some independent analyses and more easily accommodates many dynamical flat models (like Lambda CDM).

Leishmania-Specific Surface Antigens Show Sub-Genus Sequence Variation and Immune Recognition

Background: A family of hydrophilic acylated surface (HASP) proteins, containing extensive and variant amino acid repeats, is expressed at the plasma membrane in infective extracellular (metacyclic) and intracellular (amastigote) stages of Old World Leishmania species. While HASPs are antigenic in the host and can induce protective immune responses, the biological functions of these Leishmania-specific proteins remain unresolved. Previous genome analysis has suggested that parasites of the sub-genus Leishmania (Viannia) have lost HASP genes from their genomes. Methods/Principal Findings: We have used molecular and cellular methods to analyse HASP expression in New World Leishmania mexicana complex species and show that, unlike in L. major, these proteins are expressed predominantly following differentiation into amastigotes within macrophages. Further genome analysis has revealed that the L. (Viannia) species, L. (V.) braziliensis, does express HASP-like proteins of low amino acid similarity but with similar biochemical characteristics, from genes present on a region of chromosome 23 that is syntenic with the HASP/SHERP locus in Old World Leishmania species and the L. (L.) mexicana complex. A related gene is also present in Leptomonas seymouri and this may represent the ancestral copy of these Leishmania-genus specific sequences. The L. braziliensis HASP-like proteins (named the orthologous (o) HASPs) are predominantly expressed on the plasma membrane in amastigotes and are recognised by immune sera taken from 4 out of 6 leishmaniasis patients tested in an endemic region of Brazil. Analysis of the repetitive domains of the oHASPs has shown considerable genetic variation in parasite isolates taken from the same patients, suggesting that antigenic change may play a role in immune recognition of this protein family. Conclusions/Significance: These findings confirm that antigenic hydrophilic acylated proteins are expressed from genes in the same chromosomal region in species across the genus Leishmania. These proteins are surface-exposed on amastigotes (although L. (L.) major parasites also express HASPB on the metacyclic plasma membrane). The central repetitive domains of the HASPs are highly variant in their amino acid sequences, both within and between species, consistent with a role in immune recognition in the host.