Phenotypic features of carbohydrate sulfotransferase 3 (CHST3) deficiency in 24 patients: congenital dislocations and vertebral changes as principal diagnostic features.

We recently reported on the deficiency of carbohydrate sulfotransferase 3 (CHST3; chondroitin-6-sulfotransferase) in six subjects diagnosed with recessive Larsen syndrome or humero-spinal dysostosis [Hermanns et al. (2008); Am J Hum Genet 82:1368-1374]. Since then, we have identified 17 additional families with CHST3 mutations and we report here on a series of 24 patients in 23 families. The diagnostic hypothesis prior to molecular analysis had been: Larsen syndrome (15 families), humero-spinal dysostosis (four cases), chondrodysplasia with multiple dislocations (CDMD "Megarbane type"; two cases), Desbuquois syndrome (one case), and spondylo-epiphyseal dysplasia (one case). In spite of the different diagnostic labels, the clinical features in these patients were similar and included dislocation of the knees and/or hips at birth, clubfoot, elbow joint dysplasia with subluxation and limited extension, short stature, and progressive kyphosis developing in late childhood. The most useful radiographic clues were the changes of the lumbar vertebrae. Twenty-four different CHST3 mutations were identified; 16 patients had homozygous mutations. We conclude that CHST3 deficiency presents at birth with congenital dislocations of knees, hips, and elbows, and is often diagnosed initially as Larsen syndrome, humero-spinal dysostosis, or chondrodysplasia with dislocations. The incidence of CHST3 deficiency seems to be higher than assumed so far. The clinical and radiographic pattern (joint dislocations, vertebral changes, normal carpal age, lack of facial flattening, and recessive inheritance) is characteristic and distinguishes CHST3 deficiency from other disorders with congenital dislocations such as filamin B-associated dominant Larsen syndrome and Desbuquois syndrome.

CD30-positive peripheral T-cell lymphomas share molecular and phenotypic features.

Peripheral T-cell lymphoma, not otherwise specified is a heterogeneous group of aggressive neoplasms with indistinct borders. By gene expression profiling we previously reported unsupervised clusters of peripheral T-cell lymphomas, not otherwise specified correlating with CD30 expression. In this work we extended the analysis of peripheral T-cell lymphoma molecular profiles to prototypical CD30(+) peripheral T-cell lymphomas (anaplastic large cell lymphomas), and validated mRNA expression profiles at the protein level. Existing transcriptomic datasets from peripheral T-cell lymphomas, not otherwise specified and anaplastic large cell lymphomas were reanalyzed. Twenty-one markers were selected for immunohistochemical validation on 80 peripheral T-cell lymphoma samples (not otherwise specified, CD30(+) and CD30(-); anaplastic large cell lymphomas, ALK(+) and ALK(-)), and differences between subgroups were assessed. Clinical follow-up was recorded. Compared to CD30(-) tumors, CD30(+) peripheral T-cell lymphomas, not otherwise specified were significantly enriched in ALK(-) anaplastic large cell lymphoma-related genes. By immunohistochemistry, CD30(+) peripheral T-cell lymphomas, not otherwise specified differed significantly from CD30(-) samples [down-regulated expression of T-cell receptor-associated proximal tyrosine kinases (Lck, Fyn, Itk) and of proteins involved in T-cell differentiation/activation (CD69, ICOS, CD52, NFATc2); upregulation of JunB and MUM1], while overlapping with anaplastic large cell lymphomas. CD30(-) peripheral T-cell lymphomas, not otherwise specified tended to have an inferior clinical outcome compared to the CD30(+) subgroups. In conclusion, we show molecular and phenotypic features common to CD30(+) peripheral T-cell lymphomas, and significant differences between CD30(-) and CD30(+) peripheral T-cell lymphomas, not otherwise specified, suggesting that CD30 expression might delineate two biologically distinct subgroups.

Oculo-facio-cardio-dental syndrome in three succeeding generations: genotypic data and phenotypic features

