Author Archives: Bora Zivkovic

Continue reading →

Continue reading →

So, let’s see what’s new in PLoS Genetics, PLoS Computational Biology and PLoS Pathogens this week. As always, you should rate the articles, post notes and comments and send trackbacks when you blog about the papers. Here are my own picks for the week – you go and look for your own favourites:
Social Contact Networks and Disease Eradicability under Voluntary Vaccination:

Interest in infectious disease models that incorporate the effects of human behavior has been growing in recent years. However, most of these models predict that it should never be possible to eradicate a disease under voluntary vaccination, due to nonvaccinating “free riders” that emerge when vaccine coverage is high. This prediction contradicts the fact that smallpox was eradicated under a voluntary vaccination policy in many jurisdictions, and that other diseases such as polio are likewise near eradication. These previous models assumed that populations mix homogeneously. However, for some diseases, such as HIV and smallpox, individuals are more likely to get the disease from certain social contacts. Here we show that using a network model that captures this social structure can reconcile the previous theories to the empirical fact that diseases can be eradicated under voluntary vaccination. When infection is transmitted only through close contacts in the network, then an outbreak can be quickly contained using only voluntary vaccination. However, when infection can potentially be transmitted to almost everyone in the network (such as for measles), a disease outbreak can never be contained using voluntary vaccination. This latter observation may have some relevance to the Measles-Mumps-Rubella autism “vaccine scare.”

Genome-Wide Association Studies in an Isolated Founder Population from the Pacific Island of Kosrae:

Isolated populations have contributed to the discovery of loci with simple Mendelian segregation and large effects on disease risk or trait variation. We hypothesized that the use of isolated populations might also facilitate the discovery of common alleles contributing to complex traits with relatively larger effects. However, the use of association analyses to map common loci influencing trait variation in large, inbred cohorts introduces analytic challenges, as extensive relatedness between subjects violates the assumptions of independence upon which traditional association test statistics are based. We developed an analytic strategy to perform genome-wide association studies in an inbred family containing over 2,800 individuals from the island of Kosrae, Federated States of Micronesia. No alleles with large effect were observed with strong statistical support in any of the 15 traits examined, suggesting that the contribution of individual common variants to complex trait variation in Kosraens is typically not much greater than that observed in other populations. We show that the effects of many loci previously identified in Caucasian populations are indistinguishable in Caucasians and Kosraens, despite very different population genetics and environmental influences.

A Gene-Based Linkage Map for Bicyclus anynana Butterflies Allows for a Comprehensive Analysis of Synteny with the Lepidopteran Reference Genome:

Lepidopterans (butterflies and moths) are a rich and diverse order of insects, which, despite their economic impact and unusual biological properties, are relatively underrepresented in terms of genomic resources. The genome of the silkworm Bombyx mori has been fully sequenced, but comparative lepidopteran genomics has been hampered by the scarcity of information for other species. This is especially striking for butterflies, even though they have diverse and derived phenotypes (such as color vision and wing color patterns) and are considered prime models for the evolutionary and developmental analysis of ecologically relevant, complex traits. We focus on Bicyclus anynana butterflies, a laboratory system for studying the diversification of novelties and serially repeated traits. With a panel of 12 small families and a biphasic mapping approach, we first assigned 508 expressed genes to segregation groups and then ordered 297 of them within individual linkage groups. We also coarsely mapped seven color pattern loci. This is the richest gene-based map available for any butterfly species and allowed for a broad-coverage analysis of synteny with the lepidopteran reference genome. Based on 462 pairs of mapped orthologous markers in Bi. anynana and Bo. mori, we observed strong conservation of gene assignment to chromosomes, but also evidence for numerous large- and small-scale chromosomal rearrangements. With gene collections growing for a variety of target organisms, the ability to place those genes in their proper genomic context is paramount. Methods to map expressed genes and to compare maps with relevant model systems are crucial to extend genomic-level analysis outside classical model species. Maps with gene-based markers are useful for comparative genomics and to resolve mapped genomic regions to a tractable number of candidate genes, especially if there is synteny with related model species. This is discussed in relation to the identification of the loci contributing to color pattern evolution in butterflies.

