Parallel computation of a maximum-likelihood estimator of a physical map

Reconstructing a physical map of a chromosome from a genomic library presents a central computational problem in genetics. Physical map reconstruction in the presence of errors is a problem of high computational complexity that provides the motivation for parallel computing. Parallelization strategies for a maximum-likelihood estimation-based approach to physical map reconstruction are presented. The estimation procedure entails a gradient descent search for determining the optimal spacings between probes for a given probe ordering. The optimal probe ordering is determined using a stochastic optimization algorithm such as simulated annealing or microcanonical annealing. A two-level parallelization strategy is proposed wherein the gradient descent search is parallelized at the lower level and the stochastic optimization algorithm is simultaneously parallelized at the higher level. Implementation and experimental results on a distributed-memory multiprocessor cluster running the parallel virtual machine (PVM) environment are presented using simulated and real hybridization data.     

Uniformly minimum variance unbiased estimation of gene diversity

Gene diversity is an important measure of genetic variability in inbred populations. The survival of species in changing environments depends on, among other factors, the genetic variability of the population. In this communication, I have derived the uniformly minimum variance unbiased estimator of gene diversity. The proposed estimator of gene diversity does not assume that the inbreeding coefficient is known. I have also provided the approximate variance of this estimator according to Fisher's method. In addition, I have developed a numerical resampling-based method for obtaining variances and confidence intervals based on the maximum likelihood estimator and the uniformly minimum variance unbiased estimator. Efficiency in estimation of the gene diversity based on these two estimators is discussed. In accordance with the simulation results, I found that the uniformly minimum variance estimator developed in this report is more accurate for estimation of gene diversity than the maximum likelihood estimator.     

Complex segregation analysis reveals a multigene model for lung cancer

Lung cancer risk is largely attributed to tobacco exposure, but genetic predisposition also plays an etiologic role. Several studies have investigated the involvement of genetic predisposition in lung cancer aggregation in affected families, although with inconsistent results. Some studies have provided evidence for Mendelian inheritance, whereas others have suggested that environmental models are most appropriate for lung cancer aggregation in families. To examine the genetic basis of lung cancer, we performed segregation analysis on 14,378 individuals from 1,561 lung cancer case families, allowing for the effects of smoking, sex, and age. Both a Mendelian decreasing model and a Mendelian codominant model were found to be the best fitting models for susceptibility. However, when we modeled age-of-onset, all Mendelian models and the environmental model were rejected suggesting that multiple genetic factors (possibly multiple genetic loci and interactions) contribute to the age-of-onset of lung cancer. The results provide evidence that multiple genetic factors contribute to lung cancer and may act as a guide in further studies to localize susceptibility genes in lung cancer.     

Optimization of redox reactions employing whole cell biocatalysis

The ability of Saccharomyces cerevisiae to catalyse the reduction reaction of carboxylic acids into alcohols is described. Earlier reports have led to the characterization of the reduction of carbonyl groups into alcohols mediated by the enzyme alcohol dehydrogenase. We investigated the ability of this organism to catalyse the said conversion using the carboxylic acids, acetic acid and butyric acid. In the absence of any previous characterization, whole cell catalysis proved effective. The uptake of these acids from the medium was estimated using a plate assay method involving litmus-agar. The plate assay was found to be a convenient and extremely adaptable method for quantitation of acids in organic as well as aqueous medium. The comparison of existing paradigms in pure protein catalysis with whole cells catalysis proved anomalous. We report that it is solvent toxicity rather than hydrophobic index that correlates with the activity observed in non-aqueous conditions for whole cell biocatalysis. Reduction of acetic acid as well as butyric acid occurred, with efficiency of reaction with butyric acid being marginally higher. The reduction therefore occurs for both the short chain carboxylic acids used in this study. We therefore illustrate the reduction route of acids into alcohols and propose a model two-step pathway for the reaction. Process optimization may be further attempted to enhance the presently moderate reaction efficiencies. Steps made in the direction by studying the pH dependency and use of sacrificial substrate have yielded encouraging results.     

Ignoring linkage disequilibrium among tightly linked markers induces false-positive evidence of linkage for affected sib pair analysis

Most multipoint linkage programs assume linkage equilibrium among the markers being studied. The assumption is appropriate for the study of sparsely spaced markers with intermarker distances exceeding a few centimorgans, because linkage equilibrium is expected over these intervals for almost all populations. However, with recent advancements in high-throughput genotyping technology, much denser markers are available, and linkage disequilibrium (LD) may exist among the markers. Applying linkage analyses that assume linkage equilibrium to dense markers may lead to bias. Here, we demonstrated that, when some or all of the parental genotypes are missing, assuming linkage equilibrium among tightly linked markers where strong LD exists can cause apparent oversharing of multipoint identity by descent (IBD) between sib pairs and false-positive evidence for multipoint model-free linkage analysis of affected sib pair data. LD can also mimic linkage between a disease locus and multiple tightly linked markers, thus causing false-positive evidence of linkage using parametric models, particularly when heterogeneity LOD score approaches are applied. Bias can be eliminated by inclusion of parental genotype data and can be reduced when additional unaffected siblings are included in the analysis.     

Joint effects of germ-line <Emphasis Type="Italic">TP53</Emphasis> mutation, <Emphasis Type="Italic">MDM2</Emphasis> SNP309, and gender on cancer risk in family studies of Li–Fraumeni syndrome

Li–Fraumeni syndrome (LFS) is a rare familial cancer syndrome characterized by early cancer onset, diverse tumor types, and multiple primary tumors. Germ-line TP53 mutations have been identified in most LFS families. A high-frequency single-nucleotide polymorphism, SNP309 (rs2279744), in MDM2 was recently confirmed to be a modifier of cancer risk in several case-series studies: substantially earlier cancer onset was observed in SNP309 G-allele carriers than in wild-type individuals by 7–16 years. However, cancer risk analyses that jointly account for measured hereditary TP53 mutations and MDM2 SNP309 have not been systematically investigated in familial cases. Here, we determined the combined effects of measured TP53 mutations, MDM2 SNP309, and gender and their interactions simultaneously in LFS families. We used the method that is designed for extended pedigrees and structured for age-specific risk models based on Cox proportional hazards regression. We analyzed the cancer incidence in 19 extended pedigrees with germ-line TP53 mutations ascertained through the clinical LFS phenotype. The dataset consisted of 463 individuals with 129 TP53 mutation carriers. Our analyses showed that the TP53 germ-line mutation and its interaction with gender were strongly associated with familial cancer incidence and that the association between MDM2 SNP309 and increased cancer risk was modest. In contrast with several case-series studies, the interaction between MDM2 SNP309 and TP53 mutation was not statistically significant in our LFS family cohort. Our results showed that SNP309 G-alleles were associated with accelerated tumor formation in both carriers and non-carriers of germ-line TP53 mutations.     

Malsoor Virus,a Novel Bat Phlebovirus,Is Closely Related to Severe Fever with Thrombocytopenia Syndrome Virus and Heartland Virus

During a survey in the year 2010, a novel phlebovirus was isolated from the Rousettus leschenaultii species of bats in western India. The virus was identified by electron microscopy from infected Vero E6 cells. Phylogenic analysis of the complete genome showed its close relation to severe fever with thrombocytopenia syndrome (SFTS) and Heartland viruses, which makes it imperative to further study its natural ecology and potential as a novel emerging zoonotic virus.