egr-4, a target of EGFR signaling, is required for the formation of the brian primordial and head regeneration in planarians

During the regeneration of freshwater planarians, polarity and patterning programs play essential roles in determining whether a head or a tail regenerates at anterior or posterior-facing wounds. This decision is made very soon after amputation. The pivotal role of the Wnt/β-catenin and Hh signaling pathways in re-establishing anterior-posterior (AP) polarity has been well documented. However, the mechanisms that control the growth and differentiation of the blastema in accordance with its AP identity are less well understood. Previous studies have described a role of Smed-egfr-3, a planarian epidermal growth factor receptor, in blastema growth and differentiation. Here, we identify Smed-egr-4, a zinc-finger transcription factor belonging to the early growth response gene family, as a putative downstream target of Smed-egfr-3. Smed-egr-4 is mainly expressed in the central nervous system and its silencing inhibits anterior regeneration without affecting the regeneration of posterior regions. Single and combinatorial RNA interference to target different elements of the Wnt/β-catenin pathway, together with expression analysis of brain- and anterior-specific markers, revealed that Smed-egr-4: (1) is expressed in two phases - an early Smed-egfr-3-independent phase and a late Smed-egfr-3-dependent phase; (2) is necessary for the differentiation of the brain primordia in the early stages of regeneration; and (3) that it appears to antagonize the activity of the Wnt/β-catenin pathway to allow head regeneration. These results suggest that a conserved EGFR/egr pathway plays an important role in cell differentiation during planarian regeneration and indicate an association between early brain differentiation and the proper progression of head regeneration.

An update of the mechanisms of resistance to EGFR-tyrosine kinase inhibitors in breast cancer: gefitinib (Iressatrade mark)-induced changes in the expression and nucleo-cytoplasmic trafficking of HER-ligands (Review)

