The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Genome wide scan for quantitative trait loci affecting tick resistance in cattle (Bos taurus x Bos indicus)

Background : In tropical countries, losses caused by bovine tick Rhipicephalus (Boophilus) microplus infestation have a tremendous economic impact on cattle production systems. Genetic variation between Bos taurus and Bos indicus to tick resistance and molecular biology tools might allow for the identification of molecular markers linked to resistance traits that could be used as an auxiliary tool in selection programs. The objective of this work was to identify QTL associated with tick resistance/susceptibility in a bovine F2 population derived from the Gyr (Bos indicus) x Holstein (Bos taurus) cross. Results: Through a whole genome scan with microsatellite markers, we were able to map six genomic regions associated with bovine tick resistance. For most QTL, we have found that depending on the tick evaluation season (dry and rainy) different sets of genes could be involved in the resistance mechanism. We identified dry season specific QTL on BTA 2 and 10, rainy season specific QTL on BTA 5, 11 and 27. We also found a highly significant genome wide QTL for both dry and rainy seasons in the central region of BTA 23. Conclusions: The experimental F2 population derived from Gyr x Holstein cross successfully allowed the identification of six highly significant QTL associated with tick resistance in cattle. QTL located on BTA 23 might be related with the bovine histocompatibility complex. Further investigation of these QTL will help to isolate candidate genes involved with tick resistance in cattle.

Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus

We present the genome sequences of a new clinical isolate of the important human pathogen, Aspergillus fumigatus, A1163, and two closely related but rarely pathogenic species, Neosartorya fischeri NRRL181 and Aspergillus clavatus NRRL1. Comparative genomic analysis of A1163 with the recently sequenced A. fumigatus isolate Af293 has identified core, variable and up to 2% unique genes in each genome. While the core genes are 99.8% identical at the nucleotide level, identity for variable genes can be as low 40%. The most divergent loci appear to contain heterokaryon incompatibility ( het) genes associated with fungal programmed cell death such as developmental regulator rosA. Cross-species comparison has revealed that 8.5%, 13.5% and 12.6%, respectively, of A. fumigatus, N. fischeri and A. clavatus genes are species-specific. These genes are significantly smaller in size than core genes, contain fewer exons and exhibit a subtelomeric bias. Most of them cluster together in 13 chromosomal islands, which are enriched for pseudogenes, transposons and other repetitive elements. At least 20% of A. fumigatus-specific genes appear to be functional and involved in carbohydrate and chitin catabolism, transport, detoxification, secondary metabolism and other functions that may facilitate the adaptation to heterogeneous environments such as soil or a mammalian host. Contrary to what was suggested previously, their origin cannot be attributed to horizontal gene transfer ( HGT), but instead is likely to involve duplication, diversification and differential gene loss (DDL). The role of duplication in the origin of lineage-specific genes is further underlined by the discovery of genomic islands that seem to function as designated ""gene dumps'' and, perhaps, simultaneously, as "" gene factories''.

Phylogenetic relationships within the speciose family Characidae (Teleostei: Ostariophysi: Characiformes) based on multilocus analysis and extensive ingroup sampling

Background: With nearly 1,100 species, the fish family Characidae represents more than half of the species of Characiformes, and is a key component of Neotropical freshwater ecosystems. The composition, phylogeny, and classification of Characidae is currently uncertain, despite significant efforts based on analysis of morphological and molecular data. No consensus about the monophyly of this group or its position within the order Characiformes has been reached, challenged by the fact that many key studies to date have non-overlapping taxonomic representation and focus only on subsets of this diversity. Results: In the present study we propose a new definition of the family Characidae and a hypothesis of relationships for the Characiformes based on phylogenetic analysis of DNA sequences of two mitochondrial and three nuclear genes (4,680 base pairs). The sequences were obtained from 211 samples representing 166 genera distributed among all 18 recognized families in the order Characiformes, all 14 recognized subfamilies in the Characidae, plus 56 of the genera so far considered incertae sedis in the Characidae. The phylogeny obtained is robust, with most lineages significantly supported by posterior probabilities in Bayesian analysis, and high bootstrap values from maximum likelihood and parsimony analyses. Conclusion: A monophyletic assemblage strongly supported in all our phylogenetic analysis is herein defined as the Characidae and includes the characiform species lacking a supraorbital bone and with a derived position of the emergence of the hyoid artery from the anterior ceratohyal. To recognize this and several other monophyletic groups within characiforms we propose changes in the limits of several families to facilitate future studies in the Characiformes and particularly the Characidae. This work presents a new phylogenetic framework for a speciose and morphologically diverse group of freshwater fishes of significant ecological and evolutionary importance across the Neotropics and portions of Africa.

