Comprehensive approach to analyzing rare genetic variants

Recent findings suggest that rare variants play an important role in both monogenic and common diseases. Due to their rarity, however, it remains unclear how to appropriately analyze the association between such variants and disease. A common approach entails combining rare variants together based on a priori information and analyzing them as a single group. Here one must make some assumptions about what to aggregate. Instead, we propose two approaches to empirically determine the most efficient grouping of rare variants. The first considers multiple possible groupings using existing information. The second is an agnostic "step-up" approach that determines an optimal grouping of rare variants analytically and does not rely on prior information. To evaluate these approaches, we undertook a simulation study using sequence data from genes in the one-carbon folate metabolic pathway. Our results show that using prior information to group rare variants is advantageous only when information is quite accurate, but the step-up approach works well across a broad range of plausible scenarios. This agnostic approach allows one to efficiently analyze the association between rare variants and disease while avoiding assumptions required by other approaches for grouping such variants.     

Characterizing the Aedes aegypti Population in a Vietnamese Village in Preparation for a Wolbachia-Based Mosquito Control Strategy to Eliminate Dengue

Background

A life-shortening strain of the obligate intracellular bacteria Wolbachia, called wMelPop, is seen as a promising new tool for the control of Aedes aegypti. However, developing a vector control strategy based on the release of mosquitoes transinfected with wMelPop requires detailed knowledge of the demographics of the target population.

Methodology/Principal Findings

In Tri Nguyen village (611 households) on Hon Mieu Island in central Vietnam, we conducted nine quantitative entomologic surveys over 14 months to determine if Ae. aegypti populations were spatially and temporally homogenous, and to estimate population size. There was no obvious relationship between mosquito (larval, pupal or adult) abundance and temperature and rainfall, and no area of the village supported consistently high numbers of mosquitoes. In almost all surveys, key premises produced high numbers of Ae. aegypti. However, these premises were not consistent between surveys. For an intervention based on a single release of wMelPop-infected Ae. aegypti, release ratios of infected to uninfected adult mosquitoes of all age classes are estimated to be 1.8–6.7∶1 for gravid females (and similarly aged males) or teneral adults, respectively. We calculated that adult female mosquito abundance in Tri Nguyen village could range from 1.1 to 43.3 individuals of all age classes per house. Thus, an intervention could require the release of 2–78 wMelPop-infected gravid females and similarly aged males per house, or 7–290 infected teneral female and male mosquitoes per house.

Conclusions/Significance

Given the variability we encountered, this study highlights the importance of multiple entomologic surveys when evaluating the spatial structure of a vector population or estimating population size. If a single release of wMelPop-infected Ae. aegypti were to occur when wild Ae. aegypti abundance was at its maximum, a preintervention control program would be necessary to ensure that there was no net increase in mosquito numbers. However, because of the short-term temporal heterogeneity, the inconsistent spatial structure and the impact of transient key premises that we observed, the feasibility of multiple releases of smaller numbers of mosquitoes also needs to be considered. In either case, fewer wMelPop-infected mosquitoes would then need to be released, which will likely be more acceptable to householders.     

West Africa - A Safe Haven for Frogs? A Sub-Continental Assessment of the Chytrid Fungus (Batrachochytrium dendrobatidis)

A putative driver of global amphibian decline is the panzootic chytrid fungus Batrachochytrium dendrobatidis (Bd). While Bd has been documented across continental Africa, its distribution in West Africa remains ambiguous. We tested 793 West African amphibians (one caecilian and 61 anuran species) for the presence of Bd. The samples originated from seven West African countries - Bénin, Burkina Faso, Côte d''Ivoire, Ghana, Guinea, Liberia, Sierra Leone - and were collected from a variety of habitats, ranging from lowland rainforests to montane forests, montane grasslands to humid and dry lowland savannahs. The species investigated comprised various life-history strategies, but we focused particularly on aquatic and riparian species. We used diagnostic PCR to screen 656 specimen swabs and histology to analyse 137 specimen toe tips. All samples tested negative for Bd, including a widespread habitat generalist Hoplobatrachus occipitalis which is intensively traded on the West African food market and thus could be a potential dispersal agent for Bd. Continental fine-grained (30 arc seconds) environmental niche models suggest that Bd should have a broad distribution across West Africa that includes most of the regions and habitats that we surveyed. The surprising apparent absence of Bd in West Africa indicates that the Dahomey Gap may have acted as a natural barrier. Herein we highlight the importance of this Bd-free region of the African continent - especially for the long-term conservation of several threatened species depending on fast flowing forest streams (Conraua alleni (“Vulnerable”) and Petropedetes natator (“Near Threatened”)) as well as the “Critically Endangered” viviparous toad endemic to the montane grasslands of Mount Nimba (Nimbaphrynoides occidentalis).     

