Rapid Sequential Spread of Two Wolbachia Variants in Drosophila simulans

The maternally inherited intracellular bacteria Wolbachia can manipulate host reproduction in various ways that foster frequency increases within and among host populations. Manipulations involving cytoplasmic incompatibility (CI), where matings between infected males and uninfected females produce non-viable embryos, are common in arthropods and produce a reproductive advantage for infected females. CI was associated with the spread of Wolbachia variant wRi in Californian populations of Drosophila simulans, which was interpreted as a bistable wave, in which local infection frequencies tend to increase only once the infection becomes sufficiently common to offset imperfect maternal transmission and infection costs. However, maternally inherited Wolbachia are expected to evolve towards mutualism, and they are known to increase host fitness by protecting against infectious microbes or increasing fecundity. We describe the sequential spread over approximately 20 years in natural populations of D. simulans on the east coast of Australia of two Wolbachia variants (wAu and wRi), only one of which causes significant CI, with wRi displacing wAu since 2004. Wolbachia and mtDNA frequency data and analyses suggest that these dynamics, as well as the earlier spread in California, are best understood as Fisherian waves of favourable variants, in which local spread tends to occur from arbitrarily low frequencies. We discuss implications for Wolbachia-host dynamics and coevolution and for applications of Wolbachia to disease control.     

A Flexible Approach for the Analysis of Rare Variants Allowing for a Mixture of Effects on Binary or Quantitative Traits

Multiple rare variants either within or across genes have been hypothesised to collectively influence complex human traits. The increasing availability of high throughput sequencing technologies offers the opportunity to study the effect of rare variants on these traits. However, appropriate and computationally efficient analytical methods are required to account for collections of rare variants that display a combination of protective, deleterious and null effects on the trait. We have developed a novel method for the analysis of rare genetic variation in a gene, region or pathway that, by simply aggregating summary statistics at each variant, can: (i) test for the presence of a mixture of effects on a trait; (ii) be applied to both binary and quantitative traits in population-based and family-based data; (iii) adjust for covariates to allow for non-genetic risk factors and; (iv) incorporate imputed genetic variation. In addition, for preliminary identification of promising genes, the method can be applied to association summary statistics, available from meta-analysis of published data, for example, without the need for individual level genotype data. Through simulation, we show that our method is immune to the presence of bi-directional effects, with no apparent loss in power across a range of different mixtures, and can achieve greater power than existing approaches as long as summary statistics at each variant are robust. We apply our method to investigate association of type-1 diabetes with imputed rare variants within genes in the major histocompatibility complex using genotype data from the Wellcome Trust Case Control Consortium.     

Novel Staphylococcal Glycosyltransferases SdgA and SdgB Mediate Immunogenicity and Protection of Virulence-Associated Cell Wall Proteins

Infection of host tissues by Staphylococcus aureus and S. epidermidis requires an unusual family of staphylococcal adhesive proteins that contain long stretches of serine-aspartate dipeptide-repeats (SDR). The prototype member of this family is clumping factor A (ClfA), a key virulence factor that mediates adhesion to host tissues by binding to extracellular matrix proteins such as fibrinogen. However, the biological siginificance of the SDR-domain and its implication for pathogenesis remain poorly understood. Here, we identified two novel bacterial glycosyltransferases, SdgA and SdgB, which modify all SDR-proteins in these two bacterial species. Genetic and biochemical data demonstrated that these two glycosyltransferases directly bind and covalently link N-acetylglucosamine (GlcNAc) moieties to the SDR-domain in a step-wise manner, with SdgB appending the sugar residues proximal to the target Ser-Asp repeats, followed by additional modification by SdgA. GlcNAc-modification of SDR-proteins by SdgB creates an immunodominant epitope for highly opsonic human antibodies, which represent up to 1% of total human IgG. Deletion of these glycosyltransferases renders SDR-proteins vulnerable to proteolysis by human neutrophil-derived cathepsin G. Thus, SdgA and SdgB glycosylate staphylococcal SDR-proteins, which protects them against host proteolytic activity, and yet generates major eptopes for the human anti-staphylococcal antibody response, which may represent an ongoing competition between host and pathogen.     

