Population genomics from pool sequencing

Next generation sequencing of pooled samples is an effective approach for studies of variability and differentiation in populations. In this paper we provide a comprehensive set of estimators of the most common statistics in population genetics based on the frequency spectrum, namely the Watterson estimator , nucleotide pairwise diversity Π, Tajima's D, Fu and Li's D and F, Fay and Wu's H, McDonald‐Kreitman and HKA tests and , corrected for sequencing errors and ascertainment bias. In a simulation study, we show that pool and individual θ estimates are highly correlated and discuss how the performance of the statistics vary with read depth and sample size in different evolutionary scenarios. As an application, we reanalyse sequences from Drosophila mauritiana and from an evolution experiment in Drosophila melanogaster. These methods are useful for population genetic projects with limited budget, study of communities of individuals that are hard to isolate, or autopolyploid species.     

犬齿窝与人类中面部骨骼的演化

犬齿窝是包括现代人类在内的许多人族成员面部骨骼的重要性状,但在分类学上的意义仍存在争议。有学者认为该性状是一个发生于基础面形的近祖性状,除了一些例外,在灭绝和现生的大猿及人属中都存在。另有学者认为,犬齿窝是仅存在于智人及其直系祖先的衍生性状,在发育上与颧齿槽突嵴有关。这种关系并非总是成立,在智人中存在着明显的差异：弧型颧齿槽突嵴和直斜型颧齿槽突嵴与犬齿窝有时共存,有时不共存。我们由此推测,犬齿窝的发生和形态与上颌窦的前部发育有关,颧齿槽突嵴的形态与鼻窦的侧面发育有关。在人类演化的过程中,犬齿窝经历了不同的变形,比如上颌沟(如南方古猿非洲种、傍人粗壮种)、上颌小窝(如傍人粗壮种)、上颌沟(如匠人)或犬齿窝缺如(如傍人埃塞俄比亚种、傍人鲍氏种、肯尼亚扁脸人、人属鲁道夫种)。犬齿窝消失的原因各类群并不相同,如中新世和早更新世人属以及中更新世人属(如人属海德堡种/人属罗得西亚种、人属尼安德特种)。人属罗得西亚种具有弱化的犬齿窝,不具备演化为智人的可能,因此被排除在智人的演化支之外。     

Nitric oxide modulates murine yolk sac vasculogenesis and rescues glucose induced vasculopathy

Nitric oxide (NO) has been demonstrated to mediate events during ovulation, pregnancy, blastocyst invasion and preimplantation embryogenesis. However, less is known about the role of NO during postimplantation development. Therefore, in this study, we explored the effects of NO during vascular development of the murine yolk sac, which begins shortly after implantation. Establishment of the vitelline circulation is crucial for normal embryonic growth and development. Moreover, functional inactivation of the endodermal layer of the yolk sac by environmental insults or genetic manipulations during this period leads to embryonic defects/lethality, as this structure is vital for transport, metabolism and induction of vascular development. In this study, we describe the temporally/spatially regulated distribution of nitric oxide synthase (NOS) isoforms during the three stages of yolk sac vascular development (blood island formation, primary capillary plexus formation and vessel maturation/remodeling) and found NOS expression patterns were diametrically opposed. To pharmacologically manipulate vascular development, an established in vitro system of whole murine embryo culture was employed. During blood island formation, the endoderm produced NO and inhibition of NO (L-NMMA) at this stage resulted in developmental arrest at the primary plexus stage and vasculopathy. Furthermore, administration of a NO donor did not cause abnormal vascular development; however, exogenous NO correlated with increased eNOS and decreased iNOS protein levels. Additionally, a known environmental insult (high glucose) that produces reactive oxygen species (ROS) and induces vasculopathy also altered eNOS/iNOS distribution and induced NO production during yolk sac vascular development. However, administration of a NO donor rescued the high glucose induced vasculopathy, restored the eNOS/iNOS distribution and decreased ROS production. These data suggest that NO acts as an endoderm-derived factor that modulates normal yolk sac vascular development, and decreased NO bioavailability and NO-mediated sequela may underlie high glucose induced vasculopathy.     

Elevated glucose inhibits VEGF-A-mediated endocardial cushion formation: modulation by PECAM-1 and MMP-2

