MEME Suite: tools for motif discovery and searching

The MEME Suite web server provides a unified portal for online discovery and analysis of sequence motifs representing features such as DNA binding sites and protein interaction domains. The popular MEME motif discovery algorithm is now complemented by the GLAM2 algorithm which allows discovery of motifs containing gaps. Three sequence scanning algorithms—MAST, FIMO and GLAM2SCAN—allow scanning numerous DNA and protein sequence databases for motifs discovered by MEME and GLAM2. Transcription factor motifs (including those discovered using MEME) can be compared with motifs in many popular motif databases using the motif database scanning algorithm Tomtom. Transcription factor motifs can be further analyzed for putative function by association with Gene Ontology (GO) terms using the motif-GO term association tool GOMO. MEME output now contains sequence LOGOS for each discovered motif, as well as buttons to allow motifs to be conveniently submitted to the sequence and motif database scanning algorithms (MAST, FIMO and Tomtom), or to GOMO, for further analysis. GLAM2 output similarly contains buttons for further analysis using GLAM2SCAN and for rerunning GLAM2 with different parameters. All of the motif-based tools are now implemented as web services via Opal. Source code, binaries and a web server are freely available for noncommercial use at http://meme.nbcr.net

Postglacial expansion of the arctic keystone copepod calanus glacialis

