EGASP: the human ENCODE Genome Annotation Assessment Project.

BACKGROUND: We present the results of EGASP, a community experiment to assess the state-of-the-art in genome annotation within the ENCODE regions, which span 1% of the human genome sequence. The experiment had two major goals: the assessment of the accuracy of computational methods to predict protein coding genes; and the overall assessment of the completeness of the current human genome annotations as represented in the ENCODE regions. For the computational prediction assessment, eighteen groups contributed gene predictions. We evaluated these submissions against each other based on a 'reference set' of annotations generated as part of the GENCODE project. These annotations were not available to the prediction groups prior to the submission deadline, so that their predictions were blind and an external advisory committee could perform a fair assessment. RESULTS: The best methods had at least one gene transcript correctly predicted for close to 70% of the annotated genes. Nevertheless, the multiple transcript accuracy, taking into account alternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotide level, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programs relying on mRNA and protein sequences were the most accurate in reproducing the manually curated annotations. Experimental validation shows that only a very small percentage (3.2%) of the selected 221 computationally predicted exons outside of the existing annotation could be verified. CONCLUSION: This is the first such experiment in human DNA, and we have followed the standards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe the results presented here contribute to the value of ongoing large-scale annotation projects and should guide further experimental methods when being scaled up to the entire human genome sequence.

Dynamics of case-fatalilty rates of recurrent thromboembolism and major bleeding in patients treated for venous thromboembolism.

In patients with venous thromboembolism (VTE), assessment of the risk of fatal recurrent VTE and fatal bleeding during anticoagulation may help to guide intensity and duration of therapy. We aimed to provide estimates of the case-fatality rate (CFR) of recurrent VTE and major bleeding during anticoagulation in a 'real life' population, and to assess these outcomes according to the initial presentation of VTE and its etiology. The study included 41,826 patients with confirmed VTE from the RIETE registry who received different durations of anticoagulation (mean 7.8 ± 0.6 months). During 27,110 patient-years, the CFR was 12.1% (95% CI, 10.2-14.2) for recurrent VTE, and 19.7% (95% CI, 17.4-22.1) for major bleeding. During the first three months of anticoagulant therapy, the CFR of recurrent VTE was 16.1% (95% CI, 13.6-18.9), compared to 2.0% (95% CI, 0-4.2) beyond this period. The CFR of bleeding was 20.2% (95% CI, 17.5-23.1) during the first three months, compared to 18.2% (95% CI, 14.0-23.2) beyond this period. The CFR of recurrent VTE was higher in patients initially presenting with PE (18.5%; 95% CI, 15.3-22.1) than in those with DVT (6.3%; 95% CI, 4.5-8.6), and in patients with provoked VTE (16.3%; 95% CI, 13.6-19.4) than in those with unprovoked VTE (5.5%; 95% CI, 3.5-8.0). In conclusion, the CFR of recurrent VTE decreased over time during anticoagulation, while the CFR of major bleeding remained stable. The CFR of recurrent VTE was higher in patients initially presenting with PE and in those with provoked VTE.

Effects of domain morphology in phase-separation dynamics at low temperature

We present numerical results of the deterministic Ginzburg-Landau equation with a concentration-dependent diffusion coefficient, for different values of the volume fraction phi of the minority component. The morphology of the domains affects the dynamics of phase separation. The effective growth exponents, but not the scaled functions, are found to be temperature dependent.

Cation dynamics in soils with different salinity levels growing irrigated rice

Salinity levels in soils of the Outer Coastal Plain of Rio Grande do Sul, Brazil, can be high, due to excess of Na in the irrigation water, evapotranspiration and soil development from marine sediments. The cultivation of irrigated rice could be an alternative, since ion uptake as well as leaching by the establishment of a water layer could mitigate the effects of soil salinity. This study aimed to evaluate the dynamics of basic cations in the solution of Albaqualf soils with different salinity levels growing irrigated rice. The plow layer contained exchangeable Na percentages (ESP) of 5.6, 9.0, 21.2 and 32.7 %. The plant stand, dry matter, Na, K and Ca + Mg uptake at full flowering and grain yield were evaluated. The levels of Na, K, Ca + Mg and electrical conductivity (EC) in the soil solution were also measured weekly during the rice cycle at four soil depths, in the water layer and irrigation water. The Na, K and Ca + Mg uptake by rice at full flowering was used to estimate ion depletion from the layer under root influence. Soil salinity induced a reduction in the rice stand, especially in the soil with ESP of 32.7 %, resulting in lower cation uptake and very low yield at that site. As observed in the water layer and irrigation water, the Na, K, Ca + Mg and EC levels in the soil solution decreased with time at depths of 5, 10 and 20 cm, regardless of the original soil salinity, showing that cation dynamics in the plow layer was determined by leaching and root uptake, rather than by the effect of evapoconcentration of basic cations in the soil surface layer.

