Developing and applying heterogeneous phylogenetic models with XRate

Modeling sequence evolution on phylogenetic trees is a useful technique in computational biology. Especially powerful are models which take account of the heterogeneous nature of sequence evolution according to the "grammar" of the encoded gene features. However, beyond a modest level of model complexity, manual coding of models becomes prohibitively labor-intensive. We demonstrate, via a set of case studies, the new built-in model-prototyping capabilities of XRate (macros and Scheme extensions). These features allow rapid implementation of phylogenetic models which would have previously been far more labor-intensive. XRate's new capabilities for lineage-specific models, ancestral sequence reconstruction, and improved annotation output are also discussed. XRate's flexible model-specification capabilities and computational efficiency make it well-suited to developing and prototyping phylogenetic grammar models. XRate is available as part of the DART software package: http://biowiki.org/DART .Comment: 34 pages, 3 figures, glossary of XRate model terminolog

Unveiling the oldest and most massive galaxies at very high redshift

(Abridged) This work explores the existence of high redshift massive galaxies unveiled with Spitzer+IRAC, but missed by conventional selection techniques based on optical and near-infrared observations. To this end, we use the multi-wavelength imaging data available for the GOODS-South field, and select a flux-limited sample from the IRAC 3.6um image to m(AB)<23.26. We confine our study to the galaxies undetected by the optical HST+ACS imaging and close to the detection limit of the K-band image (K>23.5 AB). Our selection unveiled 20 galaxies on which we performed a detailed photometric analysis. For each galaxy, we built an SED based on optical-to-8um photometry to estimate the photo-z and to derive the main galaxies physical properties. The majority of the sample sources show degenerate/bimodal solutions for the photometric redshifts (Abridged). These can either be heavily dust-enshrouded (Av~2-4) starbursts at 210^12 Lsun, or massive post-starburst galaxies in the redshift interval 4<z<9 with stellar masses of 10^11 Msun. One galaxy, the only source in our sample with both an X-ray and a 24um detection, might be an extremely massive object at z~8 detected during a post-starburst phase with concomitant QSO activity (although a lower-z solution is not excluded). Our investigation of Spitzer-selected galaxies with very red SEDs and completely undetected in the optical reveals a potential population of massive galaxies at z>4 which appear to include significant AGN emissions. These sources may be the oldest stellar systems at z~4. These, previously unrecognized, optically obscured objects might provide an important contribution to the massive-end (M>10^11 sun) of the high-z stellar mass function and they would almost double it (Abridged).Comment: Accepted for publication in Astronomy and Astrophysic

Cluster Lenses

Clusters of galaxies are the most recently assembled, massive, bound structures in the Universe. As predicted by General Relativity, given their masses, clusters strongly deform space-time in their vicinity. Clusters act as some of the most powerful gravitational lenses in the Universe. Light rays traversing through clusters from distant sources are hence deflected, and the resulting images of these distant objects therefore appear distorted and magnified. Lensing by clusters occurs in two regimes, each with unique observational signatures. The strong lensing regime is characterized by effects readily seen by eye, namely, the production of giant arcs, multiple-images, and arclets. The weak lensing regime is characterized by small deformations in the shapes of background galaxies only detectable statistically. Cluster lenses have been exploited successfully to address several important current questions in cosmology: (i) the study of the lens(es) - understanding cluster mass distributions and issues pertaining to cluster formation and evolution, as well as constraining the nature of dark matter; (ii) the study of the lensed objects - probing the properties of the background lensed galaxy population - which is statistically at higher redshifts and of lower intrinsic luminosity thus enabling the probing of galaxy formation at the earliest times right up to the Dark Ages; and (iii) the study of the geometry of the Universe - as the strength of lensing depends on the ratios of angular diameter distances between the lens, source and observer, lens deflections are sensitive to the value of cosmological parameters and offer a powerful geometric tool to probe Dark Energy. In this review, we present the basics of cluster lensing and provide a current status report of the field.Comment: About 120 pages - Published in Open Access at: http://www.springerlink.com/content/j183018170485723/ . arXiv admin note: text overlap with arXiv:astro-ph/0504478 and arXiv:1003.3674 by other author

Care of vision and ocular health in diabetic members of a national diabetes organization: A cross-sectional study

