Mapping a quantitative trait locus via the EM algorithm and Bayesian classification.
Mapping a locus controlling a quantitative genetic trait (e.g., blood pressure) to a specific genomic region is of considerable interest. Data on the quantitative trait under consideration and several codominant genetic markers with known genomic locations are collected from members of families and statistically analyzed to draw inferences on the genomic position of the trait locus. The vector of parameters of interest comprises the pairwise recombination fractions, theta, between the putative quantitative trait locus and the marker loci. One of the major complications in estimating theta for a quantitative trait in humans is the lack of haplotype information on members of families. The purpose of this study was to devise a computationally simple and efficient method of estimation of theta in the absence of haplotype information. We have proposed a two-stage estimation procedure using the expectation-maximization (EM) algorithm. In the first stage, parameters of the QTL are estimated based on data of a sample of unrelated individuals. From estimates thus obtained, we have used a Bayes' rule to infer QTL genotypes of parents in families. Finally, in the second stage of the procedure, we have proposed an EM algorithm for obtaining the maximum likelihood estimate of theta based on data of informative families (which are identified upon inferring parental QTL genotypes performed in the first stage). We have shown, using simulated data, that the proposed procedure is cost-effective, computationally simple, and statistically efficient. As expected, analysis of data on multiple markers jointly is more efficient than the analysis based on single markers.

title: Identifying rare variants from exome scans: the GAW17 experience
abstract: Genetic Analysis Workshop 17 (GAW17) provided a platform for evaluating existing statistical genetic methods and for developing novel methods to analyze rare variants that modulate complex traits. In this article, we present an overview of the 1000 Genomes Project exome data and simulated phenotype data that were distributed to GAW17 participants for analyses, the different issues addressed by the participants, and the process of preparation of manuscripts resulting from the discussions during the workshop.

title: A novel non-parametric regression reveals linkage on chromosome 4 for the number of externalizing symptoms in sib-pairs.
abstract: In this report, we present results of a genome-wide linkage scan using as a phenotype the number of externalizing symptoms associated with alcohol use disorders. Subjects were collected by the Collaborative Study on the Genetics of Alcoholism project from families in which at least three first degree relatives were affected by alcohol dependence. We use a novel non-parametric regression method based on kernel smoothing for our analysis. We report a statistically significant linkage close to the ADH gene cluster on Chromosome 4. We also obtain evidence for epistatic interaction between a region on Chromosome 1 and one on Chromosome 15. Although alcoholism as a covariate does not have any effect on the linkage scan, it has an effect on the epistatic interaction.

title: A quantile-based method for association mapping of quantitative phenotypes: an application to rheumatoid arthritis phenotypes
abstract: Genetic association of population-based quantitative trait data has traditionally been analyzed using analysis of variance (ANOVA). However, violations of certain statistical assumptions may lead to false-positive association results. In this study, we have explored model-free alternatives to ANOVA using correlations between allele frequencies in the different quantile intervals of the quantitative trait and the quantile values. We performed genome-wide association scans on anti-cyclic citrullinated peptide and rheumatoid factor-immunoglobulin M, two quantitative traits correlated with rheumatoid arthritis, using the data provided in Genetic Analysis Workshop 16. Both the quantitative traits exhibited significant evidence of association on Chromosome 6, although not in the human leukocyte antigen region which is known to harbor a major gene predisposing to rheumatoid arthritis. We found that while a majority of the significant findings using the asymptotic thresholds of ANOVA was not validated using permutations, a relatively higher proportion of the significant findings using the asymptotic cut-offs of the correlation statistic were validated using permutations.

title: Integrating sibship data for mapping quantitative trait loci.
abstract: Sibship methods have been shown to be more powerful than traditional sib-pair methods in mapping quantative trait loci. We propose a statistical procedure which integrates data on sibships into a so-called 'contrast function', a natural extension of the classical squared sib-pair trait difference proposed by Haseman & Elston (1972). We also develop a combined mean and contrast function which provides more information on linkage compared to the contrast function. Our method is extended to multiple, epistatically interacting trait loci. Monte-Carlo simulations are included to compare the efficiencies of the proposed procedures with some currently used methods. An application of our proposed method is presented using data on alcohol dependence.

title: Linkage mapping of a complex trait in the New York population of the GAW14 simulated dataset: a multivariate phenotype approach
abstract: Multivariate phenotypes underlie complex traits. Thus, instead of using the end-point trait, it may be statistically more powerful to use a multivariate phenotype correlated to the end-point trait for detecting linkage. In this study, we develop a reverse regression method to analyze linkage of Kofendrerd Personality Disorder affection status in the New York population of the Genetic Analysis Workshop 14 (GAW14) simulated dataset. When we used the multivariate phenotype, we obtained significant evidence of linkage near four of the six putative loci in at least 25% of the replicates. On the other hand, the linkage analysis based on Kofendrerd Personality Disorder status as a phenotype produced significant findings only near two of the loci and in a smaller proportion of replicates.