There are 15 new articles in PLoS ONE today. As always, you should rate the articles, post notes and comments and send trackbacks when you blog about the papers. You can now also easily place articles on various social services (CiteULike, Connotea, Stumbleupon, Facebook and Digg) with just one click. Here are my own picks for the week – you go and look for your own favourites:
Web Queries as a Source for Syndromic Surveillance:

In the field of syndromic surveillance, various sources are exploited for outbreak detection, monitoring and prediction. This paper describes a study on queries submitted to a medical web site, with influenza as a case study. The hypothesis of the work was that queries on influenza and influenza-like illness would provide a basis for the estimation of the timing of the peak and the intensity of the yearly influenza outbreaks that would be as good as the existing laboratory and sentinel surveillance. We calculated the occurrence of various queries related to influenza from search logs submitted to a Swedish medical web site for two influenza seasons. These figures were subsequently used to generate two models, one to estimate the number of laboratory verified influenza cases and one to estimate the proportion of patients with influenza-like illness reported by selected General Practitioners in Sweden. We applied an approach designed for highly correlated data, partial least squares regression. In our work, we found that certain web queries on influenza follow the same pattern as that obtained by the two other surveillance systems for influenza epidemics, and that they have equal power for the estimation of the influenza burden in society. Web queries give a unique access to ill individuals who are not (yet) seeking care. This paper shows the potential of web queries as an accurate, cheap and labour extensive source for syndromic surveillance.

The Fetal Hypothalamus Has the Potential to Generate Cells with a Gonadotropin Releasing Hormone (GnRH) Phenotype:

Neurospheres (NS) are colonies of neural stem and precursor cells capable of differentiating into the central nervous system (CNS) cell lineages upon appropriate culture conditions: neurons, and glial cells. NS were originally derived from the embryonic and adult mouse striatum subventricular zone. More recently, experimental evidence substantiated the isolation of NS from almost any region of the CNS, including the hypothalamus. Here we report a protocol that enables to generate large quantities of NS from both fetal and adult rat hypothalami. We found that either FGF-2 or EGF were capable of inducing NS formation from fetal hypothalamic cultures, but that only FGF-2 is effective in the adult cultures. The hypothalamic-derived NS are capable of differentiating into neurons and glial cells and most notably, as demonstrated by immunocytochemical detection with a specific anti-GnRH antibody, the fetal cultures contain cells that exhibit a GnRH phenotype upon differentiation. This in vitro model should be useful to study the molecular mechanisms involved in GnRH neuronal differentiation.

The Cost of Simplifying Air Travel When Modeling Disease Spread:

Air travel plays a key role in the spread of many pathogens. Modeling the long distance spread of infectious disease in these cases requires an air travel model. Highly detailed air transportation models can be over determined and computationally problematic. We compared the predictions of a simplified air transport model with those of a model of all routes and assessed the impact of differences on models of infectious disease. Using U.S. ticket data from 2007, we compared a simplified “pipe” model, in which individuals flow in and out of the air transport system based on the number of arrivals and departures from a given airport, to a fully saturated model where all routes are modeled individually. We also compared the pipe model to a “gravity” model where the probability of travel is scaled by physical distance; the gravity model did not differ significantly from the pipe model. The pipe model roughly approximated actual air travel, but tended to overestimate the number of trips between small airports and underestimate travel between major east and west coast airports. For most routes, the maximum number of false (or missed) introductions of disease is small (<1 per day) but for a few routes this rate is greatly underestimated by the pipe model. If our interest is in large scale regional and national effects of disease, the simplified pipe model may be adequate. If we are interested in specific effects of interventions on particular air routes or the time for the disease to reach a particular location, a more complex point-to-point model will be more accurate. For many problems a hybrid model that independently models some frequently traveled routes may be the best choice. Regardless of the model used, the effect of simplifications and sensitivity to errors in parameter estimation should be analyzed.

