Survival analysis of circulating tumor cells in triple-negative breast cancer

Background. Breast cancer is the most frequently diagnosed cancer and the leading cause of cancer death among females, accounting for 23% (1.38 million) of the total new cancer cases and 14% (458,400) of the total cancer deaths in 2008. [1] Triple-negative breast cancer (TNBC) is an aggressive phenotype comprising 10–20% of all breast cancers (BCs). [2-4] TNBCs show absence of estrogen, progesterone and HER2/neu receptors on the tumor cells. Because of the absence of these receptors, TNBCs are not candidates for targeted therapies. Circulating tumor cells (CTCs) are observed in blood of breast cancer patients even at early stages (Stage I & II) of the disease. Immunological and molecular analysis can be used to detect the presence of tumor cells in the blood (Circulating tumor cells; CTCs) of many breast cancer patients. These cells may explain relapses in early stage breast cancer patients even after adequate local control. CTC detection may be useful in identifying patients at risk for disease progression, and therapies targeting CTCs may improve outcome in patients harboring them. Methods . In this study we evaluated 80 patients with TNBC who are enrolled in a larger prospective study conducted at M D Anderson Cancer Center in order to determine whether the presence of circulating tumor cells is a significant prognostic factor in relapse free and overall survival . Patients with metastatic disease at the time of presentation were excluded from the study. CTCs were assessed using CellSearch System™ (Veridex, Raritan, NJ). CTCs were defined as nucleated cells lacking the presence of CD45 but expressing cytokeratins 8, 18 or 19. The distribution of patient and tumor characteristics was analyzed using chi square test and Fisher's exact test. Log rank test and Cox regression analysis was applied to establish the association of circulating tumor cells with relapse free and overall survival. Results. The median age of the study participants was 53years. The median duration of follow-up was 40 months. Eighty-eight percent (88%) of patients were newly diagnosed (without a previous history of breast cancer), and (60%) of patients were chemo naïve (had not received chemotherapy at the time of their blood draw for CTC analysis). Tumor characteristics such as stage (P=0.40), tumor size (P=69), sentinel nodal involvement (P=0.87), axillary lymph node involvement (P=0.13), adjuvant therapy (P=0.83), and high histological grade of tumor (P=0.26) did not predict the presence of CTCs. However, CTCs predicted worse relapse free survival (1 or more CTCs log rank P value = 0.04, at 2 or more CTCs P = 0.02 and at 3 or more CTCs P < 0.0001) and overall survival (at 1 or more CTCs log rank P value = 0.08, at 2 or more CTCs P = 0.01 and at 3 or more CTCs P = 0.0001. Conclusions. The number of circulating tumor cells predicted worse relapse free survival and overall survival in TNBC patients.^

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

Retrospective analysis of Methicillin-Resistant Staphylococcus Aureus (MRSA) skin and Soft-Tissue Infections (SSTI) in a university hospital emergency department (ED)

Study Objective: Identify the most frequent risk factors of Community Acquired-MRSA (CA-MRSA) Skin and Soft-tissue Infections (SSTIs) using a case series of patients and characterize them by age, race/ethnicity, gender, abscess location, druguse and intravenous drug-user (IVDU), underlying medical conditions, homelessness, treatment resistance, sepsis, those whose last healthcare visit was within the last 12 months, and describe the susceptibility pattern from this central Texas population that have come into the University Medical Center Brackenridge (UMCB) Emergency Department (ED). ^ Methods: This study was a retrospective case-series medical record review involving a convenience sample of patients in 2007 from an urban public hospital's ED in Texas that had a SSTI that tested positive for MRSA. All positive MRSA cultures underwent susceptibility testing to determine antibiotic resistance. The demographic and clinical variables that were independently associated with MRSA were determined by univariate and multivariate analysis using logistic regression to calculate odds ratios (OR), 95% confidence intervals, and significance (p≤ 0.05). ^ Results: In 2007, there were 857 positive MRSA cultures. The demographics were: males 60% and females 40%, with the average age of 36.2 (std. dev. =13) the study population consisted of non-Hispanic white (42%), Hispanics (38%), and non-Hispanic black (18.8%). Possible risk factors addressed included using recreational drugs (not including IVDU) (27%) homelessness (13%), diabetes status (12.6%) or having an infectious disease, and IVDU (10%). The most frequent abscess location was the leg (26.6%), followed by the arm and torso (both 13.7%). Eighty-three percent of patients had one prominent susceptibility pattern that had a susceptibility rate for the following antibiotics: trimethoprim/sulfamethoxazole (TMP-SMX) and vancomycin had 100%, gentamicin 99%, clindamycin 96%, tetracycline 96%, and erythromycin 56%. ^ Conclusion: The ED is becoming an important area for disease transmission between the sterile hospital environment and the outside environment. As always, it is important to further research in the ED in an effort to better understand MRSA transmission and antibiotic resistance, as well as to keep surveillance for the introduction of new opportunistic pathogens into the population. ^

