PML tumor suppressor potentiates cell death through inhibition of the NF -kappa B survival pathway

The promyelocytic leukemia protein PML is a growth suppressor essential for induction of apoptosis by diverse apoptotic stimuli. The mechanism by which PML regulates cell death remains unclear. In this study we found that ectopic expression of PML potentiates cell death in the TNFα-resistant tumor line U2OS and significantly sensitized these cells to apoptosis induced by TNFα in a p53-independent manner. Our study demonstrated that both PML and PML/TNFα-induced cell death are associated with DNA fragmentation, activation of caspase-3, -7, -8, and degradation of DFF/ICAD. Furthermore, we found that PML-induced and PML/TNFα-induced cell death could be blocked by the caspase-8 inhibitors crmA and c-FLIP, but not by Bcl-2, the inhibitor of mitochondria-mediated apoptotic pathway. These findings indicate that this cell death event is initiated through the death receptor-dependent apoptosis pathway. Our study further showed that PML recruits NF-kappa B (NF-κB) to the PML nuclear body, blocks NF-κB binding to its cognate enhancer, and represses its transactivation function with the C-terminal region. Therefore PML inhibits the NF-κB survival pathway. Overexpression of NF-κB rescued cell death induced by PML and PML/TNFκ. These results imply that PML is a functional repressor of NF-κB. This notion was further supported by the finding that the PML−/− mouse embryo fibroblasts (MEFs) are more resistant than the wild-type MEFs to TNFκ-induced apoptosis. In conclusion, our studies convincingly demonstrated that PML potentiates cell death through inhibition of the NF-κB survival pathway. Activation of NF-κB frequently occurs during oncogenesis. Our study here suggests that a loss of PML function enhances the NF-κB survival pathway and this event may contribute to tumorigenesis. ^

Pollen and non-pollen palynomorph records of the Holocene part of the composite sediment core TMD from lake Tso Moriri (NW India) and of modern surface samples from the Tso Moriri region

The high-altitude lake Tso Moriri (32°55'46'' N, 78°19'24'' E; 4522 m a.s.l.) is situated at the margin of the ISM and westerly influences in the Trans-Himalayan region of Ladakh. Human settlements are rare and domestic and wild animals are concentrating at the alpine meadows. A set of modern surface samples and fossil pollen from deep-water TMD core was evaluated with a focus on indicator types revealing human impact, grazing activities and lake system development during the last ca. 12 cal ka BP. Furthermore, the non-pollen palynomorph (NPP) record, comprising remains of limnic algae and invertebrates as well as fungal spores and charred plant tissue fragments, were examined in order to attest palaeolimnic phases and human impact, respectively. Changes in the early and middle Holocene limnic environment are mainly influenced by regional climatic conditions and glacier-fed meltwater flow in the catchment area. The NPP record indicates low lake productivity with high influx of freshwater between ca. 11.5 and 4.5 cal ka BP which is in agreement with the regional monsoon dynamics and published climate reconstructions. Geomorphologic observations suggest that during this period of enhanced precipitation the lake had a regular outflow and contributed large amounts of water to the Sutlej River, the lower reaches of which were integral part of the Indus Civilization area. The inferred minimum fresh water input and maximum lake productivity between ca. 4.5-1.8 cal ka BP coincides with the reconstruction of greatest aridity and glaciation in the Korzong valley resulting in significantly reduced or even ceased outflow. We suggest that lowered lake levels and river discharge on a larger regional scale may have caused irrigation problems and harvest losses in the Indus valley and lowlands occupied by sedentary agricultural communities. This scenario, in turn, supports the theory that, Mature Harappan urbanism (ca. 4.5-3.9 cal ka BP) emerged in order to facilitate storage, protection, administration, and redistribution of crop yields and secondly, the eventual collapse of the Harappan Culture (ca. 3.5-3 cal ka BP) was promoted by prolonged aridity. There is no clear evidence for human impact around Tso Moriri prior to ca. 3.7 cal ka BP, with a more distinct record since ca. 2.7 cal ka BP. This suggests that the sedimentary record from Tso Moriri primarily archives the regional climate history.

