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.

Transcription factor map alignment of promoter regions

We address the problem of comparing and characterizing the promoter regions of genes with similar expression patterns. This remains a challenging problem in sequence analysis, because often the promoter regions of co-expressed genes do not show discernible sequence conservation. In our approach, thus, we have not directly compared the nucleotide sequence of promoters. Instead, we have obtained predictions of transcription factor binding sites, annotated the predicted sites with the labels of the corresponding binding factors, and aligned the resulting sequences of labels—to which we refer here as transcription factor maps (TF-maps). To obtain the global pairwise alignment of two TF-maps, we have adapted an algorithm initially developed to align restriction enzyme maps. We have optimized the parameters of the algorithm in a small, but well-curated, collection of human–mouse orthologous gene pairs. Results in this dataset, as well as in an independent much larger dataset from the CISRED database, indicate that TF-map alignments are able to uncover conserved regulatory elements, which cannot be detected by the typical sequence alignments.

Pliocene and Pleistocene diversification and multiple refugia in a Eurasian shrew (Crocidura suaveolens group).

We sequenced 998 base pairs (bp) of mitochondrial DNA cytochrome b and 799 bp of nuclear gene BRCA1 in the Lesser white-toothed shrew (Crocidura suaveolens group) over its geographic range from Portugal to Japan. The aims of the study were to identify the main clades within the group and respective refugia resulting from Pleistocene glaciations. Analyses revealed the Asian lesser white-toothed shrew (C. shantungensis) as the basal clade, followed by a major branch of C. suaveolens, subdivided sensu stricto into six clades, which split-up in the Upper Pliocene and Lower Pleistocene (1.9-0.9 Myr). The largest clade, occurring over a huge range from east Europe to Mongolia, shows evidence of population expansion after a bottleneck. West European clades originated from Iberian and Italo-Balkanic refugia. In the Near East, three clades evolved in an apparent hotspot of refugia (west Turkey, south-west and south-east of the Caucasus). Most clades include specimens of different morphotypes and the validity of many taxa in the C. suaveolens group has to be re-evaluated.

Sequence determinants of a microtubule tip localization signal (MtLS).

Microtubule plus-end-tracking proteins (+TIPs) specifically localize to the growing plus-ends of microtubules to regulate microtubule dynamics and functions. A large group of +TIPs contain a short linear motif, SXIP, which is essential for them to bind to end-binding proteins (EBs) and target microtubule ends. The SXIP sequence site thus acts as a widespread microtubule tip localization signal (MtLS). Here we have analyzed the sequence-function relationship of a canonical MtLS. Using synthetic peptide arrays on membrane supports, we identified the residue preferences at each amino acid position of the SXIP motif and its surrounding sequence with respect to EB binding. We further developed an assay based on fluorescence polarization to assess the mechanism of the EB-SXIP interaction and to correlate EB binding and microtubule tip tracking of MtLS sequences from different +TIPs. Finally, we investigated the role of phosphorylation in regulating the EB-SXIP interaction. Together, our results define the sequence determinants of a canonical MtLS and provide the experimental data for bioinformatics approaches to carry out genome-wide predictions of novel +TIPs in multiple organisms.

Protonation controls ASIC1a activity via coordinated movements in multiple domains.

Acid-sensing ion channels (ASICs) are neuronal Na(+)-conducting channels activated by extracellular acidification. ASICs are involved in pain sensation, expression of fear, and neurodegeneration after ischemic stroke. Functional ASICs are composed of three identical or homologous subunits, whose extracellular part has a handlike structure. Currently, it is unclear how protonation of residues in extracellular domains controls ASIC activity. Knowledge of these mechanisms would allow a rational development of drugs acting on ASICs. Protonation may induce conformational changes that control the position of the channel gate. We used voltage-clamp fluorometry with fluorophores attached to residues in different domains of ASIC1a to detect conformational changes. Comparison of the timing of fluorescence and current signals identified residues involved in movements that preceded desensitization and may therefore be associated with channel opening or early steps leading to desensitization. Other residues participated in movements intimately linked to desensitization and recovery from desensitization. Fluorescence signals of all mutants were detected at more alkaline pH than ionic currents. Their midpoint of pH dependence was close to that of steady-state desensitization, whereas the steepness of the pH fluorescence relationship was closer to that of current activation. A sequence of movements was observed upon acidification, and its backward movements during recovery from desensitization occurred in the reverse order, indicating that the individual steps are interdependent. Furthermore, the fluorescence signal of some labeled residues in the finger domain was strongly quenched by a Trp residue in the neighboring β-ball domain. Upon channel activation, their fluorescence intensity increased, indicating that the finger moved away from the β ball. This extensive analysis of activity-dependent conformational changes in ASICs sheds new light on the mechanisms by which protonation controls ASIC activity.

