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.

Real-Time PCR in HIV/Trypanosoma cruzi Coinfection with and without Chagas Disease Reactivation: Association with HIV Viral Load and CD4(+) Level

Background: Reactivation of chronic Chagas disease, which occurs in approximately 20% of patients coinfected with HIV/Trypanosoma cruzi (T. cruzi), is commonly characterized by severe meningoencephalitis and myocarditis. The use of quantitative molecular tests to monitor Chagas disease reactivation was analyzed. Methodology: Polymerase chain reaction (PCR) of kDNA sequences, competitive (C-) PCR and real-time quantitative (q) PCR were compared with blood cultures and xenodiagnosis in samples from 91 patients (57 patients with chronic Chagas disease and 34 with HIV/T. cruzi coinfection), of whom 5 had reactivation of Chagas disease and 29 did not. Principal Findings: qRT-PCR showed significant differences between groups; the highest parasitemia was observed in patients infected with HIV/T. cruzi with Chagas disease reactivation (median 1428.90 T. cruzi/mL), followed by patients with HIV/T. cruzi infection without reactivation (median 1.57 T. cruzi/mL) and patients with Chagas disease without HIV (median 0.00 T. cruzi/mL). Spearman's correlation coefficient showed that xenodiagnosis was correlated with blood culture, C-PCR and qRT-PCR. A stronger Spearman correlation index was found between C-PCR and qRT-PCR, the number of parasites and the HIV viral load, expressed as the number of CD4(+) cells or the CD4(+)/CD8(+) ratio. Conclusions: qRT-PCR distinguished the groups of HIV/T. cruzi coinfected patients with and without reactivation. Therefore, this new method of qRT-PCR is proposed as a tool for prospective studies to analyze the importance of parasitemia (persistent and/or increased) as a criterion for recommending pre-emptive therapy in patients with chronic Chagas disease with HIV infection or immunosuppression. As seen in this study, an increase in HIV viral load and decreases in the number of CD4(+) cells/mm(3) and the CD4(+)/CD8(+) ratio were identified as cofactors for increased parasitemia that can be used to target the introduction of early, pre-emptive therapy.

Spatio-Temporal Tracking and Phylodynamics of an Urban Dengue 3 Outbreak in Sao Paulo, Brazil

The dengue virus has a single-stranded positive-sense RNA genome of similar to 10.700 nucleotides with a single open reading frame that encodes three structural (C, prM, and E) and seven nonstructural (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) proteins. It possesses four antigenically distinct serotypes (DENV 1-4). Many phylogenetic studies address particularities of the different serotypes using convenience samples that are not conducive to a spatio-temporal analysis in a single urban setting. We describe the pattern of spread of distinct lineages of DENV-3 circulating in Sao Jose do Rio Preto, Brazil, during 2006. Blood samples from patients presenting dengue-like symptoms were collected for DENV testing. We performed M-N-PCR using primers based on NS5 for virus detection and identification. The fragments were purified from PCR mixtures and sequenced. The positive dengue cases were geo-coded. To type the sequenced samples, 52 reference sequences were aligned. The dataset generated was used for iterative phylogenetic reconstruction with the maximum likelihood criterion. The best demographic model, the rate of growth, rate of evolutionary change, and Time to Most Recent Common Ancestor (TMRCA) were estimated. The basic reproductive rate during the epidemics was estimated. We obtained sequences from 82 patients among 174 blood samples. We were able to geo-code 46 sequences. The alignment generated a 399-nucleotide-long dataset with 134 taxa. The phylogenetic analysis indicated that all samples were of DENV-3 and related to strains circulating on the isle of Martinique in 2000-2001. Sixty DENV-3 from Sao Jose do Rio Preto formed a monophyletic group (lineage 1), closely related to the remaining 22 isolates (lineage 2). We assumed that these lineages appeared before 2006 in different occasions. By transforming the inferred exponential growth rates into the basic reproductive rate, we obtained values for lineage 1 of R(0) = 1.53 and values for lineage 2 of R(0) = 1.13. Under the exponential model, TMRCA of lineage 1 dated 1 year and lineage 2 dated 3.4 years before the last sampling. The possibility of inferring the spatio-temporal dynamics from genetic data has been generally little explored, and it may shed light on DENV circulation. The use of both geographic and temporally structured phylogenetic data provided a detailed view on the spread of at least two dengue viral strains in a populated urban area.

High levels of T lymphocyte activation in Leishmania-HIV-1 co-infected individuals despite low HIV viral load

