A proof of the DBRF-MEGN method, an algorithm for deducing minimum equivalent gene networks

Background

We previously developed the DBRF-MEGN (difference-based regulation finding-minimum equivalent gene network) method, which deduces the most parsimonious signed directed graphs (SDGs) consistent with expression profiles of single-gene deletion mutants. However, until the present study, we have not presented the details of the method's algorithm or a proof of the algorithm.

Results

We describe in detail the algorithm of the DBRF-MEGN method and prove that the algorithm deduces all of the exact solutions of the most parsimonious SDGs consistent with expression profiles of gene deletion mutants.

Conclusions

The DBRF-MEGN method provides all of the exact solutions of the most parsimonious SDGs consistent with expression profiles of gene deletion mutants.     

Soluble RAGE but not endogenous secretory RAGE is associated with albuminuria in patients with type 2 diabetes

Aims

Type 2 diabetes is characterised by increased plasma concentrations of pro-inflammatory cytokines [such as tumour necrosis factor – alpha; TNF-α] and soluble forms of adhesion molecules involved in leukocyte – endothelial interactions. These molecules are synthesised as transmembrane proteins and the plasma soluble forms are generated by ectodomain cleavage from the cell surface by members of the ADAM [a disintegrin and metalloproteinase] proteinase family. We hypothesised that plasma low density lipoprotein [LDL] from subjects with Type 2 diabetes would influence in vitro monocytic ADAM and matrix metalloproteinase [MMP] gene expression differently compared to control LDL.

Methods

We examined relative mRNA expression by real time PCR in a monocytic cell line [THP-1] cultured for 4, 8 and 24 hrs with human plasma LDL derived from subjects with [n = 5] or without [n = 4] Type 2 diabetes. Gene expression for MMP-1 and 9, and ADAM – 8, 15, 17 and 28 was studied.

Results

Type 2 diabetes LDL significantly increased gene expression of MMP – 1 [p < 0.01] MMP – 9 [p < 0.001], and ADAM 17 [p < 0.05], – 28 [p < 0.01] and – 15 [p < 0.01] compared to control LDL. Type 2 diabetes LDL had disparate effects on inhibitors of MMP.

Conclusion

These data suggest that Type 2 diabetes LDL could lead to increased adhesion molecule and TNF alpha cell surface shedding, and vascular plaque instability, by promoting increased expression of ADAM and MMP genes.     

A mosaic genetic screen for novel mutations affecting Drosophila neuroblast divisions

Background

The asymmetric segregation of determinants during cell division is a fundamental mechanism for generating cell fate diversity during development. In Drosophila, neural precursors (neuroblasts) divide in a stem cell-like manner generating a larger apical neuroblast and a smaller basal ganglion mother cell. The cell fate determinant Prospero and its adapter protein Miranda are asymmetrically localized to the basal cortex of the dividing neuroblast and segregated into the GMC upon cytokinesis. Previous screens to identify components of the asymmetric division machinery have concentrated on embryonic phenotypes. However, such screens are reaching saturation and are limited in that the maternal contribution of many genes can mask the effects of zygotic loss of function, and other approaches will be necessary to identify further genes involved in neuroblast asymmetric division.

Results

We have performed a genetic screen in the third instar larval brain using the basal localization of Miranda as a marker for neuroblast asymmetry. In addition to the examination of pupal lethal mutations, we have employed the MARCM (Mosaic Analysis with a Repressible Cell Marker) system to generate postembryonic clones of mutations with an early lethal phase. We have screened a total of 2,300 mutagenized chromosomes and isolated alleles affecting cell fate, the localization of basal determinants or the orientation of the mitotic spindle. We have also identified a number of complementation groups exhibiting defects in cell cycle progression and cytokinesis, including both novel genes and new alleles of known components of these processes.

Conclusion

We have identified four mutations which affect the process of neuroblast asymmetric division. One of these, mapping to the imaginal discs arrested locus, suggests a novel role for the anaphase promoting complex/cyclosome (APC/C) in the targeting of determinants to the basal cortex. The identification and analysis of the remaining mutations will further advance our understanding of the process of asymmetric cell division. We have also isolated a number of mutations affecting cell division which will complement the functional genomics approaches to this process being employed by other laboratories. Taken together, these results demonstrate the value of mosaic screens in the identification of genes involved in neuroblast division.     

