Screening synteny blocks in pairwise genome comparisons through integer programming

<p>Abstract</p> <p>Background</p> <p>It is difficult to accurately interpret chromosomal correspondences such as true orthology and paralogy due to significant divergence of genomes from a common ancestor. Analyses are particularly problematic among lineages that have repeatedly experienced whole genome duplication (WGD) events. To compare multiple "subgenomes" derived from genome duplications, we need to relax the traditional requirements of "one-to-one" syntenic matchings of genomic regions in order to reflect "one-to-many" or more generally "many-to-many" matchings. However this relaxation may result in the identification of synteny blocks that are derived from ancient shared WGDs that are not of interest. For many downstream analyses, we need to eliminate weak, low scoring alignments from pairwise genome comparisons. Our goal is to objectively select subset of synteny blocks whose total scores are maximized while respecting the duplication history of the genomes in comparison. We call this "quota-based" screening of synteny blocks in order to appropriately fill a quota of syntenic relationships within one genome or between two genomes having WGD events.</p> <p>Results</p> <p>We have formulated the synteny block screening as an optimization problem known as "Binary Integer Programming" (BIP), which is solved using existing linear programming solvers. The computer program QUOTA-ALIGN performs this task by creating a clear objective function that maximizes the compatible set of synteny blocks under given constraints on overlaps and depths (corresponding to the duplication history in respective genomes). Such a procedure is useful for any pairwise synteny alignments, but is most useful in lineages affected by multiple WGDs, like plants or fish lineages. For example, there should be a 1:2 ploidy relationship between genome A and B if genome B had an independent WGD subsequent to the divergence of the two genomes. We show through simulations and real examples using plant genomes in the rosid superorder that the quota-based screening can eliminate ambiguous synteny blocks and focus on specific genomic evolutionary events, like the divergence of lineages (in cross-species comparisons) and the most recent WGD (in self comparisons).</p> <p>Conclusions</p> <p>The QUOTA-ALIGN algorithm screens a set of synteny blocks to retain only those compatible with a user specified ploidy relationship between two genomes. These blocks, in turn, may be used for additional downstream analyses such as identifying true orthologous regions in interspecific comparisons. There are two major contributions of QUOTA-ALIGN: 1) reducing the block screening task to a BIP problem, which is novel; 2) providing an efficient software pipeline starting from all-against-all BLAST to the screened synteny blocks with dot plot visualizations. Python codes and full documentations are publicly available <url>http://github.com/tanghaibao/quota-alignment</url>. QUOTA-ALIGN program is also integrated as a major component in SynMap <url>http://genomevolution.com/CoGe/SynMap.pl</url>, offering easier access to thousands of genomes for non-programmers.</p

Estimating EQ-5D utilities based on the Short-Form Long Term Conditions Questionnaire (LTCQ-8)

Purpose: The aim of this work was to develop a mapping algorithm for estimating EuroQoL 5 Dimension (EQ-5D) utilities from responses to the Long-Term Conditions Questionnaire (LTCQ), thus increasing LTCQ’s potential as a comprehensive outcome measure for evaluating integrated care initiatives. Methods: We combined data from three studies to give a total sample of 1334 responses. In each of the three datasets, we randomly selected 75% of the sample and combined the selected random samples to generate the estimation dataset, which consisted of 1001 patients. The unselected 25% observations from each dataset were combined to generate an internal validation dataset of 333 patients. We used direct mapping models by regressing responses to the LTCQ-8 directly onto EQ-5D-5L and EQ-5D-3L utilities as well as response (or indirect) mapping to predict the response level that patients selected for each of the five EQ-5D-5L domains. Several models were proposed and compared on mean squared error and mean absolute error. Results: A two-part model with OLS was the best performing based on the mean squared error (0.038) and mean absolute error (0.147) when estimating the EQ-5D-5L utilities. A multinomial response mapping model using LTCQ-8 responses was used to predict EQ-5D-5L responses levels. Conclusions: This study provides a mapping algorithm for estimating EQ-5D utilities from LTCQ responses. The results from this study can help broaden the applicability of the LTCQ by producing utility values for use in economic analyses

Economic Impacts of Non-Native Forest Insects in the Continental United States