Atrioventricular (AV) septal defects resulting from aberrant endocardial cushion (EC) formation are observed at increased rates in infants of diabetic mothers. EC formation occurs via an epithelial-mesenchymal transformation (EMT), involving transformation of endocardial cells into mesenchymal cells, migration, and invasion into extracellular matrix. Here, we report that elevated glucose inhibits EMT by reducing myocardial vascular endothelial growth factor A (VEGF-A). This effect is reversed with exogenous recombinant mouse VEGF-A165, whereas addition of soluble VEGF receptor-1 blocks EMT. We show that disruption of EMT is associated with persistence of platelet endothelial cell adhesion molecule-1 (PECAM-1) and decreased matrix metalloproteinase-2 (MMP-2) expression. These findings correlate with retention of a nontransformed endocardial sheet and lack of invasion. The MMP inhibitor GM6001 blocks invasion, whereas explants from PECAM-1 deficient mice exhibit MMP-2 induction and normal EMT in high glucose. PECAM-1-negative endothelial cells are highly motile and express more MMP-2 than do PECAM-1-positive endothelial cells. During EMT, loss of PECAM-1 similarly promotes single cell motility and MMP-2 expression. Our findings suggest that high glucose-induced inhibition of AV cushion morphogenesis results from decreased myocardial VEGF-A expression and is, in part, mediated by persistent endocardial cell PECAM-1 expression and failure to up-regulate MMP-2 expression.     

Pipeliner: software to evaluate the performance of bioinformatics pipelines for next‐generation resequencing

The choice of technology and bioinformatics approach is critical in obtaining accurate and reliable information from next‐generation sequencing (NGS) experiments. An increasing number of software and methodological guidelines are being published, but deciding upon which approach and experimental design to use can depend on the particularities of the species and on the aims of the study. This leaves researchers unable to produce informed decisions on these central questions. To address these issues, we developed pipeliner – a tool to evaluate, by simulation, the performance of NGS pipelines in resequencing studies. Pipeliner provides a graphical interface allowing the users to write and test their own bioinformatics pipelines with publicly available or custom software. It computes a number of statistics summarizing the performance in SNP calling, including the recovery, sensitivity and false discovery rate for heterozygous and homozygous SNP genotypes. Pipeliner can be used to answer many practical questions, for example, for a limited amount of NGS effort, how many more reliable SNPs can be detected by doubling coverage and halving sample size or what is the false discovery rate provided by different SNP calling algorithms and options. Pipeliner thus allows researchers to carefully plan their study's sampling design and compare the suitability of alternative bioinformatics approaches for their specific study systems. Pipeliner is written in C++ and is freely available from http://github.com/brunonevado/Pipeliner .     

A model of direction selectivity in the starburst amacrine cell network

Displaced starburst amacrine cells (SACs) are retinal interneurons that exhibit GABA _A receptor-mediated and Cl ^? cotransporter-mediated, directionally selective (DS) light responses in the rabbit retina. They depolarize to stimuli that move centrifugally through the receptive field surround and hyperpolarize to stimuli that move centripetally through the surround (Gavrikov et al, PNAS 100(26):16047–16052, 2003, PNAS 103(49):18793–18798, 2006). They also play a key role in the activity of DS ganglion cells (DS GC; Amthor et al, Vis Neurosci 19:495–509 2002; Euler et al, Nature 418:845–852, 2002; Fried et al, Nature 420:411– 414, 2002; Gavrikov et al, PNAS 100(26):16047–16052, 2003, PNAS 103(49):18793–18798, 2006; Lee and Zhou, Neuron 51:787–799 2006; Yoshida et al, Neuron 30:771–780, 2001). In this paper we present a model of strong DS behavior of SACs which relies on the GABA-mediated communication within a tightly interconnected network of these cells and on the glutamate signal that the SACs receive from bipolar cells (a presynaptic cell that receives input from cones). We describe how a moving light stimulus can produce a large, sustained depolarization of the SAC dendritic tips that point in the direction that the stimulus moves (i.e., centrifugal motion), but produce a minimal depolarization of the dendritic tips that point in the opposite direction (i.e., centripetal motion). This DS behavior, which is quantified based on the relative size and duration of the depolarizations evoked by stimulus motion at dendritic tips pointing in opposite directions, is robust to changes of many different parameter values and consistent with experimental data. In addition, the DS behavior is strengthened under the assumptions that the Cl^? cotransporters Na^?+?-K^?+?-Cl^?? and K^?+?-Cl^?? are located in different regions of the SAC dendritic tree (Gavrikov et al, PNAS 103(49):18793–18798, 2006) and that GABA evokes a long-lasting response (Gavrikov et al, PNAS 100(26):16047–16052, 2003, PNAS 103(49):18793–18798, 2006; Lee and Zhou, Neuron 51:787–799, 2006). A possible mechanism is discussed based on the generation of waves of local glutamate and GABA secretion, and their postsynaptic interplay as the waves travel between cell compartments.     

The chimerical genome of Isla del Coco feral pigs (Costa Rica), an isolated population since 1793 but with remarkable levels of diversity