Identification of Pyrus Single Nucleotide Polymorphisms (SNPs) and Evaluation for Genetic Mapping in European Pear and Interspecific Pyrus Hybrids

We have used new generation sequencing (NGS) technologies to identify single nucleotide polymorphism (SNP) markers from three European pear (Pyrus communis L.) cultivars and subsequently developed a subset of 1096 pear SNPs into high throughput markers by combining them with the set of 7692 apple SNPs on the IRSC apple Infinium® II 8K array. We then evaluated this apple and pear Infinium® II 9K SNP array for large-scale genotyping in pear across several species, using both pear and apple SNPs. The segregating populations employed for array validation included a segregating population of European pear (‘Old Home’בLouise Bon Jersey’) and four interspecific breeding families derived from Asian (P. pyrifolia Nakai and P. bretschneideri Rehd.) and European pear pedigrees. In total, we mapped 857 polymorphic pear markers to construct the first SNP-based genetic maps for pear, comprising 78% of the total pear SNPs included in the array. In addition, 1031 SNP markers derived from apple (13% of the total apple SNPs included in the array) were polymorphic and were mapped in one or more of the pear populations. These results are the first to demonstrate SNP transferability across the genera Malus and Pyrus. Our construction of high density SNP-based and gene-based genetic maps in pear represents an important step towards the identification of chromosomal regions associated with a range of horticultural characters, such as pest and disease resistance, orchard yield and fruit quality.     

Generation of an immortal differentiated lung type-II epithelial cell line from the adult H-2KbtsA58 transgenic mouse

Summary This paper describes a new fully differentiated Type-II alveolar epithelial cell line designated T₇, derived from transgenic H-2K^b-tsA58 mice, capable of being passaged as an immortalized cloned cell line in culture. H-2K^b-tsA58 mice harbor a temperature-sensitive (ts) mutant of the simian virus 40 (SV40) large tumor antigen (T antigen) under the control of the γ-interferon (INF)-inducible mouse major histocompatibility complex H-2K^b promoter. When cultured under permissive conditions (33°C and in the presence of γ-INF) cells isolated from H-2K^b-tsA58 mice express the large T antigen, which drives the cells to proliferate. However, upon withdrawal of the γ-INF and transfer of the cells to a higher temperature (39°C), T antigen expression is turned off, the cells stop proliferating and differentiate. The T₇ cell line is a clonal cell line originally derived from a Type-II cell-rich fraction isolated from lungs of H-2K^b-tsA58 mice. The T₇ cells form confluent monolayers, and have a polarized epithelial cell morphology with tight junctions and apical microvilli. In addition, the T₇ cells have distinct cytoplasmic lamellar bodies, which become more numerous and pronounced when the cells are grown under nonpermissive conditions. The T₇ cells synthesize and secrete phosphatidylcholine and the three surfactant proteins, SP-A, SP-B, and SP-C. The T₇ cell line is unique in that it is the first non-tumor-derived Type-II cell line capable of synthesizing and secreting the major components of surfactant. Based on the criteria studied, the T₇ cell line is phenotypically very similar to normal Type-II cells. The T₇ cell line, therefore, should prove a valuable experimental system to advance the study of the cell biology/physiology of surfactant metabolism and secretion as well as serve as a model for other studies of Type-II cell physiology.     

FlowMax: A Computational Tool for Maximum Likelihood Deconvolution of CFSE Time Courses

The immune response is a concerted dynamic multi-cellular process. Upon infection, the dynamics of lymphocyte populations are an aggregate of molecular processes that determine the activation, division, and longevity of individual cells. The timing of these single-cell processes is remarkably widely distributed with some cells undergoing their third division while others undergo their first. High cell-to-cell variability and technical noise pose challenges for interpreting popular dye-dilution experiments objectively. It remains an unresolved challenge to avoid under- or over-interpretation of such data when phenotyping gene-targeted mouse models or patient samples. Here we develop and characterize a computational methodology to parameterize a cell population model in the context of noisy dye-dilution data. To enable objective interpretation of model fits, our method estimates fit sensitivity and redundancy by stochastically sampling the solution landscape, calculating parameter sensitivities, and clustering to determine the maximum-likelihood solution ranges. Our methodology accounts for both technical and biological variability by using a cell fluorescence model as an adaptor during population model fitting, resulting in improved fit accuracy without the need for ad hoc objective functions. We have incorporated our methodology into an integrated phenotyping tool, FlowMax, and used it to analyze B cells from two NFκB knockout mice with distinct phenotypes; we not only confirm previously published findings at a fraction of the expended effort and cost, but reveal a novel phenotype of nfkb1/p105/50 in limiting the proliferative capacity of B cells following B-cell receptor stimulation. In addition to complementing experimental work, FlowMax is suitable for high throughput analysis of dye dilution studies within clinical and pharmacological screens with objective and quantitative conclusions.     