Anaerobic oxidation of methane and sulphate reduction rates, hydrocarbons, sulphate and sulphide concentrations of sediment cores from the Larsen B embayment at the eastern site of the Antarctic Peninsula

First videographic indication of an Antarctic cold seep ecosystem was recently obtained from the collapsed Larsen B ice shelf, western Weddell Sea (Domack et al., 2005). Within the framework of the R/V Polarstern expedition ANTXXIII-8, we revisited this area for geochemical, microbiological and further videographical examinations. During two dives with ROV Cherokee (MARUM, Bremen), several bivalve shell agglomerations of the seep-associated, chemosynthetic clam Calyptogena sp. were found in the trough of the Crane and Evans glacier. The absence of living clam specimens indicates that the flux of sulphide and hence the seepage activity is diminished at present. This impression was further substantiated by our geochemical observations. Concentrations of thermogenic methane were moderately elevated with 2 µM in surface sediments of a clam patch, increasing up to 9 µM at a sediment depth of about 1 m in the bottom sections of the sediment cores. This correlated with a moderate decrease in sulphate from about 28 mM at the surface down to 23.4 mM, an increase in sulphide to up to 1.43 mM and elevated rates of the anaerobic oxidation of methane (AOM) of up to 600 pmol cm**-3 d**-1 at about 1 m below the seafloor. Molecular analyses indicate that methanotrophic archaea related to ANME-3 are the most likely candidates mediating AOM in sediments of the Larsen B seep.

Nonmuscle myosin II-B is required for normal development of the mouse heart

We used targeted gene disruption in mice to ablate nonmuscle myosin heavy chain B (NMHC-B), one of the two isoforms of nonmuscle myosin II present in all vertebrate cells. Approximately 65% of the NMHC-B−/− embryos died prior to birth, and those that were born suffered from congestive heart failure and died during the first day. No abnormalities were detected in NMHC-B+/− mice. The absence of NMHC-B resulted in a significant increase in the transverse diameters of the cardiac myocytes from 7.8 ± 1.8 μm (right ventricle) and 7.8 ± 1.3 μm (left ventricle) in NMHC-B+/+ and B+/− mice to 14.7 ± 1.1 μm and 13.8 ± 2.3 μm, respectively, in NMHC-B−/− mice (in both cases, P < 0.001). The increase in size of the cardiac myocytes was seen as early as embryonic day 12.5 (4.5 ± 0.2 μm for NMHC-B+/+ and B+/− vs. 7.2 ± 0.6 μm for NMHC-B−/− mice (P < 0.01)). Six of seven NMHC-B−/− newborn mice analyzed by serial sectioning also showed structural cardiac defects, including a ventricular septal defect, an aortic root that either straddled the defect or originated from the right ventricle, and muscular obstruction to right ventricular outflow. Some of the hearts of NMHC-B−/− mice showed evidence for up-regulation of NMHC-A protein. These studies suggest that nonmuscle myosin II-B is required for normal cardiac myocyte development and that its absence results in structural defects resembling, in part, two common human congenital heart diseases, tetralogy of Fallot and double outlet right ventricle.

Molecular cloning of a rat chromosome putative recombinogenic sequence homologous to the hepatitis B virus encapsidation signal.

Previously, we reported that a 61-bp subgenomic HBV DNA sequence (designated as 15AB, nt 1855-1915) is a hot spot for genomic recombination and that a cellular protein binding to 15AB may be the putative recombinogenic protein. In the present study, we established the existence of a 15AB-like sequence in human and rat chromosomal DNA by Southern blot analysis. The 15AB-like sequence isolated from the rat chromosome demonstrated a 80.9% identity with 5'-CCAAGCTGTGCCTTGGGTGGC-3', at 1872-1892 of the hepatitis B virus genome, thought to be the essential region for recombination. Interestingly, this 15AB-like sequence also contained the pentanucleotide motifs GCTGG and CCAGC as an inverted repeat, part of the chi known hot spot for recombination in Escherichia coli. Importantly, a portion of the 15AB-like sequence is homologous (82.1%, 23/28 bp) to break point clusters of the human promyelocytic leukemia (PML) gene, characterized by a translocation [t(15;17)], and to rearranged mouse DNA for the immunoglobulin kappa light chain. Moreover, 15AB and 15AB-like sequences have striking homologies (12/15 = 80.0% and 13/15 = 86.7%, respectively) to the consensus sequence for topoisomerase II. Our present results suggest that this 15AB-like sequence in the rat genome might be a recombinogenic candidate triggering genomic instability in carcinogenesis.

