Utilization of Probabilistic Models in Short Read Assembly from Second-Generation Sequencing

With the advent of cheaper and faster DNA sequencing technologies, assembly methods have greatly changed. Instead of outputting reads that are thousands of base pairs long, new sequencers parallelize the task by producing read lengths between 35 and 400 base pairs. Reconstructing an organism’s genome from these millions of reads is a computationally expensive task. Our algorithm solves this problem by organizing and indexing the reads using n-grams, which are short, fixed-length DNA sequences of length n. These n-grams are used to efficiently locate putative read joins, thereby eliminating the need to perform an exhaustive search over all possible read pairs. Our goal was develop a novel n-gram method for the assembly of genomes from next-generation sequencers. Specifically, a probabilistic, iterative approach was utilized to determine the most likely reads to join through development of a new metric that models the probability of any two arbitrary reads being joined together. Tests were run using simulated short read data based on randomly created genomes ranging in lengths from 10,000 to 100,000 nucleotides with 16 to 20x coverage. We were able to successfully re-assemble entire genomes up to 100,000 nucleotides in length.

Towards a Variety of Family Models: Economic Contributions of Spouses and Family Power

In Western societies the increase in female employment (especially among married women) is seen as having brought about the crisis of the traditional model of the family, reinforcing the position of the "modern" model - the egalitarian family with two working spouses and a "dual-career" family. In contrast, the transitional situation in the post-communist countries during the 1990s is producing a crisis of the family with two working spouses (the basic type of the communist period) and leading to new power relations within the family. While the growth of dual-earner households in this century has implied modification of family models towards greater symmetry of responsibility for breadwinning and homemaking, there is considerable evidence that women's increased employment does not necessarily lead to a more egalitarian approach to gender roles within the family. The group set out to investigate the economic situation of families and economic power within the fame as a crucial factor in the transformation of families with two working spouses in order to reveal the specific patterns of gender contracts and power relations within the family that are emerging in response to the current political and economic transformation. They opted for a comparative approach, selecting the Czech Republic as a country where the very similar tendencies of a few years ago (almost 100% of women employed and the family as a realm of considerable private freedom where both women's and men's gender identities and the traditional distribution of family responsibilities were largely preserved) are combined with a very different experience in terms of economic inequalities during the 1990s to that of Russia. In the first stage of the study they surveyed 300 married couples (150 in each country) on the question of breadwinning. They then carried out in-depth interviews with 10 couples from each country (selected from among the educated layers of the population), focusing on the process of the social construction of gender, using breadwinning and homemaking as gender boundaries which distinguish men from women. By analysing changes in social position and the type of interpersonal interaction of spouses they distinguished two main types of family contracts: the neo-traditional "communal sharing" (with male breadwinner, traditional distribution of family chores and negotiated family power) and the modern one based on negotiated agreement. The most important pre-conditions of husband-wife agreement about breadwinning seemed to imply their overall gender ideology rather than the economic and/or family circumstances. In general, wives were more likely to express egalitarian views, supporting the blurring or even elimination of many gender boundaries. Husbands, on the other hand, more often gave responses calling for the continued maintenance of gender boundaries. The analysis showed that breadwinning is still an important gender boundary in these cultures, one that is assumed unless it is explicitly questioned and that is seen as part of what makes a man a "real man". The majority of respondents seemed to be committed to egalitarian ideology on gender roles and the distribution of family tasks, including decision making, but this is contradicted by the persistent idea of the husband as the breadwinner. This contradiction is more characteristic of the Russian situation than of the Czech. The quantitative study showed a difference in prevailing family models between the two countries, with a clearer shift towards the traditional family contract in the Russian case. The Czechs were more likely to consider their partnerships as based on negotiated agreement, while the Russians saw theirs as based on egalitarian contract, in both cases seeing this as the norm. The majority of couples said they felt satisfied with their marriage, although in both countries wives seemed to be less satisfied. There was however a difference in the issues that aroused dissatisfaction, with Czech women being more sensitive to issues such as self-realisation, personal independence, understanding and recognition in the family, and Russians to issues of love, understanding and recognition. The most disputed area for the majority of families was chores in the home, presumably because in many families both husband and wife were working hard outside the home and because a number of partners had differing views as to the ideal distribution of chores within the family. The distribution of power in the family seems to be linked to the level of well being. The analysis showed that in the dominant democratic model there is still an inverse connection between family leadership and well being: the more prominent the wife's position as head of the family is, the lower the level of family income. This may reflect both the husband's refusal to play the leading role in the family and even his rejection of any involvement in family issues in such a family. The qualitative research revealed that both men and women see the breadwinning role to be an essential part of masculine identity, a role which the female partner would take on temporarily to assist the male but not permanently since this would threaten the gender boundaries and the man's identity. At the same time, few breadwinners expressed a sense of job satisfaction and all considered their choice as imposed on them by the circumstances (i.e. having a family in difficult times). The group feel that family orientation and some loss of personal involvement in their profession is partly reflected in the fact that many of the men felt more comfortable and self-confident at home than at work. Women's work, on the other hand, was largely seen as a source of personal and self-realisation and social life. Eight out of ten of the Russian women interviewed were employed, although only two on a full-time basis, but none saw their jobs as adding substantially to the family budget. Both partners see the most important factor as the wife's wish to work or stay at home, and do not think it wise for the wife to work at the expense of her part of the "family contract", although husbands from the "egalitarian" relationships expressed more willingness to compromise. The analysis showed clearly that wives and husbands did not construct gender boundaries in isolation, with the interviews providing clear evidence of negotiation. At the same time, husbands' interpretations of their wives' employment were less susceptible to the influence of negotiation than were their gender attitudes and norms about breadwinning. One of the most interesting aspects of the spouses' negotiations was the extent to which they disagreed about what they seemed to have agreed upon. Most disagreements about the breadwinning boundaries, however, were over norms and were settled by changes in norms rather than in behavioural interpretation. Changes in norms were often a form of peace offering or were in response in changes in circumstances. The study did show, however, that many of the efforts at cooperation and compensation were more symbolic than real and the group found the plasticity of expressed gender ideology to be one of the most striking findings of their work. They conclude that the shift towards more traditional gednder distributions of incomes and domestic chores does not automatically mean the reestablishment of a patriarchal model of family power. On the contrary, it seems to be a compromise formation, relatively unstable, temporary and containing self-defeating forces as the split between the personal and professional value of work and its social value expressed in a money equivalent cannot be maintained for generations.