Rab18 Binds to Hepatitis C Virus NS5A and Promotes Interaction between Sites of Viral Replication and Lipid Droplets

Hepatitis C virus (HCV) is a single-stranded RNA virus that replicates on endoplasmic reticulum-derived membranes. HCV particle assembly is dependent on the association of core protein with cellular lipid droplets (LDs). However, it remains uncertain whether HCV assembly occurs at the LD membrane itself or at closely associated ER membranes. Furthermore, it is not known how the HCV replication complex and progeny genomes physically associate with the presumed sites of virion assembly at or near LDs. Using an unbiased proteomic strategy, we have found that Rab18 interacts with the HCV nonstructural protein NS5A. Rab18 associates with LDs and is believed to promote physical interaction between LDs and ER membranes. Active (GTP-bound) forms of Rab18 bind more strongly to NS5A than a constitutively GDP-bound mutant. NS5A colocalizes with Rab18-positive LDs in HCV-infected cells, and Rab18 appears to promote the physical association of NS5A and other replicase components with LDs. Modulation of Rab18 affects genome replication and possibly also the production of infectious virions. Our results support a model in which specific interactions between viral and cellular proteins may promote the physical interaction between membranous HCV replication foci and lipid droplets.     

Unlocking the Bottleneck in Forward Genetics Using Whole-Genome Sequencing and Identity by Descent to Isolate Causative Mutations

Forward genetics screens with N-ethyl-N-nitrosourea (ENU) provide a powerful way to illuminate gene function and generate mouse models of human disease; however, the identification of causative mutations remains a limiting step. Current strategies depend on conventional mapping, so the propagation of affected mice requires non-lethal screens; accurate tracking of phenotypes through pedigrees is complex and uncertain; out-crossing can introduce unexpected modifiers; and Sanger sequencing of candidate genes is inefficient. Here we show how these problems can be efficiently overcome using whole-genome sequencing (WGS) to detect the ENU mutations and then identify regions that are identical by descent (IBD) in multiple affected mice. In this strategy, we use a modification of the Lander-Green algorithm to isolate causative recessive and dominant mutations, even at low coverage, on a pure strain background. Analysis of the IBD regions also allows us to calculate the ENU mutation rate (1.54 mutations per Mb) and to model future strategies for genetic screens in mice. The introduction of this approach will accelerate the discovery of causal variants, permit broader and more informative lethal screens to be used, reduce animal costs, and herald a new era for ENU mutagenesis.     

SLC26A4 Targeted to the Endolymphatic Sac Rescues Hearing and Balance in Slc26a4 Mutant Mice

Mutations of SLC26A4 are a common cause of human hearing loss associated with enlargement of the vestibular aqueduct. SLC26A4 encodes pendrin, an anion exchanger expressed in a variety of epithelial cells in the cochlea, the vestibular labyrinth and the endolymphatic sac. Slc26a4 ^Δ/Δ mice are devoid of pendrin and develop a severe enlargement of the membranous labyrinth, fail to acquire hearing and balance, and thereby provide a model for the human phenotype. Here, we generated a transgenic mouse line that expresses human SLC26A4 controlled by the promoter of ATP6V1B1. Crossing this transgene into the Slc26a4 ^Δ/Δ line restored protein expression of pendrin in the endolymphatic sac without inducing detectable expression in the cochlea or the vestibular sensory organs. The transgene prevented abnormal enlargement of the membranous labyrinth, restored a normal endocochlear potential, normal pH gradients between endolymph and perilymph in the cochlea, normal otoconia formation in the vestibular labyrinth and normal sensory functions of hearing and balance. Our study demonstrates that restoration of pendrin to the endolymphatic sac is sufficient to restore normal inner ear function. This finding in conjunction with our previous report that pendrin expression is required for embryonic development but not for the maintenance of hearing opens the prospect that a spatially and temporally limited therapy will restore normal hearing in human patients carrying a variety of mutations of SLC26A4.     

An evaluation of analysis pipelines for DNA methylation profiling using the Illumina HumanMethylation450 BeadChip platform

