Inference of chromosomal inversion dynamics from Pool-Seq data in natural and laboratory populations of Drosophila melanogaster.

Sequencing of pools of individuals (Pool-Seq) represents a reliable and cost-effective approach for estimating genome-wide SNP and transposable element insertion frequencies. However, Pool-Seq does not provide direct information on haplotypes so that, for example, obtaining inversion frequencies has not been possible until now. Here, we have developed a new set of diagnostic marker SNPs for seven cosmopolitan inversions in Drosophila melanogaster that can be used to infer inversion frequencies from Pool-Seq data. We applied our novel marker set to Pool-Seq data from an experimental evolution study and from North American and Australian latitudinal clines. In the experimental evolution data, we find evidence that positive selection has driven the frequencies of In(3R)C and In(3R)Mo to increase over time. In the clinal data, we confirm the existence of frequency clines for In(2L)t, In(3L)P and In(3R)Payne in both North America and Australia and detect a previously unknown latitudinal cline for In(3R)Mo in North America. The inversion markers developed here provide a versatile and robust tool for characterizing inversion frequencies and their dynamics in Pool-Seq data from diverse D. melanogaster populations.

Molecular evidence of Pleistocene bidirectional faunal exchange between Europe and the Near East: the case of the bicoloured shrew (Crocidura leucodon, Soricidae).

We sequenced 1077 bp of the mitochondrial cytochrome b gene and 511 bp of the nuclear Apolipoprotein B gene in bicoloured shrew (Crocidura leucodon, Soricidae) populations ranging from France to Georgia. The aims of the study were to identify the main genetic clades within this species and the influence of Pleistocene climatic variations on the respective clades. The mitochondrial analyses revealed a European clade distributed from France eastwards to north-western Turkey and a Near East clade distributed from Georgia to Romania; the two clades separated during the Middle Pleistocene. We clearly identified a population expansion after a bottleneck for the European clade based on mitochondrial and nuclear sequencing data; this expansion was not observed for the eastern clade. We hypothesize that the western population was confined to a small Italo-Balkanic refugium, whereas the eastern population subsisted in several refugia along the southern coast of the Black Sea.

Isolation of a novel monkey adenovirus reveals a new phylogenetic clade in the evolutionary history of simian adenoviruses

Adenoviruses of primates include human (HAdV) and simian (SAdV) isolates classified into 8 species (Human Adenovirus A to G, and Simian Adenovirus A). In this study, a novel adenovirus was isolated from a colony of cynomolgus macaques (Macaca fascicularis) and subcultured in VERO cells. Its complete genome was purified and a region encompassing the hexon gene, the protease gene, the DNA binding protein (DBP) and the 100 kDa protein was amplified by PCR and sequenced by primer walking. Sequence analysis of these four genes showed that the new isolate had 80% identity to other primate adenoviruses and lacked recombination events. The study of the evolutionary relationships of this new monkey AdV based on the combined sequences of the four genes supported a close relationship to SAdV-3 and SAdV-6, lineages isolated from Rhesus monkeys. The clade formed by these three types is separated from the remaining clades and establishes a novel branch that is related to species HAdV-A, F and G. However, the genetic distance corresponding to the newly isolated monkey AdV considerably differs from these as to belong to a new, not yet established species. Results presented here widen our knowledge on SAdV and represents an important contribution to the understanding of the evolutionary history of primate adenoviruses.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.

Regional and cellular gene expression changes in human Huntington's disease brain.

