Trade‐offs and utility of alternative RADseq methods: Reply to Puritz et al.

Puritz et al. provide a review of several RADseq methodological approaches in response to our ‘Population Genomic Data Analysis’ workshop (Sept 2013) review (Andrews & Luikart 2014). We agree with Puritz et al. on the importance for researchers to thoroughly understand RADseq library preparation and data analysis when choosing an approach for answering their research questions. Some of us are currently using multiple RADseq protocols, and we agree that the different methods may offer advantages in different cases. Our workshop review did not intend to provide a thorough review of RADseq because the workshop covered a broad range of topics within the field of population genomics. Similarly, neither the response of Puritz et al. nor our comments here provide sufficient space to thoroughly review RADseq. Nonetheless, here we address some key points that we find unclear or potentially misleading in their evaluation of techniques.     

Identifying rare and common disease associated variants in genomic data using Parkinson’s disease as a model

Background

Genome-wide association studies have been successful in identifying common genetic variants for human diseases. However, much of the heritable variation associated with diseases such as Parkinson’s disease remains unknown suggesting that many more risk loci are yet to be identified. Rare variants have become important in disease association studies for explaining missing heritability. Methods for detecting this type of association require prior knowledge on candidate genes and combining variants within the region. These methods may suffer from power loss in situations with many neutral variants or causal variants with opposite effects.

Results

We propose a method capable of scanning genetic variants to identify the region most likely harbouring disease gene with rare and/or common causal variants. Our method assigns a score at each individual variant based on our scoring system. It uses aggregate scores to identify the region with disease association. We evaluate performance by simulation based on 1000 Genomes sequencing data and compare with three commonly used methods. We use a Parkinson’s disease case–control dataset as a model to demonstrate the application of our method.Our method has better power than CMC and WSS and similar power to SKAT-O with well-controlled type I error under simulation based on 1000 Genomes sequencing data. In real data analysis, we confirm the association of α-synuclein gene (SNCA) with Parkinson’s disease (p = 0.005). We further identify association with hyaluronan synthase 2 (HAS2, p = 0.028) and kringle containing transmembrane protein 1 (KREMEN1, p = 0.006). KREMEN1 is associated with Wnt signalling pathway which has been shown to play an important role for neurodegeneration in Parkinson’s disease.

Conclusions

Our method is time efficient and less sensitive to inclusion of neutral variants and direction effect of causal variants. It can narrow down a genomic region or a chromosome to a disease associated region. Using Parkinson’s disease as a model, our method not only confirms association for a known gene but also identifies two genes previously found by other studies. In spite of many existing methods, we conclude that our method serves as an efficient alternative for exploring genomic data containing both rare and common variants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-014-0088-9) contains supplementary material, which is available to authorized users.     

Quantification of cellular protein expression and molecular features of group 3 LEA proteins from embryos of Artemia franciscana

Late embryogenesis abundant (LEA) proteins are highly hydrophilic, low complexity proteins whose expression has been correlated with desiccation tolerance in anhydrobiotic organisms. Here, we report the identification of three new mitochondrial LEA proteins in anhydrobiotic embryos of Artemia franciscana, AfrLEA3m_47, AfrLEA3m_43, and AfrLEA3m_29. These new isoforms are recognized by antibody raised against recombinant AfrLEA3m, the original mitochondrial-targeted LEA protein previously reported from these embryos; mass spectrometry confirms all four proteins share sequence similarity. The corresponding messenger RNA (mRNA) species for the four proteins are readily amplified from total complementary DNA (cDNA) prepared from embryos. cDNA sequences of the four mRNAs are quite similar, but each has a stretch of sequence that is absent in at least one of the others, plus multiple single base pair differences. We conclude that all four mitochondrial LEA proteins are products of independent genes. Each possesses a mitochondrial targeting sequence, and indeed Western blots performed on extracts of isolated mitochondria clearly detect all four isoforms. Based on mass spectrometry and sodium dodecyl sulfate polyacrylamide gel electrophoresis migration, the cytoplasmic-localized AfrLEA2 exists primarily as a homodimer in A. franciscana. Quantification of protein expression for AfrLEA2, AfrLEA3m, AfrLEA3m_43, and AfrLEA3m_29 as a function of development shows that cellular concentrations are highest in diapause embryos and decrease during development to low levels in desiccation-intolerant nauplius larvae. When adjustment is made for mitochondria matrix volume, the effective concentrations of cytoplasmic versus mitochondrial group 3 LEA proteins are similar in vivo, and the values provide guidance for the design of in vitro functional studies with these proteins.     

Development and implementation of a multiplexed single nucleotide polymorphism genotyping tool for differentiation of ryegrass species and cultivars