<p>Abstract</p> <p>Background</p> <p>Regular examination and early treatment of diabetic retinopathy can prevent visual loss. The aim of the study was to describe the care of vision and ocular health in people with diabetes in Norway.</p> <p>Methods</p> <p>A cross-sectional questionnaire survey of a random sample (n = 1,887) of the Norwegian Diabetic Associations' (NDA) members was carried out in 2005. Questions were asked about care of vision and ocular health, history of ocular disease and visual symptoms, general medical history and diabetes management. The study was approved by the Regional Committee for Medical Research Ethics.</p> <p>Results</p> <p>The response rate was 74%. Forty-four questionnaires with incomplete data regarding gender, age or type of diabetes were excluded, leaving 1352 cases (52% females) for analysis. 451 (33%) had type 1 and 901 (67%) had type 2 diabetes, the mean duration of diabetes was respectively, 22 (sd ± 14) and 10 (sd ± 9) years. In all 1,052 (78%) had their eyes examined according to guidelines and 1,169 (87%) confirmed to have received information about regular eye examinations. One in two recalled to have received such information from their general practitioner. To have received information about the importance of eye examinations (PR 3.1, 95% CI 2.4 to 4.0), and diabetes duration > 10 years (PR 1.2, 95% CI 1.2 to 1.3), were independently associated with reporting regular eye examinations. A history of diabetic retinopathy was reported by 178 (13%) responders, of which 101 (57%) reported a history of laser treatment. Responders who had regular eye examinations reported more frequently a history of diabetic retinopathy (19% vs. 5%, p < 0.001). The frequency of retinopathy was significantly higher in responders with reported HbA1c values above treatment target (23% vs. 13%, p = 0.001). However, in responders who were not regularly examined, there was no difference in reported frequency of retinopathy with regard to HbA1c level.</p> <p>Conclusion</p> <p>Eight out of ten diabetic members of the NDA had their eyes examined according to current guidelines and the majority was well informed about the risk of vision loss due to diabetes. The results indicate that the reported history of diabetic retinopathy likely underestimates the prevalence of retinopathy.</p

Transat—A Method for Detecting the Conserved Helices of Functional RNA Structures, Including Transient, Pseudo-Knotted and Alternative Structures

The prediction of functional RNA structures has attracted increased interest, as it allows us to study the potential functional roles of many genes. RNA structure prediction methods, however, assume that there is a unique functional RNA structure and also do not predict functional features required for in vivo folding. In order to understand how functional RNA structures form in vivo, we require sophisticated experiments or reliable prediction methods. So far, there exist only a few, experimentally validated transient RNA structures. On the computational side, there exist several computer programs which aim to predict the co-transcriptional folding pathway in vivo, but these make a range of simplifying assumptions and do not capture all features known to influence RNA folding in vivo. We want to investigate if evolutionarily related RNA genes fold in a similar way in vivo. To this end, we have developed a new computational method, Transat, which detects conserved helices of high statistical significance. We introduce the method, present a comprehensive performance evaluation and show that Transat is able to predict the structural features of known reference structures including pseudo-knotted ones as well as those of known alternative structural configurations. Transat can also identify unstructured sub-sequences bound by other molecules and provides evidence for new helices which may define folding pathways, supporting the notion that homologous RNA sequence not only assume a similar reference RNA structure, but also fold similarly. Finally, we show that the structural features predicted by Transat differ from those assuming thermodynamic equilibrium. Unlike the existing methods for predicting folding pathways, our method works in a comparative way. This has the disadvantage of not being able to predict features as function of time, but has the considerable advantage of highlighting conserved features and of not requiring a detailed knowledge of the cellular environment

Efficient Algorithms for Probing the RNA Mutation Landscape

The diversity and importance of the role played by RNAs in the regulation and development of the cell are now well-known and well-documented. This broad range of functions is achieved through specific structures that have been (presumably) optimized through evolution. State-of-the-art methods, such as McCaskill's algorithm, use a statistical mechanics framework based on the computation of the partition function over the canonical ensemble of all possible secondary structures on a given sequence. Although secondary structure predictions from thermodynamics-based algorithms are not as accurate as methods employing comparative genomics, the former methods are the only available tools to investigate novel RNAs, such as the many RNAs of unknown function recently reported by the ENCODE consortium. In this paper, we generalize the McCaskill partition function algorithm to sum over the grand canonical ensemble of all secondary structures of all mutants of the given sequence. Specifically, our new program, RNAmutants, simultaneously computes for each integer k the minimum free energy structure MFE(k) and the partition function Z(k) over all secondary structures of all k-point mutants, even allowing the user to specify certain positions required not to mutate and certain positions required to base-pair or remain unpaired. This technically important extension allows us to study the resilience of an RNA molecule to pointwise mutations. By computing the mutation profile of a sequence, a novel graphical representation of the mutational tendency of nucleotide positions, we analyze the deleterious nature of mutating specific nucleotide positions or groups of positions. We have successfully applied RNAmutants to investigate deleterious mutations (mutations that radically modify the secondary structure) in the Hepatitis C virus cis-acting replication element and to evaluate the evolutionary pressure applied on different regions of the HIV trans-activation response element. In particular, we show qualitative agreement between published Hepatitis C and HIV experimental mutagenesis studies and our analysis of deleterious mutations using RNAmutants. Our work also predicts other deleterious mutations, which could be verified experimentally. Finally, we provide evidence that the 3′ UTR of the GB RNA virus C has been optimized to preserve evolutionarily conserved stem regions from a deleterious effect of pointwise mutations. We hope that there will be long-term potential applications of RNAmutants in de novo RNA design and drug design against RNA viruses. This work also suggests potential applications for large-scale exploration of the RNA sequence-structure network. Binary distributions are available at http://RNAmutants.csail.mit.edu/

Febrile seizures and mechanisms of epileptogenesis: insights from an animal model.