The proper identification of differentially methylated CpGs is central in most epigenetic studies. The Illumina HumanMethylation450 BeadChip is widely used to quantify DNA methylation; nevertheless, the design of an appropriate analysis pipeline faces severe challenges due to the convolution of biological and technical variability and the presence of a signal bias between Infinium I and II probe design types. Despite recent attempts to investigate how to analyze DNA methylation data with such an array design, it has not been possible to perform a comprehensive comparison between different bioinformatics pipelines due to the lack of appropriate data sets having both large sample size and sufficient number of technical replicates. Here we perform such a comparative analysis, targeting the problems of reducing the technical variability, eliminating the probe design bias and reducing the batch effect by exploiting two unpublished data sets, which included technical replicates and were profiled for DNA methylation either on peripheral blood, monocytes or muscle biopsies. We evaluated the performance of different analysis pipelines and demonstrated that: (1) it is critical to correct for the probe design type, since the amplitude of the measured methylation change depends on the underlying chemistry; (2) the effect of different normalization schemes is mixed, and the most effective method in our hands were quantile normalization and Beta Mixture Quantile dilation (BMIQ); (3) it is beneficial to correct for batch effects. In conclusion, our comparative analysis using a comprehensive data set suggests an efficient pipeline for proper identification of differentially methylated CpGs using the Illumina 450K arrays.     

Quantification and functional analysis of modular protein evolution in a dense phylogenetic tree

Modularity is a hallmark of molecular evolution. Whether considering gene regulation, the components of metabolic pathways or signaling cascades, the ability to reuse autonomous modules in different molecular contexts can expedite evolutionary innovation. Similarly, protein domains are the modules of proteins, and modular domain rearrangements can create diversity with seemingly few operations in turn allowing for swift changes to an organism's functional repertoire. Here, we assess the patterns and functional effects of modular rearrangements at high resolution. Using a well resolved and diverse group of pancrustaceans, we illustrate arrangement diversity within closely related organisms, estimate arrangement turnover frequency and establish, for the first time, branch-specific rate estimates for fusion, fission, domain addition and terminal loss. Our results show that roughly 16 new arrangements arise per million years and that between 64% and 81% of these can be explained by simple, single-step modular rearrangement events. We find evidence that the frequencies of fission and terminal deletion events increase over time, and that modular rearrangements impact all levels of the cellular signaling apparatus and thus may have strong adaptive potential. Novel arrangements that cannot be explained by simple modular rearrangements contain a significant amount of repeat domains that occur in complex patterns which we term “supra-repeats”. Furthermore, these arrangements are significantly longer than those with a single-step rearrangement solution, suggesting that such arrangements may result from multi-step events. In summary, our analysis provides an integrated view and initial quantification of the patterns and functional impact of modular protein evolution in a well resolved phylogenetic tree. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

The Importance of Tryptophan B28 in H2 Relaxin for RXFP2 Binding and Activation

H2 relaxin (relaxin) is a member of the insulin–relaxin superfamily and exhibits several non-reproductive functions in addition to its well-known properties as a pregnancy hormone. Over the years, the therapeutic potential of relaxin has been examined for a number of conditions. It is currently in phase III clinical trials for the treatment of acute heart failure. The 53 amino acid peptide hormone consists of two polypeptide chains (A and B) which are cross-linked by two inter-chains and one intra-A chain disulfide bridge. Although its cognate receptor is relaxin family peptide receptor (RXFP) 1, relaxin is also able to cross-react with RXFP2, for which the native ligand is INSL3. The “RXXXRXXI” motif in the B-chain of H2 relaxin is responsible for primary binding to LRR of the RXFP1 receptor (Büllesbach and Schwabe, J Biol Chem 280:14051–14056, 2005). Previous RXFP2 receptor mutation and molecular modelling studies strongly suggest that, in addition to this motif, the Trp-B28 residue in the B-chain is responsible for H2–RXFP2 interaction. To confirm this finding, here we have mutated H2 relaxin in which Trp-B28 was replaced with alanine. The synthetic relaxin analogue was then tested on cells expressing either RXFP1 or 2 to determine the affinity and potency for the respective receptors. Our results confirm that Trp-B28 in the B-chain is crucial for binding and activating RXFP2, but not for RXFP1.