Targeted disruption of the surfactant protein B gene disrupts surfactant homeostasis, causing respiratory failure in newborn mice.

Surfactant protein B (SP-B) is an 8.7-kDa, hydrophobic protein that enhances the spreading and stability of surfactant phospholipids in the alveolus. To further assess the role of SP-B in lung function, the SP-B gene was disrupted by homologous recombination in murine mouse embryonic stem cells. Mice with a single mutated SP-B allele (+/-) were unaffected, whereas homozygous SP-B -/- offspring died of respiratory failure immediately after birth. Lungs of SP-B -/- mice developed normally but remained atelectatic in spite of postnatal respiratory efforts. SP-B protein and mRNA were undetectable and tubular myelin figures were lacking in SP-B -/- mice. Type II cells of SP-B -/- mice contained no fully formed lamellar bodies. While the abundance of SP-A and SP-C mRNAs was not altered, an aberrant form of pro-SP-C, 8.5 kDa, was detected, and fully processed SP-C peptide was markedly decreased in lung homogenates of SP-B -/- mice. Ablation of the SP-B gene disrupts the routing, storage, and function of surfactant phospholipids and proteins, causing respiratory failure at birth.

Drosophila hormone receptor 38: a second partner for Drosophila USP suggests an unexpected role for nuclear receptors of the nerve growth factor-induced protein B type.

In Drosophila the response to the hormone ecdysone is mediated in part by Ultraspiracle (USP) and ecdysone receptor (EcR), which are members of the nuclear receptor superfamily. Heterodimers of these proteins bind to ecdysone response elements (EcREs) and ecdysone to modulate transcription. Herein we describe Drosophila hormone receptor 38 (DHR38) and Bombyx hormone receptor 38 (BHR38), two insect homologues of rat nerve growth factor-induced protein B (NGFI-B). Although members of the NGFI-B family are thought to function exclusively as monomers, we show that DHR38 and BHR38 in fact interact strongly with USP and that this interaction is evolutionarily conserved. DHR38 can compete in vitro against EcR for dimerization with USP and consequently disrupt EcR-USP binding to an EcRE. Moreover, transfection experiments in Schneider cells show that DHR38 can affect ecdysone-dependent transcription. This suggests that DHR38 plays a role in the ecdysone response and that more generally NGFI-B type receptors may be able to function as heterodimers with retinoid X receptor type receptors in regulating transcription.

DNA-mediated immunization to the hepatitis B surface antigen in mice: aspects of the humoral response mimic hepatitis B viral infection in humans.

Intramuscular injection of plasmid DNA expression vectors encoding the three envelope proteins of the hepatitis B virus (HBV) induced humoral responses in C57BL/6 mice specific to several antigenic determinants of the viral envelope. The first antibodies appeared within 1-2 weeks after injection of DNA and included antibodies of the IgM isotype. Over the next few weeks, an IgM to IgG class switch occurred, indicating helper T-lymphocyte activity. Peak IgG titers were reached by 4-8 weeks after a single DNA injection and were maintained for at least 6 months without further DNA injections. The antibodies to the envelope proteins reacted with group- and subtype-specific antigenic determinants of the HBV surface antigen (HBsAg). Expression vectors encoding the major (S) and middle (preS2 plus S) envelope proteins induced antibodies specific to the S protein and preS2 domain, and preS2 antibodies were prominent at early time points. In general, the expression vectors induced humoral responses in mice that mimic those observed in humans during the course of natural HBV infection.

The VLT-FLAMES Tarantula Survey: XVIII. Classifications and radial velocities of the B-type stars

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 ± 12.2 km s^-1 from 273 presumed single stars. Employing a 3σ criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (≥345 km s^-1) or small (≤65 km s^-1) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 ± 16.2 km s^-1), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12 + log (N/H) ≳ 8.5. Combined with a large rotational velocity (υb>eb> sin i = 345 ± 22 km s^-1), this is suggestive of past binary interaction for this star.