Inhibition of estrogen-responsive gene activation by the retinoid X receptor beta: evidence for multiple inhibitory pathways.

The retinoid X receptor beta (RXR beta; H-2RIIBP) forms heterodimers with various nuclear hormone receptors and binds multiple hormone response elements, including the estrogen response element (ERE). In this report, we show that endogenous RXR beta contributes to ERE binding activity in nuclear extracts of the human breast cancer cell line MCF-7. To define a possible regulatory role of RXR beta regarding estrogen-responsive transcription in breast cancer cells, RXR beta and a reporter gene driven by the vitellogenin A2 ERE were transfected into estrogen-treated MCF-7 cells. RXR beta inhibited ERE-driven reporter activity in a dose-dependent and element-specific fashion. This inhibition occurred in the absence of the RXR ligand 9-cis retinoic acid. The RXR beta-induced inhibition was specific for estrogen receptor (ER)-mediated ERE activation because inhibition was observed in ER-negative MDA-MB-231 cells only following transfection of the estrogen-activated ER. No inhibition of the basal reporter activity was observed. The inhibition was not caused by simple competition of RXR beta with the ER for ERE binding, since deletion mutants retaining DNA binding activity but lacking the N-terminal or C-terminal domain failed to inhibit reporter activity. In addition, cross-linking studies indicated the presence of an auxiliary nuclear factor present in MCF-7 cells that contributed to RXR beta binding of the ERE. Studies using known heterodimerization partners of RXR beta confirmed that RXR beta/triiodothyronine receptor alpha heterodimers avidly bind the ERE but revealed the existence of another triiodothyronine-independent pathway of ERE inhibition. These results indicate that estrogen-responsive genes may be negatively regulated by RXR beta through two distinct pathways.

The CbrA-CbrB two-component regulatory system controls the utilization of multiple carbon and nitrogen sources in Pseudomonas aeruginosa.

A novel two-component system, CbrA-CbrB, was discovered in Pseudomonas aeruginosa; cbrA and cbrB mutants of strain PAO were found to be unable to use several amino acids (such as arginine, histidine and proline), polyamines and agmatine as sole carbon and nitrogen sources. These mutants were also unable to use, or used poorly, many other carbon sources, including mannitol, glucose, pyruvate and citrate. A 7 kb EcoRI fragment carrying the cbrA and cbrB genes was cloned and sequenced. The cbrA and cbrB genes encode a sensor/histidine kinase (Mr 108 379, 983 residues) and a cognate response regulator (Mr 52 254, 478 residues) respectively. The amino-terminal half (490 residues) of CbrA appears to be a sensor membrane domain, as predicted by 12 possible transmembrane helices, whereas the carboxy-terminal part shares homology with the histidine kinases of the NtrB family. The CbrB response regulator shows similarity to the NtrC family members. Complementation and primer extension experiments indicated that cbrA and cbrB are transcribed from separate promoters. In cbrA or cbrB mutants, as well as in the allelic argR9901 and argR9902 mutants, the aot-argR operon was not induced by arginine, indicating an essential role for this two-component system in the expression of the ArgR-dependent catabolic pathways, including the aruCFGDB operon specifying the major aerobic arginine catabolic pathway. The histidine catabolic enzyme histidase was not expressed in cbrAB mutants, even in the presence of histidine. In contrast, proline dehydrogenase, responsible for proline utilization (Pru), was expressed in a cbrB mutant at a level comparable with that of the wild-type strain. When succinate or other C4-dicarboxylates were added to proline medium at 1 mM, the cbrB mutant was restored to a Pru+ phenotype. Such a succinate-dependent Pru+ property was almost abolished by 20 mM ammonia. In conclusion, the CbrA-CbrB system controls the expression of several catabolic pathways and, perhaps together with the NtrB-NtrC system, appears to ensure the intracellular carbon: nitrogen balance in P. aeruginosa.