Huntington's disease (HD) pathology is well understood at a histological level but a comprehensive molecular analysis of the effect of the disease in the human brain has not previously been available. To elucidate the molecular phenotype of HD on a genome-wide scale, we compared mRNA profiles from 44 human HD brains with those from 36 unaffected controls using microarray analysis. Four brain regions were analyzed: caudate nucleus, cerebellum, prefrontal association cortex [Brodmann's area 9 (BA9)] and motor cortex [Brodmann's area 4 (BA4)]. The greatest number and magnitude of differentially expressed mRNAs were detected in the caudate nucleus, followed by motor cortex, then cerebellum. Thus, the molecular phenotype of HD generally parallels established neuropathology. Surprisingly, no mRNA changes were detected in prefrontal association cortex, thereby revealing subtleties of pathology not previously disclosed by histological methods. To establish that the observed changes were not simply the result of cell loss, we examined mRNA levels in laser-capture microdissected neurons from Grade 1 HD caudate compared to control. These analyses confirmed changes in expression seen in tissue homogenates; we thus conclude that mRNA changes are not attributable to cell loss alone. These data from bona fide HD brains comprise an important reference for hypotheses related to HD and other neurodegenerative diseases.

Proof of concept for microarray-based detection of DNA-binding oncogenes in cell extracts.

The function of DNA-binding proteins is controlled not just by their abundance, but mainly at the level of their activity in terms of their interactions with DNA and protein targets. Moreover, the affinity of such transcription factors to their target sequences is often controlled by co-factors and/or modifications that are not easily assessed from biological samples. Here, we describe a scalable method for monitoring protein-DNA interactions on a microarray surface. This approach was designed to determine the DNA-binding activity of proteins in crude cell extracts, complementing conventional expression profiling arrays. Enzymatic labeling of DNA enables direct normalization of the protein binding to the microarray, allowing the estimation of relative binding affinities. Using DNA sequences covering a range of affinities, we show that the new microarray-based method yields binding strength estimates similar to low-throughput gel mobility-shift assays. The microarray is also of high sensitivity, as it allows the detection of a rare DNA-binding protein from breast cancer cells, the human tumor suppressor AP-2. This approach thus mediates precise and robust assessment of the activity of DNA-binding proteins and takes present DNA-binding assays to a high throughput level.

Foodborne norovirus outbreak: the role of an asymptomatic food handler

Background: In July 2005 an outbreak of acute gastroenteritis occurred on a residential summer camp in the province of Barcelona (northeast of Spain). Forty-four people were affected among residents and employees. All of them had in common a meal at lunch time on 13 July (paella, round of beef and fruit). The aim of this study was to investigate a foodborne norovirus outbreak that occurred in the residential summer camp and in which the implication of a food handler was demonstrated by laboratory tests. Methods: A retrospective cohort study was designed. Personal or telephone interview was carried out to collect demographic, clinical and microbiological data of the exposed people, as well as food consumption in the suspected lunch. Food handlers of the mentioned summer camp were interviewed. Ten stool samples were requested from symptomatic exposed residents and the three food handlers that prepared the suspected food. Stools were tested for bacteries and noroviruses. Norovirus was detected using RT-PCR and sequence analysis. Attack rate, relative risks (RR) and its 95% confidence intervals (CI) were calculated to assess the association between food consumption and disease. Results: The global attack rate of the outbreak was 55%. The main symptoms were abdominal pain (90%), nausea (85%), vomiting (70%) and diarrhoea (42.5%). The disease remitted in 24-48 hours. Norovirus was detected in seven faecal samples, one of them was from an asymptomatic food handler who had not eaten the suspected food (round of beef), but cooked and served the lunch. Analysis of the two suspected foods isolated no pathogenic bacteria and detected no viruses. Molecular analysis showed that the viral strain was the same in ill patients and in the asymptomatic food handler (genotype GII.2 Melksham-like). Conclusions: In outbreaks of foodborne disease, the search for viruses in affected patients and all food handlers, even in those that are asymptomatic, is essential. Health education of food handlers with respect to hand washing should be promoted.

Molecular cophylogenetic relationships between European bats and their ectoparasitic mites (Acari, Spinturnicidae).

Cospeciation between host-parasite species is generally thought to result in mirror-image congruent phylogenies. Incongruence can be explained by mechanisms such as host switching, duplication, failure to speciate and sorting events. To investigate the level of association in the host-parasite relationship between Spinturnicid mites and their bat hosts, we constructed the phylogenetic tree of the genus Spinturnix (Acari, Mesostigmata) and compared it to the host phylogeny. We sequenced 938bp of the mitochondrial 16S rDNA and Cytochrome Oxydase subunit I (COI) genes among eleven morphospecies of Spinturnix collected on 20 European Vespertilionid and Rhinolophid bat species. Phylogenetic reconstruction of hosts and parasites showed statistical evidence for cospeciation and suggested that their evolutionary history involved also failure to speciate events and host switches. The latter seem to be mainly promoted by similar roosting habits of the host. As currently understood, host associations of Spinturnicid mites likely results from a complex interaction between the phylogenetic history of the host and the behaviour and the ecology of both parasite and host.

Additional data for nuclear DNA give new insights into the phylogenetic position of Sorex granarius within the Sorex araneus group.

Many species contain genetic lineages that are phylogenetically intermixed with those of other species. In the Sorex araneus group, previous results based on mtDNA and Y chromosome sequence data showed an incongruent position of Sorex granarius within this group. In this study, we explored the relationship between species within the S. araneus group, aiming to resolve the particular position of S. granarius. In this context, we sequenced a total of 2447 base pairs (bp) of X-linked and nuclear genes from 47 individuals of the S. araneus group. The same taxa were also analyzed within a Bayesian framework with nine autosomal microsatellites. These analyses revealed that all markers apart from mtDNA showed similar patterns, suggesting that the problematic position of S. granarius is best explained by an incongruent behavior by mtDNA. Given their close phylogenetic relationship and their close geographic distribution, the most likely explanation for this pattern is past mtDNA introgression from S. araneus race Carlit to S. granarius.

Tenascin-W, a new marker of cancer stroma, is elevated in sera of colon and breast cancer patients.

Tenascins are extracellular matrix proteins present during the development of organisms as well as in pathological conditions. Tenascin-W, the fourth and last member of the tenascin family remains the least well-characterized one. Our study aimed to evaluate the potential significance of tenascin-W as cancer biomarker by monitoring its presence in the serum of colorectal and breast cancer patients and its expression in colorectal tumor tissues. To measure serum tenascin-W levels, a sensitive sandwich-ELISA was established. Mean tenascin-W concentration in sera of patients with nonmetastatic colorectal cancer at time of diagnosis was highly increased compared to that of healthy volunteers. A similar tendency was observed for tenascin-C in the same patient cohort. However, the increase was much more striking for tenascin-W. We also detected elevated tenascin-W levels in sera of breast cancer patients. Furthermore, we could show a prominent expression of tenascin-W in extracts from colorectal tumor tissues by immunoblot analysis, whereas tenascin-W was not detectable in the corresponding normal colon mucosa. To confirm the western blot results, we performed immunohistochemistry of frozen sections of the same patients as well as of an additional, independently chosen collection of colorectal cancer tissues. In all cases, similarly to tenascin-C, tenascin-W was detected in the tumor stroma. Our results reveal a clear association between elevated levels of tenascin-W and the presence of cancer. These results warrant further studies to evaluate the potential value of serum and tissue tenascin-W levels as diagnostic, prognostic or monitoring biomarker in colorectal, breast and possibly other solid cancers.

The fasting-induced adipose factor/angiopoietin-like protein 4 is physically associated with lipoproteins and governs plasma lipid levels and adiposity.

Proteins secreted from adipose tissue are increasingly recognized to play an important role in the regulation of glucose metabolism. However, much less is known about their effect on lipid metabolism. The fasting-induced adipose factor (FIAF/angiopoietin-like protein 4/peroxisome proliferator-activated receptor gamma angiopoietin-related protein) was previously identified as a target of hypolipidemic fibrate drugs and insulin-sensitizing thiazolidinediones. Using transgenic mice that mildly overexpress FIAF in peripheral tissues we show that FIAF is an extremely powerful regulator of lipid metabolism and adiposity. FIAF overexpression caused a 50% reduction in adipose tissue weight, partly by stimulating fatty acid oxidation and uncoupling in fat. In addition, FIAF overexpression increased plasma levels of triglycerides, free fatty acids, glycerol, total cholesterol, and high density lipoprotein (HDL)-cholesterol. Functional tests indicated that FIAF overexpression severely impaired plasma triglyceride clearance but had no effect on very low density lipoprotein production. The effects of FIAF overexpression were amplified by a high fat diet, resulting in markedly elevated plasma and liver triglycerides, plasma free fatty acids, and plasma glycerol levels, and impaired glucose tolerance in FIAF transgenic mice fed a high fat diet. Remarkably, in mice the full-length form of FIAF was physically associated with HDL, whereas truncated FIAF was associated with low density lipoprotein. In human both full-length and truncated FIAF were associated with HDL. The composite data suggest that via physical association with plasma lipoproteins, FIAF acts as a powerful signal from fat and other tissues to prevent fat storage and stimulate fat mobilization. Our data indicate that disturbances in FIAF signaling might be involved in dyslipidemia.

Molecular clock is involved in predictive circadian adjustment of renal function.

Renal excretion of water and major electrolytes exhibits a significant circadian rhythm. This functional periodicity is believed to result, at least in part, from circadian changes in secretion/reabsorption capacities of the distal nephron and collecting ducts. Here, we studied the molecular mechanisms underlying circadian rhythms in the distal nephron segments, i.e., distal convoluted tubule (DCT) and connecting tubule (CNT) and the cortical collecting duct (CCD). Temporal expression analysis performed on microdissected mouse DCT/CNT or CCD revealed a marked circadian rhythmicity in the expression of a large number of genes crucially involved in various homeostatic functions of the kidney. This analysis also revealed that both DCT/CNT and CCD possess an intrinsic circadian timing system characterized by robust oscillations in the expression of circadian core clock genes (clock, bma11, npas2, per, cry, nr1d1) and clock-controlled Par bZip transcriptional factors dbp, hlf, and tef. The clock knockout mice or mice devoid of dbp/hlf/tef (triple knockout) exhibit significant changes in renal expression of several key regulators of water or sodium balance (vasopressin V2 receptor, aquaporin-2, aquaporin-4, alphaENaC). Functionally, the loss of clock leads to a complex phenotype characterized by partial diabetes insipidus, dysregulation of sodium excretion rhythms, and a significant decrease in blood pressure. Collectively, this study uncovers a major role of molecular clock in renal function.

Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa.

PURPOSE: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. METHODS: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. RESULTS: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease. CONCLUSIONS: The MDH1 gene is not the cause of RP28-linked arRP. Our experimental strategy shows that long-range genomic PCR followed by UHTs provides an excellent system to perform a thorough screening of candidate genes for hereditary retinal degeneration.

Plastid genome characterisation in Brassica and Brassicaceae using a new set of nine SSRs.

We report a new set of nine primer pairs specifically developed for amplification of Brassica plastid SSR markers. The wide utility of these markers is demonstrated for haplotype identification and detection of polymorphism in B. napus, B. nigra, B. oleracea, B. rapa and in related genera Arabidopsis, Camelina, Raphanus and Sinapis. Eleven gene regions (ndhB-rps7 spacer, rbcL-accD spacer, rpl16 intron, rps16 intron, atpB-rbcL spacer, trnE-trnT spacer, trnL intron, trnL-trnF spacer, trnM-atpE spacer, trnR-rpoC2 spacer, ycf3-psaA spacer) were sequenced from a range of Brassica and related genera for SSR detection and primer design. Other sequences were obtained from GenBank/EMBL. Eight out of nine selected SSR loci showed polymorphism when amplified using the new primers and a combined analysis detected variation within and between Brassica species, with the number of alleles detected per locus ranging from 5 (loci MF-6, MF-1) to 11 (locus MF-7). The combined SSR data were used in a neighbour-joining analysis (SMM, D (DM) distances) to group the samples based on the presence and absence of alleles. The analysis was generally able to separate plastid types into taxon-specific groups. Multi-allelic haplotypes were plotted onto the neighbour joining tree. A total number of 28 haplotypes were detected and these differentiated 22 of the 41 accessions screened from all other accessions. None of these haplotypes was shared by more than one species and some were not characteristic of their predicted type. We interpret our results with respect to taxon differentiation, hybridisation and introgression patterns relating to the 'Triangle of U'.