Functional and molecular evolutionary analysis of the "Brevix Radix" gene family in mono-and dicotyledonous plants

ABSTRACTIn contrast to animals, plants cannot move from their place of birth and, therefore, need to adapt to their particular habitat in order to survive. Thus, plant development is remarkably plastic, making plants an ideal system for the isolation of genes that account for intraspecific natural variation and possibly environmental adaptation. However, to date, this approach mostly identified null alleles and missed mutations with subtle effects. For instance, BREVIS RADIX (BRX) has been isolated as a key regulator of root growth through a naturally occurring loss-of-function allele in the Arabidopsis thaliana accession Uk-1 and is the founding member of a highly-conserved plant-specific gene family.In this work, we show that a strong selective pressure is acting on the BRX gene family and dates back before the monocot-dicot divergence. However, functional diversification is observed mainly in dicotyledon BRX family genes and is correlated with acceleration in the evolutionary rates in the N-terminal regions. Population genetic data revealed that BRX is highly conserved across Arabidopsis accessions and presents signatures of adaptation. Interestingly, a seven amino acid deletion polymorphism in BRX sequence was found in a few accessions, which seems to be responsible for their enhanced primary root growth. Nevertheless, BRX might not only be active in the root, as suggested by its expression in the shoot. Indeed, leaves and cotyledons of brx mutants are significantly smaller than wild- type. This phenotype is a direct consequence of the absence of BRX function in the shoot rather than an indirect effect of an altered root system growth. Interestingly, cotyledons of brx plants reflect the same physiological defects as the root. Moreover, phenotypes in BRX gain-of-function plants, such as epinastic leaves and increased epidermal cell size, could be associated with an increase in leaf brassinosteroid content.Collectively, these results indicate that BRX contributes to local adaptation by ubiquitously regulating plant growth, probably through the modulation of brassinosteroid biosynthesis.RÉSUMÉContrairement à la plupart des animaux, les plantes ne peuvent se mouvoir et doivent ainsi s'adapter à leur environnement pour survivre. Pour cette raison, elles représentent un système idéal pour l'identification de gènes contribuant à la variation naturelle intra- spécifique, ainsi qu'à l'adaptation. Cependant, cette approche a, jusqu'à présent, surtout permis d'isoler des allèles nuls et non des mutations conférant des effets plus subtiles. C'est le cas du gène Β REVIS RADIX (BRX), un régulateur clé de la croissance racinaire, qui a été identifié grâce à un allèle non-fonctionnel présent dans l'accession naturelle d'Arabidopsis thaliana Uk-1. BRX et ses homologues des plantes mono- et dicotylédones forment une famille très conservée et spécifique aux plantes.Dans ce travail, nous démontrons que la famille de gènes BRX est soumise à une forte pression de sélection qui remonte avant la divergence entre mono- et dicotylédones. Cependant, une diversification fonctionnelle a été observée chez les gènes des dicotylédones et corrèle avec une accélération de la vitesse d'évolution dans leur région N- terminale. Une analyse génétique de différentes accessions naturelles d'Arabidopsis a révélé que BRX est hautement conservé et présente des signatures d'adaptation. Remarquablement, un polymorphisme de délétion de sept acides aminés a été détecté dans quelques accessions et a pour conséquence une plus forte croissance de la racine primaire. Néanmoins, il semble que le rôle de BRX ne se limite pas qu'à la racine, comme indiqué par son expression dans les parties aériennes de la plante. En effet, les mutants brx présentent des cotylédons et des feuilles significativement plus petits que le type sauvage, une conséquence directe de l'absence d'activité de BRX dans ces organes. Nous avons aussi noté que les cotylédons des mutants brx, à l'instar des racines, ont une perception altérée de l'auxine et peuvent être complémentés par l'application exogène de brassinostéroïdes. De plus, dans des plantes présentant un gain de fonction BRX, les feuilles sont épinastiques et les cellules de leur épiderme plus grandes. Ces phénotypes sont accompagnés d'une augmentation de la concentration de brassinostéroïdes dans les feuilles. Conjointement, ces résultats démontrent que BRX contribue à une adaptation locale de la plante par la régulation générale de sa croissance, probablement en modulant la biosynthèse des brassinostéroïdes.

Expression profiling of rectal tumors defines response to neoadjuvant treatment related genes.

