Common genetic risk factors in ASD and ADHD co-occurring families.

Autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD) are two major neurodevelopmental disorders that frequently co-occur. However, the genetic mechanism of the co-occurrence remains unclear. The New Jersey Language and Autism Genetics Study (NJLAGS) collected more than 100 families with at least one member affected by ASD. NJLAGS families show a high prevalence of ADHD and provide a good opportunity to study shared genetic risk factors for ASD and ADHD. The linkage study of the NJLAGS families revealed regions on chromosomes 12 and 17 that are significantly associated with ADHD. Using whole-genome sequencing data on 272 samples from 73 NJLAGS families, we identified potential risk genes for ASD and ADHD. Within the linkage regions, we identified 36 genes that are associated with ADHD using a pedigree-based gene prioritization approach. KDM6B (Lysine Demethylase 6B) is the highest-ranking gene, which is a known risk gene for neurodevelopmental disorders, including ASD and ADHD. At the whole-genome level, we identified 207 candidate genes from the analysis of both small variants and structure variants, including both known and novel genes. Using enrichment and protein-protein interaction network analyses, we identified gene ontology terms and pathways enriched for ASD and ADHD candidate genes, such as cilia function and cation channel activity. Candidate genes and pathways identified in our study improve the understanding of the genetic etiology of ASD and ADHD and will lead to new diagnostic or therapeutic interventions for ASD and ADHD in the future.

Structural Variations Contribute to the Genetic Etiology of Autism Spectrum Disorder and Language Impairments.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by restrictive interests and/or repetitive behaviors and deficits in social interaction and communication. ASD is a multifactorial disease with a complex polygenic genetic architecture. Its genetic contributing factors are not yet fully understood, especially large structural variations (SVs). In this study, we aimed to assess the contribution of SVs, including copy number variants (CNVs), insertions, deletions, duplications, and mobile element insertions, to ASD and related language impairments in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Within the cohort, ~77% of the families contain SVs that followed expected segregation or de novo patterns and passed our filtering criteria. These SVs affected 344 brain-expressed genes and can potentially contribute to the genetic etiology of the disorders. Gene Ontology and protein-protein interaction network analysis suggested several clusters of genes in different functional categories, such as neuronal development and histone modification machinery. Genes and biological processes identified in this study contribute to the understanding of ASD and related neurodevelopment disorders.

Integrative survival analysis of breast cancer with gene expression and DNA methylation data.

MOTIVATION: Integrative multi-feature fusion analysis on biomedical data has gained much attention recently. In breast cancer, existing studies have demonstrated that combining genomic mRNA data and DNA methylation data can better stratify cancer patients with distinct prognosis than using single signature. However, those existing methods are simply combining these gene features in series and have ignored the correlations between separate omics dimensions over time. RESULTS: In the present study, we propose an adaptive multi-task learning method, which combines the Cox loss task with the ordinal loss task, for survival prediction of breast cancer patients using multi-modal learning instead of performing survival analysis on each feature dataset. First, we use local maximum quasi-clique merging (lmQCM) algorithm to reduce the mRNA and methylation feature dimensions and extract cluster eigengenes respectively. Then, we add an auxiliary ordinal loss to the original Cox model to improve the ability to optimize the learning process in training and regularization. The auxiliary loss helps to reduce the vanishing gradient problem for earlier layers and helps to decrease the loss of the primary task. Meanwhile, we use an adaptive weights approach to multi-task learning which weighs multiple loss functions by considering the homoscedastic uncertainty of each task. Finally, we build an ordinal cox hazards model for survival analysis and use long short-term memory (LSTM) method to predict patients' survival risk. We use the cross-validation method and the concordance index (C-index) for assessing the prediction effect. Stringent cross-verification testing processes for the benchmark dataset and two additional datasets demonstrate that the developed approach is effective, achieving very competitive performance with existing approaches. AVAILABILITY AND IMPLEMENTATION: https://github.com/bhioswego/ML_ordCOX.

MiR-29b is associated with perinatal inflammation in extremely preterm infants.