Oculo-facio-cardio-dental (OFCD) syndrome is a rare X-linked disorder mainly manifesting in females. Patients show ocular, facial, cardiac, and dental abnormalities. OFCD syndrome is caused by heterozygous mutations in the BCOR gene, located in Xp11.4, encoding the BCL6 co-repressor. We report a Croatian family with four female members (grandmother, mother and monozygotic female twins) diagnosed with OFCD syndrome who carry the novel BCOR mutation c.4438C>T (p.R1480*). They present high intrafamilial phenotypic variability with special regard to cardiac defect and cataract that showed more severe disease expression in successive generations. Clinical and radiographic examination of the mother of the twins revealed a talon cusp involving the permanent maxillary right central incisor. This is the first known report of a talon cusp in OFCD syndrome with a novel mutation in the BCOR gene.

Interacting roles of myofibroblasts, apoptosis and fibrogenic growth factors in the pathogenesis of renal tubulo-interstitial fibrosis

The interrelationship between myofibroblasts and fibrogenic growth factors in the pathogenesis of renal fibrosis is poorly defined. A temporal and spatial analysis of myofibroblasts, their proliferation and death, and presence of transforming growth factor-beta1 (TGF-beta1) and platelet-derived growth factor-B (PDGF-B) was carried out in an established rodent model in which chronic renal scarring and fibrosis occurs after healed renal papillary necrosis (RPN), similar to that seen with analgesic nephropathy. Treated and control groups (N = 6 and 4, respectively) were compared at 2, 4, 8 and 12 weeks. A positive relationship was found between presence of tubulo-interstitial myofibroblasts and development of fibrosis. Apoptotic myofibroblasts were identified in the interstitium and their incidence peaked 2 weeks after treatment. Levels of interstitial cell apoptosis and fibrosis were negatively correlated over time (r = -0.57, p < 0.01 ), suggesting that as apoptosis progressively failed to limit myofibroblast numbers, fibrosis increased. In comparison with the diminishing apoptosis in the interstitium, the tubular epithelium had progressively increasing levels of apoptosis over time, indicative of developing atrophy of nephrons. TGF-beta1 protein expression had a close spatial and temporal association with fibrosis and myofibroblasts, whilst PDGF-B appeared to have a closer link with populations of other chronic inflammatory cells such as infiltrating lymphocytes. Peritubular myofibroblasts were often seen near apoptotic cells in the tubular epithelium, suggestive of a paracrine toxic effect of factor/s secreted by the myofibroblasts. In vitro , TGF-beta1 was found to be toxic to renal tubular epithelial cells. These findings suggest an interaction between myofibroblasts, their deletion by apoptosis, and the presence of the fibrogenic growth factor TGF-beta1 in renal fibrosis, whereby apoptotic deletion of myofibroblasts could act as a controlling factor in progression of fibrosis.

Ankylosing spondylitis: genomic and functional characterization of candidate genes and their repercussion in clinical practice