Calanus glacialis, a major contributor to zooplankton biomass in the Arctic shelf seas, is a key link between primary production and higher trophic levels that may be sensitive to climate warming. The aim of this study was to explore genetic variation in contemporary populations of this species to infer possible changes during the Quaternary period, and to assess its population structure in both space and time. Calanus glacialis was sampled in the fjords of Spitsbergen (Hornsund and Kongsfjorden) in 2003, 2004, 2006, 2009 and 2012. The sequence of a mitochondrial marker, belonging to the ND5 gene, selected for the study was 1249 base pairs long and distinguished 75 unique haplotypes among 140 individuals that formed three main clades. There was no detectable pattern in the distribution of haplotypes by geographic distance or over time. Interestingly, a Bayesian skyline plot suggested that a 1000-fold increase in population size occurred approximately 10,000 years before present, suggesting a species expansion after the Last Glacial Maximum.GAME from the National Science Centre, the Polish Ministry of Science and Higher Education Iuventus Plus [IP2014 050573]; FCT-PT [CCMAR/Multi/04326/2013]; [2011/03/B/NZ8/02876

First insights into the phylogenetic diversity of Mycobacterium tuberculosis in Nepal

BACKGROUND: Tuberculosis (TB) is a major public health problem in Nepal. Strain variation in Mycobacterium tuberculosis may influence the outcome of TB infection and disease. To date, the phylogenetic diversity of M. tuberculosis in Nepal is unknown. METHODS AND FINDINGS: We analyzed 261 M. tuberculosis isolates recovered from pulmonary TB patients recruited between August 2009 and August 2010 in Nepal. M. tuberculosis lineages were determined by single nucleotide polymorphisms (SNP) typing and spoligotyping. Drug resistance was determined by sequencing the hot spot regions of the relevant target genes. Overall, 164 (62.8%) TB patients were new, and 97 (37.2%) were previously treated. Any drug resistance was detected in 50 (19.2%) isolates, and 16 (6.1%) were multidrug-resistant. The most frequent M. tuberculosis lineage was Lineage 3 (CAS/Delhi) with 106 isolates (40.6%), followed by Lineage 2 (East-Asian lineage, includes Beijing genotype) with 84 isolates (32.2%), Lineage 4 (Euro-American lineage) with 41 (15.7%) isolates, and Lineage 1 (Indo-Oceanic lineage) with 30 isolates (11.5%). Based on spoligotyping, we found 45 different spoligotyping patterns that were previously described. The Beijing (83 isolates, 31.8%) and CAS spoligotype (52, 19.9%) were the dominant spoligotypes. A total of 36 (13.8%) isolates could not be assigned to any known spoligotyping pattern. Lineage 2 was associated with female sex (adjusted odds ratio [aOR] 2.58, 95% confidence interval [95% CI] 1.42-4.67, p = 0.002), and any drug resistance (aOR 2.79; 95% CI 1.43-5.45; p = 0.002). We found no evidence for an association of Lineage 2 with age or BCG vaccination status. CONCLUSIONS: We found a large genetic diversity of M. tuberculosis in Nepal with representation of all four major lineages. Lineages 3 and 2 were dominating. Lineage 2 was associated with clinical characteristics. This study fills an important gap on the map of the M. tuberculosis genetic diversity in the Asian reg

Revival of the magnetar PSR J1622-4950: observations with MeerKAT, Parkes, XMM-Newton, Swift, Chandra, and NuSTAR

New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x larger than during quiescence, and has been decaying exponentially on a 111+/-19 day timescale. This high-flux state, together with a radio-derived rotational ephemeris, enabled for the first time the detection of X-ray pulsations for this magnetar. At 5%, the 0.3-6 keV pulsed fraction is comparable to the smallest observed for magnetars. The overall pulsar geometry inferred from polarized radio emission appears to be broadly consistent with that determined 6-8 years earlier. However, rotating vector model fits suggest that we are now seeing radio emission from a different location in the magnetosphere than previously. This indicates a novel way in which radio emission from magnetars can differ from that of ordinary pulsars. The torque on the neutron star is varying rapidly and unsteadily, as is common for magnetars following outburst, having changed by a factor of 7 within six months of reactivation.Comment: Published in ApJ (2018 April 5); 13 pages, 4 figure

Spatio-temporal Models of Lymphangiogenesis in Wound Healing

Several studies suggest that one possible cause of impaired wound healing is failed or insufficient lymphangiogenesis, that is the formation of new lymphatic capillaries. Although many mathematical models have been developed to describe the formation of blood capillaries (angiogenesis), very few have been proposed for the regeneration of the lymphatic network. Lymphangiogenesis is a markedly different process from angiogenesis, occurring at different times and in response to different chemical stimuli. Two main hypotheses have been proposed: 1) lymphatic capillaries sprout from existing interrupted ones at the edge of the wound in analogy to the blood angiogenesis case; 2) lymphatic endothelial cells first pool in the wound region following the lymph flow and then, once sufficiently populated, start to form a network. Here we present two PDE models describing lymphangiogenesis according to these two different hypotheses. Further, we include the effect of advection due to interstitial flow and lymph flow coming from open capillaries. The variables represent different cell densities and growth factor concentrations, and where possible the parameters are estimated from biological data. The models are then solved numerically and the results are compared with the available biological literature.Comment: 29 pages, 9 Figures, 6 Tables (39 figure files in total

FONZIE: An optimized pipeline for minisatellite marker discovery and primer design from large sequence data sets

<p>Abstract</p> <p>Background</p> <p>Micro-and minisatellites are among the most powerful genetic markers known to date. They have been used as tools for a large number of applications ranging from gene mapping to phylogenetic studies and isolate typing. However, identifying micro-and minisatellite markers on large sequence data sets is often a laborious process.</p> <p>Results</p> <p>FONZIE was designed to successively 1) perform a search for markers via the external software Tandem Repeat Finder, 2) exclude user-defined specific genomic regions, 3) screen for the size and the percent matches of each relevant marker found by Tandem Repeat Finder, 4) evaluate marker specificity (i.e., occurrence of the marker as a single copy in the genome) using BLAST2.0, 5) design minisatellite primer pairs via the external software Primer3, and 6) check the specificity of each final PCR product by BLAST. A final file returns to users all the results required to amplify markers. A biological validation of the approach was performed using the whole genome sequence of the phytopathogenic fungus <it>Leptosphaeria maculans</it>, showing that more than 90% of the minisatellite primer pairs generated by the pipeline amplified a PCR product, 44.8% of which showed agarose-gel resolvable polymorphism between isolates. Segregation analyses confirmed that the polymorphic minisatellites corresponded to single-locus markers.</p> <p>Conclusion</p> <p>FONZIE is a stand-alone and user-friendly application developed to minimize tedious manual operations, reduce errors, and speed up the search for efficient minisatellite and microsatellite markers departing from whole-genome sequence data. This pipeline facilitates the integration of data and provides a set of specific primer sequences for PCR amplification of single-locus markers. FONZIE is freely downloadable at: <url>http://www.versailles-grignon.inra.fr/bioger/equipes/leptosphaeria_maculans/outils_d_analyses/fonzie</url></p