Antiproliferative Effect of Ascorbic Acid Is Associated with the Inhibition of Genes Necessary to Cell Cycle Progression:

Ascorbic acid (AA), or Vitamin C, is most well known as a nutritional supplement with antioxidant properties. Recently, we demonstrated that high concentrations of AA act on PMP22 gene expression and partially correct the Charcot-Marie-Tooth disease phenotype in a mouse model. This is due to the capacity of AA, but not other antioxidants, to down-modulate cAMP intracellular concentration by a competitive inhibition of the adenylate cyclase enzymatic activity. Because of the critical role of cAMP in intracellular signalling, we decided to explore the possibility that ascorbic acid could modulate the expression of other genes. Using human pangenomic microarrays, we found that AA inhibited the expression of two categories of genes necessary for cell cycle progression, tRNA synthetases and translation initiation factor subunits. In in vitro assays, we demonstrated that AA induced the S-phase arrest of proliferative normal and tumor cells. Highest concentrations of AA leaded to necrotic cell death. However, quiescent cells were not susceptible to AA toxicity, suggesting the blockage of protein synthesis was mainly detrimental in metabolically-active cells. Using animal models, we found that high concentrations of AA inhibited tumor progression in nude mice grafted with HT29 cells (derived from human colon carcinoma). Consistently, expression of tRNA synthetases and ieF2 appeared to be specifically decreased in tumors upon AA treatment. AA has an antiproliferative activity, at elevated concentration that could be obtained using IV injection. This activity has been observed in vitro as well in vivo and likely results from the inhibition of expression of genes involved in protein synthesis. Implications for a clinical use in anticancer therapies will be discussed.

Continue reading →

There are 16 new articles in PLoS ONE today. As always, you should rate the articles, post notes and comments and send trackbacks when you blog about the papers. You can now also easily place articles on various social services (CiteULike, Connotea, Stumbleupon, Facebook and Digg) with just one click. Here are my own picks for the week – you go and look for your own favourites:
Climate Change and Trophic Response of the Antarctic Bottom Fauna:

As Earth warms, temperate and subpolar marine species will increasingly shift their geographic ranges poleward. The endemic shelf fauna of Antarctica is especially vulnerable to climate-mediated biological invasions because cold temperatures currently exclude the durophagous (shell-breaking) predators that structure shallow-benthic communities elsewhere. We used the Eocene fossil record from Seymour Island, Antarctic Peninsula, to project specifically how global warming will reorganize the nearshore benthos of Antarctica. A long-term cooling trend, which began with a sharp temperature drop ~41 Ma (million years ago), eliminated durophagous predators–teleosts (modern bony fish), decapod crustaceans (crabs and lobsters) and almost all neoselachian elasmobranchs (modern sharks and rays)–from Antarctic nearshore waters after the Eocene. Even prior to those extinctions, durophagous predators became less active as coastal sea temperatures declined from 41 Ma to the end of the Eocene, ~33.5 Ma. In response, dense populations of suspension-feeding ophiuroids and crinoids abruptly appeared. Dense aggregations of brachiopods transcended the cooling event with no apparent change in predation pressure, nor were there changes in the frequency of shell-drilling predation on venerid bivalves. Rapid warming in the Southern Ocean is now removing the physiological barriers to shell-breaking predators, and crabs are returning to the Antarctic Peninsula. Over the coming decades to centuries, we predict a rapid reversal of the Eocene trends. Increasing predation will reduce or eliminate extant dense populations of suspension-feeding echinoderms from nearshore habitats along the Peninsula while brachiopods will continue to form large populations, and the intensity of shell-drilling predation on infaunal bivalves will not change appreciably. In time the ecological effects of global warming could spread to other portions of the Antarctic coast. The differential responses of faunal components will reduce the endemic character of Antarctic subtidal communities, homogenizing them with nearshore communities at lower latitudes.

Risks for Central Nervous System Diseases among Mobile Phone Subscribers: A Danish Retrospective Cohort Study:

