Select-divide-and-conquer method for large-scale configuration interaction

A select-divide-and-conquer variational method to approximate configuration interaction (CI) is presented. Given an orthonormal set made up of occupied orbitals (Hartree-Fock or similar) and suitable correlation orbitals (natural or localized orbitals), a large N-electron target space S is split into subspaces S0,S1,S2,...,SR. S0, of dimension d0, contains all configurations K with attributes (energy contributions, etc.) above thresholds T0={T0egy, T0etc.}; the CI coefficients in S0 remain always free to vary. S1 accommodates KS with attributes above T1≤T0. An eigenproblem of dimension d0+d1 for S0+S 1 is solved first, after which the last d1 rows and columns are contracted into a single row and column, thus freezing the last d1 CI coefficients hereinafter. The process is repeated with successive Sj(j≥2) chosen so that corresponding CI matrices fit random access memory (RAM). Davidson's eigensolver is used R times. The final energy eigenvalue (lowest or excited one) is always above the corresponding exact eigenvalue in S. Threshold values {Tj;j=0, 1, 2,...,R} regulate accuracy; for large-dimensional S, high accuracy requires S 0+S1 to be solved outside RAM. From there on, however, usually a few Davidson iterations in RAM are needed for each step, so that Hamiltonian matrix-element evaluation becomes rate determining. One μhartree accuracy is achieved for an eigenproblem of order 24 × 106, involving 1.2 × 1012 nonzero matrix elements, and 8.4×109 Slater determinants

Rates of cultural change and patterns of cultural accumulation in stochastic models of social transmission.

Cultural variation in a population is affected by the rate of occurrence of cultural innovations, whether such innovations are preferred or eschewed, how they are transmitted between individuals in the population, and the size of the population. An innovation, such as a modification in an attribute of a handaxe, may be lost or may become a property of all handaxes, which we call "fixation of the innovation." Alternatively, several innovations may attain appreciable frequencies, in which case properties of the frequency distribution-for example, of handaxe measurements-is important. Here we apply the Moran model from the stochastic theory of population genetics to study the evolution of cultural innovations. We obtain the probability that an initially rare innovation becomes fixed, and the expected time this takes. When variation in cultural traits is due to recurrent innovation, copy error, and sampling from generation to generation, we describe properties of this variation, such as the level of heterogeneity expected in the population. For all of these, we determine the effect of the mode of social transmission: conformist, where there is a tendency for each naïve newborn to copy the most popular variant; pro-novelty bias, where the newborn prefers a specific variant if it exists among those it samples; one-to-many transmission, where the variant one individual carries is copied by all newborns while that individual remains alive. We compare our findings with those predicted by prevailing theories for rates of cultural change and the distribution of cultural variation.

Patterns of plant subcellular responses to successful oomycete infections reveal differences in host cell reprogramming and endocytic trafficking.

Adapted filamentous pathogens such as the oomycetes Hyaloperonospora arabidopsidis (Hpa) and Phytophthora infestans (Pi) project specialized hyphae, the haustoria, inside living host cells for the suppression of host defence and acquisition of nutrients. Accommodation of haustoria requires reorganization of the host cell and the biogenesis of a novel host cell membrane, the extrahaustorial membrane (EHM), which envelops the haustorium separating the host cell from the pathogen. Here, we applied live-cell imaging of fluorescent-tagged proteins labelling a variety of membrane compartments and investigated the subcellular changes associated with accommodating oomycete haustoria in Arabidopsis and N. benthamiana. Plasma membrane-resident proteins differentially localized to the EHM. Likewise, secretory vesicles and endosomal compartments surrounded Hpa and Pi haustoria revealing differences between these two oomycetes, and suggesting a role for vesicle trafficking pathways for the pathogen-controlled biogenesis of the EHM. The latter is supported by enhanced susceptibility of mutants in endosome-mediated trafficking regulators. These observations point at host subcellular defences and specialization of the EHM in a pathogen-specific manner. Defence-associated haustorial encasements, a double-layered membrane that grows around mature haustoria, were frequently observed in Hpa interactions. Intriguingly, all tested plant proteins accumulated at Hpa haustorial encasements suggesting the general recruitment of default vesicle trafficking pathways to defend pathogen access. Altogether, our results show common requirements of subcellular changes associated with oomycete biotrophy, and highlight differences between two oomycete pathogens in reprogramming host cell vesicle trafficking for haustoria accommodation. This provides a framework for further dissection of the pathogen-triggered reprogramming of host subcellular changes.