The history of domestic species and of their wild ancestors is not a simple one, and feral processes can clarify key aspects of this history, including the adaptive processes triggered by new environments. Here, we provide a comprehensive genomic study of Isla del Coco (Costa Rica) feral pigs, a unique population that was allegedly founded by two individuals and has remained isolated since 1793. Using SNP arrays and genome sequencing, we show that Cocos pigs are hybrids between Asian and European pigs, as are modern international pig breeds. This conclusively shows that, as early as the 18th century, British vessels were loading crossbred pigs in Great Britain and transporting them overseas. We find that the Y chromosome has Asian origin, which has not been reported in any international pig breed. Chinese haplotypes seem to have been transmitted independently between Cocos and other pig breeds, suggesting independent introgression events and a complex pattern of admixing. Although data are compatible with a founder population of N = 2, variability levels are as high in Cocos pigs as in international pig breeds (~1.9 SNPs/kb) and higher than in European wild boars or local breeds (~1.7 SNPs/kb). Nevertheless, we also report a 10‐Mb region with a marked decrease in variability across all samples that contains four genes (CPE, H3F3C, SC4MOL and KHL2) previously identified as highly differentiated between wild and domestic pigs. This work therefore illustrates how feral population genomic studies can help to resolve the history of domestic species and associated admixture events.     

A new method to analyze boar sperm DNA fragmentation under bright-field or fluorescence microscopy

We present a new, rapid and simple method to study DNA Fragmentation Index (DFI) in sperm samples from boar under bright-field and fluorescence microscopy. Discrimination of sperm cells containing fragmented DNA relies on the extreme peripheral diffusion of their chromatin fragments, whereas those sperm nuclei without DNA fragmentation do not disperse or show very restricted spreading of DNA loops close to the flagellum. The basic methodology provided in the commercial kit Sperm-Sus-Halomax allows, in addition to a direct estimation of DFI in a sperm sample under bright field microscopy, a direct visualization of DNA breaks by incorporation of labelled nucleotides using the DNA polymerase I following the in situ nick translation assay (ISNT methodology not provided in the kit). An external control using DBD-FISH (DNA breakage detection-fluorescence in situ hybridization) on human and boar sperm samples was used in this experiment. The results obtained show (i) low levels of background DNA fragmentation (from 0.7 to 10%), (ii) no significant differences for DFI after the application of Sperm-Sus-Halomax and ISNT, with a tendency to be underestimated after using DBD-FISH and (iii) a characteristic chromatin organization in boar sperm nucleus, with a particular response to chromatin loop relaxation and preferential DNA labelling by ISNT at the proximal nuclear area, close to the flagellum. This methodology allows the routine assessment of boar sperm samples for DFI, as well as basic and clinical research on this relevant topic in any laboratory of semen analysis.     

Protein scaffolds can enhance the bistability of multisite phosphorylation systems

The phosphorylation of a substrate at multiple sites is a common protein modification that can give rise to important structural and electrostatic changes. Scaffold proteins can enhance protein phosphorylation by facilitating an interaction between a protein kinase enzyme and its target substrate. In this work we consider a simple mathematical model of a scaffold protein and show that under specific conditions, the presence of the scaffold can substantially raise the likelihood that the resulting system will exhibit bistable behavior. This phenomenon is especially pronounced when the enzymatic reactions have sufficiently large K(M), compared to the concentration of the target substrate. We also find for a closely related model that bistable systems tend to have a specific kinetic conformation. Using deficiency theory and other methods, we provide a number of necessary conditions for bistability, such as the presence of multiple phosphorylation sites and the dependence of the scaffold binding/unbinding rates on the number of phosphorylated sites.     

Evaluation in a Dog Model of Three Antimicrobial Glassy Coatings: Prevention of Bone Loss around Implants and Microbial Assessments

ObjectivesThe aim of the present study is to evaluate, in a ligature-induced peri-implantitis model, the efficacy of three antimicrobial glassy coatings in the prevention of biofilm formation, intrasulcular bacterial growth and the resulting peri-implant bone loss.MethodsMandibular premolars were bilaterally extracted from five beagle dogs. Four dental implants were inserted on each hemiarch. Eight weeks after, one control zirconia abutment and three with different bactericidal coatings (G1n-Ag, ZnO35, G3) were connected. After a plaque control period, bacterial accumulation was allowed and biofilm formation on abutments was observed by Scanning Electron Microscopy (SEM). Peri-implantitis was induced by cotton ligatures. Microbial samples and peri-implant crestal bone levels of all implant sites were obtained before, during and after the breakdown period.ResultsDuring experimental induce peri-implantitis: colony forming units counts from intrasulcular microbial samples at implants with G1n-Ag coated abutment remained close to the basal inoculum; G3 and ZnO35 coatings showed similar low counts; and anaerobic bacterias counts at control abutments exhibited a logarithmic increase by more than 2. Bone loss during passive breakdown period was no statistically significant. Additional bone loss occurred during ligature-induce breakdown: 0.71 (SD 0.48) at G3 coating, 0.57 (SD 0.36) at ZnO35 coating, 0.74 (SD 0.47) at G1n-Ag coating, and 1.29 (SD 0.45) at control abutments; and statistically significant differences (p<0.001) were found. The lowest bone loss at the end of the experiment was exhibited by implants dressing G3 coated abutments (mean 2.1; SD 0.42).SignificanceAntimicrobial glassy coatings could be a useful tool to ward off, diminish or delay peri-implantitis progression.