Oligomeric structure of ExbB and ExbB-ExbD isolated from Escherichia coli as revealed by LILBID mass spectrometry

Energy-coupled transporters in the outer membrane of Escherichia coli and other Gram-negative bacteria allow the entry of scarce substrates, toxic proteins, and bacterial viruses (phages) into the cells. The required energy is derived from the proton-motive force of the cytoplasmic membrane, which is coupled to the outer membrane via the ExbB-ExbD-TonB protein complex. Knowledge of the structure of this complex is required to elucidate the mechanisms of energy harvesting in the cytoplasmic membrane and energy transfer to the outer membrane transporters. Here we solubilized an ExbB oligomer and an ExbB-ExbD subcomplex from the cytoplasmic membrane with the detergent undecyl maltoside. Using laser-induced liquid bead ion desorption mass spectrometry (LILBID-MS), we determined at moderate desorption laser energies the oligomeric structure of ExbB to be mainly hexameric (ExbB(6)), with minor amounts of trimeric (ExbB(3)), dimeric (ExbB(2)), and monomeric (ExbB(1)) oligomers. Under the same conditions ExbB-ExbD formed a subcomplex consisting of ExbB(6)ExbD(1), with a minor amount of ExbB(5)ExbD(1). At higher desorption laser intensities, ExbB(1) and ExbD(1) and traces of ExbB(3)ExbD(1), ExbB(2)ExbD(1), ExbB(1)ExbD(1), ExbB(3), and ExbB(2) were observed. Since the ExbB(6) complex and the ExbB(6)ExbD(1) complex remained stable during solubilization and subsequent chromatographic purification on nickel-nitrilotriacetate agarose, Strep-Tactin, and Superdex 200, and during native blue gel electrophoresis, we concluded that ExbB(6) and ExbB(6)ExbD(1) are subcomplexes on which the final complex including TonB is assembled.     

A Novel Mitochondrial Genome Organization for the Blue Mussel, Mytilus Edulis

The sequence of 13.9 kilobases (kb) of the 17.1-kb mitochondrial genome of Mytilus edulis has been determined, and the arrangement of all genes has been deduced. Mytilus mitochondrial DNA (mtDNA) contains 37 genes, all of which are transcribed from the same DNA strand. The gene content of Mytilus is typically metazoan in that it includes genes for large and small ribosomal RNAs, for a complete set of transfer RNAs and for 12 proteins. The protein genes encode the cytochrome b apoenzyme, cytochrome c oxidase (CO) subunits I-III, NADH dehydrogenase (ND) subunits 1-6 and 4L, and ATP synthetase (ATPase) subunit 6. No gene for ATPase subunit 8 could be found. The reading frames for the ND1, COI, and COIII genes contain long extensions relative to those genes in other metazoan mtDNAs. There are 23 tRNA genes, one more than previously found in any metazoan mtDNA. The additional tRNA appears to specify methionine, making Mytilus mtDNA unique in having two tRNA(Met) genes. Five lengthy unassigned intergenic sequences are present, four of which vary in length from 79 to 119 nucleotides and the largest of which is 1.2 kb. The base compositions of these are unremarkable and do not differ significantly from that of the remainder of the mtDNA. The arrangement of genes in Mytilus mtDNA is remarkably unlike that found in any other known metazoan mtDNA.     

Loss of Parkin or PINK1 Function Increases Drp1-dependent Mitochondrial Fragmentation