Sequential hygienic behavior in Carniolan honey bees (Apis mellifera carnica)

We examined the sequence, order or steps of hygienic behavior (HB) from pin-killed pupae until the removal of them by the bees. We conducted our study with four colonies of Apis mellifera carnica in Germany and made four repetitions. The pin-killing method was used for evaluation of the HB of bees. The data were collected every 2 h after perforation, totaling 13 observations. Additionally, for one hygienic colony and another non-hygienic colony, individual analyses of each dead pupa were made at every observation, including all details, steps or sequences of HB. The bees recognize the cells containing dead pupae within 2 h after perforation, initially making a hole in the capping, which is the beginning of HB. Uncapping of the dead brood cell reached maximum values from 4 to 6 h after perforation; after 24 h, practically all cells were already uncapped. Another variable, called brood partially removed, was analyzed 4 h after perforation, after the cells had been perforated, which involved uncapping, followed by partial or total removal of the brood. Maximum values of brood partially removed were found 10 h after perforation, though such cells could be found up to 48 h after perforation. The most frequent sequence of events in both colonies was: capped cell -> punctured cell. brood partially removed -> empty cell. A new model of three pairs of recessive genes (uncapping u1, u2 and remover r) was proposed in order to explain the genetic control of the HB in Apis mellifera. We recommend evaluating HB 24 h after perforation and using a correction factor to compensate for control removal levels. We found a series of details of HB, which allow a study of how various factors may affect the sequence of the activities involved in HB and investigation of the genetics that controls this process.

Ecdysteroid-Dependent Expression of the Tweedle and Peroxidase Genes during Adult Cuticle Formation in the Honey Bee, Apis mellifera

Cuticle renewal is a complex biological process that depends on the cross talk between hormone levels and gene expression. This study characterized the expression of two genes encoding cuticle proteins sharing the four conserved amino acid blocks of the Tweedle family, AmelTwdl1 and AmelTwdl2, and a gene encoding a cuticle peroxidase containing the Animal haem peroxidase domain, Ampxd, in the honey bee. Gene sequencing and annotation validated the formerly predicted tweedle genes, and revealed a novel gene, Ampxd, in the honey bee genome. Expression of these genes was studied in the context of the ecdysteroid-coordinated pupal-to-adult molt, and in different tissues. Higher transcript levels were detected in the integument after the ecdysteroid peak that induces apolysis, coinciding with the synthesis and deposition of the adult exoskeleton and its early differentiation. The effect of this hormone was confirmed in vivo by tying a ligature between the thorax and abdomen of early pupae to prevent the abdominal integument from coming in contact with ecdysteroids released from the prothoracic gland. This procedure impaired the natural increase in transcript levels in the abdominal integument. Both tweedle genes were expressed at higher levels in the empty gut than in the thoracic integument and trachea of pharate adults. In contrast, Ampxd transcripts were found in higher levels in the thoracic integument and trachea than in the gut. Together, the data strongly suggest that these three genes play roles in ecdysteroid-dependent exoskeleton construction and differentiation and also point to a possible role for the two tweedle genes in the formation of the cuticle (peritrophic membrane) that internally lines the gut.

Sugarcane genes associated with sucrose content