Temporal patterns of nucleotide misincorporations and DNA fragmentation in ancient DNA.

DNA that survives in museum specimens, bones and other tissues recovered by archaeologists is invariably fragmented and chemically modified. The extent to which such modifications accumulate over time is largely unknown but could potentially be used to differentiate between endogenous old DNA and present-day DNA contaminating specimens and experiments. Here we examine mitochondrial DNA sequences from tissue remains that vary in age between 18 and 60,000 years with respect to three molecular features: fragment length, base composition at strand breaks, and apparent C to T substitutions. We find that fragment length does not decrease consistently over time and that strand breaks occur preferentially before purine residues by what may be at least two different molecular mechanisms that are not yet understood. In contrast, the frequency of apparent C to T substitutions towards the 5'-ends of molecules tends to increase over time. These nucleotide misincorporations are thus a useful tool to distinguish recent from ancient DNA sources in specimens that have not been subjected to unusual or harsh treatments.

The Paediatric Rheumatology International Trials Organisation provisional criteria for the evaluation of response to therapy in juvenile dermatomyositis.

OBJECTIVE: To develop a provisional definition for the evaluation of response to therapy in juvenile dermatomyositis (DM) based on the Paediatric Rheumatology International Trials Organisation juvenile DM core set of variables. METHODS: Thirty-seven experienced pediatric rheumatologists from 27 countries achieved consensus on 128 difficult patient profiles as clinically improved or not improved using a stepwise approach (patient's rating, statistical analysis, definition selection). Using the physicians' consensus ratings as the "gold standard measure," chi-square, sensitivity, specificity, false-positive and-negative rates, area under the receiver operating characteristic curve, and kappa agreement for candidate definitions of improvement were calculated. Definitions with kappa values >0.8 were multiplied by the face validity score to select the top definitions. RESULTS: The top definition of improvement was at least 20% improvement from baseline in 3 of 6 core set variables with no more than 1 of the remaining worsening by more than 30%, which cannot be muscle strength. The second-highest scoring definition was at least 20% improvement from baseline in 3 of 6 core set variables with no more than 2 of the remaining worsening by more than 25%, which cannot be muscle strength (definition P1 selected by the International Myositis Assessment and Clinical Studies group). The third is similar to the second with the maximum amount of worsening set to 30%. This indicates convergent validity of the process. CONCLUSION: We propose a provisional data-driven definition of improvement that reflects well the consensus rating of experienced clinicians, which incorporates clinically meaningful change in core set variables in a composite end point for the evaluation of global response to therapy in juvenile DM.

Aging effects on selective attention-related electroencephalographic patterns during face encoding.

Previous electrophysiological studies revealed that human faces elicit an early visual event-related potential (ERP) within the occipito-temporal cortex, the N170 component. Although face perception has been proposed to rely on automatic processing, the impact of selective attention on N170 remains controversial both in young and elderly individuals. Using early visual ERP and alpha power analysis, we assessed the influence of aging on selective attention to faces during delayed-recognition tasks for face and letter stimuli, examining 36 elderly and 20 young adults with preserved cognition. Face recognition performance worsened with age. Aging induced a latency delay of the N1 component for faces and letters, as well as of the face N170 component. Contrasting with letters, ignored faces elicited larger N1 and N170 components than attended faces in both age groups. This counterintuitive attention effect on face processing persisted when scenes replaced letters. In contrast with young, elderly subjects failed to suppress irrelevant letters when attending faces. Whereas attended stimuli induced a parietal alpha band desynchronization within 300-1000 ms post-stimulus with bilateral-to-right distribution for faces and left lateralization for letters, ignored and passively viewed stimuli elicited a central alpha synchronization larger on the right hemisphere. Aging delayed the latency of this alpha synchronization for both face and letter stimuli, and reduced its amplitude for ignored letters. These results suggest that due to their social relevance, human faces may cause paradoxical attention effects on early visual ERP components, but they still undergo classical top-down control as a function of endogenous selective attention. Aging does not affect the face bottom-up alerting mechanism but reduces the top-down suppression of distracting letters, possibly impinging upon face recognition, and more generally delays the top-down suppression of task-irrelevant information.

1st datings of ammonite-bearing Bajocian suites from Caucasus (Georgia) - K-Ar ages of volcanic hornblendes

Volcanic flows and tuffs interbedded with ammonite-bearing sediments directly correlatable with the stratotype section of the Bajocian stage have been dated for the first time within the Caucasus area. Three samples, each from a different section, allowed separation of well-preserved brown hornblende; these are considered reliable geochronometers in a region where subsequent volcanic activity occurred. The dated separates are V139: a volcanic layer probably near the base, V142 another layer near the top of the Lower Bajocian substage; V141 a boulder from a latest Bajocian volcanic conglomerate. From bottom to top, apparent ages at 173.5 +/- 2.6, 164.8 +/- 2.5 and 167.1 +/- 1.9 Ma (analytical uncertainty, 95% confidence level) respectively, can be calculated. The consistency of the results is obtained if the dated boulder is interpreted as derived from an underlying layer, The Bajocian-Bathonian boundary is much younger than commonly accepted and younger than 164 Ma.

Incremental growth of the Patagonian Torres del Paine laccolith over 90 k.y.

The Miocene Paine Granite in the Torres del Paine Intrusive Complex, southern Chile, is an extraordinary example of an upper crustal mafic and granitic intrusion. The granite intruded as a series of three sheets, each one underplating the previous sheet along the top of the basal Paine Mafic Complex. High-precision U/Pb geochronology on single zircons using isotope dilution-thermal ionization mass spectrometry yields distinct ages of 12.59 +/- 0.02 Ma and 12.50 +/- 0.02 Ma, respectively, for the first and last sheet of the laccolith. This age relationship is consistent with field observations. The zircon ages define a time frame of 90 +/- 40 k.y. for the emplacement of a >2000-m-thick granite laccollith. These precise U-Pb zircon ages permit identification of the pulses in a 20 k.y. range. The data obtained for the Paine Granite fill the gap between 100 k.y. and 100-1000 yr pulses described in the literature for crustal magma chambers.

Changes of beta-amyloid precursor protein splice patterns in brain cell aggregate cultures.

The splice pattern of beta-amyloid precursor protein (beta-APP) has been studied in a variety of neuronal and glial cells and in brain cell aggregate cultures by the polymerase chain reaction (PCR). The brain-typical pattern, in which beta-APP695 is the dominant form, has been found only in aggregate cultures but not in any of the other cell types including neuronal cell lines. Selective elimination of glial cells from aggregates resulted in increased quantities of beta-APP695, whereas removal of neurons led to a reduction of beta-APP695 and to an elevation of beta-APP751 and beta-APP770. This shift of splice pattern was not observed in cocultures of the neuronal cell line PC 12 with primary astrocytes combined in a variety of cellular ratios. Blood serum, which is an essential component of these cultures, tested on aggregates, did not reduce the amount of beta-APP695 or have any marked effects on splice patterns generally. From these results it is concluded that investigations on brain-typical splicing of beta-APP require primary neurons. Neuronal cell lines may be no suitable model systems. Splicing events favoring production of beta-APP695 may mark an important, very early step of amyloid formation in the brain.

Mutations in FGF17, IL17RD, DUSP6, SPRY4, and FLRT3 Are Identified in Individuals with Congenital Hypogonadotropic Hypogonadism.

Congenital hypogonadotropic hypogonadism (CHH) and its anosmia-associated form (Kallmann syndrome [KS]) are genetically heterogeneous. Among the >15 genes implicated in these conditions, mutations in FGF8 and FGFR1 account for ∼12% of cases; notably, KAL1 and HS6ST1 are also involved in FGFR1 signaling and can be mutated in CHH. We therefore hypothesized that mutations in genes encoding a broader range of modulators of the FGFR1 pathway might contribute to the genetics of CHH as causal or modifier mutations. Thus, we aimed to (1) investigate whether CHH individuals harbor mutations in members of the so-called "FGF8 synexpression" group and (2) validate the ability of a bioinformatics algorithm on the basis of protein-protein interactome data (interactome-based affiliation scoring [IBAS]) to identify high-quality candidate genes. On the basis of sequence homology, expression, and structural and functional data, seven genes were selected and sequenced in 386 unrelated CHH individuals and 155 controls. Except for FGF18 and SPRY2, all other genes were found to be mutated in CHH individuals: FGF17 (n = 3 individuals), IL17RD (n = 8), DUSP6 (n = 5), SPRY4 (n = 14), and FLRT3 (n = 3). Independently, IBAS predicted FGF17 and IL17RD as the two top candidates in the entire proteome on the basis of a statistical test of their protein-protein interaction patterns to proteins known to be altered in CHH. Most of the FGF17 and IL17RD mutations altered protein function in vitro. IL17RD mutations were found only in KS individuals and were strongly linked to hearing loss (6/8 individuals). Mutations in genes encoding components of the FGF pathway are associated with complex modes of CHH inheritance and act primarily as contributors to an oligogenic genetic architecture underlying CHH.

Ecological assembly rules in plant communities--approaches, patterns and prospects.

Understanding how communities of living organisms assemble has been a central question in ecology since the early days of the discipline. Disentangling the different processes involved in community assembly is not only interesting in itself but also crucial for an understanding of how communities will behave under future environmental scenarios. The traditional concept of assembly rules reflects the notion that species do not co-occur randomly but are restricted in their co-occurrence by interspecific competition. This concept can be redefined in a more general framework where the co-occurrence of species is a product of chance, historical patterns of speciation and migration, dispersal, abiotic environmental factors, and biotic interactions, with none of these processes being mutually exclusive. Here we present a survey and meta-analyses of 59 papers that compare observed patterns in plant communities with null models simulating random patterns of species assembly. According to the type of data under study and the different methods that are applied to detect community assembly, we distinguish four main types of approach in the published literature: species co-occurrence, niche limitation, guild proportionality and limiting similarity. Results from our meta-analyses suggest that non-random co-occurrence of plant species is not a widespread phenomenon. However, whether this finding reflects the individualistic nature of plant communities or is caused by methodological shortcomings associated with the studies considered cannot be discerned from the available metadata. We advocate that more thorough surveys be conducted using a set of standardized methods to test for the existence of assembly rules in data sets spanning larger biological and geographical scales than have been considered until now. We underpin this general advice with guidelines that should be considered in future assembly rules research. This will enable us to draw more accurate and general conclusions about the non-random aspect of assembly in plant communities.

Extent and patterns of MGMT promoter methylation in glioblastoma- and respective glioblastoma-derived spheres.

Purpose: Quantitative methylation-specific tests suggest that not all cells in a glioblastoma with detectable promoter methylation of the O6-methylguanine DNA methyltransferase (MGMT) gene carry a methylated MGMT allele. This observation may indicate cell subpopulations with distinct MGMT status, raising the question of the clinically relevant cutoff of MGMT methylation therapy. Epigenetic silencing of the MGMT gene by promoter methylation blunts repair of O6-methyl guanine and has been shown to be a predictive factor for benefit from alkylating agent therapy in glioblastoma. Experimental Design: Ten paired samples of glioblastoma and respective glioblastoma-derived spheres (GS), cultured under stem cell conditions, were analyzed for the degree and pattern of MGMT promoter methylation by methylation-specific clone sequencing, MGMT gene dosage, chromatin status, and respective effects on MGMT expression and MGMT activity. Results: In glioblastoma, MGMT-methylated alleles ranged from 10% to 90%. In contrast, methylated alleles were highly enriched (100% of clones) in respective GS, even when 2 MGMT alleles were present, with 1 exception (<50%). The CpG methylation patterns were characteristic for each glioblastoma exhibiting 25% to 90% methylated CpGs of 28 sites interrogated. Furthermore, MGMT promoter methylation was associated with a nonpermissive chromatin status in accordance with very low MGMT transcript levels and undetectable MGMT activity. Conclusions: In MGMT-methylated glioblastoma, MGMT promoter methylation is highly enriched in GS that supposedly comprise glioma-initiating cells. Thus, even a low percentage of MGMT methylation measured in a glioblastoma sample may be relevant and predict benefit from an alkylating agent therapy. Clin Cancer Res; 17(2); 255-66. (C)2010 AACR.

Extent and Patterns of MGMT Promoter Methylation in Glioblastoma- and Respective Glioblastoma-Derived Spheres.

PURPOSE: Quantitative methylation-specific tests suggest that not all cells in a glioblastoma with detectable promoter methylation of the O6-methylguanine DNA methyltransferase (MGMT) gene carry a methylated MGMT allele. This observation may indicate cell subpopulations with distinct MGMT status, raising the question of the clinically relevant cutoff of MGMT methylation therapy. Epigenetic silencing of the MGMT gene by promoter methylation blunts repair of O6-methyl guanine and has been shown to be a predictive factor for benefit from alkylating agent therapy in glioblastoma. EXPERIMENTAL DESIGN: Ten paired samples of glioblastoma and respective glioblastoma-derived spheres (GS), cultured under stem cell conditions, were analyzed for the degree and pattern of MGMT promoter methylation by methylation-specific clone sequencing, MGMT gene dosage, chromatin status, and respective effects on MGMT expression and MGMT activity. RESULTS: In glioblastoma, MGMT-methylated alleles ranged from 10% to 90%. In contrast, methylated alleles were highly enriched (100% of clones) in respective GS, even when 2 MGMT alleles were present, with 1 exception (<50%). The CpG methylation patterns were characteristic for each glioblastoma exhibiting 25% to 90% methylated CpGs of 28 sites interrogated. Furthermore, MGMT promoter methylation was associated with a nonpermissive chromatin status in accordance with very low MGMT transcript levels and undetectable MGMT activity. CONCLUSIONS: In MGMT-methylated glioblastoma, MGMT promoter methylation is highly enriched in GS that supposedly comprise glioma-initiating cells. Thus, even a low percentage of MGMT methylation measured in a glioblastoma sample may be relevant and predict benefit from an alkylating agent therapy. Clin Cancer Res; 17(2); 255-66. ©2010 AACR.

Genome-wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits.

Nonalcoholic fatty liver disease (NAFLD) clusters in families, but the only known common genetic variants influencing risk are near PNPLA3. We sought to identify additional genetic variants influencing NAFLD using genome-wide association (GWA) analysis of computed tomography (CT) measured hepatic steatosis, a non-invasive measure of NAFLD, in large population based samples. Using variance components methods, we show that CT hepatic steatosis is heritable (∼26%-27%) in family-based Amish, Family Heart, and Framingham Heart Studies (n = 880 to 3,070). By carrying out a fixed-effects meta-analysis of genome-wide association (GWA) results between CT hepatic steatosis and ∼2.4 million imputed or genotyped SNPs in 7,176 individuals from the Old Order Amish, Age, Gene/Environment Susceptibility-Reykjavik study (AGES), Family Heart, and Framingham Heart Studies, we identify variants associated at genome-wide significant levels (p<5×10(-8)) in or near PNPLA3, NCAN, and PPP1R3B. We genotype these and 42 other top CT hepatic steatosis-associated SNPs in 592 subjects with biopsy-proven NAFLD from the NASH Clinical Research Network (NASH CRN). In comparisons with 1,405 healthy controls from the Myocardial Genetics Consortium (MIGen), we observe significant associations with histologic NAFLD at variants in or near NCAN, GCKR, LYPLAL1, and PNPLA3, but not PPP1R3B. Variants at these five loci exhibit distinct patterns of association with serum lipids, as well as glycemic and anthropometric traits. We identify common genetic variants influencing CT-assessed steatosis and risk of NAFLD. Hepatic steatosis associated variants are not uniformly associated with NASH/fibrosis or result in abnormalities in serum lipids or glycemic and anthropometric traits, suggesting genetic heterogeneity in the pathways influencing these traits.

Patterns of failure and outcome in patients with carcinoma of the anal margin.

BACKGROUND: To evaluate the outcome of patients with carcinoma of anal margin in terms of recurrence, survival, and radiation toxicity. METHODS: A series of 45 consecutive patients, with anal margin carcinoma treated between 1983 and 2006 with curative intent at two institutions, was retrospectively analyzed. A surgical excision (close or positive surgical margin in 22 out of 29 patients) was realized before radiotherapy (RT). RT consisted of definitive external beam RT (EBRT) in 36 patients, brachytherapy (BT) alone in two patients, and both BT and EBRT in seven patients. The median total radiation dose was 59.4 Gy (range, 30-74 Gy). RESULTS: The 5-year locoregional control (LRC) rate was 78% [95% confidence interval (CI), 64-93%]. The 5-year disease-specific survival (DSS) and overall survival (OS) rates were respectively 86% (95% CI, 72-99%) and 55% (95% CI, 44-66%). The overall anal conservation rate was 80% for the whole series. There was no significant association between local recurrence and patient age, histological grade, tumor size, T stage, overall treatment time, RT dose, or chemotherapy. Long-term side effects were observed in 15 patients (33%). Only three patients developed grade 3-4 late toxicity (CTCAE/NCI v3.0). Significant relationship was found between dose, and complication rate (48% for dose >or=59.4 Gy versus 8% for dose < 59.4 Gy; P = 0.03). CONCLUSIONS: We conclude that definitive RT and/or BT yield a good local control and disease-specific survival comparable with published data. This study suggests that radiation dose over 59.4 Gy seems to increase treatment-related morbidity.