Two Distinct Triatoma dimidiata (Latreille, 1811) Taxa Are Found in Sympatry in Guatemala and Mexico

Approximately 10 million people are infected with Trypanosoma cruzi, the causative agent of Chagas disease, which remains the most serious parasitic disease in the Americas. Most people are infected via triatomine vectors. Transmission has been largely halted in South America in areas with predominantly domestic vectors. However, one of the main Chagas vectors in Mesoamerica, Triatoma dimidiata, poses special challenges to control due to its diversity across its large geographic range (from Mexico into northern South America), and peridomestic and sylvatic populations that repopulate houses following pesticide treatment. Recent evidence suggests T. dimidiata may be a complex of species, perhaps including cryptic species; taxonomic ambiguity which confounds control. The nuclear sequence of the internal transcribed spacer 2 (ITS2) of the ribosomal DNA and the mitochondrial cytochrome b (mt cyt b) gene were used to analyze the taxonomy of T. dimidiata from southern Mexico throughout Central America. ITS2 sequence divides T. dimidiata into four taxa. The first three are found mostly localized to specific geographic regions with some overlap: (1) southern Mexico and Guatemala (Group 2); (2) Guatemala, Honduras, El Salvador, Nicaragua, and Costa Rica (Group 1A); (3) and Panama (Group 1B). We extend ITS2 Group 1A south into Costa Rica, Group 2 into southern Guatemala and show the first information on isolates in Belize, identifying Groups 2 and 3 in that country. The fourth group (Group 3), a potential cryptic species, is dispersed across parts of Mexico, Guatemala, and Belize. We show it exists in sympatry with other groups in Peten, Guatemala, and Yucatan, Mexico. Mitochondrial cyt b data supports this putative cryptic species in sympatry with others. However, unlike the clear distinction of the remaining groups by ITS2, the remaining groups are not separated by mt cyt b. This work contributes to an understanding of the taxonomy and population subdivision of T. dimidiata, essential for designing effective control strategies

High level of conservation between genes coding for the GAMYB transcription factor in barley (Hordeum vulgare L.) and bread wheat (Triticum aestivum L.) collections

The transcription factor GAMYB is involved in gibberellin signalling in cereal aleurone cells and in plant developmental processes. Nucleotide diversity of HvGAMYB and TaGAMYB was investigated in 155 barley (Hordeum vulgare) and 42 wheat (Triticum aestivum) accessions, respectively. Polymorphisms defined 18 haplotypes in the barley collection and 1, 7 and 3 haplotypes for the A, B, and D genomes of wheat, respectively. We found that (1) Hv- and TaGAMYB genes have identical structures. (2) Both genes show a high level of nucleotide identity (>95%) in the coding sequences and the distribution of polymorphisms is similar in both collections. At the protein level the functional domain is identical in both species. (3) GAMYB genes map to a syntenic position on chromosome 3. GAMYB genes are different in both collections with respect to the Tajima D statistic and linkage disequilibrium (LD). A moderate level of LD was observed in the barley collection. In wheat, LD is absolute between polymorphic sites, mostly located in the first intron, while it decays within the gene. Differences in Tajima D values might be due to a lower selection pressure on HvGAMYB, compared to its wheat orthologue. Altogether our results provide evidence that there have been only few evolutionary changes in Hv- and TaGAMYB. This confirms the close relationship between these species and also highlights the functional importance of this transcription factor

Identification of a Mutation Associated with Fatal Foal Immunodeficiency Syndrome in the Fell and Dales Pony

