BAsE-Seq: a method for obtaining long viral haplotypes from short sequence reads.

We present a method for obtaining long haplotypes, of over 3 kb in length, using a short-read sequencer, Barcode-directed Assembly for Extra-long Sequences (BAsE-Seq). BAsE-Seq relies on transposing a template-specific barcode onto random segments of the template molecule and assembling the barcoded short reads into complete haplotypes. We applied BAsE-Seq on mixed clones of hepatitis B virus and accurately identified haplotypes occurring at frequencies greater than or equal to 0.4%, with >99.9% specificity. Applying BAsE-Seq to a clinical sample, we obtained over 9,000 viral haplotypes, which provided an unprecedented view of hepatitis B virus population structure during chronic infection. BAsE-Seq is readily applicable for monitoring quasispecies evolution in viral diseases

Single-virion sequencing of lamivudine-treated HBV populations reveal population evolution dynamics and demographic history.

BACKGROUND: Viral populations are complex, dynamic, and fast evolving. The evolution of groups of closely related viruses in a competitive environment is termed quasispecies. To fully understand the role that quasispecies play in viral evolution, characterizing the trajectories of viral genotypes in an evolving population is the key. In particular, long-range haplotype information for thousands of individual viruses is critical; yet generating this information is non-trivial. Popular deep sequencing methods generate relatively short reads that do not preserve linkage information, while third generation sequencing methods have higher error rates that make detection of low frequency mutations a bioinformatics challenge. Here we applied BAsE-Seq, an Illumina-based single-virion sequencing technology, to eight samples from four chronic hepatitis B (CHB) patients - once before antiviral treatment and once after viral rebound due to resistance. RESULTS: With single-virion sequencing, we obtained 248-8796 single-virion sequences per sample, which allowed us to find evidence for both hard and soft selective sweeps. We were able to reconstruct population demographic history that was independently verified by clinically collected data. We further verified four of the samples independently through PacBio SMRT and Illumina Pooled deep sequencing. CONCLUSIONS: Overall, we showed that single-virion sequencing yields insight into viral evolution and population dynamics in an efficient and high throughput manner. We believe that single-virion sequencing is widely applicable to the study of viral evolution in the context of drug resistance and host adaptation, allows differentiation between soft or hard selective sweeps, and may be useful in the reconstruction of intra-host viral population demographic history

A knowledge-driven GIS modeling technique for groundwater potential mapping at the Upper Langat Basin, Malaysia.

The aim of this paper is to use a knowledge-driven expert-based geographical information system (GIS) model coupling with remote-sensing-derived parameters for groundwater potential mapping in an area of the Upper Langat Basin, Malaysia. In this study, nine groundwater storage controlling parameters that affect groundwater occurrences are derived from remotely sensed imagery, available maps, and associated databases. Those parameters are: lithology, slope, lineament, land use, soil, rainfall, drainage density, elevation, and geomorphology. Then the parameter layers were integrated and modeled using a knowledge-driven GIS of weighted linear combination. The weightage and score for each parameter and their classes are based on the Malaysian groundwater expert opinion survey. The predicted groundwater potential map was classified into four distinct zones based on the classification scheme designed by Department of Minerals and Geoscience Malaysia (JMG). The results showed that about 17% of the study area falls under low-potential zone, with 66% on moderate-potential zone, 15% with high-potential zone, and only 0.45% falls under very-high-potential zone. The results obtained in this study were validated with the groundwater borehole wells data compiled by the JMG and showed 76% of prediction accuracy. In addition statistical analysis indicated that hard rock dominant of the study area is controlled by secondary porosity such as distance from lineament and density of lineament. There are high correlations between area percentage of predicted groundwater potential zones and groundwater well yield. Results obtained from this study can be useful for future planning of groundwater exploration, planning and development by related agencies in Malaysia which provide a rapid method and reduce cost as well as less time consuming. The results may be also transferable to other areas of similar hydrological characteristics

Prediction of Protein Domain with mRMR Feature Selection and Analysis