The large- and small-scale Ca IIK structure of the Milky Way from observations of Galactic and Magellanic sightlines

Aims. The large and small-scale (pc) structure of the Galactic interstellar medium can be investigated by utilising spectra of early-type stellar probes of known distances in the same region of the sky. This paper determines the variation in line strength of Ca ii at 3933.661 Å as a function of probe separation for a large sample of stars, including a number of sightlines in the Magellanic Clouds. 

Methods. FLAMES-GIRAFFE data taken with the Very Large Telescope towards early-type stars in 3 Galactic and 4 Magellanic open clusters in Ca ii are used to obtain the velocity, equivalent width, column density, and line width of interstellar Galactic calcium for a total of 657 stars, of which 443 are Magellanic Cloud sightlines. In each cluster there are between 43 and 111 stars observed. Additionally, FEROS and UVES Ca ii K and Na i D spectra of 21 Galactic and 154 Magellanic early-type stars are presented and combined with data from the literature to study the calcium column density - parallax relationship. 

Results. For the four Magellanic clusters studied with FLAMES, the strength of the Galactic interstellar Ca ii K equivalent width on transverse scales from ∼0.05-9 pc is found to vary by factors of ∼1.8-3.0, corresponding to column density variations of ∼0.3-0.5 dex in the optically-thin approximation. Using FLAMES, FEROS, and UVES archive spectra, the minimum and maximum reduced equivalent widths for Milky Way gas are found to lie in the range ∼35-125 mÅ and ∼30-160 mÅ for Ca ii K and Na i D, respectively. The range is consistent with a previously published simple model of the interstellar medium consisting of spherical cloudlets of filling factor ∼0.3, although other geometries are not ruled out. Finally, the derived functional form for parallax (π) and Ca ii column density (Nb>CaIIb>) is found to be π(mas) = 1 / (2.39 × 10^-13 × Nb>CaIIb> (cm^-2) + 0.11). Our derived parallax is ∼25 per cent lower than predicted by Megier et al. (2009, A&A, 507, 833) at a distance of ∼100 pc and ∼15 percent lower at a distance of ∼200 pc, reflecting inhomogeneity in the Ca ii distribution in the different sightlines studied.

Expression of high amounts of the CD117 molecule in a case of B-cell non-Hodgkin's lymphoma carrying the t(14:18) translocation.

The c-kit proto-oncogen (CD117) has been described to be present in normal and neoplastic hemopoietic cells including both myeloid and lymphoid lineages. Among the normal lymphoid cells CD117 expression would be restricted to a small subset of NK-cells, and to early T-cell precursors and it is not expressed by normal B-cells. Regarding chronic lymphoproliferative disorders the only data provided up to now suggests that CD117 expression is restricted to cases of Hodgkin's disease and anaplastic large-cell lymphoma. In the present paper we describe a case of a B-cell chronic lymphoproliferative disorder carrying the t(14:18) translocation as demonstrated by molecular studies, in which the flow cytometric immunophenotypic analysis of both peripheral blood and bone marrow samples revealed the expression of high amounts of the CD117 antigen in the surface of the clonal B-cell population. Further studies are necessary to explore both the functional role of c-kit expression in the neoplastic B-cells from this patient and its potential utility for the diagnosis and follow-up of patients with B-cell non-Hodgkin's lymphoma.

Evaluation of the Suchey-Brooks Method of age estimation in an Australian subpopulation using computed tomography of the pubic symphyseal surface

Despite the prominent use of the Suchey-Brooks (S-B) method of age estimation in forensic anthropological practice, it is subject to intrinsic limitations, with reports of differential inter-population error rates between geographical locations. This study assessed the accuracy of the S-B method to a contemporary adult population in Queensland, Australia and provides robust age parameters calibrated for our population. Three-dimensional surface reconstructions were generated from computed tomography scans of the pubic symphysis of male and female Caucasian individuals aged 15–70 years (n = 195) in Amira® and Rapidform®. Error was analyzed on the basis of bias, inaccuracy and percentage correct classification for left and right symphyseal surfaces. Application of transition analysis and Chi-square statistics demonstrated 63.9% and 69.7% correct age classification associated with the left symphyseal surface of Australian males and females, respectively, using the S-B method. Using Bayesian statistics, probability density distributions for each S-B phase were calculated, providing refined age parameters for our population. Mean inaccuracies of 6.77 (±2.76) and 8.28 (±4.41) years were reported for the left surfaces of males and females, respectively; with positive biases for younger individuals (<55 years) and negative biases in older individuals. Significant sexual dimorphism in the application of the S-B method was observed; and asymmetry in phase classification of the pubic symphysis was a frequent phenomenon. These results recommend that the S-B method should be applied with caution in medico-legal death investigations of Queensland skeletal remains and warrant further investigation of reliable age estimation techniques.

