High-throughput SNP genotyping in the highly heterozygous genome of Eucalyptus: assay success, polymorphism and transferability across species

Background: High-throughput SNP genotyping has become an essential requirement for molecular breeding and population genomics studies in plant species. Large scale SNP developments have been reported for several mainstream crops. A growing interest now exists to expand the speed and resolution of genetic analysis to outbred species with highly heterozygous genomes. When nucleotide diversity is high, a refined diagnosis of the target SNP sequence context is needed to convert queried SNPs into high-quality genotypes using the Golden Gate Genotyping Technology (GGGT). This issue becomes exacerbated when attempting to transfer SNPs across species, a scarcely explored topic in plants, and likely to become significant for population genomics and inter specific breeding applications in less domesticated and less funded plant genera. Results: We have successfully developed the first set of 768 SNPs assayed by the GGGT for the highly heterozygous genome of Eucalyptus from a mixed Sanger/454 database with 1,164,695 ESTs and the preliminary 4.5X draft genome sequence for E. grandis. A systematic assessment of in silico SNP filtering requirements showed that stringent constraints on the SNP surrounding sequences have a significant impact on SNP genotyping performance and polymorphism. SNP assay success was high for the 288 SNPs selected with more rigorous in silico constraints; 93% of them provided high quality genotype calls and 71% of them were polymorphic in a diverse panel of 96 individuals of five different species. SNP reliability was high across nine Eucalyptus species belonging to three sections within subgenus Symphomyrtus and still satisfactory across species of two additional subgenera, although polymorphism declined as phylogenetic distance increased. Conclusions: This study indicates that the GGGT performs well both within and across species of Eucalyptus notwithstanding its nucleotide diversity >= 2%. The development of a much larger array of informative SNPs across multiple Eucalyptus species is feasible, although strongly dependent on having a representative and sufficiently deep collection of sequences from many individuals of each target species. A higher density SNP platform will be instrumental to undertake genome-wide phylogenetic and population genomics studies and to implement molecular breeding by Genomic Selection in Eucalyptus.

Long Terminal Repeat Sequence Analysis of HTLV-2 Molecular Variants Identified in Southern Brazil

In Brazil, human T-lymphotropic virus type 2 (HTLV-2) is endemic in Amerindians and epidemic in intravenous drug users (IDUs). The long terminal repeat (LTR) is the most divergent genomic region of HTLV-2, therefore useful to characterize subtypes. Nucleotide sequence and restriction fragment length polymorphism (RFLP) analysis of LTR genomic segments of fourteen HTLV-2 strains isolated from HIV-infected patients of Londrina, Southern Brazil, were carried out. Molecular analysis disclosed that all HTLV-2 strains belonged to 2a subtype, and RFLP detected the presence of the a4, a5, and a6 subgroups according to Switzer's nomenclature. RFLP correlated with nucleotide sequence, and phylogenetic analysis clustered HTLV-2 sequences of IDUs into subgroups a5 and a6. HTLV-2 sequences from individuals of sexual risk factor clustered into the a4 subgroup. These results extend the knowledge of the genetic diversity of HTLV-2 circulating in Brazil and provide insights into HTLV-2 transmission and virus movement in this geographic area.

Sub-telomere directed gene expression during initiation of invasive aspergillosis

Aspergillus fumigatus is a common mould whose spores are a component of the normal airborne flora. Immune dysfunction permits developmental growth of inhaled spores in the human lung causing aspergillosis, a significant threat to human health in the form of allergic, and life-threatening invasive infections. The success of A. fumigatus as a pathogen is unique among close phylogenetic relatives and is poorly characterised at the molecular level. Recent genome sequencing of several Aspergillus species provides an exceptional opportunity to analyse fungal virulence attributes within a genomic and evolutionary context. To identify genes preferentially expressed during adaptation to the mammalian host niche, we generated multiple gene expression profiles from minute samplings of A. fumigatus germlings during initiation of murine infection. They reveal a highly co-ordinated A. fumigatus gene expression programme, governing metabolic and physiological adaptation, which allows the organism to prosper within the mammalian niche. As functions of phylogenetic conservation and genetic locus, 28% and 30%, respectively, of the A. fumigatus subtelomeric and lineage-specific gene repertoires are induced relative to laboratory culture, and physically clustered genes including loci directing pseurotin, gliotoxin and siderophore biosyntheses are a prominent feature. Locationally biased A. fumigatus gene expression is not prompted by in vitro iron limitation, acid, alkaline, anaerobic or oxidative stress. However, subtelomeric gene expression is favoured following ex vivo neutrophil exposure and in comparative analyses of richly and poorly nourished laboratory cultured germlings. We found remarkable concordance between the A. fumigatus host-adaptation transcriptome and those resulting from in vitro iron depletion, alkaline shift, nitrogen starvation and loss of the methyltransferase LaeA. This first transcriptional snapshot of a fungal genome during initiation of mammalian infection provides the global perspective required to direct much-needed diagnostic and therapeutic strategies and reveals genome organisation and subtelomeric diversity as potential driving forces in the evolution of pathogenicity in the genus Aspergillus.