The aim of this study was to investigate a possible link between cellular telephone use and risks for various diseases of the central nervous system (CNS). We conducted a large nationwide cohort study of 420 095 persons whose first cellular telephone subscription was between 1982 and 1995, who were followed through 2003 for hospital contacts for a diagnosis of a CNS disorder. Standardized hospitalization ratios (SHRs) were derived by dividing the number of hospital contacts in the cohort by the number expected in the Danish population. The SHRs were increased by 10-20% for migraine and vertigo. No associations were seen for amyotrophic lateral sclerosis, multiple sclerosis or epilepsy in women. SHRs decreased by 30-40% were observed for dementia (Alzheimer disease, vascular and other dementia), Parkinson disease and epilepsy among men. In analyses restricted to subscribers of 10 years or more, the SHRs remained similarly increased for migraine and vertigo and similarly decreased for Alzheimer disease and other dementia and epilepsy (in men); the other SHRs were close to unity. In conclusion, the excesses of migraine and vertigo observed in this first study on cellular telephones and CNS disease deserve further attention. An interplay of a healthy cohort effect and reversed causation bias due to prodromal symptoms impedes detection of a possible association with dementia and Parkinson disease. Identification of the factors that result in a healthy cohort might be of interest for elucidation of the etiology of these diseases.

A Ligand Channel through the G Protein Coupled Receptor Opsin:

The G protein coupled receptor rhodopsin contains a pocket within its seven-transmembrane helix (TM) structure, which bears the inactivating 11-cis-retinal bound by a protonated Schiff-base to Lys296 in TM7. Light-induced 11-cis-/all-trans-isomerization leads to the Schiff-base deprotonated active Meta II intermediate. With Meta II decay, the Schiff-base bond is hydrolyzed, all-trans-retinal is released from the pocket, and the apoprotein opsin reloaded with new 11-cis-retinal. The crystal structure of opsin in its active Ops* conformation provides the basis for computational modeling of retinal release and uptake. The ligand-free 7TM bundle of opsin opens into the hydrophobic membrane layer through openings A (between TM1 and 7), and B (between TM5 and 6), respectively. Using skeleton search and molecular docking, we find a continuous channel through the protein that connects these two openings and comprises in its central part the retinal binding pocket. The channel traverses the receptor over a distance of ca. 70 Å and is between 11.6 and 3.2 Å wide. Both openings are lined with aromatic residues, while the central part is highly polar. Four constrictions within the channel are so narrow that they must stretch to allow passage of the retinal β-ionone-ring. Constrictions are at openings A and B, respectively, and at Trp265 and Lys296 within the retinal pocket. The lysine enforces a 90° elbow-like kink in the channel which limits retinal passage. With a favorable Lys side chain conformation, 11-cis-retinal can take the turn, whereas passage of the all-trans isomer would require more global conformational changes. We discuss possible scenarios for the uptake of 11-cis- and release of all-trans-retinal. If the uptake gate of 11-cis-retinal is assigned to opening B, all-trans is likely to leave through the same gate. The unidirectional passage proposed previously requires uptake of 11-cis-retinal through A and release of photolyzed all-trans-retinal through B.

Phylogenetic Resolution and Quantifying the Phylogenetic Diversity and Dispersion of Communities:

Conservation biologists and community ecologists have increasingly begun to quantify the phylogenetic diversity and phylogenetic dispersion in species assemblages. In some instances, the phylogenetic trees used for such analyses are fully bifurcating, but in many cases the phylogenies being used contain unresolved nodes (i.e. polytomies). The lack of phylogenetic resolution in such studies, while certainly not preferred, is likely to continue particularly for those analyzing diverse communities and datasets with hundreds to thousands of taxa. Thus it is imperative that we quantify potential biases and losses of statistical power in studies that use phylogenetic trees that are not completely resolved. The present study is designed to meet both of these goals by quantifying the phylogenetic diversity and dispersion of simulated communities using resolved and gradually ‘unresolved’ phylogenies. The results show that: (i) measures of community phylogenetic diversity and dispersion are generally more sensitive to loss of resolution basally in the phylogeny and less sensitive to loss of resolution terminally; and (ii) the loss of phylogenetic resolution generally causes false negative results rather than false positives.