Multiple roles of mouse Numb in tuning developmental cell fates.

BACKGROUND: Notch signaling regulates multiple differentiation processes and cell fate decisions during both invertebrate and vertebrate development. Numb encodes an intracellular protein that was shown in Drosophila to antagonize Notch signaling at binary cell fate decisions of certain cell lineages. Although overexpression experiments suggested that Numb might also antagonize some Notch activity in vertebrates, the developmental processes in which Numb is involved remained elusive. RESULTS: We generated mice with a homozygous inactivation of Numb. These mice died before embryonic day E11.5, probably because of defects in angiogenic remodeling and placental dysfunction. Mutant embryos had an open anterior neural tube and impaired neuronal differentiation within the developing cranial central nervous system (CNS). In the developing spinal cord, the number of differentiated motoneurons was reduced. Within the peripheral nervous system (PNS), ganglia of cranial sensory neurons were formed. Trunk neural crest cells migrated and differentiated into sympathetic neurons. In contrast, a selective differentiation anomaly was observed in dorsal root ganglia, where neural crest--derived progenitor cells had migrated normally to form ganglionic structures, but failed to differentiate into sensory neurons. CONCLUSIONS: Mouse Numb is involved in multiple developmental processes and required for cell fate tuning in a variety of lineages. In the nervous system, Numb is required for the generation of a large subset of neuronal lineages. The restricted requirement of Numb during neural development in the mouse suggests that in some neuronal lineages, Notch signaling may be regulated independently of Numb.

Identification of Novel Craniofacial Regulatory Domains Located far Upstream of SOX9 and Disrupted in Pierre Robin Sequence.

Mutations in the coding sequence of SOX9 cause campomelic dysplasia (CD), a disorder of skeletal development associated with 46,XY disorders of sex development (DSDs). Translocations, deletions, and duplications within a ∼2 Mb region upstream of SOX9 can recapitulate the CD-DSD phenotype fully or partially, suggesting the existence of an unusually large cis-regulatory control region. Pierre Robin sequence (PRS) is a craniofacial disorder that is frequently an endophenotype of CD and a locus for isolated PRS at ∼1.2-1.5 Mb upstream of SOX9 has been previously reported. The craniofacial regulatory potential within this locus, and within the greater genomic domain surrounding SOX9, remains poorly defined. We report two novel deletions upstream of SOX9 in families with PRS, allowing refinement of the regions harboring candidate craniofacial regulatory elements. In parallel, ChIP-Seq for p300 binding sites in mouse craniofacial tissue led to the identification of several novel craniofacial enhancers at the SOX9 locus, which were validated in transgenic reporter mice and zebrafish. Notably, some of the functionally validated elements fall within the PRS deletions. These studies suggest that multiple noncoding elements contribute to the craniofacial regulation of SOX9 expression, and that their disruption results in PRS.

MP2RAGE Multiple Sclerosis Magnetic Resonance Imaging at 3 T.