RESUMO: Introdução: A espondilite anquilosante (EA) é uma doença inflamatória crónica caracterizada pela inflamação das articulações sacroilíacas e da coluna. A anquilose progressiva motiva uma deterioração gradual da função física e da qualidade de vida. O diagnóstico e o tratamento precoces podem contribuir para um melhor prognóstico. Neste contexto, a identificação de biomarcadores, assume-se como sendo muito útil para a prática clínica e representa hoje um grande desafio para a comunidade científica. Objetivos: Este estudo teve como objetivos: 1 - caracterizar a EA em Portugal; 2 - investigar possíveis associações entre genes, MHC e não-MHC, com a suscetibilidade e as características fenotípicas da EA; 3 - identificar genes candidatos associados a EA através da tecnologia de microarray. Material e Métodos: Foram recrutados doentes com EA, de acordo com os critérios modificados de Nova Iorque, nas consultas de Reumatologia dos diferentes hospitais participantes. Colecionaram-se dados demográficos, clínicos e radiológicos e colhidas amostras de sangue periférico. Selecionaram-se de forma aleatória, doentes HLA-B27 positivos, os quais foram tipados em termos de HLA classe I e II por PCR-rSSOP. Os haplótipos HLA estendidos foram estimados pelo algoritmo Expectation Maximization com recurso ao software Arlequin v3.11. As variantes alélicas dos genes IL23R, ERAP1 e ANKH foram estudadas através de ensaios de discriminação alélica TaqMan. A análise de associação foi realizada utilizando testes da Cochrane-Armitage e de regressão linear, tal como implementado pelo PLINK, para variáveis qualitativas e quantitativas, respetivamente. O estudo de expressão génica foi realizado por Illumina HT-12 Whole-Genome Expression BeadChips. Os genes candidatos foram validados usando qPCR-based TaqMan Low Density Arrays (TLDAs). Resultados: Foram incluídos 369 doentes (62,3% do sexo masculino, com idade média de 45,4 ± 13,2 anos, duração média da doença de 11,4 ± 10,5 anos). No momento da avaliação, 49,9% tinham doença axial, 2,4% periférica, 40,9% mista e 7,1% entesopática. A uveíte anterior aguda (33,6%) foi a manifestação extra-articular mais comum. Foram positivos para o HLA-B27, 80,3% dos doentes. Os haplótipo A*02/B*27/Cw*02/DRB1*01/DQB1*05 parece conferir suscetibilidade para a EA, e o A*02/B*27/Cw*01/DRB1*08/DQB1*04 parece conferir proteção em termos de atividade, repercussão funcional e radiológica da doença. Três variantes (2 para IL23R e 1 para ERAP1) mostraram significativa associação com a doença, confirmando a associação destes genes com a EA na população Portuguesa. O mesmo não se verificou com as variantes estudadas do ANKH. Não se verificou associação entre as variantes génicas não-MHC e as manifestações clínicas da EA. Foi identificado um perfil de expressão génica para a EA, tendo sido validados catorze genes - alguns têm um papel bem documentado em termos de inflamação, outros no metabolismo da cartilagem e do osso. Conclusões: Foi estabelecido um perfil demográfico e clínico dos doentes com EA em Portugal. A identificação de variantes génicas e de um perfil de expressão contribuem para uma melhor compreensão da sua fisiopatologia e podem ser úteis para estabelecer modelos com relevância em termos de diagnóstico, prognóstico e orientação terapêutica dos doentes. -----------ABSTRACT: Background: Ankylosing Spondylitis (AS) is a chronic inflammatory disorder characterized by inflammation in the spine and sacroiliac joints leading to progressive joint ankylosis and in progressive deterioration of physical function and quality of life. An early diagnosis and early therapy may contribute to a better prognosis. The identification of biomarkers would be helpful and represents a great challenge for the scientific community. Objectives: The present study had the following aims: 1- to characterize the pattern of AS in Portuguese patients; 2- to investigate MHC and non-MHC gene associations with susceptibility and phenotypic features of AS and; 3- to identify candidate genes associated with AS by means of whole-genome microarray. Material and Methods: AS was defined in accordance to the modified New York criteria and AS cases were recruited from hospital outcares patient clinics. Demographic and clinical data were recorded and blood samples collected. A random group of HLA-B27 positive patients and controls were selected and typed for HLA class I and II by PCR-rSSOP. The extended HLA haplotypes were estimated by Expectation Maximization Algorithm using Arlequin v3.11 software. Genotyping of IL23R, ERAP1 and ANKH allelic variants was carried out with TaqMan allelic discrimination assays. Association analysis was performed using the Cochrane-Armitage and linear regression tests as implemented in PLINK, for dichotomous and quantitative variables, respectively. Gene expression profile was carried out using Illumina HT-12 Whole-Genome Expression BeadChips and candidate genes were validated using qPCR-based TaqMan Low Density Arrays (TLDAs). Results: A total of 369 patients (62.3% male; mean age 45.4±13.2 years; mean disease duration 11.4±10.5 years), were included. Regarding clinical disease pattern, at the time of assessment, 49.9% had axial disease, 2.4% peripheral disease, 40.9% mixed disease and 7.1% isolated enthesopathic disease. Acute anterior uveitis (33.6%) was the most common extra-articular manifestation. 80.3% of AS patients were HLA-B27 positive. The haplotype A*02/B*27/Cw*02/DRB1*01/DQB1*05 seems to confer susceptibility to AS, whereas A*02/B*27/Cw*01/DRB1*08/DQB1*04 seems to provide protection in terms of disease activity, functional and radiological repercussion. Three markers (two for IL23R and one for ERAP1) showed significant single-locus disease associations. Association of these genes with AS in the Portuguese population was confirmed, whereas ANKH markers studied did not show an association with AS. No association was seen between non-MHC genes and clinical manifestations of AS. A gene expression signature for AS was established; among the fourteen validated genes, a number of them have a well-documented inflammatory role or in modulation of cartilage and bone metabolism. Conclusions: A demographic and clinical profile of patients with AS in Portugal was established. Identification of genetic variants of target genes as well as gene expression signatures could provide a better understanding of AS pathophysiology and could be useful to establish models with relevance in terms of susceptibility, prognosis, and potential therapeutic guidance.