To date, no effective method exists that predicts the response to preoperative chemoradiation (CRT) in locally advanced rectal cancer (LARC). Nevertheless, identification of patients who have a higher likelihood of responding to preoperative CRT could be crucial in decreasing treatment morbidity and avoiding expensive and time-consuming treatments. The aim of this study was to identify signatures or molecular markers related to response to pre-operative CRT in LARC. We analyzed the gene expression profiles of 26 pre-treatment biopsies of LARC (10 responders and 16 non-responders) without metastasis using Human WG CodeLink microarray platform. Two hundred and fifty seven genes were differentially over-expressed in the responder patient subgroup. Ingenuity Pathway Analysis revealed a significant ratio of differentially expressed genes related to cancer, cellular growth and proliferation pathways, and c-Myc network. We demonstrated that high Gng4, c-Myc, Pola1, and Rrm1 mRNA expression levels was a significant prognostic factor for response to treatment in LARC patients (p<0.05). Using this gene set, we were able to establish a new model for predicting the response to CRT in rectal cancer with a sensitivity of 60% and 100% specificity. Our results reflect the value of gene expression profiling to gain insight about the molecular pathways involved in the response to treatment of LARC patients. These findings could be clinically relevant and support the use of mRNA levels when aiming to identify patients who respond to CRT therapy.

Comparative analysis of human and mouse expression data illuminates tissue-specific evolutionary patterns of miRNAs.

MicroRNAs (miRNAs) constitute an important class of gene regulators. While models have been proposed to explain their appearance and expansion, the validation of these models has been difficult due to the lack of comparative studies. Here, we analyze miRNA evolutionary patterns in two mammals, human and mouse, in relation to the age of miRNA families. In this comparative framework, we confirm some predictions of previously advanced models of miRNA evolution, e.g. that miRNAs arise more frequently de novo than by duplication, or that the number of protein-coding gene targeted by miRNAs decreases with evolutionary time. We also corroborate that miRNAs display an increase in expression level with evolutionary time, however we show that this relation is largely tissue-dependent, and especially low in embryonic or nervous tissues. We identify a bias of tag-sequencing techniques regarding the assessment of breadth of expression, leading us, contrary to predictions, to find more tissue-specific expression of older miRNAs. Together, our results refine the models used so far to depict the evolution of miRNA genes. They underline the role of tissue-specific selective forces on the evolution of miRNAs, as well as the potential co-evolution patterns between miRNAs and the protein-coding genes they target.

A positive FGFR3/FOXN1 feedback loop underlies benign skin keratosis versus squamous cell carcinoma formation in humans.

Seborrheic keratoses (SKs) are common, benign epithelial tumors of the skin that do not, or very rarely, progress into malignancy, for reasons that are not understood. We investigated this by gene expression profiling of human SKs and cutaneous squamous cell carcinomas (SCCs) and found that several genes previously connected with keratinocyte tumor development were similarly modulated in SKs and SCCs, whereas the expression of others differed by only a few fold. In contrast, the tyrosine kinase receptor FGF receptor-3 (FGFR3) and the transcription factor forkhead box N1 (FOXN1) were highly expressed in SKs, and close to undetectable in SCCs. We also showed that increased FGFR3 activity was sufficient to induce FOXN1 expression, counteract the inhibitory effect of EGFR signaling on FOXN1 expression and differentiation, and induce differentiation in a FOXN1-dependent manner. Knockdown of FOXN1 expression in primary human keratinocytes cooperated with oncogenic RAS in the induction of SCC-like tumors, whereas increased FOXN1 expression triggered the SCC cells to shift to a benign SK-like tumor phenotype, which included increased FGFR3 expression. Thus,we have uncovered a positive regulatory loop between FGFR3 and FOXN1 that underlies a benign versus malignant skin tumor phenotype.

Characterization of the fasting-induced adipose factor FIAF, a novel peroxisome proliferator-activated receptor target gene.

Fasting is associated with significant changes in nutrient metabolism, many of which are governed by transcription factors that regulate the expression of rate-limiting enzymes. One factor that plays an important role in the metabolic response to fasting is the peroxisome proliferator-activated receptor alpha (PPARalpha). To gain more insight into the role of PPARalpha during fasting, and into the regulation of metabolism during fasting in general, a search for unknown PPARalpha target genes was performed. Using subtractive hybridization (SABRE) comparing liver mRNA from wild-type and PPARalpha null mice, we isolated a novel PPARalpha target gene, encoding the secreted protein FIAF (for fasting induced adipose factor), that belongs to the family of fibrinogen/angiopoietin-like proteins. FIAF is predominantly expressed in adipose tissue and is strongly up-regulated by fasting in white adipose tissue and liver. Moreover, FIAF mRNA is decreased in white adipose tissue of PPARgamma +/- mice. FIAF protein can be detected in various tissues and in blood plasma, suggesting that FIAF has an endocrine function. Its plasma abundance is increased by fasting and decreased by chronic high fat feeding. The data suggest that FIAF represents a novel endocrine signal involved in the regulation of metabolism, especially under fasting conditions.