BACKGROUND: Inflammation is strongly associated with premature birth and neonatal morbidities. Increases in infant haptoglobin, haptoglobin-related protein (Hp&HpRP), and interleukin-6 (IL-6) levels are indicators of intra-amniotic inflammation (IAI) and have been linked to poor neonatal outcomes. Inflammation causes epigenetic changes, specifically suppression of miR-29 expression. The current study sought to determine whether miR-29b levels in cord blood or neonatal venous blood are associated with IAI, identified by elevated IL-6 and Hp, and subsequent clinical morbidities in the infant. METHODS: We tested 92 cord blood samples from premature newborns and 18 venous blood samples at 36 weeks corrected gestational age. MiR-29b, Hp&HpRP, and IL-6 were measured by polymerase chain reaction and enzyme-linked immunosorbent assay, respectively. RESULTS: Decreased levels of miR-29b were observed in infants exposed to IAI with elevated Hp&HpRP and IL-6 levels and in infants delivered by spontaneous preterm birth. Lower miR-29 levels were also observed in women diagnosed with histological chorioamnionitis or funisitis and in infants with cerebral palsy. Higher levels of miR-29 were measured in infants small for gestational age and in venous samples from older infants. CONCLUSIONS: MiR-29 may be an additional biomarker of IAI and a potential therapeutic target for treating poor newborn outcomes resulting from antenatal exposure to IAI. IMPACT: Decreases in miR-29b are associated with intrauterine inflammation. Hp&HpRP increases are associated with decreased miR-29b. MiR-29b may be an additional biomarker for neonatal outcomes and a potential therapeutic target for intrauterine inflammation.

Phosphorylation of Serine 225 in Hepatitis C Virus NS5A Regulates Protein-Protein Interactions.

Hepatitis C virus (HCV) nonstructural protein 5A (NS5A) is a phosphoprotein that plays key, yet poorly defined, roles in both virus genome replication and virion assembly/release. It has been proposed that differential phosphorylation could act as a switch to regulate the various functions of NS5A; however, the mechanistic details of the role of this posttranslational modification in the virus life cycle remain obscure. We previously reported (D. Ross-Thriepland, J. Mankouri, and M. Harris, J Virol 89:3123-3135, 2015, doi:10.1128/JVI.02995-14) a role for phosphorylation at serine 225 (S225) of NS5A in the regulation of JFH-1 (genotype 2a) genome replication. A phosphoablatant (S225A) mutation resulted in a 10-fold reduction in replication and a perinuclear restricted distribution of NS5A, whereas the corresponding phosphomimetic mutation (S225D) had no phenotype. To determine the molecular mechanisms underpinning this phenotype we conducted a label-free proteomics approach to identify cellular NS5A interaction partners. This analysis revealed that the S225A mutation disrupted the interactions of NS5A with a number of cellular proteins, in particular the nucleosome assembly protein 1-like protein 1 (NAP1L1), bridging integrator 1 (Bin1, also known as amphiphysin II), and vesicle-associated membrane protein-associated protein A (VAP-A). These interactions were validated by immunoprecipitation/Western blotting, immunofluorescence, and proximity ligation assay. Importantly, small interfering RNA (siRNA)-mediated knockdown of NAP1L1, Bin1 or VAP-A impaired viral genome replication and recapitulated the perinuclear redistribution of NS5A seen in the S225A mutant. These results demonstrate that S225 phosphorylation regulates the interactions of NS5A with a defined subset of cellular proteins. Furthermore, these interactions regulate both HCV genome replication and the subcellular localization of replication complexes.IMPORTANCE Hepatitis C virus is an important human pathogen. The viral nonstructural 5A protein (NS5A) is the target for new antiviral drugs. NS5A has multiple functions during the virus life cycle, but the biochemical details of these roles remain obscure. NS5A is known to be phosphorylated by cellular protein kinases, and in this study, we set out to determine whether this modification is required for the binding of NS5A to other cellular proteins. We identified 3 such proteins and show that they interacted only with NS5A that was phosphorylated on a specific residue. Furthermore, these proteins were required for efficient virus replication and the ability of NS5A to spread throughout the cytoplasm of the cell. Our results help to define the function of NS5A and may contribute to an understanding of the mode of action of the highly potent antiviral drugs that are targeted to NS5A.

Whole-Cell Detection of C-P Bonds in Bacteria.

