Efficient large-scale protein sequence comparison and gene matching to identify orthologs and co-orthologs

Broadly, computational approaches for ortholog assignment is a three steps process: (i) identify all putative homologs between the genomes, (ii) identify gene anchors and (iii) link anchors to identify best gene matches given their order and context. In this article, we engineer two methods to improve two important aspects of this pipeline [specifically steps (ii) and (iii)]. First, computing sequence similarity data [step (i)] is a computationally intensive task for large sequence sets, creating a bottleneck in the ortholog assignment pipeline. We have designed a fast and highly scalable sort-join method (afree) based on k-mer counts to rapidly compare all pairs of sequences in a large protein sequence set to identify putative homologs. Second, availability of complex genomes containing large gene families with prevalence of complex evolutionary events, such as duplications, has made the task of assigning orthologs and co-orthologs difficult. Here, we have developed an iterative graph matching strategy where at each iteration the best gene assignments are identified resulting in a set of orthologs and co-orthologs. We find that the afree algorithm is faster than existing methods and maintains high accuracy in identifying similar genes. The iterative graph matching strategy also showed high accuracy in identifying complex gene relationships. Standalone afree available from http://vbc.med.monash.edu.au/～kmahmood/afree. EGM2, complete ortholog assignment pipeline (including afree and the iterative graph matching method) available from http://vbc.med.monash.edu.au/～kmahmood/EGM2.     

Characterization of a novel missense mutation in the prodomain of GDF5, which underlies brachydactyly type C and mild Grebe type chondrodysplasia in a large Pakistani family

All TGF-beta family members have a prodomain that is important for secretion. Lack of secretion of a TGF-beta family member GDF5 is known to underlie some skeletal abnormalities, such as brachydactyly type C that is characterized by a huge and unexplained phenotypic variability. To search for potential phenotypic modifiers regulating secretion of GDF5, we compared cells overexpressing wild type (Wt) GDF5 and GDF5 with a novel mutation in the prodomain identified in a large Pakistani family with Brachydactyly type C and mild Grebe type chondrodyslplasia (c527T>C; p.Leu176Pro). Initial in vitro expression studies revealed that the p.Leu176Pro mutant (Mut) GDF5 was not secreted outside the cells. We subsequently showed that GDF5 was capable of forming a complex with latent transforming growth factor binding proteins, LTBP1 and LTBP2. Furthermore, secretion of LTBP1 and LTBP2 was severely impaired in cells expressing the Mut-GDF5 compared to Wt-GDF5. Finally, we demonstrated that secretion of Wt-GDF5 was inhibited by the Mut-GDF5, but only when LTBP (LTBP1 or LTBP2) was co-expressed. Based on these findings, we suggest a novel model, where the dosage of secretory co-factors or stabilizing proteins like LTBP1 and LTBP2 in the microenvironment may affect the extent of GDF5 secretion and thereby function as modifiers in phenotypes caused by GDF5 mutations.     

DFT studies of hydrocarbon combustion on metal surfaces

Catalytic combustion of hydrocarbons is an important technology to produce energy. Compared to conventional flame combustion, the catalyst enables this process to operate at lower temperatures; hence, reducing the energy required for efficient combustion. The reaction and activation energies of direct combustion of hydrocarbons (CH?→?C?+?H) on a series of metal surfaces were investigated using density functional theory (DFT). The data obtained for the Ag, Au, Al, Cu, Rh, Pt, and Pd surfaces were used to investigate the validity of the Brønsted-Evans-Polanyi (BEP) and transition state scaling (TSS) relations for this reaction on these surfaces. These relations were found to be valid (R²?=?0.94 for the BEP correlation and R²?=?1.0 for the TSS correlation) and were therefore used to estimate the energetics of the combustion reaction on Ni, Co, and Fe surfaces. It was found that the estimated transition state and activation energies (E_TS?=??69.70 eV and E_a?=?1.20 eV for Ni, E_TS?=??87.93 eV and E_a?=?1.08 eV for Co and E_TS?=??92.45 eV and E_a?=?0.83 eV for Fe) are in agreement with those obtained by DFT calculations (E_TS?=??69.98 eV and E_a?=?1.23 eV for Ni, E_TS?=??87.88 eV and E_a?=?1.08 eV for Co and E_TS?=??92.57 eV and E_a?=?0.79 eV for Fe). Therefore, these relations can be used to predict energetics of this reaction on these surfaces without doing the time consuming transition state calculations. Also, the calculations show that the activation barrier for CH dissociation decreases in the order Ag ? Au ? Al ? Cu ? Pt ? Pd ? Ni?>?Co?>?Rh?>?Fe.     

Adaptive response to low dose of EMS or MMS in human peripheral blood lymphocytes

Human peripheral blood lymphocytes stimulated in vitro for 6 hr were exposed to a low (conditioning) dose of ethyl methanesulfonate (EMS; 1.5 x 10(-4) M) or methyl methanesulfonate (MMS; 1.5 x 10(-5) M). After 6 hr, the cells were treated with a high (challenging) concentration of the same agent (1.5 x 10(-3) M EMS or 1.5 x 10(-4) M MMS). The cells that received both conditioning and challenging doses became less sensitive to the induction of sister chromatid exchanges (SCEs) than those which did not receive the pretreatment with EMS or MMS. They responded with lower frequencies of SCEs. This suggests that conditioning dose of EMS or MMS has offered the lymphocytes to have decreased SCEs. This led to the realization that pre-exposure of lymphocytes to low dose can cause the induction of repair activity. This is a clear indication of the existence of adaptive response induced by alkylating agents whether it is ethylating or methylating in human lymphocytes in vitro.     