Loss-of-function mutations in the parkin gene (PARK2) and PINK1 gene (PARK6) are associated with autosomal recessive parkinsonism. PINK1 deficiency was recently linked to mitochondrial pathology in human cells and Drosophila melanogaster, which can be rescued by parkin, suggesting that both genes play a role in maintaining mitochondrial integrity. Here we demonstrate that an acute down-regulation of parkin in human SH-SY5Y cells severely affects mitochondrial morphology and function, a phenotype comparable with that induced by PINK1 deficiency. Alterations in both mitochondrial morphology and ATP production caused by either parkin or PINK1 loss of function could be rescued by the mitochondrial fusion proteins Mfn2 and OPA1 or by a dominant negative mutant of the fission protein Drp1. Both parkin and PINK1 were able to suppress mitochondrial fragmentation induced by Drp1. Moreover, in Drp1-deficient cells the parkin/PINK1 knockdown phenotype did not occur, indicating that mitochondrial alterations observed in parkin- or PINK1-deficient cells are associated with an increase in mitochondrial fission. Notably, mitochondrial fragmentation is an early phenomenon upon PINK1/parkin silencing that also occurs in primary mouse neurons and Drosophila S2 cells. We propose that the discrepant findings in adult flies can be explained by the time of phenotype analysis and suggest that in mammals different strategies may have evolved to cope with dysfunctional mitochondria.Many lines of evidence suggest that mitochondrial dysfunction plays a central role in the pathogenesis of Parkinson disease, starting from the early observation that the complex I inhibitor 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine induced acute and irreversible parkinsonism in young drug addicts (for review, see Refs. 1–3). In support of a crucial role of mitochondria in Parkinson disease, several Parkinson disease-associated gene products directly or indirectly impinge on mitochondrial integrity (for review, see Refs. 4–6). A clear link between Parkinson disease genes and mitochondria has recently emerged from studies on PINK1 (PTEN-induced putative kinase 1), a mitochondrial serine/threonine kinase, and parkin, a cytosolic E3 ubiquitin ligase. Drosophila parkin null mutants displayed reduced life span, male sterility, and locomotor defects due to apoptotic flight muscle degeneration (7). The earliest manifestation of muscle degeneration and defective spermatogenesis was mitochondrial pathology, exemplified by swollen mitochondria and disintegrated cristae. Remarkably, Drosophila PINK1 null mutants shared marked phenotypic similarities with parkin mutants, and parkin could compensate for the PINK1 loss-of-function phenotype but not vice versa, leading to the conclusion that PINK1 and parkin function in a common genetic pathway with parkin acting downstream of PINK1 (8–10). We recently demonstrated that PINK1 deficiency in cultured human cells causes alterations in mitochondrial morphology, which can be rescued by wild type parkin but not by pathogenic parkin mutants (11). We now present evidence that parkin plays an essential role in maintaining mitochondrial integrity. RNAi³-mediated knockdown of parkin increases mitochondrial fragmentation and decreases cellular ATP production. Notably, mitochondrial fragmentation induced by PINK1/parkin deficiency is observed not only in human neuroblastoma cells but also in primary mouse neurons and insect S2 cells. Alterations in mitochondrial morphology are early manifestations of parkin/PINK1 silencing that are not caused by an increase in apoptosis. The mitochondrial phenotype observed in parkin- or PINK1-deficient cells can morphologically and functionally be rescued by the increased expression of a dominant negative mutant of the fission-promoting protein Drp1. Moreover, manifestation of the PINK1/parkin knockdown phenotype is dependent on Drp1 expression, indicating that an acute loss of parkin or PINK1 function increases mitochondrial fission.     

Tropical rainforest flies carrying pathogens form stable associations with social nonhuman primates

Living in groups provides benefits but also incurs costs such as attracting disease vectors. For example, synanthropic flies associate with human settlements, and higher fly densities increase pathogen transmission. We investigated whether such associations also exist in highly mobile nonhuman primate (NHP) Groups. We studied flies in a group of wild sooty mangabeys (Cercocebus atys atys) and three communities of wild chimpanzees (Pan troglodytes verus) in Taï National Park, Côte d'Ivoire. We observed markedly higher fly densities within both mangabey and chimpanzee groups. Using a mark–recapture experiment, we showed that flies stayed with the sooty mangabey group for up to 12 days and for up to 1.3 km. We also tested mangabey‐associated flies for pathogens infecting mangabeys in this ecosystem, Bacillus cereus biovar anthracis (Bcbva), causing sylvatic anthrax, and Treponema pallidum pertenue, causing yaws. Flies contained treponemal (6/103) and Bcbva (7/103) DNA. We cultured Bcbva from all PCR‐positive flies, confirming bacterial viability and suggesting that this bacterium might be transmitted and disseminated by flies. Whole genome sequences of Bcbva isolates revealed a diversity of Bcbva, probably derived from several sources. We conclude that flies actively track mangabeys and carry infectious bacterial pathogens; these associations represent an understudied cost of sociality and potentially expose many social animals to a diversity of pathogens.