Naturally produced molecules possessing a C-P bond, such as phosphonates and phosphinates, remain vastly underexplored. Although success stories like fosfomycin have reinvigorated small molecule phosphonate discovery efforts, bioinformatic analyses predict an enormous unexplored biological reservoir of C-P bond-containing molecules, including those attached to complex macromolecules. However, high polarity, a lack of chromophores, and complex macromolecular association impede phosphonate discovery and characterization. Here we detect widespread transcriptional activation of phosphonate biosynthetic machinery across diverse bacterial phyla and describe the use of solid-state nuclear magnetic resonance to detect C-P bonds in whole cells of representative Gram-negative and Gram-positive bacterial species. These results suggest that phosphonate tailoring is more prevalent than previously recognized and set the stage for elucidating the fascinating chemistry and biology of these modifications.

Hepatitis C virus in vitro replication is efficiently inhibited by acridone Fac4.

Hepatitis C virus (HCV) affects about 170 million people worldwide. The current treatment has a high cost and variable response rates according to the virus genotype. Acridones, a group of compounds extracted from natural sources, showed potential antiviral actions against HCV. Thus, this study aimed to evaluate the effect of a panel of 14 synthetic acridones on the HCV life cycle. The compounds were screened using an Huh7.5 cell line stably harbouring the HCV genotype 2a subgenomic replicon SGR-Feo-JFH-1. Cells were incubated in the presence or absence of compounds for 72 h and cell viability and replication levels were assessed by MTT and luciferase assays, respectively. At a concentration of 5 µM the acridone Fac4 exhibited a >90 % inhibition of HCV replication with no effect on cell viability. The effects of Fac4 on virus replication, entry and release steps were evaluated in Huh7.5 cells infected with the JFH-1 isolate of HCV (HCVcc). Fac4 inhibited JFH-1 replication to approximately 70 %, while no effect was observed on virus entry. The antiviral activity of Fac4 was also observed on viral release, with almost 80 % of inhibition. No inhibitory effect was observed against genotype 3 replication. Fac4 was able to intercalate into dsRNA, however did not inhibit NS5B polymerase activity or translation driven by the HCV IRES. Although its mode of action is partly understood, Fac4 presents significant inhibition of HCV replication and can therefore be considered as a candidate for the development of a future anti-HCV treatment.

Behavioral and Molecular Genetics of Reading-Related AM and FM Detection Thresholds.

Auditory detection thresholds for certain frequencies of both amplitude modulated (AM) and frequency modulated (FM) dynamic auditory stimuli are associated with reading in typically developing and dyslexic readers. We present the first behavioral and molecular genetic characterization of these two auditory traits. Two extant extended family datasets were given reading tasks and psychoacoustic tasks to determine FM 2 Hz and AM 20 Hz sensitivity thresholds. Univariate heritabilities were significant for both AM (h 2  = 0.20) and FM (h 2  = 0.29). Bayesian posterior probability of linkage (PPL) analysis found loci for AM (12q, PPL = 81 %) and FM (10p, PPL = 32 %; 20q, PPL = 65 %). Bivariate heritability analyses revealed that FM is genetically correlated with reading, while AM was not. Bivariate PPL analysis indicates that FM loci (10p, 20q) are not also associated with reading.

Early events in the generation of autophagosomes are required for the formation of membrane structures involved in hepatitis C virus genome replication.

Hepatitis C virus (HCV) infection has been shown to induce autophagy but the mechanisms underpinning this process remain to be elucidated. Induction of autophagy requires the class III phosphatidylinositol 3-kinase, Vps34, which produces phosphatidylinositol 3-phosphate (PI3P) within the endoplasmic reticulum (ER) membrane. This recruits proteins with PI3P binding domains such as the double-FYVE-containing protein 1 (DFCP1). DFCP1 generates cup-shaped protrusions from the ER membrane, termed omegasomes, which provide a platform for the production of autophagosomes. Here we present data demonstrating that both Vps34 and DFCP1 are required for HCV genome replication, in the context of both a subgenomic replicon and virus infection, but did not affect virus entry or initial translation. Using live cell fluorescence microscopy we demonstrated that early during HCV infection the nascent viral genome replication complexes (identified by using non-structural protein NS5A as a marker) transiently colocalize with DFCP1-positive punctae (omegasomes), before the two structures move apart from each other. This observation is reminiscent of the transient association of LC3 and DFCP1 during omegasome formation, and therefore we propose that omegasomes are utilized by HCV to generate the double-membrane vesicles which are the hallmark of HCV replication complexes.