Background -: Sucrose content is a highly desirable trait in sugarcane as the worldwide demand for cost-effective biofuels surges. Sugarcane cultivars differ in their capacity to accumulate sucrose and breeding programs routinely perform crosses to identify genotypes able to produce more sucrose. Sucrose content in the mature internodes reach around 20% of the culms dry weight. Genotypes in the populations reflect their genetic program and may display contrasting growth, development, and physiology, all of which affect carbohydrate metabolism. Few studies have profiled gene expression related to sugarcane's sugar content. The identification of signal transduction components and transcription factors that might regulate sugar accumulation is highly desirable if we are to improve this characteristic of sugarcane plants. Results -: We have evaluated thirty genotypes that have different Brix (sugar) levels and identified genes differentially expressed in internodes using cDNA microarrays. These genes were compared to existing gene expression data for sugarcane plants subjected to diverse stress and hormone treatments. The comparisons revealed a strong overlap between the drought and sucrose-content datasets and a limited overlap with ABA signaling. Genes associated with sucrose content were extensively validated by qRT-PCR, which highlighted several protein kinases and transcription factors that are likely to be regulators of sucrose accumulation. The data also indicate that aquaporins, as well as lignin biosynthesis and cell wall metabolism genes, are strongly related to sucrose accumulation. Moreover, sucrose-associated genes were shown to be directly responsive to short term sucrose stimuli, confirming their role in sugar-related pathways. Conclusion -: Gene expression analysis of sugarcane populations contrasting for sucrose content indicated a possible overlap with drought and cell wall metabolism processes and suggested signaling and transcriptional regulators to be used as molecular markers in breeding programs. Transgenic research is necessary to further clarify the role of the genes and define targets useful for sugarcane improvement programs based on transgenic plants.

A novel mutation in the glycogen synthase 2 gene in a child with glycogen storage disease type 0

Background: Glycogen storage disease type 0 is an autosomal recessive disease presenting in infancy or early childhood and characterized by ketotic hypoglycemia after prolonged fasting and postprandial hyperglycemia and hyperlactatemia. Sixteen different mutations have been identified to date in the gene which encodes hepatic glycogen synthase, resulting in reduction of glycogen storage in the liver. Case Presentation: Biochemical evaluation as well as direct sequencing of exons and exon-intron boundary regions of the GYS2 gene were performed in a patient presenting fasting hypoglycemia and postprandial hyperglycemia and her parents. The patient was found to be compound heterozygous for one previously reported nonsense mutation (c. 736 C>T; R243X) and a novel frameshift mutation (966_967delGA/insC) which introduces a stop codon 21 aminoacids downstream from the site of the mutation that presumably leads to loss of 51% of the COOH-terminal part of the protein. The glycemia and lactatemia of the parents after an oral glucose tolerance test were evaluated to investigate a possible impact of the carrier status on the metabolic profile. The mother, who presented a positive family history of type 2 diabetes, was classified as glucose intolerant and the father, who did not exhibit metabolic changes after the glucose overload, had an antecedent history of hypoglycemia after moderate alcohol ingestion. Conclusion: The current results expand the spectrum of known mutations in GYS2 and suggest that haploinsufficiency could explain metabolic abnormalities in heterozygous carriers in presence of predisposing conditions.

Characterization of Dengue Virus Type 2: New Insights on the 2010 Brazilian Epidemic

Dengue viruses (DENV) serotypes 1, 2, and 3 have been causing yearly outbreaks in Brazil. In this study, we report the reintroduction of DENV2 in the coast of Sao Paulo State. Partial envelope viral genes were sequenced from eighteen patients with dengue fever during the 2010 epidemic. Phylogenetic analysis showed this strain belongs to the American/Asian genotype and was closely related to the virus that circulated in Rio de Janeiro in 2007 and 2008. The phylogeny also showed no clustering by clinical presentation, suggesting that the disease severity could not be explained by distinct variants or genotypes. The time of the most recent common ancestor of American/Asian genotype and the Sao Paulo and Rio de Janeiro (SP/RJ) monophyletic cluster was estimated to be around 40 and 10 years, respectively. Since this virus was first identified in Brazil in 2007, we suggest that it was already circulating in the country before causing the first documented outbreak. This is the first description of the 2010 outbreak in the State of Sao Paulo, Brazil, and should contribute to efforts to control and monitor the spread of DENVs in endemic areas.

TCF7L2 variant genotypes and type 2 diabetes risk in Brazil: significant association, but not a significant tool for risk stratification in the general population