Observation and identification of metabolites emerging during postmortem decomposition of brain tissue by means of in situ 1H-magnetic resonance spectroscopy

Postmortem decomposition of brain tissue was investigated by (1)H-magnetic resonance spectroscopy (MRS) in a sheep head model and selected human cases. Aiming at the eventual estimation of postmortem intervals in forensic medicine, this study focuses on the characterization and identification of newly observed metabolites. In situ single-voxel (1)H-MRS at 1.5 T was complemented by multidimensional homo- and heteronuclear high-resolution NMR spectroscopy of an extract of sheep brain tissue. The inclusion of spectra of model solutions in the program LC Model confirmed the assignments in situ. The first postmortem phase was characterized mainly by changes in the concentrations of metabolites usually observed in vivo and by the appearance of previously reported decay products. About 3 days postmortem, new metabolites, including free trimethylammonium, propionate, butyrate, and iso-butyrate, started to appear in situ. Since the observed metabolites and the time course is comparable in sheep and human brain tissue, the model system seems to be appropriate.

Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximately 141 thousand years ago (Kya), an exit out-of-Africa approximately 51 Kya, and a recent colonization of the Americas approximately 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

Effects of elevated atmospheric CO2 and O3 on wood density, antomical proerties and decomposition of Northern Hardwoods

Anthropogenic activities continue to drive atmospheric CO2 and O3 concentrations to levels higher than during the pre-industrial era. Accumulating evidence indicates that both elevated CO2 and elevated O3 could modify the quantity and biochemistry of woody plant biomass. Anatomical properties of woody plants are largely influenced by the activity of the cambium and the growth characteristics of wood cells, which are in turn influenced by a range of environmental factors. Hence, alterations in the concentrations of atmospheric CO2 and / or O3 could also impact wood anatomical properties. Many fungi derive their metabolic resources for growth from plant litter, including woody tissue, and therefore modifications in the quantity, biochemistry and anatomical properties of woody plants in response to elevated CO2 and / or O3 could impact the community of wood-decaying fungi and rates of wood decomposition. Consequently carbon and nutrient cycling and productivity of terrestrial ecosystem could also be impacted. Alterations in wood structure and biochemistry of woody plants could also impact wood density and subsequently impact wood quality. This dissertation examined the long term effects of elevated CO2 and / or O3 on wood anatomical properties, wood density, wood-decaying fungi and wood decomposition of northern hardwood tree species at the Aspen Free-Air CO2 and O3 Enrichment (Aspen FACE) project, near Rhinelander, WI, USA. Anatomical properties of wood varied significantly with species and aspen genotypes and radial position within the stem. Elevated CO2 did not have significant effects on wood anatomical properties in trembling aspen, paper birch or sugar maple, except for marginally increasing (P < 0.1) the number of vessels per square millimeter. Elevated O3 marginally or significantly altered vessel lumen diameter, cell wall area and vessel lumen area proportions depending on species and radial position. In line with the modifications in the anatomical properties, elevated CO2 and O3, alone, significantly modified wood density but effects were species and / or genotype specific. However, the effects of elevated CO2 and O3, alone, on wood anatomical properties and density were ameliorated when in combination. Wood species had a much greater impact on the wood-decaying fungal community and initial wood decomposition rate than did growth or decomposition of wood in elevated CO2 and / or O3. Polyporales, Agaricales, and Russulales were the dominant orders of fungi isolated. Based on the current results, future higher levels of CO2 and O3 may have moderate effects on wood quality of northern hardwoods, but for utilization purposes these may not be considered significant. However, wood-decaying fungal community composition and decomposition of northern hardwoods may be altered via shifts in species and / or genotype composition under future higher levels of CO2 and O3.