New biphenyl iminium salt catalysts for highly enantioselective asymmetric epoxidation: role of additional substitution and dihedral angle.

New biaryl iminium salt catalysts for enantioselective alkene epoxidation containing additional substitution in the heterocyclic ring are reported. The effects upon conformation and enantioselectivity of this additional substitution, and the influence of dihedral angle in these systems, has been investigated using a synthetic approach supported by density functional theory. Enantioselectivities of up to 97% ee were observed.

An integrative analysis of regional gene expression profiles in the human brain.

Studies of the brain's transcriptome have become prominent in recent years, resulting in an accumulation of datasets with somewhat distinct attributes. These datasets, which are often analyzed only in isolation, also are often collected with divergent goals, which are reflected in their sampling properties. While many researchers have been interested in sampling gene expression in one or a few brain areas in a large number of subjects, recent efforts from the Allen Institute for Brain Sciences and others have focused instead on dense neuroanatomical sampling, necessarily limiting the number of individual donor brains studied. The purpose of the present work is to develop methods that draw on the complementary strengths of these two types of datasets for study of the human brain, and to characterize the anatomical specificity of gene expression profiles and gene co-expression networks derived from human brains using different specific technologies. The approach is applied using two publicly accessible datasets: (1) the high anatomical resolution Allen Human Brain Atlas (AHBA, Hawrylycz et al., 2012) and (2) a relatively large sample size, but comparatively coarse neuroanatomical dataset described previously by Gibbs et al. (2010). We found a relatively high degree of correspondence in differentially expressed genes and regional gene expression profiles across the two datasets. Gene co-expression networks defined in individual brain regions were less congruent, but also showed modest anatomical specificity. Using gene modules derived from the Gibbs dataset and from curated gene lists, we demonstrated varying degrees of anatomical specificity based on two classes of methods, one focused on network modularity and the other focused on enrichment of expression levels. Two approaches to assessing the statistical significance of a gene set's modularity in a given brain region were studied, which provide complementary information about the anatomical specificity of a gene network of interest. Overall, the present work demonstrates the feasibility of cross-dataset analysis of human brain microarray studies, and offers a new approach to annotating gene lists in a neuroanatomical context.

Flavour preferences in youth versus adults: a review.

OBJECTIVE: To understand the available evidence of how children and adults differ in their preferences for flavours that may be used in tobacco products. DATA SOURCES: A total of 474 articles published between 1931 and August 2015 were retrieved through searches conducted in PubMed, EMBASE, Web of Science and PsycINFO. STUDY SELECTION AND EXTRACTION: A 2-phase relevancy review process resulted in the identification of 59 articles and information was extracted by 2 independent reviewers. DATA SYNTHESIS: Findings were grouped by taste and smell preferences, which are important components of overall flavour. For taste, evidence is summarised in the following categories: sweet, salty, sour, bitter, umami and fat; within each of them, findings are organised by age categories. For smell, evidence is summarised as follows: fruit/herbal/spices, tobacco and coffee and other odours. Major findings from this search indicated that sweet preference in children and adolescents was higher than in adults. Examples of preferred food-related tastes and odours for young people included cherry, candy, strawberry, orange, apple and cinnamon. Currently, all these are used to flavour cigars, cartridges for electronic cigarettes, hookah (waterpipe) and smokeless tobacco products. CONCLUSIONS: Infants and children exhibited elevated sweet and salty preference relative to adults. Age-related changes in bitter, sour, umami and fat taste were not clear and more research would be useful. 'Sweet' food odours were highly preferred by children. Tobacco products in flavours preferred by young people may impact tobacco use and initiation, while flavours preferred by adults may impact product switching or dual use.

Molecular Genetic Evidence for Shared Etiology of Autism and Prodigy.