The role of Drosophila Merlin in spermatogenesis

Background

Drosophila Merlin, the homolog of the human Neurofibromatosis 2 (NF2) gene, is important for the regulation of cell proliferation and receptor endocytosis. Male flies carrying a Mer ³ allele, a missense mutation (Met¹⁷⁷→Ile) in the Merlin gene, are viable but sterile; however, the cause of sterility is unknown.

Results

Testis examination reveals that hemizygous Mer ³ mutant males have small seminal vesicles that contain only a few immotile sperm. By cytological and electron microscopy analyses of the Mer ³, Mer ⁴ (Gln¹⁷⁰→stop), and control testes at various stages of spermatogenesis, we show that Merlin mutations affect meiotic cytokinesis of spermatocytes, cyst polarization and nuclear shaping during spermatid elongation, and spermatid individualization. We also demonstrate that the lethality and sterility phenotype of the Mer ⁴ mutant is rescued by the introduction of a wild-type Merlin gene. Immunostaining demonstrates that the Merlin protein is redistributed to the area associated with the microtubules of the central spindle in telophase and its staining is less in the region of the contractile ring during meiotic cytokinesis. At the onion stage, Merlin is concentrated in the Nebenkern of spermatids, and this mitochondrial localization is maintained throughout sperm formation. Also, Merlin exhibits punctate staining in the acrosomal region of mature sperm.

Conclusion

Merlin mutations affect spermatogenesis at multiple stages. The Merlin protein is dynamically redistributed during meiosis of spermatocytes and is concentrated in the Nebenkern of spermatids. Our results demonstrated for the first time the mitochondrial localization of Merlin and suggest that Merlin may play a role in mitochondria formation and function during spermatogenesis.     

Genome skimming identifies polymorphism in tern populations and species

Background

Large-scale protein structure alignment, an indispensable tool to structural bioinformatics, poses a tremendous challenge on computational resources. To ensure structure alignment accuracy and efficiency, efforts have been made to parallelize traditional alignment algorithms in grid environments. However, these solutions are costly and of limited accessibility. Others trade alignment quality for speedup by using high-level characteristics of structure fragments for structure comparisons.

Findings

We present ppsAlign, a parallel protein structure Alignment framework designed and optimized to exploit the parallelism of Graphics Processing Units (GPUs). As a general-purpose GPU platform, ppsAlign could take many concurrent methods, such as TM-align and Fr-TM-align, into the parallelized algorithm design. We evaluated ppsAlign on an NVIDIA Tesla C2050 GPU card, and compared it with existing software solutions running on an AMD dual-core CPU. We observed a 36-fold speedup over TM-align, a 65-fold speedup over Fr-TM-align, and a 40-fold speedup over MAMMOTH.

Conclusions

ppsAlign is a high-performance protein structure alignment tool designed to tackle the computational complexity issues from protein structural data. The solution presented in this paper allows large-scale structure comparisons to be performed using massive parallel computing power of GPU.     

Genome-scale identification of Caenorhabditis elegans regulatory elements by tiling-array mapping of DNase I hypersensitive sites

Background

High throughput techniques have generated a huge set of biological data, which are deposited in various databases. Efficient exploitation of these databases is often hampered by a lack of appropriate tools, which allow easy and reliable identification of genes that miss functional characterization but are correlated with specific biological conditions (e.g. organotypic expression).

Results

We have developed a simple algorithm (DGSA = Database-dependent Gene Selection and Analysis) to identify genes with unknown functions involved in organ development concentrating on the heart. Using our approach, we identified a large number of yet uncharacterized genes, which are expressed during heart development. An initial functional characterization of genes by loss-of-function analysis employing morpholino injections into zebrafish embryos disclosed severe developmental defects indicating a decisive function of selected genes for developmental processes.

Conclusion

We conclude that DGSA is a versatile tool for database mining allowing efficient selection of uncharacterized genes for functional analysis.     

Fine-scale genetic breaks driven by historical range dynamics and ongoing density-barrier effects in the estuarine seaweed Fucus ceranoides L.

Background

In the laboratory, the Drosophila melanogaster heat shock protein Hsp90 can buffer the phenotypic effects of genetic variation. Laboratory experiments either manipulate Hsp90 activity pharmacologically, or they induce mutations with strong effects in the gene Hsp83, the single-copy fly gene encoding Hsp90. It is unknown whether observations from such laboratory experiments are relevant in the wild.