Temporal lobe epilepsy (TLE) is the most prevalent type of human epilepsy, yet the causes for its development, and the processes involved, are not known. Most individuals with TLE do not have a family history, suggesting that this limbic epilepsy is a consequence of acquired rather than genetic causes. Among suspected etiologies, febrile seizures have frequently been cited. This is due to the fact that retrospective analyses of adults with TLE have demonstrated a high prevalence (20--&gt;60%) of a history of prolonged febrile seizures during early childhood, suggesting an etiological role for these seizures in the development of TLE. Specifically, neuronal damage induced by febrile seizures has been suggested as a mechanism for the development of mesial temporal sclerosis, the pathological hallmark of TLE. However, the statistical correlation between febrile seizures and TLE does not necessarily indicate a causal relationship. For example, preexisting (genetic or acquired) 'causes' that result independently in febrile seizures and in TLE would also result in tight statistical correlation. For obvious reasons, complex febrile seizures cannot be induced in the human, and studies of their mechanisms and of their consequences on brain molecules and circuits are severely limited. Therefore, an animal model was designed to study these seizures. The model reproduces the fundamental key elements of the human condition: the age specificity, the physiological temperatures seen in fevers of children, the length of the seizures and their lack of immediate morbidity. Neuroanatomical, molecular and functional methods have been used in this model to determine the consequences of prolonged febrile seizures on the survival and integrity of neurons, and on hyperexcitability in the hippocampal-limbic network. Experimental prolonged febrile seizures did not lead to death of any of the seizure-vulnerable populations in hippocampus, and the rate of neurogenesis was also unchanged. Neuronal function was altered sufficiently to promote synaptic reorganization of granule cells, and transient and long-term alterations in the expression of specific genes were observed. The contribution of these consequences of febrile seizures to the epileptogenic process is discussed

Dosage-Sensitive Function of RETINOBLASTOMA RELATED and Convergent Epigenetic Control Are Required during the Arabidopsis Life Cycle

The plant life cycle alternates between two distinct multi-cellular generations, the reduced gametophytes and the dominant sporophyte. Little is known about how generation-specific cell fate, differentiation, and development are controlled by the core regulators of the cell cycle. In Arabidopsis, RETINOBLASTOMA RELATED (RBR), an evolutionarily ancient cell cycle regulator, controls cell proliferation, differentiation, and regulation of a subset of Polycomb Repressive Complex 2 (PRC2) genes and METHYLTRANSFERASE 1 (MET1) in the male and female gametophytes, as well as cell fate establishment in the male gametophyte. Here we demonstrate that RBR is also essential for cell fate determination in the female gametophyte, as revealed by loss of cell-specific marker expression in all the gametophytic cells that lack RBR. Maintenance of genome integrity also requires RBR, because diploid plants heterozygous for rbr (rbr/RBR) produce an abnormal portion of triploid offspring, likely due to gametic genome duplication. While the sporophyte of the diploid mutant plants phenocopied wild type due to the haplosufficiency of RBR, genetic analysis of tetraploid plants triplex for rbr (rbr/rbr/rbr/RBR) revealed that RBR has a dosage-dependent pleiotropic effect on sporophytic development, trichome differentiation, and regulation of PRC2 subunit genes CURLY LEAF (CLF) and VERNALIZATION 2 (VRN2), and MET1 in leaves. There were, however, no obvious cell cycle and cell proliferation defects in these plant tissues, suggesting that a single functional RBR copy in tetraploids is capable of maintaining normal cell division but is not sufficient for distinct differentiation and developmental processes. Conversely, in leaves of mutants in sporophytic PRC2 subunits, trichome differentiation was also affected and expression of RBR and MET1 was reduced, providing evidence for a RBR-PRC2-MET1 regulatory feedback loop involved in sporophyte development. Together, dosage-sensitive RBR function and its genetic interaction with PRC2 genes and MET1 must have been recruited during plant evolution to control distinct generation-specific cell fate, differentiation, and development

Coordinated Regulation of Intestinal Functions in C. elegans by LIN-35/Rb and SLR-2

LIN-35 is the sole C. elegans representative of the pocket protein family, which includes the mammalian Retinoblastoma protein pRb and its paralogs p107 and p130. In addition to having a well-established and central role in cell cycle regulation, pocket proteins have been increasingly implicated in the control of critical and diverse developmental and cellular processes. To gain a greater understanding of the roles of pocket proteins during development, we have characterized a synthetic genetic interaction between lin-35 and slr-2, which we show encodes a C2H2-type Zn-finger protein. Whereas animals harboring single mutations in lin-35 or slr-2 are viable and fertile, lin-35; slr-2 double mutants arrest uniformly in early larval development without obvious morphological defects. Using a combination of approaches including transcriptome profiling, mosaic analysis, starvation assays, and expression analysis, we demonstrate that both LIN-35 and SLR-2 act in the intestine to regulate the expression of many genes required for normal nutrient utilization. These findings represent a novel role for pRb family members in the maintenance of organ function. Our studies also shed light on the mechanistic basis of genetic redundancy among transcriptional regulators and suggest that synthetic interactions may result from the synergistic misregulation of one or more common targets