Error and uncertainty of adult age estimation of the pubic symphysis in an Australian sub-population using computed tomography

After attending this presentation, attendees will gain awareness of: (1) the error and uncertainty associated with the application of the Suchey-Brooks (S-B) method of age estimation of the pubic symphysis to a contemporary Australian population; (2) the implications of sexual dimorphism and bilateral asymmetry of the pubic symphysis through preliminary geometric morphometric assessment; and (3) the value of three-dimensional (3D) autopsy data acquisition for creating forensic anthropological standards. This presentation will impact the forensic science community by demonstrating that, in the absence of demographically sound skeletal collections, post-mortem autopsy data provides an exciting platform for the construction of large contemporary ‘virtual osteological libraries’ for which forensic anthropological research can be conducted on Australian individuals. More specifically, this study assesses the applicability and accuracy of the S-B method to a contemporary adult population in Queensland, Australia, and using a geometric morphometric approach, provides an insight to the age-related degeneration of the pubic symphysis. Despite the prominent use of the Suchey-Brooks (1990) method of age estimation in forensic anthropological practice, it is subject to intrinsic limitations, with reports of differential inter-population error rates between geographical locations1-4. Australian forensic anthropology is constrained by a paucity of population specific standards due to a lack of repositories of documented skeletons. Consequently, in Australian casework proceedings, standards constructed from predominately American reference samples are applied to establish a biological profile. In the global era of terrorism and natural disasters, more specific population standards are required to improve the efficiency of medico-legal death investigation in Queensland. The sample comprises multi-slice computed tomography (MSCT) scans of the pubic symphysis (slice thickness: 0.5mm, overlap: 0.1mm) on 195 individuals of caucasian ethnicity aged 15-70 years. Volume rendering reconstruction of the symphyseal surface was conducted in Amira® (v.4.1) and quantitative analyses in Rapidform® XOS. The sample was divided into ten-year age sub-sets (eg. 15-24) with a final sub-set of 65-70 years. Error with respect to the method’s assigned means were analysed on the basis of bias (directionality of error), inaccuracy (magnitude of error) and percentage correct classification of left and right symphyseal surfaces. Morphometric variables including surface area, circumference, maximum height and width of the symphyseal surface and micro-architectural assessment of cortical and trabecular bone composition were quantified using novel automated engineering software capabilities. The results of this study demonstrated correct age classification utilizing the mean and standard deviations of each phase of the S-B method of 80.02% and 86.18% in Australian males and females, respectively. Application of the S-B method resulted in positive biases and mean inaccuracies of 7.24 (±6.56) years for individuals less than 55 years of age, compared to negative biases and mean inaccuracies of 5.89 (±3.90) years for individuals greater than 55 years of age. Statistically significant differences between chronological and S-B mean age were demonstrated in 83.33% and 50% of the six age subsets in males and females, respectively. Asymmetry of the pubic symphysis was a frequent phenomenon with 53.33% of the Queensland population exhibiting statistically significant (χ2 - p<0.01) differential phase classification of left and right surfaces of the same individual. Directionality was found in bilateral asymmetry, with the right symphyseal faces being slightly older on average and providing more accurate estimates using the S-B method5. Morphometric analysis verified these findings, with the left surface exhibiting significantly greater circumference and surface area than the right (p<0.05). Morphometric analysis demonstrated an increase in maximum height and width of the surface with age, with most significant changes (p<0.05) occurring between the 25-34 and 55-64 year age subsets. These differences may be attributed to hormonal components linked to menopause in females and a reduction in testosterone in males. Micro-architectural analysis demonstrated degradation of cortical composition with age, with differential bone resorption between the medial, ventral and dorsal surfaces of the pubic symphysis. This study recommends that the S-B method be applied with caution in medico-legal death investigations of unknown skeletal remains in Queensland. Age estimation will always be accompanied by error; therefore this study demonstrates the potential for quantitative morphometric modelling of age related changes of the pubic symphysis as a tool for methodological refinement, providing a rigor and robust assessment to remove the subjectivity associated with current pelvic aging methods.

