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.

Frequency of LCT-13910C > T single nucleotide polymorphism associated with adult-type hypolactasia/lactase persistence among Brazilians of different ethnic groups

Background: Adult-type hypolactasia, the physiological decline of lactase some time after weaning, was previously associated with the LCT -13910C>T polymorphism worldwide except in Africa. Lactase non-persistence is the most common phenotype in humans, except in northwestern Europe with its long history of pastoralism and milking. We had previously shown association of LCT -13910C>T polymorphism with adult-type hypolactasia in Brazilians; thus, we assessed its frequency among different Brazilian ethnic groups. Methods: We investigated the ethnicity-related frequency of this polymorphism in 567 Brazilians [mean age, 42.1 +/- 16.8 years; 157 (27.7%) men]; 399 (70.4%) White, 50 (8.8%) Black, 65 (11.5%) Brown, and 53 (9.3%) Japanese-Brazilian. DNA was extracted from leukocytes; LCT -13910C>T polymorphism was analyzed by PCR-restriction fragment length polymorphism. Results: Prevalence of the CC genotype associated with hypolactasia was similar (57%) among White and Brown groups; however, prevalence was higher among Blacks (80%) and those of Japanese descent (100%). Only 2 (4%) Blacks had TT genotype, and 8 (16%) had the CT genotype. Assuming an association between CC genotype and hypolactasia, and CT and TT genotypes with lactase persistence, 356 (62.8%) individuals had hypolactasia and 211 (37.2%) had lactase persistence. The White and Brown groups had the same hypolactasia prevalence (similar to 57%); nevertheless, was 80% among Black individuals and 100% among Japanese-Brazilians (P < 0.01). Conclusion: The lactase persistence allele, LCT -13910T, was found in about 43% of both White and Brown and 20% of the Black Brazilians, but was absent among all Japanese Brazilians studied.

Variations in Human Herpesvirus Type 8 Seroprevalence in Native Americans, South America

To determine the epidemiology of human herpesvirus type 8 (HHV-8) among non-Amazonian native populations, we conducted a cross-sectional study in Brazil, Bolivia, and Paraguay. Our data show striking ethnic and geographic variations in the distribution of HHV-8 seroprevalences in Amazonian (77%) and non-Amazonian native populations (range 0%-83%).

Age of Child, More than HPV Type, Is Associated with Clinical Course in Recurrent Respiratory Papillomatosis