Results

We here study naturally occurring mutations in Hsp83, and their effects on fitness and phenotypic buffering in flies derived from wild populations. We examined more than 4500 flies from 42 Drosophila populations distributed world-wide for insertions or deletions of mobile DNA in or near the Hsp83 gene. The insertions we observed occur at low population frequencies, and reduce Hsp83 gene expression. In competition experiments, mutant flies performed much more poorly than wild-type flies. Mutant flies were also significantly less fecund and shorter-lived than wild-type flies, as well as less well buffered against cryptic deleterious variation, as we show through inbreeding experiments. Specifically, in Hsp83 mutant flies female fecundity dropped to much lower levels after inbreeding than in wild-type flies. At even slightly elevated temperatures, inbred mutant Hsp83 populations went extinct, whereas inbred wild-type populations persisted.

Conclusions

Our work shows that Hsp90, a regulator of the stress response and of signaling, helps buffer deleterious variation in fruit flies derived from wild population, and that its buffering role becomes even more important under heat stress.     

Coccidioides posadasii infection alters the expression of pulmonary surfactant proteins (SP)-A and SP-D

Background

Patients with asthma demonstrate circadian variations in the airway inflammation and lung function. Pinealectomy reduces the total inflammatory cell number in the asthmatic rat lung. We hypothesize that melatonin, a circadian rhythm regulator, may modulate the circadian inflammatory variations in asthma by stimulating the chemotaxins expression in the lung epithelial cell.

Methods

Lung epithelial cells (A549) were stimulated with melatonin in the presence or absence of TNF-α(100 ng/ml). RANTES (Regulated on Activation Normal T-cells Expressed and Secreted) and eotaxin expression were measured using ELISA and real-time RT-PCR, eosinophil chemotactic activity (ECA) released by A549 was measured by eosinophil chemotaxis assay.

Results

TNF-α increased the expression of RANTES (307.84 ± 33.56 versus 207.64 ± 31.27 pg/ml of control, p = 0.025) and eotaxin (108.97 ± 10.87 versus 54.00 ± 5.29 pg/ml of control, p = 0.041). Melatonin(10^-10 to 10^-6M) alone didn't change the expression of RNATES (204.97 ± 32.56 pg/ml) and eotaxin (55.28 ± 6.71 pg/ml). However, In the presence of TNF-α (100 ng/ml), melatonin promoted RANTES (410.88 ± 52.03, 483.60 ± 55.37, 559.92 ± 75.70, 688.42 ± 95.32, 766.39 ± 101.53 pg/ml, treated with 10^-10, 10^-9, 10^-8, 10^-7,10^-6M melatonin, respectively) and eotaxin (151.95 ± 13.88, 238.79 ± 16.81, 361.62 ± 36.91, 393.66 ± 44.89, 494.34 ± 100.95 pg/ml, treated with 10^-10, 10^-9, 10^-8, 10^-7, 10^-6M melatonin, respectively) expression in a dose dependent manner in A549 cells (compared with TNF-α alone, P < 0.05). The increased release of RANTES and eotaxin in A549 cells by above treatment were further confirmed by both real-time RT-PCR and the ECA assay.

Conclusion

Taken together, our results suggested that melatonin might synergize with pro-inflammatory cytokines to modulate the asthma airway inflammation through promoting the expression of chemotaxins in lung epithelial cell.     

Large-scale preparation of active caspase-3 in E. coli by designing its thrombin-activatable precursors

Background

Combination of CHD (chromo-helicase-DNA binding protein)-specific polymerase chain reaction (PCR) with electrophoresis (PCR/electrophoresis) is the most common avian molecular sexing technique but it is lab-intensive and gel-required. Gender determination often fails when the difference in length between the PCR products of CHD-Z and CHD-W genes is too short to be resolved.

Results