Background: Concomitant infections may influence HIV progression by causing chronic activation leading to decline in T-cell function. In the Americas, visceral (AVL) and tegumentary leishmaniasis (ATL) have emerged as important opportunistic infections in HIV-AIDS patients and both of those diseases have been implicated as potentially important co-factors in disease progression. We investigated whether leishmaniasis increases lymphocyte activation in HIV-1 co-infected patients. This might contribute to impaired cellular immune function. Methods: To address this issue we analyzed CD4(+) T absolute counts and the proportion of CD8(+) T cells expressing CD38 in Leishmania/HIV co-infected patients that recovered after anti-leishmanial therapy. Results: We found that, despite clinical remission of leishmaniasis, AVL co-infected patients presented a more severe immunossupression as suggested by CD4(+) T cell counts under 200 cells/mm(3), differing from ATL/HIV-AIDS cases that tends to show higher lymphocytes levels (over 350 cells/mm(3)). Furthermore, five out of nine, AVL/HIV-AIDS presented low CD4(+) T cell counts in spite of low or undetectable viral load. Expression of CD38 on CD8(+) T lymphocytes was significantly higher in AVL or ATL/HIV-AIDS cases compared to HIV/AIDS patients without leishmaniasis or healthy subjects. Conclusions: Leishmania infection can increase the degree of immune system activation in individuals concomitantly infected with HIV. In addition, AVL/HIV-AIDS patients can present low CD4(+) T cell counts and higher proportion of activated T lymphocytes even when HIV viral load is suppressed under HAART. This fact can cause a misinterpretation of these laboratorial markers in co-infected patients.

Near full-length genome analysis of low prevalent human immunodeficiency virus type 1 subclade F1 in Sao Paulo, Brazil

Background: The genetic diversity of the human immunodeficiency virus type 1 (HIV-1) is critical to lay the groundwork for the design of successful drugs or vaccine. In this study we aimed to characterize and define the molecular prevalence of HIV-1 subclade F1 currently circulating in Sao Paulo, Brazil. Methods: A total of 36 samples were selected from 888 adult patients residing in Sao Paulo who had previously been diagnosed in two independent studies in our laboratory as being infected with subclade F1 based on pol subgenomic fragment sequencing. Proviral DNA was amplified from the purified genomic DNA of all 36 blood samples by 5 fragments overlapping PCR followed by direct sequencing. Sequence data were obtained from the 5 fragments of pure subclade F1 and phylogenetic trees were constructed and compared with previously published sequences. Subclades F1 that exhibited mosaic structure with other subtypes were omitted from any further analysis Results: Our methods of fragment amplification and sequencing confirmed that only 5 sequences inferred from pol region as subclade F1 also holds true for the genome as a whole and, thus, estimated the true prevalence at 0.56%. The results also showed a single phylogenetic cluster of the Brazilian subclade F1 along with non-Brazilian South American isolates in both subgenomic and the full-length genomes analysis with an overall intrasubtype nucleotide divergence of 6.9%. The nucleotide differences within the South American and Central African F1 strains, in the C2-C3 env, were 8.5% and 12.3%, respectively. Conclusion: All together, our findings showed a surprisingly low prevalence rate of subclade F1 in Brazil and suggest that these isolates originated in Central Africa and subsequently introduced to South America.

An assessment of population structure in eight breeds of cattle using a whole genome SNP panel

Background: Analyses of population structure and breed diversity have provided insight into the origin and evolution of cattle. Previously, these studies have used a low density of microsatellite markers, however, with the large number of single nucleotide polymorphism markers that are now available, it is possible to perform genome wide population genetic analyses in cattle. In this study, we used a high-density panel of SNP markers to examine population structure and diversity among eight cattle breeds sampled from Bos indicus and Bos taurus. Results: Two thousand six hundred and forty one single nucleotide polymorphisms ( SNPs) spanning all of the bovine autosomal genome were genotyped in Angus, Brahman, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black, Limousin and Nelore cattle. Population structure was examined using the linkage model in the program STRUCTURE and Fst estimates were used to construct a neighbor-joining tree to represent the phylogenetic relationship among these breeds. Conclusion: The whole-genome SNP panel identified several levels of population substructure in the set of examined cattle breeds. The greatest level of genetic differentiation was detected between the Bos taurus and Bos indicus breeds. When the Bos indicus breeds were excluded from the analysis, genetic differences among beef versus dairy and European versus Asian breeds were detected among the Bos taurus breeds. Exploration of the number of SNP loci required to differentiate between breeds showed that for 100 SNP loci, individuals could only be correctly clustered into breeds 50% of the time, thus a large number of SNP markers are required to replace the 30 microsatellite markers that are currently commonly used in genetic diversity studies.