Evaluation of the Suchey-Brooks method for aging Australian Caucasian Populations using multislice Computed Tomography reconstructions of the Pubic Symphyseal Surface

Establishing age-at-death for skeletal remains is a vital component of forensic anthropology. The Suchey-Brooks (S-B) method of age estimation has been widely utilised since 1986 and relies on a visual assessment of the pubic symphyseal surface in comparison to a series of casts. Inter-population studies (Kimmerle et al., 2005; Djuric et al., 2007; Sakaue, 2006) demonstrate limitations of the S-B method, however, no assessment of this technique specific to Australian populations has been published. Aim: This investigation assessed the accuracy and applicability of the S-B method to an adult Australian Caucasian population by highlighting error rates associated with this technique. Methods: Computed tomography (CT) and contact scans of the S-B casts were performed; each geometrically modelled surface was extracted and quantified for reference purposes. A Queensland skeletal database for Caucasian remains aged 15 – 70 years was initiated at the Queensland Health Forensic and Scientific Services – Forensic Pathology Mortuary (n=350). Three-dimensional reconstruction of the bone surface using innovative volume visualisation protocols in Amira® and Rapidform® platforms was performed. Samples were allocated into 11 sub-sets of 5-year age intervals and changes associated with the surface geometry were quantified in relation to age, gender and asymmetry. Results: Preliminary results indicate that computational analysis was successfully applied to model morphological surface changes. Significant differences in observed versus actual ages were noted. Furthermore, initial morphological assessment demonstrates significant bilateral asymmetry of the pubic symphysis, which is unaccounted for in the S-B method. These results propose refinements to the S-B method, when applied to Australian casework. Conclusion: This investigation promises to transform anthropological analysis to be more quantitative and less invasive using CT imaging. The overarching goal contributes to improving skeletal identification and medico-legal death investigation in the coronial process by narrowing the range of age-at-death estimation in a biological profile.

Trace-element modeling of the magmatic evolution of rare-earth-rich carbonatite from the Miaoya deposit, Central China

Carbonatites are known to contain the highest concentrations of rare-earth elements (REE) among all igneous rocks. The REE distribution of carbonatites is commonly believed to be controlled by that of the rock forming Ca minerals (i.e., calcite, dolomite, and ankerite) and apatite because of their high modal content and tolerance for the substitution of Ca by light REE (LREE). Contrary to this conjecture, calcite from the Miaoya carbonatite (China), analyzed in situ by laser-ablation inductively-coupled-plasma mass-spectrometry, is characterized by low REE contents (100–260 ppm) and relatively !at chondrite-normalized REE distribution patterns [average (La/Yb)CN=1.6]. The carbonatite contains abundant REE-rich minerals, including monazite and !uorapatite, both precipitated earlier than the REE-poor calcite, and REE-fluorocarbonates that postdated the calcite. Hydrothermal REE-bearing !uorite and barite veins are not observed at Miaoya. The textural and analytical evidence indicates that the initially high concentrations of REE and P in the carbonatitic magma facilitated early precipitation of REE-rich phosphates. Subsequent crystallization of REE-poor calcite led to enrichment of the residual liquid in REE, particularly LREE. This implies that REE are generally incompatible with respect to calcite and the calcite/melt partition coefficients for heavy REE (HREE) are significantly greater than those for LREE. Precipitation of REE-fluorocarbonates late in the evolutionary history resulted in depletion of the residual liquid in LREE, as manifested by the development of HREE-enriched late-stage calcite [(La/Yb)CN=0.7] in syenites associated with the carbonatite. The observed variations of REE distribution between calcite and whole rocks are interpreted to arise from multistage fractional crystallization (phosphates!calcite!REE-!uorocarbonates) from an initially REE-rich carbonatitic liquid.