Patenting genes and gene sequences: Analysis of public health implications and the law

At issue is whether or not isolated DNA is patent eligible under the U.S. Patent Law and the implications of that determination on public health. The U.S. Patent and Trademark Office has issued patents on DNA since the 1980s, and scientists and researchers have proceeded under that milieu since that time. Today, genetic research and testing related to the human breast cancer genes BRCA1 and BRCA2 is conducted within the framework of seven patents that were issued to Myriad Genetics and the University of Utah Research Foundation between 1997 and 2000. In 2009, suit was filed on behalf of multiple researchers, professional associations and others to invalidate fifteen of the claims underlying those patents. The Court of Appeals for the Federal Circuit, which hears patent cases, has invalidated claims for analyzing and comparing isolated DNA but has upheld claims to isolated DNA. The specific issue of whether isolated DNA is patent eligible is now before the Supreme Court, which is expected to decide the case by year's end. In this work, a systematic review was performed to determine the effects of DNA patents on various stakeholders and, ultimately, on public health; and to provide a legal analysis of the patent eligibility of isolated DNA and the likely outcome of the Supreme Court's decision. ^ A literature review was conducted to: first, identify principle stakeholders with an interest in patent eligibility of the isolated DNA sequences BRCA1 and BRCA2; and second, determine the effect of the case on those stakeholders. Published reports that addressed gene patents, the Myriad litigation, and implications of gene patents on stakeholders were included. Next, an in-depth legal analysis of the patent eligibility of isolated DNA and methods for analyzing it was performed pursuant to accepted methods of legal research and analysis based on legal briefs, federal law and jurisprudence, scholarly works and standard practice legal analysis. ^ Biotechnology, biomedical and clinical research, access to health care, and personalized medicine were identified as the principle stakeholders and interests herein. Many experts believe that the patent eligibility of isolated DNA will not greatly affect the biotechnology industry insofar as genetic testing is concerned; unlike for therapeutics, genetic testing does not require tremendous resources or lead time. The actual impact on biomedical researchers is uncertain, with greater impact expected for researchers whose work is intended for commercial purposes (versus basic science). The impact on access to health care has been surprisingly difficult to assess; while invalidating gene patents might be expected to decrease the cost of genetic testing and improve access to more laboratories and physicians' offices that provide the test, a 2010 study on the actual impact was inconclusive. As for personalized medicine, many experts believe that the availability of personalized medicine is ultimately a public policy issue for Congress, not the courts. ^ Based on the legal analysis performed in this work, this writer believes the Supreme Court is likely to invalidate patents on isolated DNA whose sequences are found in nature, because these gene sequences are a basic tool of scientific and technologic work and patents on isolated DNA would unduly inhibit their future use. Patents on complementary DNA (cDNA) are expected to stand, however, based on the human intervention required to craft cDNA and the product's distinction from the DNA found in nature. ^ In the end, the solution as to how to address gene patents may lie not in jurisprudence but in a fundamental change in business practices to provide expanded licenses to better address the interests of the several stakeholders. ^

Analysis of p53 function and regulation