OBJECTIVES: Lesion detection and characterization in multiple sclerosis (MS) are an essential part of its clinical diagnosis and an important research field. In this pilot study, we applied the recently introduced two inversion-contrast magnetization-prepared rapid gradient echo sequence (MP2RAGE) to patients with early-stage MS.¦MATERIALS AND METHODS: The MP2RAGE is a 3-dimensional (3D) magnetization-prepared rapid gradient echo derivative providing homogeneous T1 weighting and simultaneous T1 mapping. The MP2RAGE performance was compared with that of 2 clinical routine sequences (2D fluid-attenuated inversion recovery [FLAIR] and 3D magnetization-prepared rapid gradient echo [MP-RAGE]) and 2 state-of-the art clinical research sequences (the 3D FLAIR-SPACE [sampling perfection with application-optimized contrasts by using different flip-angle evolutions], a fluid-attenuated variable flip-angle fast spin echo technique, and the 3D double-inversion recovery SPACE). A cohort of 10 early-stage female MS patients (age, 31.6 ± 4.7 years; disease duration, 3.8 ± 1.9 years; median expanded disability status scale score, 1.75) and 10 age- and gender-matched controls were enrolled after approval of the local institutional review board was obtained. Multiple sclerosis lesions were identified and assigned to brain locations and tissue types by two experienced physicians in all 5 contrasts. Subsequently, lesions were manually delineated for comparison and statistical analysis of lesion count, volume and quantitative measures.¦RESULTS AND CONCLUSIONS: The results show that the 3D T1-weighted high-resolution MP2RAGE contrast provides a sensitive means for MS lesion assessment. The additional quantitative T1 relaxation time maps obtained with the MP2RAGE provide further potential diagnostic and prognostic information that could help (a) to better discriminate lesion subtypes and (b) to stage and predict the activity and the evolution of MS. Results also indicate that the T2-weighted double-inversion recovery and FLAIR-SPACE contrasts are attractive complements to the MP2RAGE for lesion detection.

Effect of bicarbonate and lactate buffer on glucose and lactate metabolism during hemodiafiltration in patients with multiple organ failure.

OBJECTIVE: To compare the effects of sodium bicarbonate and lactate for continuous veno-venous hemodiafiltration (CVVHDF) in critically ill patients. DESIGN AND SETTINGS: Prospective crossed-over controlled trial in the surgical and medical ICUs of a university hospital. PATIENTS: Eight patients with multiple organ dysfunction syndrome (MODS) requiring CVVHDF. INTERVENTION: Each patient received the two buffers in a randomized sequence over two consecutive days. MEASUREMENTS AND RESULTS: The following variables were determined: acid-base parameters, lactate production and utilization ((13)C lactate infusion), glucose turnover (6,6(2)H(2)-glucose), gas exchange (indirect calorimetry). No side effect was observed during lactate administration. Baseline arterial acid-base variables were equal with the two buffers. Arterial lactate (2.9 versus 1.5 mmol/l), glycemia (+18%) and glucose turnover (+23%) were higher in the lactate period. Bicarbonate and glucose losses in CVVHDF were substantial, but not lactate elimination. Infusing (13)C lactate increased plasma lactate levels equally with the two buffers. Lactate clearance (7.8+/-0.8 vs 7.5+/-0.8 ml/kg per min in the bicarbonate and lactate periods) and endogenous production rates (14.0+/-2.6 vs 13.6+/-2.6 mmol/kg per min) were similar. (13)C lactate was used as a metabolic substrate, as shown by (13)CO(2) excretion. Glycemia and metabolic rate increased significantly and similarly during the two periods during lactate infusion. CONCLUSION: Lactate was rapidly cleared from the blood of critically ill patients without acute liver failure requiring CVVHDF, being transformed into glucose or oxidized. Lactate did not exert undesirable effects, except moderate hyperglycemia, and achieved comparable effects on acid-base balance to bicarbonate.

Genome sequence of the pea aphid Acyrthosiphon pisum

Aphids are important agricultural pests and also biological models for studies of insect-plant interactions, symbiosis, virus vectoring, and the developmental causes of extreme phenotypic plasticity. Here we present the 464 Mb draft genome assembly of the pea aphid Acyrthosiphon pisum. This first published whole genome sequence of a basal hemimetabolous insect provides an outgroup to the multiple published genomes of holometabolous insects. Pea aphids are host-plant specialists, they can reproduce both sexually and asexually, and they have coevolved with an obligate bacterial symbiont. Here we highlight findings from whole genome analysis that may be related to these unusual biological features. These findings include discovery of extensive gene duplication in more than 2000 gene families as well as loss of evolutionarily conserved genes. Gene family expansions relative to other published genomes include genes involved in chromatin modification, miRNA synthesis, and sugar transport. Gene losses include genes central to the IMD immune pathway, selenoprotein utilization, purine salvage, and the entire urea cycle. The pea aphid genome reveals that only a limited number of genes have been acquired from bacteria; thus the reduced gene count of Buchnera does not reflect gene transfer to the host genome. The inventory of metabolic genes in the pea aphid genome suggests that there is extensive metabolite exchange between the aphid and Buchnera, including sharing of amino acid biosynthesis between the aphid and Buchnera. The pea aphid genome provides a foundation for post-genomic studies of fundamental biological questions and applied agricultural problems.