Background: RRP is a devastating disease in which papillomas in the airway cause hoarseness and breathing difficulty. The disease is caused by human papillomavirus (HPV) 6 or 11 and is very variable. Patients undergo multiple surgeries to maintain a patent airway and in order to communicate vocally. Several small studies have been published in which most have noted that HPV 11 is associated with a more aggressive course. Methodology/Principal Findings: Papilloma biopsies were taken from patients undergoing surgical treatment of RRP and were subjected to HPV typing. 118 patients with juvenile-onset RRP with at least 1 year of clinical data and infected with a single HPV type were analyzed. HPV 11 was encountered in 40% of the patients. By our definition, most of the patients in the sample (81%) had run an aggressive course. The odds of a patient with HPV 11 running an aggressive course were 3.9 times higher than that of patients with HPV 6 (Fisher's exact p = 0.017). However, clinical course was more closely associated with age of the patient (at diagnosis and at the time of the current surgery) than with HPV type. Patients with HPV 11 were diagnosed at a younger age (2.4y) than were those with HPV 6 (3.4y) (p = 0.014). Both by multiple linear regression and by multiple logistic regression HPV type was only weakly associated with metrics of disease course when simultaneously accounting for age. Conclusions/Significance Abstract: The course of RRP is variable and a quarter of the variability can be accounted for by the age of the patient. HPV 11 is more closely associated with a younger age at diagnosis than it is associated with an aggressive clinical course. These data suggest that there are factors other than HPV type and age of the patient that determine disease course.

Genetic Diversity on the Integrase Region of the pol Gene among HIV Type 1-Infected Patients Naive for Integrase Inhibitors in Sao Paulo City, Brazil

The presence of mutations associated with integrase inhibitor (INI) resistance among INI-naive patients may play an important clinical role in the use of those drugs Samples from 76 HIV-1-infected subjects naive to INIs were submitted to direct sequencing. No differences were found between naive (25%) subjects and subjects on HAART (75%). No primary mutation associated with raltegravir or elvitegravir resistance was found. However, 78% of sequences showed at least one accessory mutation associated with resistance. The analysis of the 76 IN sequences showed a high polymorphic level on this region among Brazilian HIV-1-infected subjects, including a high prevalence of aa substitutions related to INI resistance. The impact of these findings remains unclear and further studies are necessary to address these questions.

Abnormal collagen deposition in synovia after collagen type V immunization in rabbits

Sinovitis in Scleroderma (SSc) is rare, usually aggressive and fully resembles rheumatoid arthritis. Experimental models of SSc have been used in an attempt to understand its pathogenesis. Previous studies done in our laboratory had already revealed the presence of a synovial remodeling process in rabbits immunized with collagen V. To validate the importance of collagen type V and to explore the quantitative relationship between this factor and synovia remodeling as well as the relationship between collagen type V and other collagens, we studied the synovial tissue in immunized rabbits. Rabbits (N= 10) were immunized with collagen V plus Freund's adjuvant and compared with animals inoculated with adjuvant only (N= 10). Synovial tissues were submitted to histological analysis, immunolocalization to collagen I, III and V and biochemical analysis by eletrophoresis, immunoblot and densitometric method. The synovial tissue presented an intense remodeling process with deposits of collagen types I, III and V after 75 and 120 days of immunization, mainly distributed around the vessels and interstitium of synovial extracellular matrix. Densitometric analysis confirmed the increased synthesis of collagen I, III and V chains (407.69 +/- 80.31; 24.46 +/- 2.58; 70.51 +/- 7.66, respectively) in immunized rabbits when compared with animals from control group (164.91 +/- 15.67; 12.89 +/- 1.05; 32 +/- 3.57) (p<0.0001). We conclude that synovial remodeling observed in the experimental model can reflect the articular compromise present in patients with scleroderma. Certainly, this experimental model induced by collagen V immunization will bring new insights in to pathogenic mechanisms and allow the testing of new therapeutic strategies to ameliorate the prognosis for scleroderma patients.

Genotypic Characteristics of HIV Type 1 Based on gp120 Hypervariable Region 3 of Isolates from Southern Brazil

The aim of this study was to investigate HIV-1 molecular diversity and the epidemiological profile of HIV-1-infected patients from Ribeirao Preto, Brazil. A nested PCR followed by sequencing of a 302-base pair fragment of the env gene (C2-V3 region) was performed in samples from HIV-1-positive patients. A total of 45 sequences were aligned with final manual adjustments. The phylogenetic analyses showed a higher prevalence of HIV-1 subtype B in the studied population (97.8%) with only one sample yielding an F1 subtype. The viral genotyping prediction showed that CCR5 tropism was the most prevalent in the studied cohort. Geno2pheno analysis showed that R5 and CXCR4 prediction were 69% and 31%, respectively. There was no statistical significance, either in viral load or in CD4(+) T cell count when R5 and X4 prediction groups were compared. Moreover, the GPGR tetramer was the most common V3 loop core motif identified in the HIV-1 strains studied (34.1%) followed by GWGR, identified in 18.1% of the samples. The high level of B subtype in this Brazilian population reinforces the nature of the HIV epidemic in Brazil, and corroborates previous data obtained in the Brazilian HIV-infected population.

Antiretroviral resistance among HIV type 1-infected women first exposed to antiretrovirals during pregnancy: Plasma versus PBMCs

Resistance-associated mutations (RAMs) in plasma samples from HIV-1-infected women who received antiretroviral (ARV) prophylaxis during pregnancy was assessed and correlated with the detection of RAMs in peripheral blood mononuclear cells (PMBCs). The study population was composed of HIV-1-infected women enrolled in a prospective cohort study in Latin America and the Caribbean (NISDI Perinatal Study) as of March 1, 2005, who were diagnosed with HIV-1 infection during the current pregnancy, who received ARVs during pregnancy for prevention of mother-to-child transmission of HIV-1, and who were followed through at least the 6-12 week postpartum visit. Plasma samples collected at enrollment during pregnancy and at 6-12 weeks postpartum were assayed for RAMs. Plasma results were compared to previously described PBMC results from the same study population. Of 819 enrolled subjects, 197 met the eligibility criteria. Nucleic acid amplification was accomplished in 123 plasma samples at enrollment or 6-12 weeks postpartum, and RAMs were detected in 22 (17.9%; 95% CI: 11.7-25.9%). Previous analyses had demonstrated detection of RAMs in PBMCs in 19 (16.1%). There was high concordance between RAMs detected in plasma and PBMC samples, with only eight discordant pairs. The prevalence of RAMs among these pregnant, HIV-1-infected women is high (>15%). Rates of detection of RAMs in plasma and PBMC samples were similar.

The role of major histocompatibility complex Alleles in the susceptibility of Brazilian patients to develop the myogenic type of Graves' orbitopathy

Objective: To assess the frequency of the genetic markers HLA-DRB1 and DQB1 in patients with Graves' orbitopathy ( GO) with and without extraocular muscle involvement. Design: The frequencies of class II HLA-DRB1 and DQB1 allele groups were determined for 81 Brazilian patients with GO and 161 normal subjects. The patients were divided into myogenic and nonmyogenic groups based on the clinical characteristics of the orbitopathy and quantitative computed tomography analysis of the extraocular muscle ( EOM) dimensions. Main outcome: Compared to the frequency obtained for samples of normal subjects of the Brazilian population, HLA-DRB1*16 (p(c)= 0.008) was overrepresented in myogenic and HLA-DRB1*03 (p(c)= 0.02) in nonmyogenic patients. Conclusions: The association between the HLA-DRB1* 16 and the myogenic subtype of GO suggests that EOM involvement in GO may be genetically predisposed.

Primary Hyperparathyroidism as the First Clinical Manifestation of Multiple Endocrine Neoplasia Type 2A in a 5-Year-Old Child

Background: Primary hyperparathyroidism occurs in only 10%-30% of patients with multiple endocrine neoplasia type 2A (MEN2A), rarely as the sole clinical manifestation, and is usually diagnosed after the third decade of life. Summary: A5-year-old girl was referred for prophylactic thyroidectomy as she carried the p.C634R RET mutation. She was clinically asymptomatic, with a normally palpable thyroid and with the cervical region free of lymphadenopathy or other nodules. Preoperative tests revealed hypercalcemia associated with elevation of parathyroid hormone (PTH) (calcium = 11.2mg/dL, calcium ion = 1.48mmol/L, phosphorus = 4.0 mg/dL, alkaline phosphatase = 625U/L, parathyroid hormone (PTH) PTH = 998 pg/mL). A thyroid ultrasound was normal and parathyroid scintigraphy with (99m)Tc-Sestamibi revealed an area of radioconcentration in the upper half of the left thyroid lobe suggesting hyperfunctioning parathyroid tissue. She underwent total thyroidectomy and parathyroidectomy and developed hypocalcemia. The anatomopathological examination showed no histopathological changes in the thyroid tissue and an adenoma of the parathyroid gland, confirming the diagnosis of hyperparathyroidism. Conclusions: Primary hyperparathyroidism can be a precocious manifestation of MEN2A. This case report highlights that asymptomatic hypercalcemia should be scrutinized in children related to patients with MEN2A who carry a mutation in the RET proto-oncogene, especially mutations in the codon 634, before the currently recommended age of 8 years.

Determination and Molecular Analysis of the Complete Genome Sequence of Two Wild-Type Rabies Viruses Isolated from a Haematophagous Bat and a Frugivorous Bat in Brazil

The complete genome sequences of two Brazilian wild-type rabies viruses (RABV), a BR-DR1 isolate from a haematophagous bat (Desmodus rotundus) and a BR-AL1 isolate from a frugivorous bat (Artibeus lituratus), were determined. The genomes of the BR-DR1 and RR-AL1 had 11,923 and 11,922 nt, respectively, and both encoded the five standard genes of rhabdoviruses. The complete nucleotide sequence identity between the BR-DR1 and BR-AL1 isolates was 97%. The BR-DR1 and BR-AL1 isolates had some conserved functional sites revealed by the fixed isolates, whereas both isolates had unique amino acid substitutions in the antigenic region IV of the nucleocapsid gene. Therefore, it is speculated that both isolates were nearly identical in virologic character. According to our phylogenetic analysis based on the complete genomes, both isolates belonged to genotype 1, and to the previously defined ""vampire bat-related RABV lineage"" which consisted of mainly D. rotundus- and A. lituratus- isolates; however, a branch pattern with high bootstrap values suggested that BR-DR1 was more closely related to the 9001FRA isolate, which was collected from a dog bitten by a bat in French Guiana, than to BR-AL1. This result suggests that the vampire bat-related RABV lineage includes Brazilian vampire bat and Brazilian frugivorous bat RABV and is further divided into Brazilian vampire bat and Brazilian frugivorous bat RABV sub-lineages. The phylogenetic analysis based on the complete genomes was valuable in discriminating among very closely related isolates.

Low-Level Laser Intensity Application in Masseter Muscle for Treatment Purposes

Objective: This study evaluated with histochemical analysis how the number of laser applications can affect the masseter muscle. Background: In dentistry today, the laser is used in patients with temporomandibular disorders (TMDs), mainly for radiating pain in the masticatory muscles, whose origins may be associated with malocclusion, although the laser effects are not well understood on the cellular level. Materials and Methods: Thirty mice (HRS/J lineage) were randomly distributed into groups according to the number of laser applications (three, six, and 10). For each group of laser applications (experimental, n = 5), it was considered the control group (n = 5), which was not irradiated. All animals inhaled halothane (2-bromo-2-chloro-1, 1, 1-trifluoroethane, minimum 99%, Sigma Aldrich, India) before each laser irradiation performed on the left masseter muscle region, on alternate days with 20 J/cm(2), 40mW, for 20 sec. The muscle samples were collected for histochemical analysis with succinate dehydrogenase (SDH) enzyme 72 h after the last application. Results: (a) A decrease in area of light fibers type (35.91% +/- 6.9%; 32.08% +/- 6.3%, and 27.88% +/- 6.3%), according to the increase of laser applications (p < 0.05); (b) significant increase (p < 0.05) in the area of intermediate fibers, with an increase of laser application (11.08% +/- 3.9%; 16.52% +/- 5.7%, and 15.96% +/- 3.9%), although the increase with 10 applications was small; (c) area increase of dark fibers in the group with three laser applications (0.16% +/- 0.3%) (p < 0.05), and in groups with six and 10 laser applications, respectively (9.68% +/- 6.0% and 9.60% +/- 4.0%). Conclusions: The SDH enzyme activity revealed that the number of laser applications increases the metabolic pattern of the muscle fibers. A minimal difference in metabolic activity between six and 10 applications of a laser suggests that further analyses should be done to confirm that six applications are enough to produce the same clinical effects, thereby contributing data to professionals from different fields in regard to the cost-benefit ratio of this therapy.

Transiting exoplanets from the CoRoT space mission XIV. CoRoT-11b: a transiting massive ""hot-Jupiter"" in a prograde orbit around a rapidly rotating F-type star

The CoRoT exoplanet science team announces the discovery of CoRoT-11b, a fairly massive hot-Jupiter transiting a V = 12.9 mag F6 dwarf star (M(*) = 1.27 +/- 0.05 M(circle dot), R(*) = 1.37 +/- 0.03 R(circle dot), T(eff) = 6440 +/- 120 K), with an orbital period of P = 2.994329 +/- 0.000011 days and semi-major axis a = 0.0436 +/- 0.005 AU. The detection of part of the radial velocity anomaly caused by the Rossiter-McLaughlin effect shows that the transit-like events detected by CoRoT are caused by a planet-sized transiting object in a prograde orbit. The relatively high projected rotational velocity of the star (upsilon sin i(star) = 40 +/- 5 km s(-1)) places CoRoT-11 among the most rapidly rotating planet host stars discovered so far. With a planetary mass of M(p) = 2.33 +/- 0.34 M(Jup) and radius R(p) = 1.43 +/- 0.03 R(Jup), the resulting mean density of CoRoT-11b (rho(p) = 0.99 +/- 0.15 g/cm(3)) can be explained with a model for an inflated hydrogen-planet with a solar composition and a high level of energy dissipation in its interior.

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.

Pulsations in the late-type Be star HD 50 209 detected by CoRoT

Context. The presence of pulsations in late-type Be stars is still a matter of controversy. It constitutes an important issue to establish the relationship between non-radial pulsations and the mass-loss mechanism in Be stars. Aims. To contribute to this discussion, we analyse the photometric time series of the B8IVe star HD 50 209 observed by the CoRoT mission in the seismology field. Methods. We use standard Fourier techniques and linear and non-linear least squares fitting methods to analyse the CoRoT light curve. In addition, we applied detailed modelling of high-resolution spectra to obtain the fundamental physical parameters of the star. Results. We have found four frequencies which correspond to gravity modes with azimuthal order m = 0,-1,-2,-3 with the same pulsational frequency in the co-rotating frame. We also found a rotational period with a frequency of 0.679 cd(-1) (7.754 mu Hz). Conclusions. HD 50 209 is a pulsating Be star as expected from its position in the HR diagram, close to the SPB instability strip.