Systems biology approaches for the design of novel Saccharomyces cerevisiae winemaking strains for enhanced flavour compounds synthesis

Tese de Doutoramento em Biologia Ambiental e Molecular

The family haloplasmataceae

Haloplasmataceae is a family within the order Haloplasmatales, which currently includes one single genus and species: Haloplasma contractile. This family has unusual phenotypic features the most noticeable being a unique morphology and cellular contractility cycle and a distinct phylogenetic position between the Firmicutes and the Tenericutes (Mollicutes). Members of the Haloplasmataceae have been isolated from the upper sediments of a deep-sea anoxic brine in the Red Sea, but cultivation-independent studies have found related sequences in a wide range of biotopes including other extreme environments, contaminated soils and marine sediments, as well as intestinal samples. The isolation and description of new representatives of this family might therefore result in significant changes to the current description.

Kinetics of T cell-activation molecules in response to Mycobacterium tuberculosis antigens

The phenotypic features acquired subsequent to antigen-specific stimulation in vitro were evaluated by means of the kinetic expressions of CD69 and CD25 activation molecules on T lymphocytes and assayed by flow cytometry in response to PPD, Ag85B, and ferritin in PPD-positive healthy control individuals. In response to PHA, CD69 staining on both CD4+ and CD8+ T cells became initially marked after 4 h, peaked at 24 h, and quickly decreased after 120 h. For CD25, a latter expression was detected around 8 h, having increased after 96 h. As expected, the response rate to the mycobacterial antigens was much lower than that to the mitogen. Positive staining was high after 96 h for CD25 and after 24 h for CD69. CD69 expression was significantly enhanced (p < 0.05) on CD8+ as compared to CD4+ T cells. High levels were also found between 96-120 h. Regarding Ag85B, CD25+ cells were mostly CD4+ instead of CD8+ T cells. Moreover, in response to ferritin, a lower CD25 expression was noted. The present data will allow further characterization of the immune response to new mycobacterial-specific antigens and their evaluation for possible inclusion in developing new diagnostic techniques for tuberculosis as well in a new vaccine to prevent the disease.

Genome size rather than content might affect call properties in toads of three ploidy levels (Anura: Bufonidae: Bufo viridis subgroup).

In vertebrates, genome size has been shown to correlate with nuclear and cell sizes, and influences phenotypic features, such as brain complexity. In three different anuran families, advertisement calls of polyploids exhibit longer notes and intervals than diploids, and difference in cellular dimensions have been hypothesized to cause these modifications. We investigated this phenomenon in green toads (Bufo viridis subgroup) of three ploidy levels, in a different call type (release calls) that may evolve independently from advertisement calls, examining 1205 calls, from ten species, subspecies, and hybrid forms. Significant differences between pulse rates of six diploid and four polyploid (3n, 4n) green toad forms across a range of temperatures from 7 to 27 °C were found. Laboratory data supported differences in pulse rates of triploids vs. tetraploids, but failed to reach significance when including field recordings. This study supports the idea that genome size, irrespective of call type, phylogenetic context, and geographical background, might affect call properties in anurans and suggests a common principle governing this relationship. The nuclear-cell size ratio, affected by genome size, seems the most plausible explanation. However, we cannot rule out hypotheses under which call-influencing genes from an unexamined diploid ancestral species might also affect call properties in the hybrid-origin polyploids.

High expression of co-stimulatory and adhesion molecules are observed on eosinophils during human Schistosoma mansoni infection

Herein we have focused attention on major phenotypic features of peripheral blood eosinophils from chronic Schistosoma mansoni-infected patients. For this purpose, detailed immunophenotypic profiles of a range of cell surface markers were performed, including activation markers (CD23/CD69/CD25/HLA-DR), co-stimulatory molecules (CD28/CD80/CD86), chemokine receptors (CXCR1/CXCR2/CCR3/CCR5) besides L-selectin-CD62L and adhesion molecules (CD18/CD54). Our major findings pointed out increased frequency of CD23+-cells, besides decreased percentages of CD69+-eosinophils, suggesting a chronic activation status with low frequency of early activated eosinophils in chronic S. mansoni-infected patients (INT) in comparison to non-infected individuals (NI). Moreover, a dichotomic expression of beta-chemokine receptors was observed during human schistosomiasis mansoni with higher CCR5 and lower levels of CCR3 observed between groups. Enhanced expression of co-stimulatory receptors (CD28/CD86) and adhesion molecules (CD54/CD18), besides striking lower frequency of L-selectin+ were reported for eosinophils from INT group as compared to NI. Interestingly, the frequency of CD62L+-eosinophils and a range of cell activation related molecules pointed out an opposite pattern of association in NI and INT, where only INT patients that display lower frequency of CD62L+-eosinophils (first CD62L tertile) kept the unusual relationship between the expression of L-selectin and the CD23 activation marker. These findings suggest that distinct dynamic of activation markers expressed by eosinophils may occur during chronic S. mansoni infection.