How Much Does It Cost? Optimization of Costs in Sequence Analysis of Social Science Data

One major methodological problem in analysis of sequence data is the determination of costs from which distances between sequences are derived. Although this problem is currently not optimally dealt with in the social sciences, it has some similarity with problems that have been solved in bioinformatics for three decades. In this article, the authors propose an optimization of substitution and deletion/insertion costs based on computational methods. The authors provide an empirical way of determining costs for cases, frequent in the social sciences, in which theory does not clearly promote one cost scheme over another. Using three distinct data sets, the authors tested the distances and cluster solutions produced by the new cost scheme in comparison with solutions based on cost schemes associated with other research strategies. The proposed method performs well compared with other cost-setting strategies, while it alleviates the justification problem of cost schemes.

Fast and simple epidemiological typing of Pseudomonas aeruginosa using the double-locus sequence typing (DLST) method.

Although the molecular typing of Pseudomonas aeruginosa is important to understand the local epidemiology of this opportunistic pathogen, it remains challenging. Our aim was to develop a simple typing method based on the sequencing of two highly variable loci. Single-strand sequencing of three highly variable loci (ms172, ms217, and oprD) was performed on a collection of 282 isolates recovered between 1994 and 2007 (from patients and the environment). As expected, the resolution of each locus alone [number of types (NT) = 35-64; index of discrimination (ID) = 0.816-0.964] was lower than the combination of two loci (NT = 78-97; ID = 0.966-0.971). As each pairwise combination of loci gave similar results, we selected the most robust combination with ms172 [reverse; R] and ms217 [R] to constitute the double-locus sequence typing (DLST) scheme for P. aeruginosa. This combination gave: (i) a complete genotype for 276/282 isolates (typability of 98%), (ii) 86 different types, and (iii) an ID of 0.968. Analysis of multiple isolates from the same patients or taps showed that DLST genotypes are generally stable over a period of several months. The high typability, discriminatory power, and ease of use of the proposed DLST scheme makes it a method of choice for local epidemiological analyses of P. aeruginosa. Moreover, the possibility to give unambiguous definition of types allowed to develop an Internet database ( http://www.dlst.org ) accessible by all.

Improved multiple-locus variable-number tandem-repeat assay for Staphylococcus aureus genotyping, providing a highly informative technique together with strong phylogenetic value.

We describe an improved multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) scheme for genotyping Staphylococcus aureus. We compare its performance to those of multilocus sequence typing (MLST) and spa typing in a survey of 309 strains. This collection includes 87 epidemic methicillin-resistant S. aureus (MRSA) strains of the Harmony collection, 75 clinical strains representing the major MLST clonal complexes (CCs) (50 methicillin-sensitive S. aureus [MSSA] and 25 MRSA), 135 nasal carriage strains (133 MSSA and 2 MRSA), and 13 published S. aureus genome sequences. The results show excellent concordance between the techniques' results and demonstrate that the discriminatory power of MLVA is higher than those of both MLST and spa typing. Two hundred forty-two genotypes are discriminated with 14 VNTR loci (diversity index, 0.9965; 95% confidence interval, 0.9947 to 0.9984). Using a cutoff value of 45%, 21 clusters are observed, corresponding to the CCs previously defined by MLST. The variability of the different tandem repeats allows epidemiological studies, as well as follow-up of the evolution of CCs and the identification of potential ancestors. The 14 loci can conveniently be analyzed in two steps, based upon a first-line simplified assay comprising a subset of 10 loci (panel 1) and a second subset of 4 loci (panel 2) that provides higher resolution when needed. In conclusion, the MLVA scheme proposed here, in combination with available on-line genotyping databases (including http://mlva.u-psud.fr/), multiplexing, and automatic sizing, can provide a basis for almost-real-time large-scale population monitoring of S. aureus.