p53 is a tumor suppressor gene that is the most frequent target inactivated in cancers. Overexpression of wild-type p53 in rat embryo fibroblasts suppresses foci formation by other cooperating oncogenes. Introduction of wild-type p53 into cells that lack p53 arrests them at the G1/S boundary and reverses the transformed phenotype of some cells. The function of p53 in normal cells is illustrated by the ability of p53 to arrest cells at G1 phase of the cell cycle upon exposure to DNA-damaging agents including UV-irradiation and biosynthesis inhibitors.^ Since the amino acid sequence of p53 suggested that it may function as a transcription factor, we used GAL4 fusion assays to test that possibility. We found that wild-type p53 could specifically activate transcription when anchored by the GAL4 DNA binding domain. Mutant p53s, which have lost the ability to suppress foci formation by other oncogenes, were not able to activate transcription in this assay. Thus, we established a direct correlation between the tumor suppression and transactivation functions of p53.^ Having learned that p53 was a transcriptional activator, we next sought targets of p53 activation. Because many transcription factors regulate their own expression, we tested whether p53 had this autoregulatory property. Transient expression of wild-type p53 in cells increased the levels of endogenous p53 mRNA. Cotransfection of p53 together with a reporter bearing the p53 promoter confirmed that wild-type p53 specifically activates its own promoter. Deletion analysis from both the 5$\sp\prime$ and 3$\sp\prime$ ends of the promoter minimized the region responsible for p53 autoregulation to 45 bp. Methylation interference identified nucleotides involved in protein-DNA interaction. Mutations within this protected site specifically eliminated the response of the promoter to p53. In addition, multiple copies of this element confer responsiveness to wild-type p53 expression. Thus, we identified a F53 responsive element within the p53 promoter.^ The presence of a consensus NF-$\kappa$B site in the p53 promoter suggested that NF-KB may regulate p53 expression. Gel-shift experiments showed that both the p50 homodimer and the p50/p65 heterodimer bind to the p53 promoter. In addition, the p65 subunit of NF-$\kappa$B activates the p53 promoter in transient transfection experiments. TNF $\alpha$, a natural NF-$\kappa$B inducer, also activates the p53 promoter. Both p65 activation and TNF $\alpha$ induction require an intact NF-$\kappa$B site in the p53 promoter. Since NF-$\kappa$B activation occurs as a response to stress and p53 arrests cells in G1/S, where DNA repair occurs, activation of p53 by NF-$\kappa$B could be a mechanism by which cells recover from stress.^ In conclusion, we provided the first data that wild-type p53 functions as a transcriptional activator, whereas mutant p53 cannot. The correlation between growth suppression and transcriptional activation by p53 implies a pathway of tumor suppression. We have analyzed upstream components of the pathway by the identification of both p53 and NF-$\kappa$B as regulators of the p53 promoter. ^

A comparative analysis of hospital -based physician monitoring systems

The purpose of the multiple case-study was to determine how hospital subsystems (such as physician monitoring and credentialing; quality assurance; risk management; and peer review) were supporting the monitoring of physicians? Three large metropolitan hospitals in Texas were studied and designated as hospitals #1, #2, and #3. Realizing that hospital subsystems are a unique entity and part of a larger system, conclusions were made on the premises of a quality control system, in relation to the tools of government (particularly the Health Care Quality Improvement Act (HCQIA)), and in relation to itself as a tool of a hospital.^ Three major analytical assessments were performed. First, the subsystems were analyzed as to their "completeness"; secondly, the subsystems were analyzed for "performance"; and thirdly, the subsystems were analyzed in reference to the interaction of completeness and performance.^ The physician credentialing and monitoring and the peer review subsystems as quality control systems were most complete, efficient, and effective in hospitals #1 and #3. The HCQIA did not seem to be an influencing factor in the completeness of the subsystem in hospital #1. The quality assurance and risk management subsystem in hospital #2 was not representative of completeness and performance and the HCQIA was not an influencing factor in the completeness of the Q.A. or R.M. systems in any hospital. The efficiency (computerization) of the physician credentialing, quality assurance and peer review subsystems in hospitals #1 and #3 seemed to contribute to their effectiveness (system-wide effect).^ The results indicated that the more complete, effective, and efficient subsystems were characterized by (1) all defined activities being met, (2) the HCQIA being an influencing factor, (3) a decentralized administrative structure, (4) computerization an important element, and (5) staff was sophisticated in subsystem operations. However, other variables were identified which deserve further research as to their effect on completeness and performance of subsystems. They include (1) medical staff affiliations, (2) system funding levels, (3) the system's administrative structure, and (4) the physician staff "cultural" characteristics. Perhaps by understanding other influencing factors, health care administrators may plan subsystems that will be compatible with legislative requirements and administrative objectives. ^