Background: Genetic polymorphisms of the TCF7L2 gene are strongly associated with large increments in type 2 diabetes risk in different populations worldwide. In this study, we aimed to confirm the effect of the TCF7L2 polymorphism rs7903146 on diabetes risk in a Brazilian population and to assess the use of this genetic marker in improving diabetes risk prediction in the general population. Methods: We genotyped the single nucleotide polymorphisms (SNP) rs7903146 of the TCF7L2 gene in 560 patients with known coronary disease enrolled in the MASS II (Medicine, Angioplasty, or Surgery Study) Trial and in 1,449 residents of Vitoria, in Southeast Brazil. The associations of this gene variant to diabetes risk and metabolic characteristics in these two different populations were analyzed. To access the potential benefit of using this marker for diabetes risk prediction in the general population we analyzed the impact of this genetic variant on a validated diabetes risk prediction tool based on clinical characteristics developed for the Brazilian general population. Results: SNP rs7903146 of the TCF7L2 gene was significantly associated with type 2 diabetes in the MASS-II population (OR = 1.57 per T allele, p = 0.0032), confirming, in the Brazilian population, previous reports of the literature. Addition of this polymorphism to an established clinical risk prediction score did not increased model accuracy (both area under ROC curve equal to 0.776). Conclusion: TCF7L2 rs7903146 T allele is associated with a 1.57 increased risk for type 2 diabetes in a Brazilian cohort of patients with known coronary heart disease. However, the inclusion of this polymorphism in a risk prediction tool developed for the general population resulted in no improvement of performance. This is the first study, to our knowledge, that has confirmed this recent association in a South American population and adds to the great consistency of this finding in studies around the world. Finally, confirming the biological association of a genetic marker does not guarantee improvement on already established screening tools based solely on demographic variables.

An insight into the sialotranscriptome of the brown dog tick, Rhipicephalus sanguineus

Background: Rhipicephalus sanguineus, known as the brown dog tick, is a common ectoparasite of domestic dogs and can be found worldwide. R. sanguineus is recognized as the primary vector of the etiological agent of canine monocytic ehrlichiosis and canine babesiosis. Here we present the first description of a R. sanguineus salivary gland transcriptome by the production and analysis of 2,034 expressed sequence tags (EST) from two cDNA libraries, one consctructed using mRNA from dissected salivary glands from female ticks fed for 3-5 days (early to mid library, RsSGL1) and the another from ticks fed for 5 days (mid library, RsSGL2), identifying 1,024 clusters of related sequences. Results: Based on sequence similarities to nine different databases, we identified transcripts of genes that were further categorized according to function. The category of putative housekeeping genes contained similar to 56% of the sequences and had on average 2.49 ESTs per cluster, the secreted protein category contained 26.6% of the ESTs and had 2.47 EST's/clusters, while 15.3% of the ESTs, mostly singletons, were not classifiable, and were annotated as ""unknown function"". The secreted category included genes that coded for lipocalins, proteases inhibitors, disintegrins, metalloproteases, immunomodulatory and antiinflammatory proteins, as Evasins and Da-p36, as well as basic-tail and 18.3 kDa proteins, cement proteins, mucins, defensins and antimicrobial peptides. Comparison of the abundance of ESTs from similar contigs of the two salivary gland cDNA libraries allowed the identification of differentially expressed genes, such as genes coding for Evasins and a thrombin inhibitor, which were over expressed in the RsSGL1 (early to mid library) versus RsSGL2 (mid library), indicating their role in inhibition of inflammation at the tick feeding site from the very beginning of the blood meal. Conversely, sequences related to cement (64P), which function has been correlated with tick attachment, was largely expressed in the mid library. Conclusions: Our survey provided an insight into the R. sanguineus sialotranscriptome, which can assist the discovery of new targets for anti-tick vaccines, as well as help to identify pharmacologically active proteins.

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.

Suppression subtractive hybridization analysis reveals expression of conserved and novel genes in male accessory glands of the ant Leptothorax gredleri