Comparison between molecular dynamics abd Monte Carlo simulations of an ordering process in a binary alloy

Ordering in a binary alloy is studied by means of a molecular-dynamics (MD) algorithm which allows to reach the domain growth regime. Results are compared with Monte Carlo simulations using a realistic vacancy-atom (MC-VA) mechanism. At low temperatures fast growth with a dynamical exponent x>1/2 is found for MD and MC-VA. The study of a nonequilibrium ordering process with the two methods shows the importance of the nonhomogeneity of the excitations in the system for determining its macroscopic kinetics.

Soil phosphorus dynamics and availability and irrigated coffee yield

Research data have demonstrated that the P demand of coffee (Coffea arabica L.) is similar to that of short-cycle crops. In this context, the objective of this study was to evaluate the influence of annual P fertilization on the soil P status by the quantification of labile, moderately labile, low-labile, and total P fractions, associating them to coffee yield. The experiment was installed in a typical dystrophic Red Latosol (Oxisol) cultivated with irrigated coffee annually fertilized with triple superphosphate at rates of 0, 50, 100, 200, and 400 kg ha-1 P2O5. Phosphorus fractions were determined in two soil layers: 0-10 and 10-20 cm. The P leaf contents and coffee yield in 2008 were also evaluated. The irrigated coffee responded to phosphate fertilization in the production phase with gains of up to 138 % in coffee yield by the application of 400 kg ha-1 P2O5. Coffee leaf P contents increased with P applications and stabilized around 1.98 g kg-1, at rates of 270 kg ha-1 P2O5 and higher. Soil P application caused, in general, an increase in bioavailable P fractions, which constitute the main soil P reservoir.

Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.

Spanning avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics: Field dependence and geometrical properties

Spanning avalanches in the 3D Gaussian Random Field Ising Model (3D-GRFIM) with metastable dynamics at T=0 have been studied. Statistical analysis of the field values for which avalanches occur has enabled a Finite-Size Scaling (FSS) study of the avalanche density to be performed. Furthermore, a direct measurement of the geometrical properties of the avalanches has confirmed an earlier hypothesis that several types of spanning avalanches with two different fractal dimensions coexist at the critical point. We finally compare the phase diagram of the 3D-GRFIM with metastable dynamics with the same model in equilibrium at T=0.

Hysteresis and avalanches in the T=0 random field ising model with 2-spin flip dynamics

We study the nonequilibrium behavior of the three-dimensional Gaussian random-field Ising model at T=0 in the presence of a uniform external field using a two-spin-flip dynamics. The deterministic, history-dependent evolution of the system is compared with the one obtained with the standard one-spin-flip dynamics used in previous studies of the model. The change in the dynamics yields a significant suppression of coercivity, but the distribution of avalanches (in number and size) stays remarkably similar, except for the largest ones that are responsible for the jump in the saturation magnetization curve at low disorder in the thermodynamic limit. By performing a finite-size scaling study, we find strong evidence that the change in the dynamics does not modify the universality class of the disorder-induced phase transition.

Diluted three-dimensional random field Ising model at zero temperature with metastable dynamics.

The influence of vacancy concentration on the behavior of the three-dimensional random field Ising model with metastable dynamics is studied. We have focused our analysis on the number of spanning avalanches which allows us a clear determination of the critical line where the hysteresis loops change from continuous to discontinuous. By a detailed finite-size scaling analysis we determine the phase diagram and numerically estimate the critical exponents along the whole critical line. Finally, we discuss the origin of the curvature of the critical line at high vacancy concentration.