The cDNA cloning of murine HIP /L29 and the functional analysis of HIP/L29 domains

Human heparin/heparan sulfate interacting protein/L29 (HIP/L29) is a heparin/heparan sulfate (Hp/HS) binding protein found in many adult human tissues. Potential functions of this protein are promotion of embryo adhesion, modulation of blood coagulation, and control of cell growth. While these activities are diverse, the ability of human HIP/L29 to interact with Hp/HS at the cell surface may be a unifying mechanism of action since Hp/HS influences all of these processes. A murine ortholog has been identified that has 78.8% homology over the entire sequence and identity over the N-terminal 64 amino acids when compared to human HIP/L29. Northern, Western, and immunohistochemical analysis shows that murine HIP/L29 mRNA and protein are expressed in a tissue specific manner. Murine HIP/L29 is enriched in the membrane fraction of NmuMG cells where it is eluted with high salt, suggesting that it is a peripheral membrane protein. The ability of murine HIP/L29 to bind Hp is verified by studies using native and recombinant forms of murine HIP/L29. A synthetic peptide (HIP peptide-2) derived from the identical N-terminal region of HIP/L29 proteins was tested for the ability to bind Hp and support cell adhesion. This peptide was chosen because it conforms to a proposed consensus sequence for Hp/HS binding peptides. HIP peptide-2 binds Hp in a dose-dependent, saturable, and selective manner and supports Hp-dependent cell adhesion. However, a scrambled form of this peptide displayed similar activities indicating a lack of peptide sequence specificity required for activity. Lastly, an unbiased approach was used to identify sequences within human and mouse HIP/L29 proteins necessary for Hp/HS binding. A panel of recombinant proteins was made that collectively are deficient in every human HIP/L29 domain. The activities of these deletion mutants and recombinant murine HIP/L29 were compared to the activity of recombinant human HIP/L29 in a number of assays designed to look at differences in the ability to bind Hp/HS. These studies suggest that each domain within human HIP/L29 is important for binding to Hp/HS and divergences in the C-terminus of human and mouse HIP/L29 account for a decrease in murine HIP/L29 affinity for Hp/HS. It is apparent that multiple domains within human and mouse HIP/L29 contribute to the function of Hp/HS binding. The interaction of multiple HIP/L29 domains with Hp/HS will influence the biological activity of HIP/L29 proteins. ^

Functional analysis of GNA -1, a Galphai superfamily member found in the filamentous fungus Neurospora crassa

