An Introduction to Programming for Bioscientists: A Python-based Primer

Computing has revolutionized the biological sciences over the past several decades, such that virtually all contemporary research in the biosciences utilizes computer programs. The computational advances have come on many fronts, spurred by fundamental developments in hardware, software, and algorithms. These advances have influenced, and even engendered, a phenomenal array of bioscience fields, including molecular evolution and bioinformatics; genome-, proteome-, transcriptome- and metabolome-wide experimental studies; structural genomics; and atomistic simulations of cellular-scale molecular assemblies as large as ribosomes and intact viruses. In short, much of post-genomic biology is increasingly becoming a form of computational biology. The ability to design and write computer programs is among the most indispensable skills that a modern researcher can cultivate. Python has become a popular programming language in the biosciences, largely because (i) its straightforward semantics and clean syntax make it a readily accessible first language; (ii) it is expressive and well-suited to object-oriented programming, as well as other modern paradigms; and (iii) the many available libraries and third-party toolkits extend the functionality of the core language into virtually every biological domain (sequence and structure analyses, phylogenomics, workflow management systems, etc.). This primer offers a basic introduction to coding, via Python, and it includes concrete examples and exercises to illustrate the language's usage and capabilities; the main text culminates with a final project in structural bioinformatics. A suite of Supplemental Chapters is also provided. Starting with basic concepts, such as that of a 'variable', the Chapters methodically advance the reader to the point of writing a graphical user interface to compute the Hamming distance between two DNA sequences.Comment: 65 pages total, including 45 pages text, 3 figures, 4 tables, numerous exercises, and 19 pages of Supporting Information; currently in press at PLOS Computational Biolog

Frequency-tunable metamaterials using broadside-coupled split ring resonators

We present frequency tunable metamaterial designs at terahertz (THz) frequencies using broadside-coupled split ring resonator (BC-SRR) arrays. Frequency tuning, arising from changes in near field coupling, is obtained by in-plane horizontal or vertical displacements of the two SRR layers. For electrical excitation, the resonance frequency continuously redshifts as a function of displacement. The maximum frequency shift occurs for displacement of half a unit cell, with vertical displacement resulting in a shift of 663 GHz (51% of f0) and horizontal displacement yielding a shift of 270 GHz (20% of f0). We also discuss the significant differences in tuning that arise for electrical excitation in comparison to magnetic excitation of BC-SRRs

Conodonts of Middle Moscovian Age from the Kongul Formation (Bolkardag Unit), Northwest of Hadim, Central Taurus, Turkey

A rich shallow water conodont fauna with Idiognathodus incurvus DUNN, I. cf. suberectus (DUNN), Idioprioniodus sp., and Neognathodus columbiensis STIBANE was found in the lower part of the Kongul Formation of the Bolkardag Unit, about 12 km NW of Hadim town, central Taurus, Turkey. The conodonts indicate a middle to late Atokan age (on the American scale) for this part of the formation. On the Russian scale, the fauna corresponds to the middle Moscovian. According to previous correlations between the Russian and North American scales, a (middle) late Atokan (middle Moscovian age) is assumed for this fauna. Variscan deep-water sequences are not present in the study area

Surface preparation of powder metallurgical tool steels by means of wire electrical discharge machining

The surface of two types of powder metallurgical (PM) tool steels (i.e., with and without nitrogen) was prepared using wire electrical discharge machining (WEDM). From each grade of tool steel, seven surfaces corresponding to one to seven passes of WEDM were prepared. The WEDM process was carried out using a brass wire as electrode and deionized water as dielectric. After eachWEDM pass the surface of the tool steels was thoroughly examined. Surface residual stresses were measured by the X-ray diffraction (XRD) technique. The measured stresses were found to be of tensile nature. The surface roughness of the WEDM specimens was measured using interference microscopy. The surface roughness as well as the residual stress measurements indicated an insignificant improvement of these parameters after four passes of WEDM. In addition, the formed recast layer was characterized by means of scanning electron microscopy (SEM), XRD, and X-ray photoelectron spectroscopy (XPS). The characterization investigation clearly shows diffusion of copper and zinc from the wire electrode into the work material, even after the final WEDM step. Finally, the importance of eliminating excessive WEDM steps is thoroughly discussed

LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.We acknowledge “Red Española de Surpercomputación” (RES) for awarding us access to the MareNostrum III machine based in Barcelona, Spain (Ref. FI-2015-2-0026 and FI-2015-3-0011). We also acknowledge PRACE for awarding us access to Fermi and Marconi Supercomputers at Cineca, Italy (Ref. 2015133120). Oriol Lehmkuhl acknowledges a PDJ 2014 Grant by AGAUR (Generalitat de Catalunya). Ugo Piomelli acknowledges the support of the Natural Sciences and Engineering Research Council (NSERC) of Canada under the Discovery Grant Programme (Grant No. RGPIN-2016-04391). Ricard Borrell acknowledges a Juan de la Cierva postdoctoral grant (IJCI-2014-21034). Ivette Rodriguez, Oriol Lehmkuhl, Ricard Borrell and Assensi Oliva acknowledge Ministerio de Economía y Competitividad, Secretaría de Estado de Investigación, Desarrollo e Innovación, Spain (ref. ENE2014-60577-R).Peer ReviewedPostprint (author's final draft

Monocytes and neutrophils expressing myeloperoxidase occur in fibrous caps and thrombi in unstable coronary plaques

<p>Abstract</p> <p>Background</p> <p>Myeloperoxidase (MPO) -containing macrophages and neutrophils have been described at sites of plaque rupture. The presence of these cells in precursor lesions to acute rupture (thin cap atheroma, or vulnerable plaque) and within thrombi adjacent to ruptures has not been described, nor an association with iron-containing macrophages within unstable plaques.</p> <p>Methods</p> <p>We studied 61 acute ruptures, 15 organizing ruptures, 31 thin cap fibroatheromas, and 28 fibroatheromas from 72 sudden coronary death victims by immunohistochemical and histochemical techniques. Inflammatory cells were typed with anti-CD68 (macrophages), anti-BP-30 (neutrophil bactericidal glycoprotein), and anti-MPO. Iron was localized by Mallory's Prussian blue stain. In selected plaques alpha smooth muscle actin (DAKO, Carpinteria, CA, clone M0851) was performed.</p> <p>Results</p> <p>MPO positive cells were present in 79% of ruptured caps, 28% of thin cap fibroatheroma, and no fibroatheromas; neutrophils were present in 72% of ruptures, 8% of thin cap fibroatheromas, and no fibroatheromas. Iron containing foam cells were present in the caps of 93% of acute ruptures, of 85% of organizing ruptures, 20% of thin cap atheromas, and 10% of fibroatheromas. MPO positive cells were more frequent in occlusive than non-occlusive thrombi adjacent to ruptures (p = .006) and were more numerous in diabetics compared to non-diabetics (p = .002)</p> <p>Conclusion</p> <p>Unstable fibrous caps are more likely to contain MPO-positive cells, neutrophils, and iron-containing macrophages than fibrous caps of stable fibroatheromas. MPO-positive cells in thrombi adjacent to disrupted plaques are associated with occlusive thrombi and are more numerous in diabetic patients.</p

Times of Minima of Some Eclipsing Binaries

We present new times of minima in the light curves of some eclipsing binarie

Exome-Wide Association Study of Competitive Performance in Elite Athletes

The aim of the study was to identify genetic variants associated with personal best scores in Turkish track and field athletes and to compare allelic frequencies between sprint/power and endurance athletes and controls using a whole-exome sequencing (WES) approach, followed by replication studies in independent cohorts. The discovery phase involved 60 elite Turkish athletes (31 sprint/power and 29 endurance) and 20 ethnically matched controls. The replication phase involved 1132 individuals (115 elite Russian sprinters, 373 elite Russian endurance athletes (of which 75 athletes were with VO2max measurements), 209 controls, 148 Russian and 287 Finnish individuals with muscle fiber composition and cross-sectional area (CSA) data). None of the single nucleotide polymorphisms (SNPs) reached an exome-wide significance level (p < 2.3 × 10−7) in genotype–phenotype and case–control studies of Turkish athletes. However, of the 53 nominally (p < 0.05) associated SNPs, four functional variants were replicated. The SIRT1 rs41299232 G allele was significantly over-represented in Turkish (p = 0.047) and Russian (p = 0.018) endurance athletes compared to sprint/power athletes and was associated with increased VO2max (p = 0.037) and a greater proportion of slow-twitch muscle fibers (p = 0.035). The NUP210 rs2280084 A allele was significantly over-represented in Turkish (p = 0.044) and Russian (p = 0.012) endurance athletes compared to sprint/power athletes. The TRPM2 rs1785440 G allele was significantly over-represented in Turkish endurance athletes compared to sprint/power athletes (p = 0.034) and was associated with increased VO2max (p = 0.008). The AGRN rs4074992 C allele was significantly over-represented in Turkish sprint/power athletes compared to endurance athletes (p = 0.037) and was associated with a greater CSA of fast-twitch muscle fibers (p = 0.024). In conclusion, we present the first WES study of athletes showing that this approach can be used to identify novel genetic markers associated with exercise- and sport-related phenotypes

A global-scale screening of non-native aquatic organisms to identify potentially invasive species under current and future climate conditions

The threat posed by invasive non-native species worldwide requires a global approach to identify which introduced species are likely to pose an elevated risk of impact to native species and ecosystems. To inform policy, stakeholders and management decisions on global threats to aquatic ecosystems, 195 assessors representing 120 risk assessment areas across all six inhabited continents screened 819 non-native species from 15 groups of aquatic organisms (freshwater, brackish, marine plants and animals) using the Aquatic Species Invasiveness Screening Kit. This multi-lingual decision-support tool for the risk screening of aquatic organisms provides assessors with risk scores for a species under current and future climate change conditions that, following a statistically based calibration, permits the accurate classification of species into high-, medium-and low-risk categories under current and predicted climate conditions. The 1730 screenings undertaken encompassed wide geographical areas (regions, political entities, parts thereof, water bodies, river basins, lake drainage basins, and marine regions), which permitted thresholds to be identified for almost all aquatic organismal groups screened as well as for tropical, temperate and continental climate classes, and for tropical and temperate marine ecoregions. In total, 33 species were identified as posing a 'very high risk' of being or becoming invasive, and the scores of several of these species under current climate increased under future climate conditions, primarily due to their wide thermal tolerances. The risk thresholds determined for taxonomic groups and climate zones provide a basis against which area-specific or climate-based calibrated thresholds may be interpreted. In turn, the risk rankings help decision-makers identify which species require an immediate 'rapid' management action (e.g. eradication, control) to avoid or mitigate adverse impacts, which require a full risk assessment, and which are to be restricted or banned with regard to importation and/or sale as ornamental or aquarium/fishery enhancement. Decision support tools AS-ISK Hazard identification Non-native species Risk analysis Climate changepublishedVersio