Cost-effectiveness analysis of the first-line therapies for nicotine dependence.

BACKGROUND: Nicotine dependence is the major obstacle for smokers who want to quit. Guidelines have identified five effective first-line therapies, four nicotine replacement therapies (NRTs)--gum, patch, nasal spray and inhaler--and bupropion. Studying the extent to which these various treatments are cost-effective requires additional research. OBJECTIVES: To determine cost-effectiveness (CE) ratios of pharmacotherapies for nicotine dependence provided by general practitioners (GPs) during routine visits as an adjunct to cessation counselling. METHODS: We used a Markov model to generate two cohorts of one-pack-a-day smokers: (1) the reference cohort received only cessation counselling from a GP during routine office visits; (2) the second cohort received the same counselling plus an offer to use a pharmacological treatment to help them quit smoking. The effectiveness of adjunctive therapy was expressed in terms of the resultant differential in mortality rate between the two cohorts. Data on the effectiveness of therapies came from meta-analyses, and we used odds ratio for quitting as the measure of effectiveness. The costs of pharmacotherapies were based on the cost of the additional time spent by GPs offering, prescribing and following-up treatment, and on the retail prices of the therapies. We used the third-party-payer perspective. Results are expressed as the incremental cost per life-year saved. RESULTS: The cost per life-year saved for only counselling ranged from Euro 385 to Euro 622 for men and from Euro 468 to Euro 796 for women. The CE ratios for the five pharmacological treatments varied from Euro 1768 to Euro 6879 for men, and from Euro 2146 to Euro 8799 for women. Significant variations in CE ratios among the five treatments were primarily due to differences in retail prices. The most cost-effective treatments were bupropion and the patch, and, then, in descending order, the spray, the inhaler and, lastly, gum. Differences in CE between men and women across treatments were due to the shape of their respective mortality curve. The lowest CE ratio in men was for the 45- to 49-year-old group and for women in the 50- to 54-year-old group. Sensitivity analysis showed that changes in treatment efficacy produced effects only for less-well proven treatments (spray, inhaler, and bupropion) and revealed a strong influence of the discount rate and natural quit rate on the CE of pharmacological treatments. CONCLUSION: The CE of first-line treatments for nicotine dependence varied widely with age and sex and was sensitive to the assumption for the natural quit rate. Bupropion and the nicotine patch were the two most cost-effective treatments.

Analysis of genes encoding penicillin-binding proteins in clinical isolates of Acinetobacter baumannii.

