Optimized mixed Markov models for motif identification

BACKGROUND: Identifying functional elements, such as transcriptional factor binding sites, is a fundamental step in reconstructing gene regulatory networks and remains a challenging issue, largely due to limited availability of training samples. RESULTS: We introduce a novel and flexible model, the Optimized Mixture Markov model (OMiMa), and related methods to allow adjustment of model complexity for different motifs. In comparison with other leading methods, OMiMa can incorporate more than the NNSplice's pairwise dependencies; OMiMa avoids model over-fitting better than the Permuted Variable Length Markov Model (PVLMM); and OMiMa requires smaller training samples than the Maximum Entropy Model (MEM). Testing on both simulated and actual data (regulatory cis-elements and splice sites), we found OMiMa's performance superior to the other leading methods in terms of prediction accuracy, required size of training data or computational time. Our OMiMa system, to our knowledge, is the only motif finding tool that incorporates automatic selection of the best model. OMiMa is freely available at [1]. CONCLUSION: Our optimized mixture of Markov models represents an alternative to the existing methods for modeling dependent structures within a biological motif. Our model is conceptually simple and effective, and can improve prediction accuracy and/or computational speed over other leading methods

Korarchaeota Diversity, Biogeography, and Abundance in Yellowstone and Great Basin Hot Springs and Ecological Niche Modeling Based on Machine Learning

Over 100 hot spring sediment samples were collected from 28 sites in 12 areas/regions, while recording as many coincident geochemical properties as feasible (>60 analytes). PCR was used to screen samples for Korarchaeota 16S rRNA genes. Over 500 Korarchaeota 16S rRNA genes were screened by RFLP analysis and 90 were sequenced, resulting in identification of novel Korarchaeota phylotypes and exclusive geographical variants. Korarchaeota diversity was low, as in other terrestrial geothermal systems, suggesting a marine origin for Korarchaeota with subsequent niche-invasion into terrestrial systems. Korarchaeota endemism is consistent with endemism of other terrestrial thermophiles and supports the existence of dispersal barriers. Korarchaeota were found predominantly in >55°C springs at pH 4.7–8.5 at concentrations up to 6.6×106 16S rRNA gene copies g−1 wet sediment. In Yellowstone National Park (YNP), Korarchaeota were most abundant in springs with a pH range of 5.7 to 7.0. High sulfate concentrations suggest these fluids are influenced by contributions from hydrothermal vapors that may be neutralized to some extent by mixing with water from deep geothermal sources or meteoric water. In the Great Basin (GB), Korarchaeota were most abundant at spring sources of pH<7.2 with high particulate C content and high alkalinity, which are likely to be buffered by the carbonic acid system. It is therefore likely that at least two different geological mechanisms in YNP and GB springs create the neutral to mildly acidic pH that is optimal for Korarchaeota. A classification support vector machine (C-SVM) trained on single analytes, two analyte combinations, or vectors from non-metric multidimensional scaling models was able to predict springs as Korarchaeota-optimal or sub-optimal habitats with accuracies up to 95%. To our knowledge, this is the most extensive analysis of the geochemical habitat of any high-level microbial taxon and the first application of a C-SVM to microbial ecology

Historical nectar assessment reveals the fall and rise of floral resources in Britain

There is considerable concern over declines in insect pollinator communities and potential impacts on the pollination of crops and wildflowers. Among the multiple pressures facing pollinators, decreasing floral resources due to habitat loss and degradation has been suggested as a key contributing factor. However, a lack of quantitative data has hampered testing for historical changes in floral resources. Here we show that overall floral rewards can be estimated at a national scale by combining vegetation surveys and direct nectar measurements. We find evidence for substantial losses in nectar resources in England and Wales between the 1930s and 1970s; however, total nectar provision in Great Britain as a whole had stabilized by 1978, and increased from 1998 to 2007. These findings concur with trends in pollinator diversity, which declined in the mid-twentieth century but stabilized more recently. The diversity of nectar sources declined from 1978 to 1990 and thereafter in some habitats, with four plant species accounting for over 50% of national nectar provision in 2007. Calcareous grassland, broadleaved woodland and neutral grassland are the habitats that produce the greatest amount of nectar per unit area from the most diverse sources, whereas arable land is the poorest with respect to amount of nectar per unit area and diversity of nectar sources. Although agri-environment schemes add resources to arable landscapes, their national contribution is low. Owing to their large area, improved grasslands could add substantially to national nectar provision if they were managed to increase floral resource provision. This national-scale assessment of floral resource provision affords new insights into the links between plant and pollinator declines, and offers considerable opportunities for conservation

DEAD-Box Protein Ddx46 Is Required for the Development of the Digestive Organs and Brain in Zebrafish

Spatially and temporally controlled gene expression, including transcription, several mRNA processing steps, and the export of mature mRNA to the cytoplasm, is essential for developmental processes. It is well known that RNA helicases of the DExD/H-box protein family are involved in these gene expression processes, including transcription, pre-mRNA splicing, and rRNA biogenesis. Although one DExD/H-box protein, Prp5, a homologue of vertebrate Ddx46, has been shown to play important roles in pre-mRNA splicing in yeast, the in vivo function of Ddx46 remains to be fully elucidated in metazoans. In this study, we isolated zebrafish morendo (mor), a mutant that shows developmental defects in the digestive organs and brain, and found that it encodes Ddx46. The Ddx46 transcript is maternally supplied, and as development proceeds in zebrafish larvae, its ubiquitous expression gradually becomes restricted to those organs. The results of whole-mount in situ hybridization showed that the expression of various molecular markers in these organs is considerably reduced in the Ddx46 mutant. Furthermore, splicing status analysis with RT-PCR revealed unspliced forms of mRNAs in the digestive organ and brain tissues of the Ddx46 mutant, suggesting that Ddx46 may be required for pre-mRNA splicing during zebrafish development. Therefore, our results suggest a model in which zebrafish Ddx46 is required for the development of the digestive organs and brain, possibly through the control of pre-mRNA splicing