Here, we are the first to introduce a PCR-melting curve analysis (PCR/MCA) to identify the gender of birds by genomic DNA, which is gel-free, quick, and inexpensive. Spilornis cheela hoya (S. c. hoya) and Pycnonotus sinensis (P. sinensis) were used to illustrate this novel molecular sexing technique. The difference in the length of CHD genes in S. c. hoya and P. sinensis is 13-, and 52-bp, respectively. Using Griffiths' P2/P8 primers, molecular sexing failed both in PCR/electrophoresis of S. c. hoya and in PCR/MCA of S. c. hoya and P. sinensis. In contrast, we redesigned sex-specific primers to yield 185- and 112-bp PCR products for the CHD-Z and CHD-W genes of S. c. hoya, respectively, using PCR/MCA. Using this specific primer set, at least 13 samples of S. c. hoya were examined simultaneously and the Tm peaks of CHD-Z and CHD-W PCR products were distinguished.

Conclusion

In this study, we introduced a high-throughput avian molecular sexing technique and successfully applied it to two species. This new method holds a great potential for use in high throughput sexing of other avian species, as well.     

The COG database: an updated version includes eukaryotes

Background

The availability of multiple, essentially complete genome sequences of prokaryotes and eukaryotes spurred both the demand and the opportunity for the construction of an evolutionary classification of genes from these genomes. Such a classification system based on orthologous relationships between genes appears to be a natural framework for comparative genomics and should facilitate both functional annotation of genomes and large-scale evolutionary studies.

Results

We describe here a major update of the previously developed system for delineation of Clusters of Orthologous Groups of proteins (COGs) from the sequenced genomes of prokaryotes and unicellular eukaryotes and the construction of clusters of predicted orthologs for 7 eukaryotic genomes, which we named KOGs after eukaryotic orthologous groups. The COG collection currently consists of 138,458 proteins, which form 4873 COGs and comprise 75% of the 185,505 (predicted) proteins encoded in 66 genomes of unicellular organisms. The eukaryotic orthologous groups (KOGs) include proteins from 7 eukaryotic genomes: three animals (the nematode Caenorhabditis elegans, the fruit fly Drosophila melanogaster and Homo sapiens), one plant, Arabidopsis thaliana, two fungi (Saccharomyces cerevisiae and Schizosaccharomyces pombe), and the intracellular microsporidian parasite Encephalitozoon cuniculi. The current KOG set consists of 4852 clusters of orthologs, which include 59,838 proteins, or ~54% of the analyzed eukaryotic 110,655 gene products. Compared to the coverage of the prokaryotic genomes with COGs, a considerably smaller fraction of eukaryotic genes could be included into the KOGs; addition of new eukaryotic genomes is expected to result in substantial increase in the coverage of eukaryotic genomes with KOGs. Examination of the phyletic patterns of KOGs reveals a conserved core represented in all analyzed species and consisting of ~20% of the KOG set. This conserved portion of the KOG set is much greater than the ubiquitous portion of the COG set (~1% of the COGs). In part, this difference is probably due to the small number of included eukaryotic genomes, but it could also reflect the relative compactness of eukaryotes as a clade and the greater evolutionary stability of eukaryotic genomes.

Conclusion

The updated collection of orthologous protein sets for prokaryotes and eukaryotes is expected to be a useful platform for functional annotation of newly sequenced genomes, including those of complex eukaryotes, and genome-wide evolutionary studies.     

The role of polyamines in protein-dependent hypoxic tolerance of Drosophila

Background

Chronic hypoxia is a major component of ischemic diseases such as stroke or myocardial infarction. Drosophila is more tolerant to hypoxia than most mammalian species. It is considered as a useful model organism to identify new mechanisms of hypoxic tolerance. The hypoxic tolerance of flies has previously been reported to be enhanced by low protein diets. This study analyses the mechanisms involved.

Results

Feeding adult Drosophila on a yeast diet dramatically reduced their longevities under chronic hypoxic conditions (5% O₂). Mean and maximum longevities became close to the values observed for starving flies. The action of dietary yeast was mimicked by a whole casein hydrolysate and by anyone of the 20 natural amino acids that compose proteins. It was mimicked by amino acid intermediates of the urea cycle such as L-citrulline and L-ornithine, and by polyamines (putrescine, spermidine and spermine). α-difluoromethylornithine, a specific inhibitor of ornithine decarboxylase, partially protected hypoxic flies from amino acid toxicity but not from polyamine toxicity. N¹-guanyl-1,7 diaminoheptane, a specific inhibitor of eIF5A hypusination, partially relieved the toxicities of both amino acids and polyamines.

Conclusion