Integration of computational models and experimental characterization to study internal frost damage in cementitious materials

The objective of this doctoral research is to investigate the internal frost damage due to crystallization pore pressure in porous cement-based materials by developing computational and experimental characterization tools. As an essential component of the U.S. infrastructure system, the durability of concrete has significant impact on maintenance costs. In cold climates, freeze-thaw damage is a major issue affecting the durability of concrete. The deleterious effects of the freeze-thaw cycle depend on the microscale characteristics of concrete such as the pore sizes and the pore distribution, as well as the environmental conditions. Recent theories attribute internal frost damage of concrete is caused by crystallization pore pressure in the cold environment. The pore structures have significant impact on freeze-thaw durability of cement/concrete samples. The scanning electron microscope (SEM) and transmission X-ray microscopy (TXM) techniques were applied to characterize freeze-thaw damage within pore structure. In the microscale pore system, the crystallization pressures at sub-cooling temperatures were calculated using interface energy balance with thermodynamic analysis. The multi-phase Extended Finite Element Modeling (XFEM) and bilinear Cohesive Zone Modeling (CZM) were developed to simulate the internal frost damage of heterogeneous cement-based material samples. The fracture simulation with these two techniques were validated by comparing the predicted fracture behavior with the captured damage from compact tension (CT) and single-edge notched beam (SEB) bending tests. The study applied the developed computational tools to simulate the internal frost damage caused by ice crystallization with the two dimensional (2-D) SEM and three dimensional (3-D) reconstructed SEM and TXM digital samples. The pore pressure calculated from thermodynamic analysis was input for model simulation. The 2-D and 3-D bilinear CZM predicted the crack initiation and propagation within cement paste microstructure. The favorably predicted crack paths in concrete/cement samples indicate the developed bilinear CZM techniques have the ability to capture crack nucleation and propagation in cement-based material samples with multiphase and associated interface. By comparing the computational prediction with the actual damaged samples, it also indicates that the ice crystallization pressure is the main mechanism for the internal frost damage in cementitious materials.

Using P-GENESIS for Parallel Simulation of GENESIS Models

P-GENESIS is an extension to the GENESIS neural simulator that allows users to take advantage of parallel machines to speed up the simulation of their network models or concurrently simulate multiple models. P-GENESIS adds several commands to the GENESIS script language that let a script running on one processor execute remote procedure calls on other processors, and that let a script synchronize its execution with the scripts running on other processors. We present here some brief comments on the mechanisms underlying parallel script execution. We also offer advice on parallelizing parameter searches, partitioning network models, and selecting suitable parallel hardware on which to run P-GENESIS.

Topographic time-frequency decomposition of the EEG

Frequency-transformed EEG resting data has been widely used to describe normal and abnormal brain functional states as function of the spectral power in different frequency bands. This has yielded a series of clinically relevant findings. However, by transforming the EEG into the frequency domain, the initially excellent time resolution of time-domain EEG is lost. The topographic time-frequency decomposition is a novel computerized EEG analysis method that combines previously available techniques from time-domain spatial EEG analysis and time-frequency decomposition of single-channel time series. It yields a new, physiologically and statistically plausible topographic time-frequency representation of human multichannel EEG. The original EEG is accounted by the coefficients of a large set of user defined EEG like time-series, which are optimized for maximal spatial smoothness and minimal norm. These coefficients are then reduced to a small number of model scalp field configurations, which vary in intensity as a function of time and frequency. The result is thus a small number of EEG field configurations, each with a corresponding time-frequency (Wigner) plot. The method has several advantages: It does not assume that the data is composed of orthogonal elements, it does not assume stationarity, it produces topographical maps and it allows to include user-defined, specific EEG elements, such as spike and wave patterns. After a formal introduction of the method, several examples are given, which include artificial data and multichannel EEG during different physiological and pathological conditions.