Influence of the driving mechanism on the response of systems with athermal dynamics: The example of the random-field Ising model

We investigate the influence of the driving mechanism on the hysteretic response of systems with athermal dynamics. In the framework of local mean-field theory at finite temperature (but neglecting thermally activated processes), we compare the rate-independent hysteresis loops obtained in the random field Ising model when controlling either the external magnetic field H or the extensive magnetization M. Two distinct behaviors are observed, depending on disorder strength. At large disorder, the H-driven and M-driven protocols yield identical hysteresis loops in the thermodynamic limit. At low disorder, when the H-driven magnetization curve is discontinuous (due to the presence of a macroscopic avalanche), the M-driven loop is reentrant while the induced field exhibits strong intermittent fluctuations and is only weakly self-averaging. The relevance of these results to the experimental observations in ferromagnetic materials, shape memory alloys, and other disordered systems is discussed.

Statistics of transverse mode turn-on dynamics in VCSEL's

The turn-on process of a multimode VCSEL is investigated from a statistical point of view. Special attention is paid to quantities such as time jitter and bit error rate. The single-mode performance of VCSEL¿s during current modulation is compared to that of edge-emitting lasers.

Large scale assessment of regulatory evolution and transcriptome complexity in mammals

AbstractIn addition to genetic changes affecting the function of gene products, changes in gene expression have been suggested to underlie many or even most of the phenotypic differences among mammals. However, detailed gene expression comparisons were, until recently, restricted to closely related species, owing to technological limitations. Thus, we took advantage of the latest technologies (RNA-Seq) to generate extensive qualitative and quantitative transcriptome data for a unique collection of somatic and germline tissues from representatives of all major mammalian lineages (placental mammals, marsupials and monotremes) and birds, the evolutionary outgroup.In the first major project of my thesis, we performed global comparative analyses of gene expression levels based on these data. Our analyses provided fundamental insights into the dynamics of transcriptome change during mammalian evolution (e.g., the rate of expression change across species, tissues and chromosomes) and allowed the exploration of the functional relevance and phenotypic implications of transcription changes at a genome-wide scale (e.g., we identified numerous potentially selectively driven expression switches).In a second project of my thesis, which was also based on the unique transcriptome data generated in the context of the first project we focused on the evolution of alternative splicing in mammals. Alternative splicing contributes to transcriptome complexity by generating several transcript isoforms from a single gene, which can, thus, perform various functions. To complete the global comparative analysis of gene expression changes, we explored patterns of alternative splicing evolution. This work uncovered several general and unexpected patterns of alternative splicing evolution (e.g., we found that alternative splicing evolves extremely rapidly) as well as a large number of conserved alternative isoforms that may be crucial for the functioning of mammalian organs.Finally, the third and final project of my PhD consisted in analyzing in detail the unique functional and evolutionary properties of the testis by exploring the extent of its transcriptome complexity. This organ was previously shown to evolve rapidly both at the phenotypic and molecular level, apparently because of the specific pressures that act on this organ and are associated with its reproductive function. Moreover, my analyses of the amniote tissue transcriptome data described above, revealed strikingly widespread transcriptional activity of both functional and nonfunctional genomic elements in the testis compared to the other organs. To elucidate the cellular source and mechanisms underlying this promiscuous transcription in the testis, we generated deep coverage RNA-Seq data for all major testis cell types as well as epigenetic data (DNA and histone methylation) using the mouse as model system. The integration of these complete dataset revealed that meiotic and especially post-meiotic germ cells are the major contributors to the widespread functional and nonfunctional transcriptome complexity of the testis, and that this "promiscuous" spermatogenic transcription is resulting, at least partially, from an overall transcriptionally permissive chromatin state. We hypothesize that this particular open state of the chromatin results from the extensive chromatin remodeling that occurs during spermatogenesis which ultimately leads to the replacement of histones by protamines in the mature spermatozoa. Our results have important functional and evolutionary implications (e.g., regarding new gene birth and testicular gene expression evolution).Generally, these three large-scale projects of my thesis provide complete and massive datasets that constitute valuables resources for further functional and evolutionary analyses of mammalian genomes.