Adipose Tissue-Derived Stem Cells from Humans and Mice Differ in Proliferative Capacity and Genome Stability in Long-Term Cultures

Adipose tissue-derived stem cells (ASCs) are among the more attractive adult stem cell options for potential therapeutic applications. Here, we studied and compared the basic biological characteristics of ASCs isolated from humans (hASCs) and mice (mASCs) and maintained in identical culture conditions, which must be examined prior to considering further potential clinical applications. hASCs and mASCs were compared for immunophenotype, differentiation potential, cell growth characteristics, senescence, nuclear morphology, and DNA content. Although both strains of ASCs displayed a similar immunophenotype, the percentage of CD73(+) cells was markedly lower and CD31(+) was higher in mASC than in hASC cultures. The mean population doubling time was 98.08 +/- 6.15 h for hASCs and 52.58 +/- 3.74 h for mASCs. The frequency of nuclear aberrations was noticeably lower in hASCs than in mASCs regardless of the passage number. Moreover, as the cells went through several in vitro passages, mASCs showed changes in DNA content and cell cycle kinetics (frequency of hypodiploid, G0/G1, G2/M, and hyperdiploid cells), whereas all of these parameters remained constant in hASCs. Collectively, these results suggest that mASCs display higher proliferative capacity and are more unstable than hASCs in long-term cultures. These results underscore the need to consider specificities among model systems that may influence outcomes when designing potential human applications.

Integration and visualization of systems biology data in context of the genome

Background: High-density tiling arrays and new sequencing technologies are generating rapidly increasing volumes of transcriptome and protein-DNA interaction data. Visualization and exploration of this data is critical to understanding the regulatory logic encoded in the genome by which the cell dynamically affects its physiology and interacts with its environment. Results: The Gaggle Genome Browser is a cross-platform desktop program for interactively visualizing high-throughput data in the context of the genome. Important features include dynamic panning and zooming, keyword search and open interoperability through the Gaggle framework. Users may bookmark locations on the genome with descriptive annotations and share these bookmarks with other users. The program handles large sets of user-generated data using an in-process database and leverages the facilities of SQL and the R environment for importing and manipulating data. A key aspect of the Gaggle Genome Browser is interoperability. By connecting to the Gaggle framework, the genome browser joins a suite of interconnected bioinformatics tools for analysis and visualization with connectivity to major public repositories of sequences, interactions and pathways. To this flexible environment for exploring and combining data, the Gaggle Genome Browser adds the ability to visualize diverse types of data in relation to its coordinates on the genome. Conclusions: Genomic coordinates function as a common key by which disparate biological data types can be related to one another. In the Gaggle Genome Browser, heterogeneous data are joined by their location on the genome to create information-rich visualizations yielding insight into genome organization, transcription and its regulation and, ultimately, a better understanding of the mechanisms that enable the cell to dynamically respond to its environment.

Differential gene expression profiling in mucus glands of honey bee (Apis mellifera) drones during sexual maturation

The mating sign that each drone leaves when mating with a queen essentially consists of mucus gland proteins. We employed a Representational Difference Analysis (RDA) methodology to identify genes that are differentially expressed in mucus glands during sexual maturation of drones. The RDA library for mucus glands of newly emerged drones was more complex than that of 8 day-old drones, with matches to 20 predicted genes. Another 26 reads matched to the Apis genome but not to any predicted gene. Since these ESTs were located within ORFs they may represent novel honey bee genes, possibly fast evolving mucus gland proteins. In the RDA library for mucus glands of 8 day-old drones, most reads corresponded to a capsid protein of deformed wing virus, indicating high viral loads in these glands. The expression of two genes encoding venom allergens, acid phosphatase-1 and hyaluronidase, in drone mucus glands argues for their homology with the female venom glands, both associated with the reproductive system.

Genome Survey Sequence Analysis and Identification of Homologs of Major Surface Protease (gp63) Genes in Trypanosoma rangeli