The Fell and Dales are rare native UK pony breeds at risk due to falling numbers, in-breeding, and inherited disease. Specifically, the lethal Mendelian recessive disease Foal Immunodeficiency Syndrome (FIS), which manifests as B-lymphocyte immunodeficiency and progressive anemia, is a substantial threat. A significant percentage (∼10%) of the Fell ponies born each year dies from FIS, compromising the long-term survival of this breed. Moreover, the likely spread of FIS into other breeds is of major concern. Indeed, FIS was identified in the Dales pony, a related breed, during the course of this work. Using a stepwise approach comprising linkage and homozygosity mapping followed by haplotype analysis, we mapped the mutation using 14 FIS–affected, 17 obligate carriers, and 10 adults of unknown carrier status to a ∼1 Mb region (29.8 – 30.8 Mb) on chromosome (ECA) 26. A subsequent genome-wide association study identified two SNPs on ECA26 that showed genome-wide significance after Bonferroni correction for multiple testing: BIEC2-692674 at 29.804 Mb and BIEC2-693138 at 32.19 Mb. The associated region spanned 2.6 Mb from ∼29.6 Mb to 32.2 Mb on ECA26. Re-sequencing of this region identified a mutation in the sodium/myo-inositol cotransporter gene (SLC5A3); this causes a P446L substitution in the protein. This gene plays a crucial role in the regulatory response to osmotic stress that is essential in many tissues including lymphoid tissues and during early embryonic development. We propose that the amino acid substitution we identify here alters the function of SLC5A3, leading to erythropoiesis failure and compromise of the immune system. FIS is of significant biological interest as it is unique and is caused by a gene not previously associated with a mammalian disease. Having identified the associated gene, we are now able to eradicate FIS from equine populations by informed selective breeding

Interaction of plant growth regulators and reactive oxygen species to regulate petal senescence in wallflowers (Erysimum linifolium)

Background In many species floral senescence is coordinated by ethylene. Endogenous levels rise, and exogenous application accelerates senescence. Furthermore, floral senescence is often associated with increased reactive oxygen species, and is delayed by exogenously applied cytokinin. However, how these processes are linked remains largely unresolved. Erysimum linifolium (wallflower) provides an excellent model for understanding these interactions due to its easily staged flowers and close taxonomic relationship to Arabidopsis. This has facilitated microarray analysis of gene expression during petal senescence and provided gene markers for following the effects of treatments on different regulatory pathways. Results In detached Erysimum linifolium (wallflower) flowers ethylene production peaks in open flowers. Furthermore senescence is delayed by treatments with the ethylene signalling inhibitor silver thiosulphate, and accelerated with ethylene released by 2-chloroethylphosphonic acid. Both treatments with exogenous cytokinin, or 6-methyl purine (which is an inhibitor of cytokinin oxidase), delay petal senescence. However, treatment with cytokinin also increases ethylene biosynthesis. Despite the similar effects on senescence, transcript abundance of gene markers is affected differentially by the treatments. A significant rise in transcript abundance of WLS73 (a putative aminocyclopropanecarboxylate oxidase) was abolished by cytokinin or 6-methyl purine treatments. In contrast, WFSAG12 transcript (a senescence marker) continued to accumulate significantly, albeit at a reduced rate. Silver thiosulphate suppressed the increase in transcript abundance both of WFSAG12 and WLS73. Activity of reactive oxygen species scavenging enzymes changed during senescence. Treatments that increased cytokinin levels, or inhibited ethylene action, reduced accumulation of hydrogen peroxide. Furthermore, although auxin levels rose with senescence, treatments that delayed early senescence did not affect transcript abundance of WPS46, an auxin-induced gene. Conclusions A model for the interaction between cytokinins, ethylene, reactive oxygen species and auxin in the regulation of floral senescence in wallflowers is proposed. The combined increase in ethylene and reduction in cytokinin triggers the initiation of senescence and these two plant growth regulators directly or indirectly result in increased reactive oxygen species levels. A fall in conjugated auxin and/or the total auxin pool eventually triggers abscission