title: Full genomic analysis of a simian SA11-like G3P[2] rotavirus strain isolated from an asymptomatic infant: identification of novel VP1, VP6 and NSP4 genotypes.
abstract: We report here the full genomic analysis of a simian SA11-like G3P[2] group A rotavirus (GAR) strain, B10, isolated from an asymptomatic infant in Kenya in 1987. By nucleotide sequence identities and phylogenetic analyses, the VP7-VP4-VP2-VP3-NSP1-NSP2-NSP3-NSP5 genes of strain B10 exhibited maximum genetic relatedness to those of the different isolates of simian strain SA11, and were assigned to the G3-P[2]-C5-M5-A5-N5-T5-H5 genotypes, respectively. On the other hand, the VP1, VP6 and NSP4 genes of strain B10 did not belong to any of the established GAR genotypes, and therefore, were assigned to new genotype numbers R8, I16 and E13, respectively, by the Rotavirus Classification Working Group. These observations suggested that strain B10 might have originated from reassortment event/s involving simian SA11-like strains and GAR strains from unknown animal host species (possibly other wild animals) preceding transmission to humans. Alternatively, considering the lack of data on simian GARs, it might be also possible that the VP1, VP6 and NSP4 genes of strain B10 are those of unknown simian strains, and that strain B10 might be a typical simian strain that was directly transmitted to humans. Therefore, either hypothesis pointed towards a rare instance of possible direct transmission of GARs from an animal host (possibly a monkey or some other wild animal) to humans. This was corroborated by the presence of different species of wild animals including non-human primates, and unhygienic conditions at the sampling site. To our knowledge, the present study is the first report on the detection of a simian SA11-like G3P[2] GAR strain in humans.

title: Application of rbcL based molecular diversity analysis to algae in wastewater treatment plants.
abstract: The molecular diversity of algae in the final clarifier or denitrification filter outfall from three wastewater treatment plants (WWTPs) with activated sludge based treatment was analyzed using the rbcL gene as a phylogenetic marker. The rbcL gene encodes the large subunit of the CO(2) fixing enzyme RuBisCO. Among algae identified at the WWTPs were diatoms, green algae, cyanobacteria, Eustigmatophyceae, and unknown heterokonts. A high level of diversity was observed within WWTPs with 19-24 unique rbcL sequences detected at each plant. Algae composition also varied between treatment plants. Our results show that the rbcL gene can be used as a phylogenetic marker for algae diversity analysis in a wastewater treatment context.

title: Full genomic analyses of two human G2P[4] rotavirus strains detected in 2005: identification of a caprine-like VP3 gene.
abstract: Although G2P[4] rotaviruses are common causes of infantile diarrhoea, to date only the full genomes of the prototype (strain DS-1) and another old strain, TB-Chen, have been analysed. We report here the full genomic analyses of two Bangladeshi G2P[4] strains, MMC6 and MMC88, detected in 2005. Both the strains exhibited a DS-1-like genotype constellation. Excluding the VP4 and VP7 genes, and except for VP3 of MMC88, the MMC strains were genetically more closely related to the contemporary G2P[4] and several non-G2P[4] human strains than the prototype G2P[4] strain. However, by phylogenetic analyses, the VP2, VP3 (except MMC88), NSP1 and NSP3-5 genes of these strains appeared to share a common origin with those of the prototype strain, whilst their VP1, VP6 and NSP2 genes clustered near a caprine strain. The VP3 gene of MMC88 exhibited maximum relatedness to a local caprine strain, representing the first reported human G2P[4] strain with a gene of animal origin.

title: Ultrafast fluorescence resonance energy transfer in the micelle and the gel phase of a PEO-PPO-PEO triblock copolymer: excitation wavelength dependence.
abstract: Fluorescence resonance energy transfer (FRET) from coumarin 480 (C480) to rhodamine 6G (R6G) is studied in the micelle and the gel phase of a triblock copolymer, (PEO)20-(PPO)70-(PEO)20 (Pluronic P123 (P123)) by picosecond and femtosecond emission spectroscopy. The time constants of FRET were obtained from the rise time of the acceptor (R6G) emission. In a P123 micelle, FRET occurs in multiple time scales: 2.5, 100, and 1700 ps. In the gel phase, three rise components are observed: 3, 150, and 2600 ps. According to a simple Förster model, the ultrafast (2.5 and 3 ps) components of FRET correspond to donor-acceptor distance RDA=13 +/- 2 A. The ultrafast FRET occurs between a donor and an acceptor residing at close contact at the corona (PEO) region of a P123 micelle. With increase in the excitation wavelength (lambdaex) from 375 to 435 nm, the relative contribution of the ultrafast component of FRET ( approximately 3 ps) increases from 13% to 100% in P123 micelle and from 1% to 100% in P123 gel. It is suggested that at lambdaex = 435 nm, mainly the highly polar peripheral region is probed where FRET is very fast due to close proximity of the donor and the acceptor. The 100 and 150 ps components correspond to RDA = 25 +/- 2 A and are ascribed to FRET from C480 deep inside the micelle to an acceptor (R6G) in the peripheral region. The very long component of FRET (1700 ps in micelle and 2600 ps component in gel) may arise from diffusion of the donor from outside the micelle to the interior followed by fast FRET.

title: Idiopathic intracranial hypertension in a child after hemodialysis.
abstract: Idiopathic intracranial hypertension remains a disease of unknown etiology. Epidemiology, associations, features, and prognosis in children have changed over the years. The main clinical complaint of pediatric idiopathic intracranial hypertension is headache with visual disturbance. We report on a patient with childhood idiopathic intracranial hypertension, who presented without headache after hemodialysis for acute renal failure. Idiopathic intracranial hypertension without headache is believed to have a poor prognosis. However, this child demonstrated a favorable outcome with medical therapy. The child was followed with serial optical coherence tomography of retinal nerve fiber layer thickness around the optic nerve head. In the absence of a reliable procedure to evaluate the therapeutic efficacy in the follow-up of patients with idiopathic intracranial hypertension, the importance of optical coherence tomography, a noninvasive, objective, and reproducible procedure, is highlighted as a useful adjunct in the management of idiopathic intracranial hypertension.