Reliable estimates of the impacts and costs of biological invasions are critical to developing credible management, trade and regulatory policies. Worldwide, forests and urban trees provide important ecosystem services as well as economic and social benefits, but are threatened by non-native insects. More than 450 non-native forest insects are established in the United States but estimates of broad-scale economic impacts associated with these species are largely unavailable. We developed a novel modeling approach that maximizes the use of available data, accounts for multiple sources of uncertainty, and provides cost estimates for three major feeding guilds of non-native forest insects. For each guild, we calculated the economic damages for five cost categories and we estimated the probability of future introductions of damaging pests. We found that costs are largely borne by homeowners and municipal governments. Wood- and phloem-boring insects are anticipated to cause the largest economic impacts by annually inducing nearly $1.7 billion in local government expenditures and approximately$830 million in lost residential property values. Given observations of new species, there is a 32% chance that another highly destructive borer species will invade the U.S. in the next 10 years. Our damage estimates provide a crucial but previously missing component of cost-benefit analyses to evaluate policies and management options intended to reduce species introductions. The modeling approach we developed is highly flexible and could be similarly employed to estimate damages in other countries or natural resource sectors

Comparison of Influenza and SIV Specific CD8 T Cell Responses in Macaques

Macaques are a potentially useful non-human primate model to compare memory T-cell immunity to acute virus pathogens such as influenza virus and effector T-cell responses to chronic viral pathogens such as SIV. However, immunological reagents to study influenza CD8+ T-cell responses in the macaque model are limited. We recently developed an influenza-SIV vaccination model of pigtail macaques (Macaca nemestrina) and used this to study both influenza-specific and SIV-specific CD8+ T-cells in 39 pigtail macaques expressing the common Mane-A*10+ (Mane-A01*084) MHC-I allele. To perform comparative studies between influenza and SIV responses a common influenza nucleoprotein-specific CD8+ T-cell response was mapped to a minimal epitope (termed RA9), MHC-restricted to Mane-A*10 and an MHC tetramer developed to study this response. Influenza-specific memory CD8+ T-cell response maintained a highly functional profile in terms of multitude of effector molecule expression (CD107a, IFN-γ, TNF-α, MIP-1β and IL-2) and showed high avidity even in the setting of SIV infection. In contrast, within weeks following active SIV infection, SIV-specific CD8+ effector T-cells expressed fewer cytokines/degranulation markers and had a lower avidity compared to influenza specific CD8+ T-cells. Further, the influenza specific memory CD8 T-cell response retained stable expression of the exhaustion marker programmed death-marker-1 (PD-1) and co-stimulatory molecule CD28 following infection with SIV. This contrasted with the effector SIV-specific CD8+ T-cells following SIV infection which expressed significantly higher amounts of PD-1 and lower amounts of CD28. Our results suggest that strategies to maintain a more functional CD8+ T-cell response, profile may assist in controlling HIV disease

Individuals with chronic low back pain have greater difficulty in engaging in positive lifestyle behaviours than those without back pain: An assessment of health literacy

Background: Despite the large volume of research dedicated to understanding chronic low back pain (CLBP), patient outcomes remain modest while healthcare costs continue to rise, creating a major public health burden. Health literacy - the ability to seek, understand and utilise health information - has been identified as an important factor in the course of other chronic conditions and may be important in the aetiology of CLBP. Many of the currently available health literacy measurement tools are limited since they measure narrow aspects of health literacy. The Health Literacy Measurement Scale (HeLMS) was developed recently to measure broader elements of health literacy. The aim of this study was to measure broad elements of health literacy among individuals with CLBP and without LBP using the HeLMS.Methods: Thirty-six community-dwelling adults with CLBP and 44 with no history of LBP responded to the HeLMS. Individuals were recruited as part of a larger community-based spinal health study in Western Australia. Scores for the eight domains of the HeLMS as well as individual item responses were compared between the groups.Results: HeLMS scores were similar between individuals with and without CLBP for seven of the eight health literacy domains (p &gt; 0.05). However, compared to individuals with no history of LBP, those with CLBP had a significantly lower score in the domain &lsquo;Patient attitudes towards their health&rsquo; (mean difference [95% CI]: 0.46 [0.11- 0.82]) and significantly lower scores for each of the individual items within this domain (p &lt; 0.05). Moderate effect sizes ranged from d = 0.47-0.65.Conclusions: Although no differences were identified in HeLMS scores between the groups for seven of the health literacy domains, adults with CLBP reported greater difficulty in engaging in general positive health behaviours. This aspect of health literacy suggests that self-management support initiatives may benefit individuals with CLBP.<br /

Integrated transcriptomic and proteomic analysis of the global response of Wolbachia to doxycycline-induced stress