There is limited information on the role of penicillin-binding proteins (PBPs) in the resistance of Acinetobacter baumannii to β-lactams. This study presents an analysis of the allelic variations of PBP genes in A. baumannii isolates. Twenty-six A. baumannii clinical isolates (susceptible or resistant to carbapenems) from three teaching hospitals in Spain were included. The antimicrobial susceptibility profile, clonal pattern, and genomic species identification were also evaluated. Based on the six complete genomes of A. baumannii, the PBP genes were identified, and primers were designed for each gene. The nucleotide sequences of the genes identified that encode PBPs and the corresponding amino acid sequences were compared with those of ATCC 17978. Seven PBP genes and one monofunctional transglycosylase (MGT) gene were identified in the six genomes, encoding (i) four high-molecular-mass proteins (two of class A, PBP1a [ponA] and PBP1b [mrcB], and two of class B, PBP2 [pbpA or mrdA] and PBP3 [ftsI]), (ii) three low-molecular-mass proteins (two of type 5, PBP5/6 [dacC] and PBP6b [dacD], and one of type 7 (PBP7/8 [pbpG]), and (iii) a monofunctional enzyme (MtgA [mtgA]). Hot spot mutation regions were observed, although most of the allelic changes found translated into silent mutations. The amino acid consensus sequences corresponding to the PBP genes in the genomes and the clinical isolates were highly conserved. The changes found in amino acid sequences were associated with concrete clonal patterns but were not directly related to susceptibility or resistance to β-lactams. An insertion sequence disrupting the gene encoding PBP6b was identified in an endemic carbapenem-resistant clone in one of the participant hospitals.

Nucleotide sequence of the 5' noncoding region and part of the gag gene of mouse mammary tumor virus; identification of the 5' splicing site for subgenomic mRNAs.

We have determined the sequence of the first 1371 nucleotides at the 5' end of the genome of mouse mammary tumor virus using molecularly cloned proviral DNA of the GR virus strain. The most likely initiation codon used for the gag gene of mouse mammary tumor virus is the first one, located 312 nucleotides from the 5' end of the viral RNA. The 5' splicing site for the subgenomic mRNA's is located approximately 288 nucleotides downstream from the 5' end of the viral RNA. From the DNA sequence the amino acid sequence of the N-terminal half of the gag precursor protein, including p10 and p21, was deduced (353 amino acids).

Purification and characterization of two enantioselective alpha-ketoglutarate-dependent dioxygenases, RdpA and SdpA, from Sphingomonas herbicidovorans MH.

Alpha-ketoglutarate-dependent (R)-dichlorprop dioxygenase (RdpA) and alpha-ketoglutarate-dependent (S)-dichlorprop dioxygenase (SdpA), which are involved in the degradation of phenoxyalkanoic acid herbicides in Sphingomonas herbicidovorans MH, were expressed and purified as His6-tagged fusion proteins from Escherichia coli BL21(DE3)(pLysS). RdpA and SdpA belong to subgroup II of the alpha-ketoglutarate-dependent dioxygenases and share the specific motif HXDX(24)TX(131)HX(10)R. Amino acids His-111, Asp-113, and His-270 and amino acids His-102, Asp-104, and His 257 comprise the 2-His-1-carboxylate facial triads and were predicted to be involved in iron binding in RdpA and SdpA, respectively. RdpA exclusively transformed the (R) enantiomers of mecoprop [2-(4-chloro-2-methylphenoxy)propanoic acid] and dichlorprop [2-(2,4-dichlorophenoxy)propanoic acid], whereas SdpA was specific for the (S) enantiomers. The apparent Km values were 99 microM for (R)-mecoprop, 164 microM for (R)-dichlorprop, and 3 microM for alpha-ketoglutarate for RdpA and 132 microM for (S)-mecoprop, 495 microM for (S)-dichlorprop, and 20 microM for alpha-ketoglutarate for SdpA. Both enzymes had high apparent Km values for oxygen; these values were 159 microM for SdpA and >230 microM for RdpA, whose activity was linearly dependent on oxygen at the concentration range measured. Both enzymes had narrow cosubstrate specificity; only 2-oxoadipate was able to replace alpha-ketoglutarate, and the rates were substantially diminished. Ferrous iron was necessary for activity of the enzymes, and other divalent cations could not replace it. Although the results of growth experiments suggest that strain MH harbors a specific 2,4-dichlorophenoxyacetic acid-converting enzyme, tfdA-, tfdAalpha-, or cadAB-like genes were not discovered in a screening analysis in which heterologous hybridization and PCR were used.

Comprehensive analysis of Arabidopsis expression level polymorphisms with simple inheritance.

In Arabidopsis thaliana, gene expression level polymorphisms (ELPs) between natural accessions that exhibit simple, single locus inheritance are promising quantitative trait locus (QTL) candidates to explain phenotypic variability. It is assumed that such ELPs overwhelmingly represent regulatory element polymorphisms. However, comprehensive genome-wide analyses linking expression level, regulatory sequence and gene structure variation are missing, preventing definite verification of this assumption. Here, we analyzed ELPs observed between the Eil-0 and Lc-0 accessions. Compared with non-variable controls, 5' regulatory sequence variation in the corresponding genes is indeed increased. However, approximately 42% of all the ELP genes also carry major transcription unit deletions in one parent as revealed by genome tiling arrays, representing a &gt;4-fold enrichment over controls. Within the subset of ELPs with simple inheritance, this proportion is even higher and deletions are generally more severe. Similar results were obtained from analyses of the Bay-0 and Sha accessions, using alternative technical approaches. Collectively, our results suggest that drastic structural changes are a major cause for ELPs with simple inheritance, corroborating experimentally observed indel preponderance in cloned Arabidopsis QTL.

Functional analysis of Ectodysplasin-A mutations causing selective tooth agenesis.

Mutations of the Ectodysplasin-A (EDA) gene are generally associated with the syndrome hypohidrotic ectodermal dysplasia (MIM 305100), but they can also manifest as selective, non-syndromic tooth agenesis (MIM300606). We have performed an in vitro functional analysis of six selective tooth agenesis-causing EDA mutations (one novel and five known) that are located in the C-terminal tumor necrosis factor homology domain of the protein. Our study reveals that expression, receptor binding or signaling capability of the mutant EDA1 proteins is only impaired in contrast to syndrome-causing mutations, which we have previously shown to abolish EDA1 expression, receptor binding or signaling. Our results support a model in which the development of the human dentition, especially of anterior teeth, requires the highest level of EDA-receptor signaling, whereas other ectodermal appendages, including posterior teeth, have less stringent requirements and form normally in response to EDA mutations with reduced activity.

The region of mouse mammary tumor virus DNA containing the long terminal repeat includes a long coding sequence and signals for hormonally regulated transcription.

Starting from a biologically active recombinant DNA clone of exogenous unintegrated GR mouse mammary tumor virus, we have generated three subclones of PstI fragments of 1.45, 1.1, and 2.0 kb in the plasmid vector PBR322. The nucleotide sequence has been determined for the clone of 1.45 kb which includes almost the complete region of the long terminal repeat (LTR) plus an adjacent stretch of unique sequence DNA. A short region of the 2.0 kb clone, containing the beginning of the LTR, has also been sequenced. Starting with the A of an initiation codon outside the LTR, we detected an open reading frame of 960 nucleotides, potentially coding for a protein of 320 amino acids (36K). Two hundred nucleotides downstream from the termination codon, and approximately 25 nucleotides upstream from the presumptive initiation site of viral RNA synthesis, we found a promoter-like sequence. The sequence AGTAAA was detected approximately 15-20 nucleotides upstream from the 3' end of virion RNA and probably serves as a polyadenylation signal. The 1.45 kb PstI fragment has been transfected into Ltk- cells together with a plasmid containing the thymidine kinase gene of herpes simplex virus. The virus-specific RNA synthesis detected in a Tk+ cell clone was strongly stimulated by the addition of dexamethasone.

Usefulness of double locus sequence typing (DLST) for regional and international epidemiological surveillance of methicilin-resistant Staphylococcus aureus.

Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of nosocomial infections worldwide. To differentiate reliably among S. aureus isolates, we recently developed double locus sequence typing (DLST) based on the analysis of partial sequences of clfB and spa genes. In the present study, we evaluated the usefulness of DLST for epidemiological investigations of MRSA by routinely typing 1242 strains isolated in Western Switzerland. Additionally, particular local and international collections were typed by pulsed field gel electrophoresis (PFGE) and DLST to check the compatibility of DLST with the results obtained by PFGE, and for international comparisons. Using DLST, we identified the major MRSA clones of Western Switzerland, and demonstrated the close relationship between local and international clones. The congruence of 88% between the major PFGE and DLST clones indicated that our results obtained by DLST were compatible with earlier results obtained by PFGE. DLST could thus easily be incorporated in a routine surveillance procedure. In addition, the unambiguous definition of DLST types makes this method more suitable than PFGE for long-term epidemiological surveillance. Finally, the comparison of the results obtained by DLST, multilocus sequence typing, PFGE, Staphylococcal cassette chromosome mec typing and the detection of Panton-Valentine leukocidin genes indicated that no typing scheme should be used on its own. It is only the combination of data from different methods that gives the best chance of describing precisely the epidemiology and phylogeny of MRSA.

Complex organization of CTF/NF-I, C/EBP, and HNF3 binding sites within the promoter of the liver-specific vitellogenin gene.

Vitellogenin genes are expressed specifically in the liver of female oviparous vertebrates under the strict control of estrogen. To explain this tissue-specific expression, we performed a detailed analysis of the Xenopus laevis vitellogenin gene B1 promoter by DNase I footprinting and gel mobility-shift assays. We characterized five binding sites for the ubiquitous factor CTF/NF-I. Two of these sites are close to the TATA-box, whereas the others are located on both sides of the estrogen responsive unit formed by two imperfect estrogen response elements. Moreover two liver-enriched factors, C/EBP and HNF3, were found to interact with multiple closely spaced proximal promoter elements in the first 100 base pairs upstream of the TATA-box. To confirm the physiological significance of this in vitro analysis, in vivo DNase I footprinting experiments were carried out using the ligation-mediated polymerase chain reaction technique. The various cis-elements characterized in vitro as binding sites for known transcription factors and more particularly for liver-enriched transcription factors are efficiently recognized in vivo as well, suggesting that they play an important role in the control of the liver-specific vitellogenin gene B1 expression.