Background: During mating, insect males eject accessory gland proteins (Acps) into the female genital tract. These substances are known to affect female post-mating behavior and physiology. In addition, they may harm the female, e. g., in reducing its lifespan. This is interpreted as a consequence of sexual antagonistic co-evolution. Whereas sexual conflict abounds in non-social species, the peculiar life history of social insects (ants, bees, wasps) with lifelong pair-bonding and no re-mating aligns the reproductive interests of the sexes. Harming the female during mating would negatively affect male fitness and sexual antagonism is therefore not expected. Indeed, mating appears to increase female longevity in at least one ant species. Acps are presumed to play a role in this phenomenon, but the underlying mechanisms are unknown. In this study, we investigated genes, which are preferentially expressed in male accessory glands of the ant Leptothorax gredleri, to determine which proteins might be transferred in the seminal fluid. Results: By a suppression subtractive hybridization protocol we obtained 20 unique sequences (USs). Twelve had mutual best matches with genes predicted for Apis mellifera and Nasonia vitripennis. Functional information (Gene Ontology) was available only for seven of these, including intracellular signaling, energy-dependent transport and metabolic enzyme activities. The remaining eight USs did not match sequences from other species. Six genes were further analyzed by quantitative RT-PCR in three life cycle stages of male ants. A gene with carboxy-lyase activity and one of unpredicted function were significantly overexpressed in accessory glands of sexually mature males. Conclusions: Our study is the first one to investigate differential gene expression in ants in a context related to mating. Our findings indicate that male accessory glands of L. gredleri express a series of genes that are unique to this species, possibly representing novel genes, in addition to conserved ones for which functions can be predicted. Identifying differentially expressed genes might help to better understand molecular mechanisms involved in reproductive processes in eusocial Hymenoptera. While the novel genes could account for rapidly evolving ones driven by intra-sexual conflict between males, conserved genes imply that rather beneficial traits might get fixed by a process described as inter-sexual cooperation between males and females.

Interethnic Differences in the Distribution of Matrix Metalloproteinases Genetic Polymorphisms Are Consistent with Interethnic Differences in Disease Prevalence

Interethnic differences exist in disease prevalence, especially with regard to cancer and cardiovascular diseases, which involve altered expression or activity of matrix metalloproteinases (MMPs). The hypothesis being tested in this study is that interethnic differences exist between blacks and whites with regard to the distribution of genetic variants of MMP polymorphisms and haplotypes. We examined the distribution of polymorphisms of MMP-2 and MMP-9 genes in 177 black and 140 white subjects. We studied the following polymorphisms: the C(-1306)T in the promoter of the MMP-2 gene, the C(-1562)T and a microsatellite -90(CA)(14-24) in the promoter, and the Q279R in exon 6 of the MMP-9 gene. We have also compared our results with those from Hapmap or Seattle SNPs Projects and estimated the haplotype frequency in these two ethnic groups. The ""C'' allele for the C(-1306)T polymorphism was more common in blacks (91.5%) than in whites (80.4%; p<0.0001). The ""T'' allele for the C(-1562)T polymorphism was more common in blacks (15.0%) than in whites (8.9%; p=0.0279), as well as the alleles with >21 repeats for the -90(CA)(14-24) were more common in blacks than in whites (61.9% in blacks and 49.3% in whites; p=0.0017). We found no interethnic differences for the Q279R polymorphism. Moreover, two haplotypes that combine ""detrimental'' alleles were found at higher frequencies in blacks than in whites (31% vs. 16.4%, respectively; p<0.05). The interethnic differences being reported here replicate those previously found with smaller number of subjects in the Hapmap or Seattle SNPs data and may help explain the higher prevalence of cancer and cardiovascular diseases in blacks compared with whites. Our findings suggest a proportional significance of these polymorphisms in each ethnic group.

Endothelial Nitric Oxide Synthase Haplotypes Associated with Hypertension Do Not Predispose to Cardiac Hypertrophy

Left ventricular hypertrophy (LVH) is a complication that may result from chronic hypertension. While nitric oxide (NO) deficiency has been associated with LVH, inconsistent results have been reported with regards to the association of endothelial NO synthase (eNOS) polymorphisms and LVH in hypertensive patients. This study aims to assess whether eNOS haplotypes are associated with LVH in hypertensive patients. This study included 101 healthy controls and 173 hypertensive patients submitted to echocardiography examination. Genotypes for three eNOS polymorphisms were determined: a single-nucleotide polymorphism in the promoter region (T-786C) and in exon 7 (Glu298Asp), and variable number of tandem repeats in intron 4. We found no significant association between eNOS genotypes and hypertension or with LVH (all p>0.05). However, while we found two eNOS haplotypes associated with variable risk of hypertension (all p<0.05), we found no significant associations between eNOS haplotypes and LVH (all p>0.05), even after adjustment in multiple linear regression analysis. These findings suggest that eNOS haplotypes that have been associated with variable susceptibility to hypertension were not associated with LVH in hypertensive patients. Further studies are necessary to examine whether other genes downstream may interact with eNOS polymorphisms and predispose to LVH in hypertensive patients.