Probing Evolutionary Patterns in Neotropical Birds through DNA Barcodes:

The Neotropical avifauna is more diverse than that of any other biogeographic region, but our understanding of patterns of regional divergence is limited. Critical examination of this issue is currently constrained by the limited genetic information available. This study begins to address this gap by assembling a library of mitochondrial COI sequences, or DNA barcodes, for Argentinian birds and comparing their patterns of genetic diversity to those of North American birds. Five hundred Argentinian species were examined, making this the first major examination of DNA barcodes for South American birds. Our results indicate that most southern Neotropical bird species show deep sequence divergence from their nearest-neighbour, corroborating that the high diversity of this fauna is not based on an elevated incidence of young species radiations. Although species ages appear similar in temperate North and South American avifaunas, patterns of regional divergence are more complex in the Neotropics, suggesting that the high diversity of the Neotropical avifauna has been fueled by greater opportunities for regional divergence. Deep genetic splits were observed in at least 21 species, though distribution patterns of these lineages were variable. The lack of shared polymorphisms in species, even in species with less than 0.5M years of reproductive isolation, further suggests that selective sweeps could regularly excise ancestral mitochondrial polymorphisms. These findings confirm the efficacy of species delimitation in birds via DNA barcodes, even when tested on a global scale. Further, they demonstrate how large libraries of a standardized gene region provide insight into evolutionary processes.

Ticks Associated with Macquarie Island Penguins Carry Arboviruses from Four Genera:

Macquarie Island, a small subantarctic island, is home to rockhopper, royal and king penguins, which are often infested with the globally distributed seabird tick, Ixodes uriae. A flavivirus, an orbivirus, a phlebovirus, and a nairovirus were isolated from these ticks and partial sequences obtained. The flavivirus was nearly identical to Gadgets Gully virus, isolated some 30 year previously, illustrating the remarkable genetic stability of this virus. The nearest relative to the orbivirus (for which we propose the name Sandy Bay virus) was the Scottish Broadhaven virus, and provided only the second available sequences from the Great Island orbivirus serogroup. The phlebovirus (for which we propose the name Catch-me-cave virus) and the previously isolated Precarious Point virus were distinct but related, with both showing homology with the Finnish Uukuniemi virus. These penguin viruses provided the second and third available sequences for the Uukuniemi group of phleboviruses. The nairovirus (for which we propose the name Finch Creek virus) was shown to be related to the North American Tillamook virus, the Asian Hazara virus and Nairobi sheep disease virus. Macquarie Island penguins thus harbour arboviruses from at least four of the seven arbovirus-containing genera, with related viruses often found in the northern hemisphere.

Evolution of Genome Size and Complexity in Pinus:

Genome evolution in the gymnosperm lineage of seed plants has given rise to many of the most complex and largest plant genomes, however the elements involved are poorly understood. Gymny is a previously undescribed retrotransposon family in Pinus that is related to Athila elements in Arabidopsis. Gymny elements are dispersed throughout the modern Pinus genome and occupy a physical space at least the size of the Arabidopsis thaliana genome. In contrast to previously described retroelements in Pinus, the Gymny family was amplified or introduced after the divergence of pine and spruce (Picea). If retrotransposon expansions are responsible for genome size differences within the Pinaceae, as they are in angiosperms, then they have yet to be identified. In contrast, molecular divergence of Gymny retrotransposons together with other families of retrotransposons can account for the large genome complexity of pines along with protein-coding genic DNA, as revealed by massively parallel DNA sequence analysis of Cot fractionated genomic DNA. Most of the enormous genome complexity of pines can be explained by divergence of retrotransposons, however the elements responsible for genome size variation are yet to be identified. Genomic resources for Pinus including those reported here should assist in further defining whether and how the roles of retrotransposons differ in the evolution of angiosperm and gymnosperm genomes.

Continue reading →