Perennial ryegrass (Lolium perenne L.) and Italian ryegrass (Lolium multiflorum Lam.) are important temperate forage grasses which are closely related, generating fertile interspecific hybrids. All groups are represented by multiple cultivars in the commercial pasture seeds market. Due to the close taxonomic relationship between the two species, differentiation based on morphophysiological criteria is not always readily achievable. In addition, an obligate outbreeding reproductive habit produces high levels of individual heterozygosity and intrapopulation diversity, which presents problems for discrimination between cultivars. Molecular genetic marker polymorphism provides an effective means of addressing these challenges. An iterative process of resequencing from loci distributed across the perennial ryegrass genome was used to identify single nucleotide polymorphism (SNP) markers, which were then validated and formatted in a highly multiplexed (384-plex) assay system. SNP genotyping was then performed across samples of 48–192 individuals from a total of 27 ryegrass cultivars (19 of perennial ryegrass, seven of Italian ryegrass and one hybrid cultivar). SNP markers from perennial ryegrass exhibited a high level of transfer to Italian ryegrass. Data analysis permitted quantification of intra- and inter-species diversity, as well as discrimination between cultivars within each species, including diploid and autotetraploid cultivars of perennial ryegrass. Lower levels of SNP-based diversity were detected in Italian ryegrass than in perennial ryegrass. A neighbour-joining tree based on genetic distance analysis located a hybrid cultivar to an intermediate position between the two species-specific cultivar groups. The resulting catalogue of ryegrass cultivars will provide support for the processes of cultivar accreditation and quality assurance.     

Impact of temperature on Na2Sr(PO4)F:Eu3+ phosphor

In this article, we report the synthesis of Na₂Sr_1‐x(PO₄)F:Eu_x phosphor via a combustion method. The influence of different annealing temperatures on the photoluminescence properties was investigated. The phosphor was excited at both 254 and 393 nm. Na₂Sr_1‐x(PO₄)F:Eu_x³⁺ phosphors emit strong orange and red color at 593 and 612 nm, respectively, under both excitation wavelengths. Na₂Sr_1‐x(PO₄)F:Eu_x³⁺ phosphors annealed at 1050°C showed stronger emission intensity compared with 600, 900 and 1200°C. Moreover, Na₂Sr_1‐x(PO₄)F:Eu_x³⁺ phosphor was found to be more intense when compared with commercial Y₂O₃:Eu³⁺ phosphor. Copyright © 2013 John Wiley & Sons, Ltd.     

Conformational analyses of mycothiol, a critical intracellular glycothiol in Mycobacteria

Intracellular thiols are essential biomolecules, which play several critical roles in living organisms including controlling intracellular redox potential and acting as cofactors for several vital detoxification enzymes including S-transferases and formaldehyde dehydrogenases. The tripeptide gamma-L-glutamyl-L-cysteinylglycine, more commonly known as glutathione, is well known as the major intracellular thiol in eukaryotes and in some bacteria. However, glutathione is absent in the Actinomycetales bacteria such as Mycobacteria and Streptomyces and is believed to be replaced by 1-D-myo-inosityl-2-(N-acetyl-L-cysteinyl)amido-2-deoxy-alpha-D-glucopyranoside, mycothiol, in these organisms. Although much is known about the chemistry and biochemistry of glutathione, currently much less is known concerning mycothiol and its properties. The structure of mycothiol is composed of a glycoside linkage between myo-inositol and D-glucosamine with an N-acetyl-L-cysteine linked to the 2'-amino group of the d-glucosamine moiety. Mycothiol is currently of intense interest due to its essential role in the cellular physiology of Mycobacteria, such as Mycobacterium tuberculosis, and its possible role in antimycobacterial drug resistance. A detailed investigation of its chemistry is therefore essential in ameliorating our knowledge of this key glycothiol, and in shedding additional light on its biochemical role in these pathogenic organisms. This report presents a detailed conformational analysis of mycothiol utilizing a variety of force fields and stochastic search protocols. Cluster analyses of energetically low lying conformations have indicated the presence of several key conformations that are populated in the gas phase and with implicit water solvation. These conformations are compared to recent NMR studies on a derivative of mycothiol. This information should be an important contribution to our basic understanding of the chemistry of this glycothiol and critical in the design of novel inhibitors of pathogen enzymes that require it.     

Mutation of a conserved threonine in the third transmembrane helix of alpha- and beta-connexins creates a dominant-negative closed gap junction channel

Single site mutations in connexins have provided insights about the influence specific amino acids have on gap junction synthesis, assembly, trafficking, and functionality. We have discovered a single point mutation that eliminates functionality without interfering with gap junction formation. The mutation occurs at a threonine residue located near the cytoplasmic end of the third transmembrane helix. This threonine is strictly conserved among members of the alpha- and beta-connexin subgroups but not the gamma-subgroup. In HeLa cells, connexin43 and connexin26 mutants are synthesized, traffic to the plasma membrane, and make gap junctions with the same overall appearance as wild type. We have isolated connexin26T135A gap junctions both from HeLa cells and baculovirus-infected insect Sf9 cells. By using cryoelectron microscopy and correlation averaging, difference images revealed a small but significant size change within the pore region and a slight rearrangement of the subunits between mutant and wild-type connexons expressed in Sf9 cells. Purified, detergent-solubilized mutant connexons contain both hexameric and partially disassembled structures, although wild-type connexons are almost all hexameric, suggesting that the three-dimensional mutant connexon is unstable. Mammalian cells expressing gap junction plaques composed of either connexin43T154A or connexin26T135A showed an absence of dye coupling. When expressed in Xenopus oocytes, these mutants, as well as a cysteine substitution mutant of connexin50 (connexin50T157C), failed to produce electrical coupling in homotypic and heteromeric pairings with wild type in a dominant-negative effect. This mutant may be useful as a tool for knocking down or knocking out connexin function in vitro or in vivo.     

Erratum

Neurotransmitter receptor trafficking and the regulation of synaptic strength. Traffic 2001:2(7):437–448.