Intrinsic resistance to the epidermal growth factor receptor (EGFR; HER1) tyrosine kinase inhibitor (TKI) gefitinib, and more generally to EGFR TKIs, is a common phenomenon in breast cancer. The availability of molecular criteria for predicting sensitivity to EGFR-TKIs is, therefore, the most relevant issue for their correct use and for planning future research. Though it appears that in non-small-cell lung cancer (NSCLC) response to gefitinib is directly related to the occurrence of specific mutations in the EGFR TK domain, breast cancer patients cannot be selected for treatment with gefitinib on the same basis as such EGFR mutations have beenreported neither in primary breast carcinomas nor in several breast cancer cell lines. Alternatively, there is a generalagreement on the hypothesis that the occurrence of molecular alterations that activate transduction pathways downstreamof EGFR (i.e., MEK1/MEK2 - ERK1/2 MAPK and PI-3'K - AKT growth/survival signaling cascades) significantly affect the response to EGFR TKIs in breast carcinomas. However,there are no studies so far addressing a role of EGF-related ligands as intrinsic breast cancer cell modulators of EGFR TKIefficacy. We recently monitored gene expression profiles andsub-cellular localization of HER-1/-2/-3/-4 related ligands (i.e., EGF, amphiregulin, transforming growth factor-α, ß-cellulin,epiregulin and neuregulins) prior to and after gefitinib treatment in a panel of human breast cancer cell lines. First, gefitinibinduced changes in the endogenous levels of EGF-related ligands correlated with the natural degree of breast cancer cellsensitivity to gefitinib. While breast cancer cells intrinsically resistant to gefitinib (IC50 ≥15 μM) markedly up-regulated(up to 600 times) the expression of genes codifying for HERspecific ligands, a significant down-regulation (up to 106 times)of HER ligand gene transcription was found in breast cancer cells intrinsically sensitive to gefitinib (IC50 ≤1 μM). Second,loss of HER1 function differentially regulated the nuclear trafficking of HER-related ligands. While gefitinib treatment induced an active import and nuclear accumulation of the HER ligand NRG in intrinsically gefitinib-resistant breastcancer cells, an active export and nuclear loss of NRG was observed in intrinsically gefitinib-sensitive breast cancer cells.In summary, through in vitro and pharmacodynamic studies we have learned that, besides mutations in the HER1 gene,oncogenic changes downstream of HER1 are the key players regulating gefitinib efficacy in breast cancer cells. It now appears that pharmacological inhibition of HER1 functionalso leads to striking changes in both the gene expression and the nucleo-cytoplasmic trafficking of HER-specific ligands,and that this response correlates with the intrinsic degree of breast cancer sensitivity to the EGFR TKI gefitinib. Therelevance of this previously unrecognized intracrine feedback to gefitinib warrants further studies as cancer cells could bypassthe antiproliferative effects of HER1-targeted therapeutics without a need for the overexpression and/or activation of other HER family members and/or the activation of HER-driven downstream signaling cascades

Valor pronóstico de la densidad microvascular y de la expresión del VEGF, EGFR y HIF-1 en pacientes con cáncer de cérvix localmente avanzado tratados con quimioradiación

Introducción: Colombia cuenta con poca información sobre el comportamiento del cáncer, no obstante, el carcinoma de cuello uterino representa la segunda causa de muerte por la enfermedad entre las mujeres de nuestro entorno. El patrón epidemiológico de la enfermedad es preocupante porque los estados localmente avanzados constituyen el estado más frecuente al momento del diagnóstico y la mortalidad siendo bastante alta a pesar de la presencia de un programa de cribado organizado. Objetivo: Describir el valor pronóstico de la densidad microvascular (DMV) y de la expresión proteica de varios genes relacionados con la supervivencia y proliferación del cáncer de cérvix localmente avanzado en un grupo de mujeres tratadas con quimioradiación y braquiterapia intracavitaria. Se estimaron la tasa de respuesta global (TRG), la supervivencia libre de progresión (SLP) y la supervivencia global (SG). Resultados: Se incluyeron 61 mujeres con una edad media de 52 ± 10 años; todas tenían diagnóstico de cáncer de cérvix localmente avanzado (IIA 2.3%/IIB 47.5%/IIIA 4.9%/IIIB 37.7%/IVA 3.3%/no definido 3.3%), con un volumen tumoral promedio de 6.4cm (DE ± 1.8cm) e infección por VPH en 46% de los casos; 58 sujetos (95%) tenían un patrón escamoso, dos fueron adenocarcinomas y &50% presentaba neoplasias moderada o pobremente diferenciadas. Todas fueron tratadas con quimioradiación (interrupción transitoria en teleterapia por toxicidad y otras causas en 19% y 21.4%, respectivamente/media de ciclos de platino concomitante 4.8 series ± 1.0) y braquiterapia (77% completaron el tratamiento intracavitario). La mediana para la SLP y global fue de 6.6 meses (r, 4.0-9.1) y 30 meses (r, 11-48), respectivamente. Ninguna de las variables tuvo un efecto positivo sobre la SLP, mientras el análisis multivariado demostró que los niveles de expresión del VEGF (P=0.026), EGFR (P=0.030), y el volumen tumoral menor de 6 cm (P=0.02) influyeron positivamente sobre éste desenlace. Conclusión: Existe una influencia positiva sobre el pronóstico, de la tipificación en el cáncer de cérvix localmente avanzado tratado con quimioradiación basada en platino.

Relació i valor pronòstic de la metilació del promotor MGMT , amplificació del EGFR i mutació de l’ EGFRvIII en una sèrie de pacients amb glioblastoma.

Introducció. La metilació de l’enzim MGMT és un factor predictiu de resposta en pacients amb glioblastoma tractats amb alquilants, però el seu valor pronòstic s’està estudiant . Tanmateix, en el cas de l’amplificació de l ‘EGFR i la mutació EGFRvIII no se s’ha demostrat cap mena de valor pronòstic ni predictiu. Objetius. Determinar la distribució de la metilació de l’enzim MGMT , amplificació de l’ EGFR i mutació EGFR vIII , l’asociació entre ells i avaluar la supervivència global i supervivència lliure de progresió en relació amb l’estat de metilació del MGMT , amplificació de l ‘EGFR mutació EGFR vIII . Material i mètodes. S’hi han determinat les variacions genètiques en una població de 70 pacients amb glioblastoma i s’han relacionat, a més, amb la supervivènciia global i temps a la progressió. Resultats. La metilació del MGMT , amplificació de l’ EGFR i detecció de ’l EGFR vIII foren detectats en un 38,6 % , 50% y 25,7% DELS pacients , respectivament. No es va demostrar pas cap relació entre l’estat de metilació del MGMT i l’estat de l’ EGFR , a diferència de la relació entre l’amplificació de l’ EGFR i la mutació EGFRvIII. En l’anàlisi multivariant de supervivència global i temps a la progressió, els factors edat , radioteràpia i metilació del MGMT foren significatius. . Conclusions. L’estat de metilació del MGMT és, doncs, un factor pronòstic de supervivència en pacients diagnosticats de glioblastoma tractats amb radioteràpia i quimioteràpia basada en alquilants , a diferencia de l’amplificació de l’ EGFR i mutació EGFRvIII.

“Unmixing” Tissue Gene Expression Signatures from Tumor Biopsies

Emergent molecular measurement methods, such as DNA microarray, qRTPCR, andmany others, offer tremendous promise for the personalized treatment of cancer. Thesetechnologies measure the amount of specific proteins, RNA, DNA or other moleculartargets from tumor specimens with the goal of “fingerprinting” individual cancers. Tumorspecimens are heterogeneous; an individual specimen typically contains unknownamounts of multiple tissues types. Thus, the measured molecular concentrations resultfrom an unknown mixture of tissue types, and must be normalized to account for thecomposition of the mixture.For example, a breast tumor biopsy may contain normal, dysplastic and cancerousepithelial cells, as well as stromal components (fatty and connective tissue) and bloodand lymphatic vessels. Our diagnostic interest focuses solely on the dysplastic andcancerous epithelial cells. The remaining tissue components serve to “contaminate”the signal of interest. The proportion of each of the tissue components changes asa function of patient characteristics (e.g., age), and varies spatially across the tumorregion. Because each of the tissue components produces a different molecular signature,and the amount of each tissue type is specimen dependent, we must estimate the tissuecomposition of the specimen, and adjust the molecular signal for this composition.Using the idea of a chemical mass balance, we consider the total measured concentrationsto be a weighted sum of the individual tissue signatures, where weightsare determined by the relative amounts of the different tissue types. We develop acompositional source apportionment model to estimate the relative amounts of tissuecomponents in a tumor specimen. We then use these estimates to infer the tissuespecificconcentrations of key molecular targets for sub-typing individual tumors. Weanticipate these specific measurements will greatly improve our ability to discriminatebetween different classes of tumors, and allow more precise matching of each patient tothe appropriate treatment

Diversity of the bacterial community in the surface soil of a pear orchard based on 16S rRNA gene analysis

A cultivation-independent approach based on polymerase chain reaction (PCR)-amplified partial small subunit rRNA genes was used to characterize bacterial populations in the surface soil of a commercial pear orchard consisting of different pear cultivars during two consecutive growing seasons. Pyrus communis L. cvs Blanquilla, Conference, and Williams are among the most widely cultivated cultivars in Europe and account for the majority of pear production in Northeastern Spain. To assess the heterogeneity of the community structure in response to environmental variables and tree phenology, bacterial populations were examined using PCR-denaturing gradient gel electrophoresis (DGGE) followed by cluster analysis of the 16S ribosomal DNA profiles by means of the unweighted pair group method with arithmetic means. Similarity analysis of the band patterns failed to identify characteristic fingerprints associated with the pear cultivars. Both environmentally and biologically based principal-component analyses showed that the microbial communities changed significantly throughout the year depending on temperature and, to a lesser extent, on tree phenology and rainfall. Prominent DGGE bands were excised and sequenced to gain insight into the identities of the predominant bacterial populations. Most DGGE band sequences were related to bacterial phyla, such as Bacteroidetes, Cyanobacteria, Acidobacteria, Proteobacteria, Nitrospirae, and Gemmatimonadetes, previously associated with typical agronomic crop environments

CAV2 vector gene transfer into central nervious System

Estudi realitzat a partir d’una estada al Institut de Génétique Moléculaire de Montpellier, França, entre 2010 i 2012. En aquest projecte s’ha avaluat les avantatges dels vectors adenovirals canins tipus 2 (CAV2) com a vectors de transferència gènica al sistema nerviós central (SNC) en un model primat no-humà i en un model caní del síndrome de Sly (mucopolisacaridosis tipus 7, MPS VII), malaltia monogènica que cursa amb neurodegeneració. En una primera part del projecte s’ha avaluat la biodistribució, l’eficàcia i la durada de l’expressió del transgen en un model primat no humà, (Microcebus murinus). Com ha vector s’ha utilitzat un CAV2 de primera generació que expressa la proteïna verda fluorescent (CAVGFP). Els resultats aportats en aquesta memòria demostren que en primats no humans, com en d’altres espècies testades anteriorment per l’equip de l’EJ Kremer, la injecció intracerebral de CAV2 resulta en una extensa transducció del SNC, siguent les neurones i els precursors neuronals les cèl•lules preferencialment transduïdes. Els vectors canins, servint-se de vesícules intracel•lulars són transportats, majoritàriament, des de les sinapsis cap al soma neuronal, aquest transport intracel•lular permet una extensa transducció del SNC a partir d’una única injecció intracerebral dels vectors virals. En una segona part d’aquest projecte s’ha avaluat l’ús terapèutic dels CAV2. S’ha injectat un vector helper-dependent que expressa el gen la b-glucuronidasa i el gen de la proteïna verda fluorescent (HD-RIGIE), en el SNC del model caní del síndrome de Sly (MPS VII). La biodistribució i la eficàcia terapèutica han estat avaluades. Els nivells d’activitat enzimàtica en animals malalts injectats amb el vector terapèutic va arribar a valors similars als dels animals no afectes. A més a més s’ha observat una reducció en la quantitat dels GAGs acumulats en les cèl•lules dels animals malalts tractats amb el vector terapèutic, demostrant la potencialitat terapèutica dels CAV2 per a malalties que afecten al SNC. Els resultats aportats en aquest treball ens permeten dir que els CAV2 són unes bones eines terapèutiques per al tractament de malalties que afecten al SNC.

Whole genome analysis of p38 SAPK-mediated gene expression upon stress

Background: Cells have the ability to respond and adapt to environmental changes through activation of stress-activated protein kinases (SAPKs). Although p38 SAPK signalling is known to participate in the regulation of gene expression little is known on the molecular mechanisms used by this SAPK to regulate stress-responsive genes and the overall set of genes regulated by p38 in response to different stimuli.Results: Here, we report a whole genome expression analyses on mouse embryonic fibroblasts (MEFs) treated with three different p38 SAPK activating-stimuli, namely osmostress, the cytokine TNFα and the protein synthesis inhibitor anisomycin. We have found that the activation kinetics of p38α SAPK in response to these insults is different and also leads to a complex gene pattern response specific for a given stress with a restricted set of overlapping genes. In addition, we have analysed the contribution of p38α the major p38 family member present in MEFs, to the overall stress-induced transcriptional response by using both a chemical inhibitor (SB203580) and p38α deficient (p38α-/-) MEFs. We show here that p38 SAPK dependency ranged between 60% and 88% depending on the treatments and that there is a very good overlap between the inhibitor treatment and the ko cells. Furthermore, we have found that the dependency of SAPK varies depending on the time the cells are subjected to osmostress. Conclusions: Our genome-wide transcriptional analyses shows a selective response to specific stimuli and a restricted common response of up to 20% of the stress up-regulated early genes that involves an important set of transcription factors, which might be critical for either cell adaptation or preparation for continuous extra-cellular changes. Interestingly, up to 85% of the up-regulated genes are under the transcriptional control of p38 SAPK. Thus, activation of p38 SAPK is critical to elicit the early gene expression program required for cell adaptation to stress.

Human synthetic lethal inference as potential anti-cancer target gene detection

Background: Two genes are called synthetic lethal (SL) if mutation of either alone is not lethal, but mutation of both leads to death or a significant decrease in organism's fitness. The detection of SL gene pairs constitutes a promising alternative for anti-cancer therapy. As cancer cells exhibit a large number of mutations, the identification of these mutated genes' SL partners may provide specific anti-cancer drug candidates, with minor perturbations to the healthy cells. Since existent SL data is mainly restricted to yeast screenings, the road towards human SL candidates is limited to inference methods. Results: In the present work, we use phylogenetic analysis and database manipulation (BioGRID for interactions, Ensembl and NCBI for homology, Gene Ontology for GO attributes) in order to reconstruct the phylogenetically-inferred SL gene network for human. In addition, available data on cancer mutated genes (COSMIC and Cancer Gene Census databases) as well as on existent approved drugs (DrugBank database) supports our selection of cancer-therapy candidates.Conclusions: Our work provides a complementary alternative to the current methods for drug discovering and gene target identification in anti-cancer research. Novel SL screening analysis and the use of highly curated databases would contribute to improve the results of this methodology.

Fusion of the human gene for the polyubiquitination co-effector UEV-1 with Kua, a newly identified gene

UEV proteins are enzymatically inactive variants of the E2 ubiquitin-conjugating enzymes that regulate noncanonical elongation of ubiquitin chains. In Saccharomyces cerevisiae, UEV is part of the RAD6-mediated error-free DNA repair pathway. In mammalian cells, UEV proteins can modulate c-FOS transcription and the G2-M transition of the cell cycle. Here we show that the UEV genes from phylogenetically distant organisms present a remarkable conservation in their exon–intron structure. We also show that the human UEV1 gene is fused with the previously unknown gene Kua. In Caenorhabditis elegans and Drosophila melanogaster, Kua and UEV are in separated loci, and are expressed as independent transcripts and proteins. In humans, Kua and UEV1 are adjacent genes, expressed either as separate transcripts encoding independent Kua and UEV1 proteins, or as a hybrid Kua–UEV transcript, encoding a two-domain protein. Kua proteins represent a novel class of conserved proteins with juxtamembrane histidine-rich motifs. Experiments with epitope-tagged proteins show that UEV1A is a nuclear protein, whereas both Kua and Kua–UEV localize to cytoplasmic structures, indicating that the Kua domain determines the cytoplasmic localization of Kua–UEV. Therefore, the addition of a Kua domain to UEV in the fused Kua–UEV protein confers new biological properties to this regulator of variant polyubiquitination.[Kua cDNAs isolated by RT-PCR and described in this paper have been deposited in the GenBank data library under accession nos. AF1155120 (H. sapiens) and AF152361 (D. melanogaster). Genomic clones containing UEV genes: S. cerevisiae, YGL087c (accession no. Z72609); S. pombe, c338 (accession no. AL023781); P. falciparum, MAL3P2 (accession no. AL034558); A. thaliana, F26F24 (accession no. AC005292); C. elegans, F39B2 (accession no. Z92834); D. melanogaster, AC014908; and H. sapiens, 1185N5 (accession no. AL034423). Accession numbers for Kua cDNAs in GenBank dbEST: M. musculus, AA7853; T. cruzi, AI612534. Other Kua-containing sequences: A. thaliana genomic clones F10M23 (accession no. AL035440), F19K23 (accession no. AC000375), and T20K9 (accession no. AC004786).

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.

Improving gene annotation using peptide mass spectrometry

Annotation of protein-coding genes is a key goal of genome sequencing projects. In spite of tremendous recent advances in computational gene finding, comprehensive annotation remains a challenge. Peptide mass spectrometry is a powerful tool for researching the dynamic proteome and suggests an attractive approach to discover and validate protein-coding genes. We present algorithms to construct and efficiently search spectra against a genomic database, with no prior knowledge of encoded proteins. By searching a corpus of 18.5 million tandem mass spectra (MS/MS) from human proteomic samples, we validate 39,000 exons and 11,000 introns at the level of translation. We present translation-level evidence for novel or extended exons in 16 genes, confirm translation of 224 hypothetical proteins, and discover or confirm over 40 alternative splicing events. Polymorphisms are efficiently encoded in our database, allowing us to observe variant alleles for 308 coding SNPs. Finally, we demonstrate the use of mass spectrometry to improve automated gene prediction, adding 800 correct exons to our predictions using a simple rescoring strategy. Our results demonstrate that proteomic profiling should play a role in any genome sequencing project.

Comparative gene prediction in human and mouse

The completion of the sequencing of the mouse genome promises to help predict human genes with greater accuracy. While current ab initio gene prediction programs are remarkably sensitive (i.e., they predict at least a fragment of most genes), their specificity is often low, predicting a large number of false-positive genes in the human genome. Sequence conservation at the protein level with the mouse genome can help eliminate some of those false positives. Here we describe SGP2, a gene prediction program that combines ab initio gene prediction with TBLASTX searches between two genome sequences to provide both sensitive and specific gene predictions. The accuracy of SGP2 when used to predict genes by comparing the human and mouse genomes is assessed on a number of data sets, including single-gene data sets, the highly curated human chromosome 22 predictions, and entire genome predictions from ENSEMBL. Results indicate that SGP2 outperforms purely ab initio gene prediction methods. Results also indicate that SGP2 works about as well with 3x shotgun data as it does with fully assembled genomes. SGP2 provides a high enough specificity that its predictions can be experimentally verified at a reasonable cost. SGP2 was used to generate a complete set of gene predictions on both the human and mouse by comparing the genomes of these two species. Our results suggest that another few thousand human and mouse genes currently not in ENSEMBL are worth verifying experimentally.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Altered brain-derived neurotrophic factor blood levels and gene variability are associated with anorexia and bulimia

Murine models and association studies in eating disorder (ED) patients have shown a role for the brain-derived neurotrophic factor (BDNF) in eating behavior. Some studies have shown association of BDNF -270C/T single-nucleotide polymorphism (SNP) with bulimia nervosa (BN), while BDNF Val66Met variant has been shown to be associated with both BN and anorexia nervosa (AN). To further test the role of this neurotrophin in humans, we screened 36 SNPs in the BDNF gene and tested for their association with ED and plasma BDNF levels as a quantitative trait. We performed a family-based association study in 106 ED nuclear families and analyzed BDNF blood levels in 110 ED patients and in 50 sib pairs discordant for ED. The rs7124442T/rs11030102C/rs11030119G haplotype was found associated with high BDNF levels (mean BDNF TCG haplotype carriers = 43.6 ng/ml vs. mean others 23.0 ng/ml, P = 0.016) and BN (Z = 2.64; P recessive = 0.008), and the rs7934165A/270T haplotype was associated with AN (Z =-2.64; P additive = 0.008). The comparison of BDNF levels in 50 ED discordant sib pairs showed elevated plasma BDNF levels for the ED group (mean controls = 41.0 vs. mean ED = 52.7; P = 0.004). Our data strongly suggest that altered BDNF levels modulated by BDNF gene variability are associated with the susceptibility to ED, providing physiological evidence that BDNF plays a role in the development of AN and BN, and strongly arguing for its involvement in eating behavior and body weight regulation.