Molecular and clinical characterization of 25 individuals with exonic deletions of NRXN1 and comprehensive review of the literature.

This study aimed to elucidate the observed variable phenotypic expressivity associated with NRXN1 (Neurexin 1) haploinsufficiency by analyses of the largest cohort of patients with NRXN1 exonic deletions described to date and by comprehensively reviewing all comparable copy number variants in all disease cohorts that have been published in the peer reviewed literature (30 separate papers in all). Assessment of the clinical details in 25 previously undescribed individuals with NRXN1 exonic deletions demonstrated recurrent phenotypic features consisting of moderate to severe intellectual disability (91%), severe language delay (81%), autism spectrum disorder (65%), seizures (43%), and hypotonia (38%). These showed considerable overlap with previously reported NRXN1-deletion associated phenotypes in terms of both spectrum and frequency. However, we did not find evidence for an association between deletions involving the β-isoform of neurexin-1 and increased head size, as was recently published in four cases with a deletion involving the C-terminus of NRXN1. We identified additional rare copy number variants in 20% of cases. This study supports a pathogenic role for heterozygous exonic deletions of NRXN1 in neurodevelopmental disorders. The additional rare copy number variants identified may act as possible phenotypic modifiers as suggested in a recent digenic model of neurodevelopmental disorders. © 2013 Wiley Periodicals, Inc.

Complement Factor B Polymorphism and the Phenotype of Early Age-related Macular Degeneration.

Abstract Purpose: Age-related macular degeneration (AMD) has been associated with a number of polymorphisms in genes in the complement pathway. We examined the potential genotype-phenotype correlation of complement factor B (CFB) (R32Q) polymorphisms in Caucasian patients with AMD. Methods: Data from a Central European cohort of 349 patients with early AMD in at least one eye were analyzed for potential associations of the CFB (R32Q/rs641153) polymorphism with phenotypic features of early AMD. Early AMD was classified according to the International Classification and Grading System into predominant drusen size, largest drusen, drusen covered surface, central or ring-like location, peripheral drusen, and pigmentary changes. The potential association with single nucleotide polymorphisms on CFB (R32Q/rs641153) was evaluated for all patients, corrected for age, sex, and the polymorphisms of CFH (Y402H) and ARMS2 (A69S). Results: CFB (R32Q) polymorphisms showed a significant association with smaller drusen size (largest drusen ≤250 µm, p = 0.021, predominant drusen ≤125 µm, p = 0.016), with smaller surface covered by drusen (≤10%; p = 0.02), and with more frequent occurrence of peripheral drusen (p = 0.007). No association was found for pigmentary changes. Conclusions: The CFB (R32Q) polymorphism was associated with AMD characterized by small drusen only, and appeared to be protective of large drusen (OR 0.48/0.45) and of larger drusen covered area (OR 0.34). Furthermore, peripheral drusen were more frequently found (OR 2.27). This result supports the role of complement components and their polymorphisms in drusen formation and may enable a better understanding of AMD pathogenesis.

Aarskog-Scott syndrome: clinical update and report of nine novel mutations of the FGD1 gene.