Genomes of the Most Dangerous Epidemic Bacteria Have a Virulence Repertoire Characterized by Fewer Genes but More Toxin-Antitoxin Modules

We conducted a comparative genomic study based on a neutral approach to identify genome specificities associated with the virulence capacity of pathogenic bacteria. We also determined whether virulence is dictated by rules, or if it is the result of individual evolutionary histories. We systematically compared the genomes of the 12 most dangerous pandemic bacteria for humans ("bad bugs") to their closest non-epidemic related species ("controls").We found several significantly different features in the "bad bugs", one of which was a smaller genome that likely resulted from a degraded recombination and repair system. The 10 Cluster of Orthologous Group (COG) functional categories revealed a significantly smaller number of genes in the "bad bugs", which lacked mostly transcription, signal transduction mechanisms, cell motility, energy production and conversion, and metabolic and regulatory functions. A few genes were identified as virulence factors, including secretion system proteins. Five "bad bugs" showed a greater number of poly (A) tails compared to the controls, whereas an elevated number of poly (A) tails was found to be strongly correlated to a low GC% content. The "bad bugs" had fewer tandem repeat sequences compared to controls. Moreover, the results obtained from a principal component analysis (PCA) showed that the "bad bugs" had surprisingly more toxin-antitoxin modules than did the controls.We conclude that pathogenic capacity is not the result of "virulence factors" but is the outcome of a virulent gene repertoire resulting from reduced genome repertoires. Toxin-antitoxin systems could participate in the virulence repertoire, but they may have developed independently of selfish evolution

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

Gravitational Waves and Gamma-Rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A

On 2017 August 17, the gravitational-wave event GW170817 was observed by the Advanced LIGO and Virgo detectors, and the gamma-ray burst (GRB) GRB 170817A was observed independently by the Fermi Gamma-ray Burst Monitor, and the Anti-Coincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory. The probability of the near-simultaneous temporal and spatial observation of GRB 170817A and GW170817 occurring by chance is $5.0\times {10}^{-8}$. We therefore confirm binary neutron star mergers as a progenitor of short GRBs. The association of GW170817 and GRB 170817A provides new insight into fundamental physics and the origin of short GRBs. We use the observed time delay of $(+1.74\pm 0.05)\,{\rm{s}}$ between GRB 170817A and GW170817 to: (i) constrain the difference between the speed of gravity and the speed of light to be between $-3\times {10}^{-15}$ and $+7\times {10}^{-16}$ times the speed of light, (ii) place new bounds on the violation of Lorentz invariance, (iii) present a new test of the equivalence principle by constraining the Shapiro delay between gravitational and electromagnetic radiation. We also use the time delay to constrain the size and bulk Lorentz factor of the region emitting the gamma-rays. GRB 170817A is the closest short GRB with a known distance, but is between 2 and 6 orders of magnitude less energetic than other bursts with measured redshift. A new generation of gamma-ray detectors, and subthreshold searches in existing detectors, will be essential to detect similar short bursts at greater distances. Finally, we predict a joint detection rate for the Fermi Gamma-ray Burst Monitor and the Advanced LIGO and Virgo detectors of 0.1-1.4 per year during the 2018-2019 observing run and 0.3-1.7 per year at design sensitivity

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline, and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams

Assessment of population genetic structure in the arbovirus vector midge, Culicoides brevitarsis (Diptera Ceratopogonidae), using multi-locus DNA microsatellites

Bluetongue virus (BTV) is a major pathogen of ruminants that is transmitted by biting midges (Culicoides spp.). Australian BTV serotypes have origins in Asia and are distributed across the continent into two distinct episystems, one in the north and another in the east. Culicoides brevitarsis is the major vector of BTV in Australia and is distributed across the entire geographic range of the virus. Here, we describe the isolation and use of DNA microsatellites and gauge their ability to determine population genetic connectivity of C. brevitarsis within Australia and with countries to the north. Eleven DNA microsatellite markers were isolated using a novel genomic enrichment method and identified as useful for genetic analyses of sampled populations in Australia, northern Papua New Guinea (PNG) and Timor-Leste. Significant (P < 0.05) population genetic subdivision was observed between all paired regions, though the highest levels of genetic sub-division involved pair-wise tests with PNG (PNG vs. Australia (F-ST = 0.120) and PNG vs. Timor-Leste (F-ST = 0.095)). Analysis of multi-locus allelic distributions using STRUCTURE identified a most probable two-cluster population model, which separated PNG specimens from a cluster containing specimens from Timor-Leste and Australia. The source of incursions of this species in Australia is more likely to be Timor-Leste than PNG. Future incursions of BTV positive C. brevitarsis into Australia may be genetically identified to their source populations using these microsatellite loci. The vector's panmictic genetic structure within Australia cannot explain the differential geographic distribution of BTV serotypes