Economic activity in the South-Asian population in Britain: the impact of ethnicity,religion, and class

This paper expands the existing literature on ethnicity and economic activity in Britain by studying the impact of religion and class. It argues that while the class location of the different South-Asian groups is important in determining their labour market outcomes, it does not operate independently from ethnicity; rather it is highly influenced by ethnicity in the process of determining the labour market participation of these groups. We use data obtained from the 2001 UK Census on Indian, Pakistani, and Bangladeshi men and women aged between twenty and twenty nine. Our findings confirm that class structure of the South-Asian groups is highly ethnicized, in that the ethno-religious background and class are interwoven to the extent that the separation between them is not easy, if not impossible.     

In-vitro and In-vivo management of Meloidogyne incognita (Kofoid and White) Chitwood and Rhizoctonia bataticola (Taub.) Butler in cotton using organic’s

Root-knot nematodes Meloidogyne incognita (Kofoid and White) Chitwood and Rhizoctonia bataticola (Taub.) Butler, fungus, are very dangerous root damaging pathogens. Present study was planned to establish a chemical control of these root deteriorating pathogens under lab conditions as well as in field. Maximum death rate of nematode juveniles and minimum numbers of nematode eggs hatched were recorded in plates treated with Cadusafos (Rugby® 100G) @12 g/100 ml and Cartap® (4% G) @9g/100 ml. Chemical treatment of Rhizoctonia bataticola with Trifloxystrobin + Tebuconazole (Nativo®) @0.2 g/100 ml and Mancozeb + Matalaxyl (Axiom) @0.25 g/100 ml significantly controlled the mycelial growth in plates. The best treatments tested in laboratory were applied in field as protective and curative treatments. Results proved that chemical control of root-knot nematode and root rot fungi by tested chemicals at recommended time and dose is a significant management technique under field conditions.     

Prediction of Protein-Protein Interactions Using Protein Signature Profiling

Protein domains are conserved and functionally independent structures that play an important role in interactions among related proteins. Domain-domain inter- actions have been recently used to predict protein-protein interactions (PPI). In general, the interaction probability of a pair of domains is scored using a trained scoring function. Satisfying a threshold, the protein pairs carrying those domains are regarded as "interacting". In this study, the signature contents of proteins were utilized to predict PPI pairs in Saccharomyces cerevisiae, Caenorhabditis ele- gans, and Homo sapiens. Similarity between protein signature patterns was scored and PPI predictions were drawn based on the binary similarity scoring function. Results show that the true positive rate of prediction by the proposed approach is approximately 32% higher than that using the maximum likelihood estimation method when compared with a test set, resulting in 22% increase in the area un- der the receiver operating characteristic (ROC) curve. When proteins containing one or two signatures were removed, the sensitivity of the predicted PPI pairs in- creased significantly. The predicted PPI pairs are on average 11 times more likely to interact than the random selection at a confidence level of 0.95, and on aver- age 4 times better than those predicted by either phylogenetic profiling or gene expression profiling.     

Nitrogen removal from secondary effluent by using integrated constructed wetland system

The treatment capacity of an integrated constructed wetland system (CWS) that was designed to reduce nitrogen (N) from secondary effluent was explored. The integrated CWS consisted of vertical-flow constructed wetland, floating bed and sand filter. The vertical-flow wetland was filled with gravel, steel slag and peat from the bottom to the top. Vetiver zizanioides was selected to grow in the vertical-flow constructed wetland and Coix lacrymajobi L. was grown in the floating bed. The results showed that the integrated CWS displayed superior removal efficiency for nitrate nitrogen (NO₃⁻-N), ammonia nitrogen (NH₄⁺-N), nitrite nitrogen (NO₂⁻-N), and total nitrogen (TN). The average NO₃⁻-N, NO₂⁻-N, NH₄⁺-N and TN removal efficiencies of the integrated CWS were 98.83%, 95.60%, 98.05% and 92.41%, respectively, during the whole experimental operation. The integrated CWS may have a good potential for removing N from secondary effluent.     

Fgf3 and Fgf8 are required together for formation of the otic placode and vesicle

Fgf3 has long been implicated in otic placode induction and early development of the otocyst; however, the results of experiments in mouse and chick embryos to determine its function have proved to be conflicting. In this study, we determined fgf3 expression in relation to otic development in the zebrafish and used antisense morpholino oligonucleotides to inhibit Fgf3 translation. Successful knockdown of Fgf3 protein was demonstrated and this resulted in a reduction of otocyst size together with reduction in expression of early markers of the otic placode. fgf3 is co-expressed with fgf8 in the hindbrain prior to otic induction and, strikingly, when Fgf3 morpholinos were co-injected together with Fgf8 morpholinos, a significant number of embryos failed to form otocysts. These effects were made manifest at early stages of otic development by an absence of early placode markers (pax2.1 and dlx3) but were not accompanied by effects on cell division or death. The temporal requirement for Fgf signalling was established as being between 60% epiboly and tailbud stages using the Fgf receptor inhibitor SU5402. However, the earliest molecular event in induction of the otic territory, pax8 expression, did not require Fgf signalling, indicating an inductive event upstream of signalling by Fgf3 and Fgf8. We propose that Fgf3 and Fgf8 are required together for formation of the otic placode and act during the earliest stages of its induction.