Dietary amino acids reduced the longevity of chronically hypoxic flies fed on a sucrose diet. Pharmacological evidence suggests that the synthesis of polyamines and the hypusination of eIF5A contributed to the life-shortening effect of dietary amino acids.     

6-OHDA generated ROS induces DNA damage and p53- and PUMA-dependent cell death

Background

Parkinson's disease (PD) is characterized by the selective loss of dopaminergic neurons in the substantia nigra (SN), resulting in tremor, rigidity, and bradykinesia. Although the etiology is unknown, insight into the disease process comes from the dopamine (DA) derivative, 6-hydroxydopamine (6-OHDA), which produces PD-like symptoms. Studies show that 6-OHDA activates stress pathways, such as the unfolded protein response (UPR), triggers mitochondrial release of cytochrome-c, and activates caspases, such as caspase-3. Because the BH3-only protein, Puma (p53-upregulated mediator of apoptosis), is activated in response to UPR, it is thought to be a link between cell stress and apoptosis.

Results

To test the hypothesis that Puma serves such a role in 6-OHDA-mediated cell death, we compared the response of dopaminergic neurons from wild-type and Puma-null mice to 6-OHDA. Results indicate that Puma is required for 6-OHDA-induced cell death, in primary dissociated midbrain cultures as well as in vivo. In these cultures, 6-OHDA-induced DNA damage and p53 were required for 6-OHDA-induced cell death. In contrast, while 6-OHDA led to upregulation of UPR markers, loss of ATF3 did not protect against 6-OHDA.

Conclusions

Together, our results indicate that 6-OHDA-induced upregulation of Puma and cell death are independent of UPR. Instead, p53 and DNA damage repair pathways mediate 6-OHDA-induced toxicity.     

The Smc complexes in DNA damage response

Background

Caenorhabditis elegans sarcomeres have been studied extensively utilizing both forward and reverse genetic techniques to provide insight into muscle development and the mechanisms behind muscle contraction. A previous genetic screen investigating early muscle development produced 13 independent mutant genes exhibiting a Pat (paralyzed and arrested elongation at the two-fold length of embryonic development) muscle phenotype. This study reports the identification and characterization of one of those genes, pat-9.

Results

Positional cloning, reverse genetics, and plasmid rescue experiments were used to identify the predicted C. elegans gene T27B1.2 (recently named ztf-19) as the pat-9 gene. Analysis of pat-9 showed it is expressed early in development and within body wall muscle lineages, consistent with a role in muscle development and producing a Pat phenotype. However, unlike most of the other known Pat gene family members, which encode structural components of muscle attachment sites, PAT-9 is an exclusively nuclear protein. Analysis of the predicted PAT-9 amino acid sequence identified one putative nuclear localization domain and three C2H2 zinc finger domains. Both immunocytochemistry and PAT-9::GFP fusion expression confirm that PAT-9 is primarily a nuclear protein and chromatin immunoprecipitation (ChIP) experiments showed that PAT-9 is present on certain gene promoters.

Conclusions

We have shown that the T27B1.2 gene is pat-9. Considering the Pat-9 mutant phenotype shows severely disrupted muscle attachment sites despite PAT-9 being a nuclear zinc finger protein and not a structural component of muscle attachment sites, we propose that PAT-9 likely functions in the regulation of gene expression for some necessary structural or regulatory component(s) of the muscle attachment sites.     

Genome-wide acceleration of protein evolution in flies (Diptera)

Background

The rate of molecular evolution varies widely between proteins, both within and among lineages. To what extent is this variation influenced by genome-wide, lineage-specific effects? To answer this question, we assess the rate variation between insect lineages for a large number of orthologous genes.

Results

When compared to the beetle Tribolium castaneum, we find that the stem lineage of flies and mosquitoes (Diptera) has experienced on average a 3-fold increase in the rate of evolution. Pairwise gene comparisons between Drosophila and Tribolium show a high correlation between evolutionary rates of orthologous proteins.

Conclusion

Gene specific divergence rates remain roughly constant over long evolutionary times, modulated by genome-wide, lineage-specific effects. Among the insects analysed so far, it appears that the Tribolium genes show the lowest rates of divergence. This has the practical consequence that homology searches for human genes yield significantly better matches in Tribolium than in Drosophila. We therefore suggest that Tribolium is better suited for comparisons between phyla than the widely employed dipterans.