The bacterium Wolbachia (order Rickettsiales), representing perhaps the most abundant vertically transmitted microbe worldwide, infects arthropods and filarial nematodes. In arthropods, Wolbachia can induce reproductive alterations and interfere with the transmission of several arthropod-borne pathogens. In addition, Wolbachia is an obligate mutualist of the filarial parasites that cause lymphatic filariasis and onchocerciasis in the tropics. Targeting Wolbachia with tetracycline antibiotics leads to sterilisation and ultimately death of adult filariae. However, several weeks of treatment are required, restricting the implementation of this control strategy. To date, the response of Wolbachia to stress has not been investigated, and almost nothing is known about global regulation of gene expression in this organism. We exposed an arthropod Wolbachia strain to doxycycline in vitro, and analysed differential expression by directional RNA-seq and label-free, quantitative proteomics. We found that Wolbachia responded not only by modulating expression of the translation machinery, but also by upregulating nucleotide synthesis and energy metabolism, while downregulating outer membrane proteins. Moreover, Wolbachia increased the expression of a key component of the twin-arginine translocase (tatA) and a phosphate ABC transporter ATPase (PstB); the latter is associated with decreased susceptibility to antimicrobials in free-living bacteria. Finally, the downregulation of 6S RNA during translational inhibition suggests that this small RNA is involved in growth rate control. Despite its highly reduced genome, Wolbachia shows a surprising ability to regulate gene expression during exposure to a potent stressor. Our findings have general relevance for the chemotherapy of obligate intracellular bacteria and the mechanistic basis of persistence in the Rickettsiales

Evolutionary genomics of a cold-adapted diatom: Fragilariopsis cylindrus

The Southern Ocean houses a diverse and productive community of organisms1, 2. Unicellular eukaryotic diatoms are the main primary producers in this environment, where photosynthesis is limited by low concentrations of dissolved iron and large seasonal fluctuations in light, temperature and the extent of sea ice3, 4, 5, 6, 7. How diatoms have adapted to this extreme environment is largely unknown. Here we present insights into the genome evolution of a cold-adapted diatom from the Southern Ocean, Fragilariopsis cylindrus8, 9, based on a comparison with temperate diatoms. We find that approximately 24.7 per cent of the diploid F. cylindrus genome consists of genetic loci with alleles that are highly divergent (15.1 megabases of the total genome size of 61.1 megabases). These divergent alleles were differentially expressed across environmental conditions, including darkness, low iron, freezing, elevated temperature and increased CO2. Alleles with the largest ratio of non-synonymous to synonymous nucleotide substitutions also show the most pronounced condition-dependent expression, suggesting a correlation between diversifying selection and allelic differentiation. Divergent alleles may be involved in adaptation to environmental fluctuations in the Southern Ocean

Six RNA Viruses and Forty-One Hosts: Viral Small RNAs and Modulation of Small RNA Repertoires in Vertebrate and Invertebrate Systems

We have used multiplexed high-throughput sequencing to characterize changes in small RNA populations that occur during viral infection in animal cells. Small RNA-based mechanisms such as RNA interference (RNAi) have been shown in plant and invertebrate systems to play a key role in host responses to viral infection. Although homologs of the key RNAi effector pathways are present in mammalian cells, and can launch an RNAi-mediated degradation of experimentally targeted mRNAs, any role for such responses in mammalian host-virus interactions remains to be characterized. Six different viruses were examined in 41 experimentally susceptible and resistant host systems. We identified virus-derived small RNAs (vsRNAs) from all six viruses, with total abundance varying from “vanishingly rare” (less than 0.1% of cellular small RNA) to highly abundant (comparable to abundant micro-RNAs “miRNAs”). In addition to the appearance of vsRNAs during infection, we saw a number of specific changes in host miRNA profiles. For several infection models investigated in more detail, the RNAi and Interferon pathways modulated the abundance of vsRNAs. We also found evidence for populations of vsRNAs that exist as duplexed siRNAs with zero to three nucleotide 3′ overhangs. Using populations of cells carrying a Hepatitis C replicon, we observed strand-selective loading of siRNAs onto Argonaute complexes. These experiments define vsRNAs as one possible component of the interplay between animal viruses and their hosts

Формирование эмоциональной культуры как компонента инновационной культуры студентов

Homozygosity has long been associated with rare, often devastating, Mendelian disorders1 and Darwin was one of the first to recognise that inbreeding reduces evolutionary fitness2. However, the effect of the more distant parental relatedness common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity, ROH), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power3,4. Here we use ROH to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts and find statistically significant associations between summed runs of homozygosity (SROH) and four complex traits: height, forced expiratory lung volume in 1 second (FEV1), general cognitive ability (g) and educational attainment (nominal p<1 × 10−300, 2.1 × 10−6, 2.5 × 10−10, 1.8 × 10−10). In each case increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing convincing evidence for the first time that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples5,6, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein (LDL) cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection7, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been