Heterotrimeric GTP-binding proteins, G proteins, are integral components of eukaryotic signaling systems linking extracellular signals to intracellular responses. Through coupling to seven-transmembrane helix receptors, G proteins convey primary signaling events into multi-leveled cascades of intracellular activity by regulating downstream enzymes, collectively called effectors. The effector enzymes regulated by G proteins include adenylyl cyclase, cAMP phosphodiesterase, phospolipase C-β, mitogen-activated protein kinases, and ion channels. ^ Neurospora crassa is a multicellular, filamentous fungus that is capable of both asexual and sexual reproduction by elaboration of specialized, developmentally controlled structures that give rise to either asexual or sexual spores, respectively. N. crassa possesses at least three heterotrimeric Gα proteins (GNA-1–3) and one Gβ subunit (GNB-1). GNA-1 was the first microbial protein that could be classified in the Gαi superfamily based on its amino acid identity and demonstration that it is a substrate for ADP-ribosylation by pertussis toxin. ^ Experiments were designed to identify the signal transduction pathways and the effector enzymes regulated by GNA-1. Targeted gene-replacement of gna-1 revealed that GNA-1 controls multiple developmental pathways including both asexual and sexual reproduction, maintenance of growth, and resistance to osmotic stress. The Gαi and Gαz members of the Gαi superfamily negatively regulate adenylyl cyclase activity in mammalian cells; therefore, adenylyl cyclase and cAMP levels were measured in Δgna-1 strains and also in strains that were deleted for both gna-1 and gna-2, a second Gα in N. crassa shown to have overlapping functions with GNA-1. Direct measurements of adenylyl cyclase activity revealed that GNA-1, but not GNA-2, was responsible for GTP-stimulated adenylyl cyclase activity in N. crassa. Furthermore, anti-GNA-1 IgG could specifically inhibit GTP-stimulated adenylyl cyclase activity in wild-type strain extracts. These studies also provided evidence that N. crassa possesses feedback mechanisms that control steady-state cAMP levels through indirect regulation of cAMP-phosphodiesterase activity; mutations in gna-1 and gna-2 were additive in their effect on lowering cAMP-phosphodiesterase activity under growth conditions where steady-state cAMP levels were normal but GTP-stimulated adenylyl cyclase activity was reduced 90% in comparison to control strains. ^ Genetic and biochemical epistasis experiments utilizing a Δ gna-1 cr-1 mutant suggest that GNA-1 is essential for female fertility in a cAMP-independent pathway. Furthermore, deletion of gna-1 in a cr-1 background exacerbated many of the defects already observed in the cr-1 strain including more severe growth restriction and developmental defects. However, deletion of gna-1 had no effect on the increased thermotolerance of cr-1, which has been attributed to loss of cAMP. cr-1 possesses GNA-1 protein, and crude membrane fractions from this strain reconstituted GTP-stimulated adenylyl cyclase activity in Δgna-1 membrane fractions. These studies provide direct evidence for the involvement of Gα proteins in the regulation of adenylyl cyclase activity in eukaryotic microbes. ^

Functional analysis of SOX9 in chondrogenesis by targeted gene inactivation

Sox9 is a Sry-related HMG-domain containing transcription factor. Lines of evidence suggest that Sox9 has a potential role in skeletal development. During mouse development, Sox9 is predominantly expressed in all chondroprogenitors and differentiated chondrocytes, throughout the deposition of cartilage matrix. Mutations in one allele of SOX9 in humans result in campomelic dysplasia (CD), a skeletal dysplasia. syndrome characterized by the bowing of long bones. Moreover, Sox9 binds to and activates chondrocyte-specific enhancers in Col2a1 and Col11a2 genes. To further investigate the function of Sox9 in chondrogenesis, we analyzed chimeras derived from Sox9 heterozygous and homozygous null embryonic stem (ES) cells. In mouse chimeras, Sox9 −/− cells were excluded from all cartilages and did not express chondrocyte-specific genes. The segregation occurred during mesenchymal condensation. No cartilages developed in teratocarcinomas derived from Sox9 −/− ES cells. Mice heterozygous for the Sox9 mutation died neonatally and exhibited skeletal abnormalities resembling those of the CD patients. The Sox9 +/− mutants had a cleft palate and hypoplasia of scapula, pelvis and other skeletal structures derived by endochondral ossification. Bending of the radius, ulna and tibia cartilage was prominent at embryonic day 14.5 (E14.5). At E12.5 many pre-cartilaginous condensations were already defective. Advanced ossification was observed and the hypertrophic zone was enlarged in the growth plates, suggesting that Sox9 also regulates hypertrophic chondrocyte differentiation. Our results identify Sox9 as the first essential regulator of chondrocyte differentiation, which plays multiple roles in chondrogenesis. ^

Separating ITCZ- and ENSO-related rainfall changes in the Galápagos over the last 3 kyr using D/H ratios of multiple lipid biomarkers