Child prodigies are rare individuals with an exceptional working memory and unique attentional skills that may facilitate the attainment of professional skill levels at an age well before what is observed in the general population. Some characteristics of prodigy have been observed to be quantitatively similar to those observed in autism spectrum disorder (ASD), suggesting possible shared etiology, though objectively validated prodigies are so rare that evidence has been sparse. We performed a family-based genome-wide linkage analysis on 5 nuclear and extended families to search for genetic loci that influence the presence of both prodigy and ASD, assuming that the two traits have the same genetic etiology in the analysis model in order to find shared loci. A shared locus on chromosome 1p31-q21 reached genome-wide significance with two extended family-based linkage methods consisting of the Bayesian PPL method and the LOD score maximized over the trait parameters (i.e., MOD), yielding a simulation-based empirical significance of p = 0.000742 and p = 0.000133, respectively. Within linkage regions, we performed association analysis and assessed if copy number variants could account for the linkage signal. No evidence of specificity for either the prodigy or the ASD trait was observed. This finding suggests that a locus on chromosome 1 increases the likelihood of both prodigy and autism in these families.

Mapping autism risk loci using genetic linkage and chromosomal rearrangements.

Autism spectrum disorders (ASDs) are common, heritable neurodevelopmental conditions. The genetic architecture of ASDs is complex, requiring large samples to overcome heterogeneity. Here we broaden coverage and sample size relative to other studies of ASDs by using Affymetrix 10K SNP arrays and 1,181 [corrected] families with at least two affected individuals, performing the largest linkage scan to date while also analyzing copy number variation in these families. Linkage and copy number variation analyses implicate chromosome 11p12-p13 and neurexins, respectively, among other candidate loci. Neurexins team with previously implicated neuroligins for glutamatergic synaptogenesis, highlighting glutamate-related genes as promising candidates for contributing to ASDs.

Catechol-O-methyltransferase Val158met polymorphism interacts with early experience to predict executive functions in early childhood.

Numerous studies demonstrate that the Methionine variant of the catechol-O-methyltransferase Val158Met polymorphism, which confers less efficient catabolism of catecholamines, is associated with increased focal activation of prefrontal cortex (PFC) and higher levels of executive function abilities. By and large, however, studies of COMT Val158Met have been conducted with adult samples and do not account for the context in which development is occurring. Effects of early adversity on stress response physiology and the inverted U shape relating catecholamine levels to neural activity in PFC indicate the need to take into account early experience when considering relations between genes such as COMT and executive cognitive ability. Consistent with this neurobiology, we find in a prospective longitudinal sample of children and families (N = 1292) that COMT Val158Met interacts with early experience to predict executive function abilities in early childhood. Specifically, the Valine variant of the COMT Val158Met polymorphism, which confers more rather than less efficient catabolism of catecholamines is associated with higher executive function abilities at child ages 48 and 60 months and with faster growth of executive function for children experiencing early adversity, as indexed by cumulative risk factors in the home at child ages 7, 15, 24, and 36 months. Findings indicate the importance of the early environment for the relation between catecholamine genes and developmental outcomes and demonstrate that the genetic moderation of environmental risk is detectable in early childhood.

Addressing the unmet need for visualizing conditional random fields in biological data.

BACKGROUND: The biological world is replete with phenomena that appear to be ideally modeled and analyzed by one archetypal statistical framework - the Graphical Probabilistic Model (GPM). The structure of GPMs is a uniquely good match for biological problems that range from aligning sequences to modeling the genome-to-phenome relationship. The fundamental questions that GPMs address involve making decisions based on a complex web of interacting factors. Unfortunately, while GPMs ideally fit many questions in biology, they are not an easy solution to apply. Building a GPM is not a simple task for an end user. Moreover, applying GPMs is also impeded by the insidious fact that the "complex web of interacting factors" inherent to a problem might be easy to define and also intractable to compute upon. DISCUSSION: We propose that the visualization sciences can contribute to many domains of the bio-sciences, by developing tools to address archetypal representation and user interaction issues in GPMs, and in particular a variety of GPM called a Conditional Random Field(CRF). CRFs bring additional power, and additional complexity, because the CRF dependency network can be conditioned on the query data. CONCLUSIONS: In this manuscript we examine the shared features of several biological problems that are amenable to modeling with CRFs, highlight the challenges that existing visualization and visual analytics paradigms induce for these data, and document an experimental solution called StickWRLD which, while leaving room for improvement, has been successfully applied in several biological research projects. Software and tutorials are available at http://www.stickwrld.org/.

Who is afraid of math? Two sources of genetic variance for mathematical anxiety.