The domains are the structural and functional units of proteins. With the avalanche of protein sequences generated in the postgenomic age, it is highly desired to develop effective methods for predicting the protein domains according to the sequences information alone, so as to facilitate the structure prediction of proteins and speed up their functional annotation. However, although many efforts have been made in this regard, prediction of protein domains from the sequence information still remains a challenging and elusive problem. Here, a new method was developed by combing the techniques of RF (random forest), mRMR (maximum relevance minimum redundancy), and IFS (incremental feature selection), as well as by incorporating the features of physicochemical and biochemical properties, sequence conservation, residual disorder, secondary structure, and solvent accessibility. The overall success rate achieved by the new method on an independent dataset was around 73%, which was about 28–40% higher than those by the existing method on the same benchmark dataset. Furthermore, it was revealed by an in-depth analysis that the features of evolution, codon diversity, electrostatic charge, and disorder played more important roles than the others in predicting protein domains, quite consistent with experimental observations. It is anticipated that the new method may become a high-throughput tool in annotating protein domains, or may, at the very least, play a complementary role to the existing domain prediction methods, and that the findings about the key features with high impacts to the domain prediction might provide useful insights or clues for further experimental investigations in this area. Finally, it has not escaped our notice that the current approach can also be utilized to study protein signal peptides, B-cell epitopes, HIV protease cleavage sites, among many other important topics in protein science and biomedicine

Impact of family structure on long-term survivors of osteosarcoma.

GOALS OF WORK: Long-term outcomes of osteosarcoma have dramatically improved with the use of modern combination therapies. Such aggressive treatments, however, entail chronic complications. In the present study, we assessed the functional, psychological, and familial status of long-term survivors of osteosarcoma treated at our institution. MATERIALS AND METHODS: Fifteen long-term survivors of osteosarcoma were evaluated for functional and psychological sequelae. Functional assessment was based on a method described by Enneking et al. Psychological assessment was based on General Health Questionnaire 28, Inventory Scale for Traumatic Neurosis, and Family System Test. MAIN RESULTS: Ten patients showed mild functional impairments; only five patients were handicapped more seriously. Depressive symptoms were diagnosed in four patients. A total of six patients revealed unbalanced family structures, including three of the four patients with depressive symptoms, all four patients with symptoms of posttraumatic stress disorder, and five of seven patients who showed poor emotional acceptance. CONCLUSIONS: Osteosarcoma survivors will generally recover good functional performance. Only a minority of them remain seriously impaired. One third of the patients present depressive symptoms and posttraumatic stress disorder. Poor coping is closely associated with unbalanced family structures. Therefore, the psychological and familial situation of patients with newly diagnosed osteosarcoma should be carefully assessed

Estimating global injuries morbidity and mortality: methods and data used in the Global Burden of Disease 2017 study

BACKGROUND: While there is a long history of measuring death and disability from injuries, modern research methods must account for the wide spectrum of disability that can occur in an injury, and must provide estimates with sufficient demographic, geographical and temporal detail to be useful for policy makers. The Global Burden of Disease (GBD) 2017 study used methods to provide highly detailed estimates of global injury burden that meet these criteria. METHODS: In this study, we report and discuss the methods used in GBD 2017 for injury morbidity and mortality burden estimation. In summary, these methods included estimating cause-specific mortality for every cause of injury, and then estimating incidence for every cause of injury. Non-fatal disability for each cause is then calculated based on the probabilities of suffering from different types of bodily injury experienced. RESULTS: GBD 2017 produced morbidity and mortality estimates for 38 causes of injury. Estimates were produced in terms of incidence, prevalence, years lived with disability, cause-specific mortality, years of life lost and disability-adjusted life-years for a 28-year period for 22 age groups, 195 countries and both sexes. CONCLUSIONS: GBD 2017 demonstrated a complex and sophisticated series of analytical steps using the largest known database of morbidity and mortality data on injuries. GBD 2017 results should be used to help inform injury prevention policy making and resource allocation. We also identify important avenues for improving injury burden estimation in the future

Global migration of influenza A viruses in swine

The complex and unresolved evolutionary origins of the 2009 H1N1 influenza pandemic exposed major gaps in our knowledge of the global spatial ecology and evolution of influenza A viruses in swine (swIAVs). Here we undertake an expansive phylogenetic analysis of swIAV sequence data and demonstrate that the global live swine trade strongly predicts the spatial dissemination of swIAVs, with Europe and North America acting as sources of viruses in Asian countries. In contrast, China has the world's largest swine population but is not a major exporter of live swine, and is not an important source of swIAVs in neighbouring Asian countries or globally. A meta-population simulation model incorporating trade data predicts that the global ecology of swIAVs is more complex than previously thought, and the United States and China's large swine populations are unlikely to be representative of swIAV diversity in their respective geographic regions, requiring independent surveillance efforts throughout Latin America and Asia.status: publishe