We present a 3000-yr rainfall reconstruction from the Galápagos Islands that is based on paired biomarker records from the sediment of El Junco Lake. Located in the eastern equatorial Pacific, the climate of the Galápagos Islands is governed by movements of the Intertropical Convergence Zone (ITCZ) and the El Niño-Southern Oscillation (ENSO). We use a novel method for reconstructing past ENSO- and ITCZ-related rainfall changes through analysis of molecular and isotopic biomarker records representing several types of plants and algae that grow under differing climatic conditions. We propose that ?D values of dinosterol, a sterol produced by dinoflagellates, record changes in mean rainfall in El Junco Lake, while dD values of C34 botryococcene, a hydrocarbon unique to the green alga Botryococcus braunii, record changes in rainfall associated with moderate-to-strong El Niño events. We use these proxies to infer changes in mean rainfall and El Niño-related rainfall over the past 3000 yr. During periods in which the inferred change in El Niño-related rainfall opposed the change in mean rainfall, we infer changes in the amount of ITCZ-related rainfall. Simulations with an idealized isotope hydrology model of El Junco Lake help illustrate the interpretation of these proxy reconstructions. Opposing changes in El Niño- and ITCZ-related rainfall appear to account for several of the largest inferred hydrologic changes in El Junco Lake. We propose that these reconstructions can be used to infer changes in frequency and/or intensity of El Niño events and changes in the position of the ITCZ in the eastern equatorial Pacific over the past 3000 yr. Comparison with El Junco Lake sediment grain size records indicates general agreement of inferred rainfall changes over the late Holocene.

Calcareous dinoflagellate cysts analysis of sediment cores from the Equatorial South Atlantic Ocean

An overview is presented of the current state of knowledge on paleo-ecological aspects of calcareous dinoflagellate resting cysts. Apart from literature-based information, a discussion of new results is also provided from Equatorial Atlantic surface plankton samples, surface sediment samples and Late Quaternary sediments from two gravity cores. With the aid of redundancy analysis statistics, variations in the calcareous cyst content of both cores are correlated to variations in total organic carbon (TOC). On a global scale, the calcareous cyst distribution in bottom sediments varies with latitude and inshore-offshore gradients. In the Equatorial Atlantic Ocean, enhanced calcareous cyst production can be observed in regions and time intervals with stratified, oligotrophic conditions in the upper water masses.

Geospatial analysis of groundwater depth for manual well drilling suitability in the Zinder Region, Niger

Manual and low-tech well drilling techniques have potential to assist in reaching the United Nations' millennium development goal for water in sub-Saharan Africa. This study used publicly available geospatial data in a regression tree analysis to predict groundwater depth in the Zinder region of Niger to identify suitable areas for manual well drilling. Regression trees were developed and tested on a database for 3681 wells in the Zinder region. A tree with 17 terminal leaves provided a range of ground water depth estimates that were appropriate for manual drilling, though much of the tree's complexity was associated with depths that were beyond manual methods. A natural log transformation of groundwater depth was tested to see if rescaling dataset variance would result in finer distinctions for regions of shallow groundwater. The RMSE for a log-transformed tree with only 10 terminal leaves was almost half that of the untransformed 17 leaf tree for groundwater depths less than 10 m. This analysis indicated important groundwater relationships for commonly available maps of geology, soils, elevation, and enhanced vegetation index from the MODIS satellite imaging system.

Sedimentological and palynological analysis of sediment cores from the northeastern Vietnamese Mekong River Delta

Three radiocarbon-dated sediment cores from the northeastern Vietnamese Mekong River Delta have been analysed with a multiproxy approach (grain size, pollen and spores, macro-charcoal, carbon content) to unravel the palaeoenvironmental history of the region since the mid Holocene. During the mid-Holocene sea-level highstand a diverse, zoned and widespread mangrove belt (dominated by Rhizophora) covered the extended tidal flats. The subsequent regression and coeval delta progradation led to the rapid development of a back-mangrove community dominated by Ceriops and Bruguiera but also represented locally by e.g. Kandelia, Excoecaria and Phoenix. Along rivers this community seems to have endured even when the adjoining floodplain had already shifted to freshwater vegetation. Generally this freshwater vegetation has a strong swamp signature but locally Arecaceae, Fabaceae, Moraceae/Urticaceae and Myrsinaceae are important and mirror the geomorphological diversity of the delta plain. The macro-charcoal record implies that natural burning of vegetation occurred throughout the records, however, the occurrence of the highest amounts of macro-charcoal particles is linked with modern human activity.