Mutations in the FGD1 gene have been shown to cause Aarskog-Scott syndrome (AAS), or facio-digito-genital dysplasia (OMIM#305400), an X-linked disorder characterized by distinctive genital and skeletal developmental abnormalities with a broad spectrum of clinical phenotypes. To date, 20 distinct mutations have been reported, but little phenotypic data are available on patients with molecularly confirmed AAS. In the present study, we report on our experience of screening for mutations in the FGD1 gene in a cohort of 60 European patients with a clinically suspected diagnosis of AAS. We identified nine novel mutations in 11 patients (detection rate of 18.33%), including three missense mutations (p.R402Q; p.S558W; p.K748E), four truncating mutations (p.Y530X; p.R656X; c.806delC; c.1620delC), one in-frame deletion (c.2020_2022delGAG) and the first reported splice site mutation (c.1935+3A>C). A recurrent mutation (p.R656X) was detected in three independent families. We did not find any evidence for phenotype-genotype correlations between type and position of mutations and clinical features. In addition to the well-established phenotypic features of AAS, other clinical features are also reported and discussed.

MYC-IG rearrangements are negative predictors of survival in DLBCL patients treated with immunochemotherapy: a GELA/LYSA study.

Diffuse large B-cell lymphoma (DLBCL) with MYC rearrangement (MYC-R) carries an unfavorable outcome. We explored the prognostic value of the MYC translocation partner gene in a series of MYC-R de novo DLBCL patients enrolled in first-line prospective clinical trials (Groupe d'Etudes des Lymphomes de l'Adulte/Lymphoma Study Association) and treated with rituximab-anthracycline-based chemotherapy. A total of 774 DLBCL cases characterized for cell of origin by the Hans classifier were analyzed using fluorescence in situ hybridization with BCL2, BCL6, MYC, immunoglobulin (IG)K, and IGL break-apart and IGH/MYC, IGK/MYC, and IGL/MYC fusion probes. MYC-R was observed in 51/574 (8.9%) evaluable DLBCL cases. MYC-R cases were predominantly of the germinal center B-cell-like subtype 37/51 (74%) with no distinctive morphologic and phenotypic features. Nineteen cases were MYC single-hit and 32 cases were MYC double-hit (MYC plus BCL2 and/or BCL6) DLBCL. MYC translocation partner was an IG gene in 24 cases (MYC-IG) and a non-IG gene (MYC-non-IG) in 26 of 50 evaluable cases. Noteworthy, MYC-IG patients had shorter overall survival (OS) (P = .0002) compared with MYC-negative patients, whereas no survival difference was observed between MYC-non-IG and MYC-negative patients. In multivariate analyses, MYC-IG predicted poor progression-free survival (P = .0051) and OS (P = .0006) independently from the International Prognostic Index and the Hans classifier. In conclusion, we show in this prospective randomized trial that the adverse prognostic impact of MYC-R is correlated to the MYC-IG translocation partner gene in DLBCL patients treated with immunochemotherapy. These results may have an important impact on the clinical management of DLBCL patients with MYC-R who should be routinely characterized according to MYC partner gene. These trials are individually registered at www.clinicaltrials.gov as #NCT00144807, #NCT01087424, #NCT00169143, #NCT00144755, #NCT00140660, #NCT00140595, and #NCT00135499.

Genome-wide analysis of adaptive molecular evolution in the carnivorous plant Utricularia gibba

The genome of the bladderwort Utricularia gibba provides an unparalleled opportunity to uncover the adaptive landscape of an aquatic carnivorous plant with unique phenotypic features such as absence of roots, development of water-filled suction bladders, and a highly ramified branching pattern. Despite its tiny size, the U. gibba genome accommodates approximately as many genes as other plant genomes. To examine the relationship between the compactness of its genome and gene turnover, we compared the U. gibba genome with that of four other eudicot species, defining a total of 17,324 gene families (orthogroups). These families were further classified as either 1) lineage-specific expanded/contracted or 2) stable in size. The U. gibba-expanded families are generically related to three main phenotypic features: 1) trap physiology, 2) key plant morphogenetic/developmental pathways, and 3) response to environmental stimuli, including adaptations to life in aquatic environments. Further scans for signatures of protein functional specialization permitted identification of seven candidate genes with amino acid changes putatively fixed by positive Darwinian selection in the U. gibba lineage. The Arabidopsis orthologs of these genes (AXR, UMAMIT41, IGS, TAR2, SOL1, DEG9, and DEG10) are involved in diverse plant biological functions potentially relevant for U. gibba phenotypic diversification, including 1) auxin metabolism and signal transduction, 2) flowering induction and floral meristem transition, 3) root development, and 4) peptidases. Taken together, our results suggest numerous candidate genes and gene families as interesting targets for further experimental confirmation of their functional and adaptive roles in the U. gibba's unique lifestyle and highly specialized body plan.