There is some interesting stuff published in PLoS Medicine and PLoS Biology today:
What Should Be Done To Tackle Ghostwriting in the Medical Literature?:

Background to the debate: Ghostwriting occurs when someone makes substantial contributions to a manuscript without attribution or disclosure. It is considered bad publication practice in the medical sciences, and some argue it is scientific misconduct. At its extreme, medical ghostwriting involves pharmaceutical companies hiring professional writers to produce papers promoting their products but hiding those contributions and instead naming academic physicians or scientists as the authors. To improve transparency, many editors’ associations and journals allow professional medical writers to contribute to the writing of papers without being listed as authors provided their role is acknowledged. This debate examines how best to tackle ghostwriting in the medical literature from the perspectives of a researcher, an editor, and the professional medical writer.

Mutation Patterns in the Human Genome: More Variable Than Expected:

The development, survival, and reproduction of an organism depend on the genetic information that is carried in its genome, yet the transmission of genetic information is not perfectly accurate: new mutations occur at each generation. These mutations are the primary cause of the genetic diversity on which natural selection can operate, and hence are the sine qua non of evolution. A better knowledge of mutation processes is crucial for investigating the causes of genetic diseases or cancer and for understanding evolutionary processes. This knowledge is also important for different practical reasons. First, comparative sequence analysis is widely used to find functional elements within genomes. The basic principle of this approach is that functional elements are affected by natural selection, and hence can be recognized because they evolve either slower or faster than expected given the local mutation rate. Hence, to be able to annotate genomic sequences, it is necessary to have a good knowledge of the underlying pattern of mutation. Moreover, this knowledge is also essential for ensuring the accuracy of the methods that analyze sequence divergence to determine the phylogeny of species or the demographical history of populations. Finally, the study of mutational processes also provides valuable information about genome function in processes such as replication, repair, transcription, and recombination. During the last few years, several important factors affecting mutation rates have been uncovered. However, a paper in this issue of PLoS Biology [1] reveals an unexpected additional layer of complexity in the determinants of mutation rates.

Cryptic Variation in the Human Mutation Rate:

Understanding the process of mutation is important, not only mechanistically, but also because it has implications for the analysis of sequence evolution and population genetic inference. The mutation rate is known to differ between sites within the human genome. The most dramatic example of this is when a C is followed by G; both the C and G nucleotides have a rate of mutation that is between 10- and 20-fold higher than the rate at other sites. In addition, is it known that the mutation rate may be influenced by the nucleotides flanking the site. Here we show that there is also very substantial variation in the mutation rate that is not associated with the flanking nucleotides, or the CpG effect. Although this variation does not depend upon the adjacent nucleotides, there are nonrandom patterns of nucleotides surrounding sites that appear to be hypermutable, suggesting there are complex context effects that influence the mutation rate.

Gene Regulatory Network Interactions in Sea Urchin Endomesoderm Induction:

In recent years, “gene regulatory networks” (GRNs) have provided integrated views of gene interactions that control biological processes. One of the earliest networks to be activated in the developing zygotes is the one controlling endomesoderm development. In the sea urchin, this network includes several subnetworks that function in adjacent tiers of cells that form the endoderm and mesoderm of the developing embryo. Although classic embryological manipulations have shown that the precursors of the embryonic skeleton induce endomesoderm fate in adjacent cells, the GRNs regulating this interaction are not understood. To investigate these networks, we ectopically activated a GRN that operates in skeletogenic precursors and characterized the responding GRN in neighboring cells, which adopt an endomesoderm fate. By testing the responsiveness of every core factor in the responding GRN, which allowed us to identify a subset that executes the response to the induction, we demonstrated that the signaling molecule, ActivinB, is an essential component of this induction and that its function is physiologically relevant: it is required during normal embryonic development to activate the same GRN that responds to signals from skeletogenic precursors. Furthermore, the network response to ActivinB signaling reveals greater complexity in an additional uncharacterized inductive signal emitted by skeletogenic precursors. Our results thus highlight how interacting GRNs can be used to understand a fundamental signaling process.

Continue reading →

Continue reading →