ProbFAST: Probabilistic Functional Analysis System Tool

Background: The post-genomic era has brought new challenges regarding the understanding of the organization and function of the human genome. Many of these challenges are centered on the meaning of differential gene regulation under distinct biological conditions and can be performed by analyzing the Multiple Differential Expression (MDE) of genes associated with normal and abnormal biological processes. Currently MDE analyses are limited to usual methods of differential expression initially designed for paired analysis. Results: We proposed a web platform named ProbFAST for MDE analysis which uses Bayesian inference to identify key genes that are intuitively prioritized by means of probabilities. A simulated study revealed that our method gives a better performance when compared to other approaches and when applied to public expression data, we demonstrated its flexibility to obtain relevant genes biologically associated with normal and abnormal biological processes. Conclusions: ProbFAST is a free accessible web-based application that enables MDE analysis on a global scale. It offers an efficient methodological approach for MDE analysis of a set of genes that are turned on and off related to functional information during the evolution of a tumor or tissue differentiation. ProbFAST server can be accessed at http://gdm.fmrp.usp.br/probfast.

Gene Expression Profiles in Radiation Workers Occupationally Exposed to Ionizing Radiation

Ionizing radiation OR) imposes risks to human health and the environment. IR at low doses and low (lose rates has the potency to initiate carcinogenesis. Genotoxic environmental agents such as IR trigger a cascade of signal transduction pathways for cellular protection. In this study, using cDNA microarray technique, we monitored the gene expression profiles in lymphocytes derived from radiation-ex posed individuals (radiation workers). Physical dosimetry records on these patients indicated that the absorbed dose ranged from 0.696 to 39.088 mSv. Gene expression analysis revealed statistically significant transcriptional changes in a total of 78 genes (21 up-regulated and 57 clown-regulated) involved in several biological processes such as ubiquitin cycle (UHRF2 and PIAS1), DNA repair (LIG3, XPA, ERCC5, RAD52, DCLRE1C), cell cycle regulation/proliferation (RHOA, CABLES2, TGFB2, IL16), and stress response (GSTP1, PPP2R5A, DUSP22). Some of the genes that showed altered expression profiles in this study call be used as biomarkers for monitoring the chronic low level exposure in humans. Additionally, alterations in gene expression patterns observed in chronically exposed radiation workers reinforces the need for defining the effective radiation dose that causes immediate genetic damage as well as the long-term effects on genomic instability, including cancer.

Ring chromosome instability evaluation in six patients with autosomal rings

Ring chromosomes are often associated with abnormal phenotypes due to loss of genomic material and also because of ring instability at mitosis after sister chromatid exchange events. We investigated ring chromosome instability in six patients with ring chromosomes 4, 14, 15, and 18 by examining 48- and 72-h lymphocyte cultures at the first, second and subsequent cell divisions after bromodeoxyuridine incorporation. Although most cells from all patients showed only one monocentric ring chromosome, ring chromosome loss and secondary aberrations were observed both in 48-and 72-h lymphocyte cultures and in metaphase cells of the different cell generations. We found no clear-cut correlation between ring size and ring instability; we also did not find differences between apparently complete rings and rings with genetic material loss. The cytogenetic findings revealed secondary aberrations in all ring chromosome patients. We concluded that cells with ring chromosome instability can multiply and survive in vivo, and that they can influence the patient's phenotype.

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.

Methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism and high plasma homocysteine in chronic hepatitis C (CHC) infected patients from the Northeast of Brazil

Background/Aim: Hyperhomocysteinemia due to Methylenetetrahydrofolate Reductase (MTHFR) gene, in particular the C677T (Ala222Val) polymorphism were recently associated to steatosis and fibrosis. We analyzed the frequency of MTHFR gene in a cross-sectional study of patients affected by Chronic Hepatitis C (CHC) from Northeast of Brazil. Method: One hundred seven-four untreated patients with CHC were genotyped for the C677T MTHFR. Genomic DNA was extracted from peripheral blood cells and the C677T MTHFR polymorphism was identified by PCR-RFLP. The homocysteine (Hcy) levels were determined by chemiluminescence method. All patients were negative for markers of Wilson's disease, hemochromatosis and autoimmune diseases and have current and past daily alcohol intake less than 100 g/week. Results: Among subjects infected with CHC genotype non-1 the frequency of MTHFR genotypes TT was 9.8% versus 4.4% genotype 1 (p = 0.01). Nevertheless, association was found between the MTHFR genotype TT x CT/CC polymorphism and the degree of steatosis and fibrosis in both hepatitis C genotype (p < 0.05). A significant difference was found on plasma Hcy levels in patients with steatosis regardless of HCV genotype (p = 0.03). Conclusion: Our results indicate that plasma Hcy levels is highly prevalent in subjects with chronic hepatits C with steatosis regardless of HCV genotype and vitamin deficiency. The presence of genotype TT of MTHFR C677T polymorphism was more common in CHC genotype non-1 infected patient regardless of histopathological classification and genotype TT+CT frequencies were significant in the presence of fibrosis grade 1+2 and of steatosis in CHC infected patients from the northeast of Brazil regardless of HCV genotype. The genetic susceptibility of MTHFR C677T polymorphism should be confirmed in a large population.

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.

Genetic Background of Patients from a University Medical Center in Manhattan: Implications for Personalized Medicine

Background: The rapid progress currently being made in genomic science has created interest in potential clinical applications; however, formal translational research has been limited thus far. Studies of population genetics have demonstrated substantial variation in allele frequencies and haplotype structure at loci of medical relevance and the genetic background of patient cohorts may often be complex. Methods and Findings: To describe the heterogeneity in an unselected clinical sample we used the Affymetrix 6.0 gene array chip to genotype self-identified European Americans (N = 326), African Americans (N = 324) and Hispanics (N = 327) from the medical practice of Mount Sinai Medical Center in Manhattan, NY. Additional data from US minority groups and Brazil were used for external comparison. Substantial variation in ancestral origin was observed for both African Americans and Hispanics; data from the latter group overlapped with both Mexican Americans and Brazilians in the external data sets. A pooled analysis of the African Americans and Hispanics from NY demonstrated a broad continuum of ancestral origin making classification by race/ethnicity uninformative. Selected loci harboring variants associated with medical traits and drug response confirmed substantial within-and between-group heterogeneity. Conclusion: As a consequence of these complementary levels of heterogeneity group labels offered no guidance at the individual level. These findings demonstrate the complexity involved in clinical translation of the results from genome-wide association studies and suggest that in the genomic era conventional racial/ethnic labels are of little value.

Cytogenetic findings in an epithelioid sarcoma with angiomatoid features. A case report

Epithelioid sarcoma is a rare, aggressive soft tissue tumor of unknown histogenesis showing predominantly epithelioid cytomorphology. We conducted a conventional and molecular cytogenetic study of a 27-year-old male with epithelioid sarcoma with angiomatoid features. Cytogenetic analysis of epithelioid sarcoma metaphase spreads by GTG-banding revealed a diploid chromosome complement with structural and numerical aberrations. Comparative genomic hybridization analysis demonstrated the amplification of 3p24 similar to pter, 4p15.2-p16 and 18q23, while chromosome losses involved 3p13-p14, 3q24-q26.1, 9q21, and 11q21. Fluorescence in situ hybridization assessment showed normal hybridization patterns for the C-MYC and CCND1 loci; CCND1 RNA overexpression was detected by real-time polymerase chain reaction analysis. Genetic evaluation of this rare condition may be useful in determining if epithelioid sarcoma is associated with a distinct genetic background.

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.

A scientific note on the isolation and characterization of microsatellite loci of Frieseomelitta varia (Hymenoptera, Apidae, Meliponini)

FAPESP[BIOTA 2004/15801-0]

Endothelial Nitric Oxide Synthase Polymorphisms and Haplotypes in Amerindians

Interethnic disparities in the distribution of endothelial nitric oxide synthase (eNOS) polymorphisms may affect nitric oxide (NO)-mediated effects of and responses to drugs. While there are differences between black and white subjects there is no information regarding the distribution of eNOS gene alleles and haplotypes in Amerindians. We studied three clinically relevant eNOS polymorphisms (T(-786) C in the promoter, a variable number of tandem repeats in intron 4, and the Glu298Asp in exon 7) and eNOS haplotypes in 170 Amerindians from three tribes of the Brazilian Amazon. The results were compared with previous findings for black and white Brazilians. The Asp298, C(-786), and 4a alleles were much less common in Amerindians (5.0%, 3.2%, and 4.1%, respectively) than in blacks (15.1%, 19.5%, and 32.0%, respectively) or whites (32.8%, 41.9%, and 17.9%, respectively) (p<0.001). The haplotype including the most common alleles for each polymorphism was much more common in Amerindians (89%) than in blacks (45%) or whites (41%). Our findings are consistent with a lower genetic diversity in Amerindians compared with blacks and whites. These striking differences may be of major relevance for case-control association studies focusing on eNOS gene polymorphisms and may explain, at least in part, differences in the responses to cardiovascular drugs.

Concurrent sequence variation of TP53 and TP73 genes in anaplastic astrocytoma

Disruption or loss of tumor suppressor gene TP53 is implicated in the development or progression of almost all different types of human malignancies. Other members of the p53 family have been identified. One member, p73, not only shares a high degree of similarity with p53 in its primary sequence, but also has similar functions. Like p53, p73 can bind to DNA and activate transcription. Using PCR-SSCP and gene sequencing, we analyzed the TP53 and TP73 genes in a case of a grade III anaplastic astrocytoma that progressed to glioblastoma. We found a deletion of AAG at position 595-597 of TP53 (exon 6), resulting in the deletion of Glu 199 in the protein and a genomic polymorphism of TP73, identified as an A-to-G change, at position E8/+15 at intron 8 (IVS8-15A>G). The mutation found at exon 6 of the gene TP53 could be associated with the rapid tumoral progression found in this case, since the mutated p53 may inactivate the wild-type p53 and the p73 alpha protein, which was conserved here, leading to an increase in cellular instability.