BACKGROUND: Emerging work suggests that academic achievement may be influenced by the management of affect as well as through efficient information processing of task demands. In particular, mathematical anxiety has attracted recent attention because of its damaging psychological effects and potential associations with mathematical problem solving and achievement. This study investigated the genetic and environmental factors contributing to the observed differences in the anxiety people feel when confronted with mathematical tasks. In addition, the genetic and environmental mechanisms that link mathematical anxiety with math cognition and general anxiety were also explored. METHODS: Univariate and multivariate quantitative genetic models were conducted in a sample of 514 12-year-old twin siblings. RESULTS: Genetic factors accounted for roughly 40% of the variation in mathematical anxiety, with the remaining being accounted for by child-specific environmental factors. Multivariate genetic analyses suggested that mathematical anxiety was influenced by the genetic and nonfamilial environmental risk factors associated with general anxiety and additional independent genetic influences associated with math-based problem solving. CONCLUSIONS: The development of mathematical anxiety may involve not only exposure to negative experiences with mathematics, but also likely involves genetic risks related to both anxiety and math cognition. These results suggest that integrating cognitive and affective domains may be particularly important for mathematics and may extend to other areas of academic achievement.

Gene-environment interplay in child psychology and psychiatry--challenges and ways forward.

This special issue in the Journal of Child Psychology and Psychiatry presents several invited articles examining gene-environment interplay in child development and psychopathology. Models of gene-environment interplay have been exhaustively discussed in the literature, including an important contribution by Rutter, Moffitt and Caspi (2006) published in this journal.

Gene × smoking interactions on human brain gene expression: finding common mechanisms in adolescents and adults.

BACKGROUND: Numerous studies have examined gene × environment interactions (G × E) in cognitive and behavioral domains. However, these studies have been limited in that they have not been able to directly assess differential patterns of gene expression in the human brain. Here, we assessed G × E interactions using two publically available datasets to assess if DNA variation is associated with post-mortem brain gene expression changes based on smoking behavior, a biobehavioral construct that is part of a complex system of genetic and environmental influences. METHODS: We conducted an expression quantitative trait locus (eQTL) study on two independent human brain gene expression datasets assessing G × E for selected psychiatric genes and smoking status. We employed linear regression to model the significance of the Gene × Smoking interaction term, followed by meta-analysis across datasets. RESULTS: Overall, we observed that the effect of DNA variation on gene expression is moderated by smoking status. Expression of 16 genes was significantly associated with single nucleotide polymorphisms that demonstrated G × E effects. The strongest finding (p = 1.9 × 10⁻¹¹) was neurexin 3-alpha (NRXN3), a synaptic cell-cell adhesion molecule involved in maintenance of neural connections (such as the maintenance of smoking behavior). Other significant G × E associations include four glutamate genes. CONCLUSIONS: This is one of the first studies to demonstrate G × E effects within the human brain. In particular, this study implicated NRXN3 in the maintenance of smoking. The effect of smoking on NRXN3 expression and downstream behavior is different based upon SNP genotype, indicating that DNA profiles based on SNPs could be useful in understanding the effects of smoking behaviors. These results suggest that better measurement of psychiatric conditions, and the environment in post-mortem brain studies may yield an important avenue for understanding the biological mechanisms of G × E interactions in psychiatry.

Evaluation of a bayesian model integration-based method for censored data.

OBJECTIVE: Non-random missing data can adversely affect family-based linkage detection through loss of power and possible introduction of bias depending on how censoring is modeled. We examined the statistical properties of a previously proposed quantitative trait threshold (QTT) model developed for when censored data can be reasonably inferred to be beyond an unknown threshold. METHODS: The QTT model is a Bayesian model integration approach implemented in the PPL framework that requires neither specification of the threshold nor imputation of the missing data. This model was evaluated under a range of simulated data sets and compared to other methods with missing data imputed. RESULTS: Across the simulated conditions, the addition of a threshold parameter did not change the PPL's properties relative to quantitative trait analysis on non-censored data except for a slight reduction in the average PPL as a reflection of the lowered information content due to censoring. This remained the case for non-normally distributed data and extreme sampling of pedigrees. CONCLUSIONS: Overall, the QTT model showed the smallest loss of linkage information relative to alternative approaches and therefore provides a unique analysis tool that obviates the need for ad hoc imputation of censored data in gene mapping studies.