In this study, 222 genome survey sequences were generated for Trypanosoma rangeli strain P07 isolated from an opossum (Didelphis albiventris) in Minas Gerais State, Brazil. T. rangeli sequences were compared by BLASTX (Basic Local Alignment Search Tool X) analysis with the assembled contigs of Leishmania braziliensis, Leishmania infantum, Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi. Results revealed that 82% (182/222) of the sequences were associated with predicted proteins described, whereas 18% (40/222) of the sequences did not show significant identity with sequences deposited in databases, suggesting that they may represent T. rangeli-specific sequences. Among the 182 predicted sequences, 179 (80.6%) had the highest similarity with T. cruzi, 2 (0.9%) with T. brucei, and 1 (0.5%) with L. braziliensis. Computer analysis permitted the identification of members of various gene families described for trypanosomatids in the genome of T. rangeli, such as trans-sialidases, mucin-associated surface proteins, and major surface proteases (MSP or gp63). This is the first report identifying sequences of the MSP family in T. rangeli. Multiple sequence alignments showed that the predicted MSP of T. rangeli presented the typical characteristics of metalloproteases, such as the presence of the HEXXH motif, which corresponds to a region previously associated with the catalytic site of the enzyme, and various cysteine and proline residues, which are conserved among MSPs of different trypanosomatid species. Reverse transcriptase-polymerase chain reaction analysis revealed the presence of MSP transcripts in epimastigote forms of T. rangeli.

Viral mutation and its influence in the time evolution of the immunizations

In this paper, we study the behavior of immune memory against antigenic mutation. Using a dynamic model proposed by one of the authors in a previous study (A. de Castro [Phys. J. Appl. Phys. 33, 147 (2006) and Simul. Mod. Pract. Theory. 15, 831 (2007)]), we have performed simulations of several inoculations, where in each virtual sample the viral population undergoes mutations. Our results suggest that the sustainability of the immunizations is dependent on viral variability and that the memory lifetimes are not random, what contradicts what was suggested by Tarlinton et al. [Curr. Opin. Immunol. 20, 162 (2008)]. We show that what may cause an apparent random behavior of the immune memory is the antigenic variability.

Detection of Human Bocavirus mRNA in Respiratory Secretions Correlates with High Viral Load and Concurrent Diarrhea

Human bocavirus (HBoV) is a parvovirus recently identified in association with acute respiratory infections (ARI). Despite its worldwide occurrence, little is known on the pathogenesis of HBoV infections. In addition, few systematic studies of HBoV in ARI have been conducted in Latin America. Therefore, in order to test whether active viral replication of human bocavirus is associated with respiratory diseases and to understand the clinical impact of this virus in patients with these diseases, we performed a 3-year retrospective hospital-based study of HBoV in outpatients and inpatients with symptoms of Acute Respiratory Infections (ARI) in Brazil. Nasopharyngeal aspirates (NPAs) from 1015 patients with respiratory symptoms were tested for HBoV DNA by PCR. All samples positive for HBoV were tested by PCR for all other respiratory viruses, had HBoV viral loads determined by quantitative real time PCR and, when possible, were tested by RT-PCR for HBoV VP1 mRNA, as evidence of active viral replication. HBoV was detected in 4.8% of patients, with annual rates of 10.0%, 3.0% and 3.0% in 2005, 2006 and 2007, respectively. The range of respiratory symptoms was similar between HBoV-positive and HBoV-negative ARI patients. However, a higher rate of diarrhea was observed in HBoV-positive patients. High HBoV viral loads (> 10(8) copies/mL) and diarrhea were significantly more frequent in patients with exclusive infection by HBoV and in patients with detection of HBoV VP1 mRNA than in patients with viral co-infection, detected in 72.9% of patients with HBoV. In summary, our data demonstrated that active HBoV replication was detected in a small percentage of patients with ARI and was correlated with concurrent diarrhea and lack of other viral co-infections.

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.

Genetically Diverse Coronaviruses in Captive Bird Populations in a Brazilian Zoological Park

This study aimed to investigate the occurrence of coronaviruses (CoVs) in captive birds placed inside a zoological park in Brazil. The role of captive birds in the epidemiology of CoVs in the tropics is poorly understood. A total of 25 (n = 25) different species were tested for viral RNA using individual fecal samples collected from healthy birds. Reverse transcription-polymerase chain reaction targeting the 30 untranslated region was used to detect CoV RNA, and positive samples were submitted for sequence analysis. The phylogenetic search revealed nine mutations in the black shouldered peafowl (Pavus cristatus) CoV sequence, which clustered separately from samples previously described in England. This is the first report on the detection of the CoV genome in captive birds in Brazil.

A Genome-Wide Association Study of Upper Aerodigestive Tract Cancers Conducted within the INHANCE Consortium

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p <= 5 x 10(-7)). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1 x 10(-8)) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p = 2 x 10(-8)) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5 x 10(-8); rs1229984-ADH1B, p = 7 x 10(-9); and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.