Assessment of transcript reconstruction methods for RNA-seq.

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.

Nonshivering thermogenesis capacity associated to mitochondrial DNA haplotypes and gender in the greater white-toothed shrew, Crocidura russula.

A selection gradient was recently suggested as one possible cause for a clinal distribution of mitochondrial DNA (mtDNA) haplotypes along an altitudinal transect in the greater white-toothed shrew, Crocidura russula (Ehinger et al. 2002). One mtDNA haplotype (H1) rare in lowland, became widespread when approaching the altitudinal margin of the distribution. As H1 differs from the main lowland haplotype by several nonsynonymous mutations (including on ATP6), and as mitochondria play a crucial role in metabolism and thermogenesis, distribution patterns might stem from differences in the thermogenic capacity of different mtDNA haplotypes. In order to test this hypothesis, we measured the nonshivering thermogenesis (NST) associated with different mtDNA haplotypes. Sixty-two shrews, half of which had the H1 haplotype, were acclimated in November at semioutdoor conditions and measured for NST throughout winter. Our results showed the crucial role of NST for winter survival in C. russula. The individuals that survived winter displayed a higher significant increase in NST during acclimation, associated with a significant gain in body mass, presumably from brown fat accumulation. The NST capacity (ratio of NST to basal metabolic rate) was exceptionally high for such a small species. NST was significantly affected by a gender x haplotype interaction after winter-acclimation: females bearing the H1 haplotype displayed a better thermogenesis at the onset of the breeding season, while the reverse was true for males. Altogether, our results suggest a sexually antagonistic cyto-nuclear selection on thermogenesis.

The TPTE gene family: cellular expression, subcellular localization and alternative splicing.

The human TPTE (Transmembrane Phosphatase with TEnsin homology) gene family encodes a PTEN-related tyrosine phosphatase with four potential transmembrane domains. Chromosomal mapping revealed multiple copies of the TPTE gene on chromosomes 13, 15, 21, 22 and Y. Human chromosomes 13 and 21 copies encode two functional proteins, TPIP (TPTE and PTEN homologous Inositol lipid Phosphatase) and TPTE, respectively, whereas only one copy of the gene exists in the mouse genome. In the present study, we show that TPTE and TPIP proteins are expressed in secondary spermatocytes and/or prespermatids. In addition, we report the existence of several novel alternatively spliced isoforms of these two proteins with variable number of transmembrane domains. The latter has no influence on the subcellular localization of these different peptides as shown by co-immunofluorescence experiments. Finally, we identify another expressed TPTE copy, mapping to human chromosome 22, whose transcription appears to be under the control of the LTR of human endogenous retrovirus RTVL-H3.

A large deletion in the adRP gene PRPF31: evidence that haploinsufficiency is the cause of disease.

PURPOSE: To report a large deletion that encompasses more than 90% of PRPF31 gene and two other neighboring genes in their entirety in an adRP pedigree that appears to show only the typical clinical features of retinitis pigmentosa. METHODS: To identify PRPF31 mutation in a dominant RP family (ADRP2) previously linked to the RP11 locus, the 14 exons of PRPF31 were screened for mutations by direct sequencing. To investigate the possibility of a large deletion, microsatellite markers near PRPF31 gene were analyzed by non-denaturing PAGE. RESULTS: Initial screening of PRPF31 gene in the ADRP2 family did not reveal an obvious mutation. A large deletion was however suspected due to lack of heterozygosity for nearly all PRPF31 intragenic single nucleotide polymorphysm (SNPs). In order to estimate the size of the deletion, SNPs and microsatellite markers spanning and flanking PRPF31 were analyzed in the entire ADRP2 family. Haplotype analysis with the above markers suggested a deletion of approximately 30 kb that included the putative promoter region of a novel gene OSCAR, the entire genomic content of genes NDUFA3, TFPT and more than 90% of PRPF31 gene. Sequence analysis of the region flanking the potential deletion showed a high presence of Alu elements implicating Alu mediated recombination as the mechanism responsible for this event. CONCLUSIONS: This mutation provides evidence that haploinsufficiency rather than aberrant function of mutated proteins is the cause of disease in these adRP patients with mutations in PRPF31 gene.

Altered expression of the transcription factor Mef2c during retinal degeneration in Rpe65-/- mice.

Purpose. To investigate the role of the myocyte enhancer factor 2 (Mef2) transcription factor family in retinal diseases, Mef2c expression was assessed during retinal degeneration in the Rpe65(-/-) mouse model of Leber's congenital amaurosis (LCA). Mef2c-dependent expression of photoreceptor-specific genes was further addressed. Methods. Expression of Mef2 members was analyzed by oligonucleotide microarray, quantitative PCR (qPCR) and in situ hybridization. Mef2c-dependent transcriptional activity was assayed by luciferase assay in HEK293T cells. Results. Mef2c was the only Mef2 member markedly downregulated during retinal degeneration in Rpe65(-/-) mice. Mef2c mRNA level was decreased by more than 2 fold at 2 and 4 months and by 3.5 fold at 6 months in retinas of Rpe65(-/-) mice. Downregulation of Mef2c at the protein level was confirmed in Rpe65(-/-) retinas. The decrease in Mef2c mRNA levels in the developing Rpe65(-/-) retinas, from post-natal day (P)13 onward, was concomitant with the decreased expression of the rod-specific transcription factors Nrl and Nr2e3. Nrl was further shown to drive Mef2c transcriptional activity, supporting a physiological role for Mef2c in the retina. In addition, Mef2c appeared to act as a transcriptional repressor of its own expression, as well as those of the retina-specific retinal G-protein coupled receptor (Rgr), rhodopsin and M-opsin genes. Conclusions. These findings highlight the early altered regulation of the rod-specific transcriptional network in Rpe65-related disease. They further indicate that Mef2c may act as a novel transcription factor involved in the development and the maintenance of photoreceptor cells.

Temporal patterns of nucleotide misincorporations and DNA fragmentation in ancient DNA.

DNA that survives in museum specimens, bones and other tissues recovered by archaeologists is invariably fragmented and chemically modified. The extent to which such modifications accumulate over time is largely unknown but could potentially be used to differentiate between endogenous old DNA and present-day DNA contaminating specimens and experiments. Here we examine mitochondrial DNA sequences from tissue remains that vary in age between 18 and 60,000 years with respect to three molecular features: fragment length, base composition at strand breaks, and apparent C to T substitutions. We find that fragment length does not decrease consistently over time and that strand breaks occur preferentially before purine residues by what may be at least two different molecular mechanisms that are not yet understood. In contrast, the frequency of apparent C to T substitutions towards the 5'-ends of molecules tends to increase over time. These nucleotide misincorporations are thus a useful tool to distinguish recent from ancient DNA sources in specimens that have not been subjected to unusual or harsh treatments.

Loss-of-Function Mutations in PTPN11 Cause Metachondromatosis, but Not Ollier Disease or Maffucci Syndrome.

Metachondromatosis (MC) is a rare, autosomal dominant, incompletely penetrant combined exostosis and enchondromatosis tumor syndrome. MC is clinically distinct from other multiple exostosis or multiple enchondromatosis syndromes and is unlinked to EXT1 and EXT2, the genes responsible for autosomal dominant multiple osteochondromas (MO). To identify a gene for MC, we performed linkage analysis with high-density SNP arrays in a single family, used a targeted array to capture exons and promoter sequences from the linked interval in 16 participants from 11 MC families, and sequenced the captured DNA using high-throughput parallel sequencing technologies. DNA capture and parallel sequencing identified heterozygous putative loss-of-function mutations in PTPN11 in 4 of the 11 families. Sanger sequence analysis of PTPN11 coding regions in a total of 17 MC families identified mutations in 10 of them (5 frameshift, 2 nonsense, and 3 splice-site mutations). Copy number analysis of sequencing reads from a second targeted capture that included the entire PTPN11 gene identified an additional family with a 15 kb deletion spanning exon 7 of PTPN11. Microdissected MC lesions from two patients with PTPN11 mutations demonstrated loss-of-heterozygosity for the wild-type allele. We next sequenced PTPN11 in DNA samples from 54 patients with the multiple enchondromatosis disorders Ollier disease or Maffucci syndrome, but found no coding sequence PTPN11 mutations. We conclude that heterozygous loss-of-function mutations in PTPN11 are a frequent cause of MC, that lesions in patients with MC appear to arise following a "second hit," that MC may be locus heterogeneous since 1 familial and 5 sporadically occurring cases lacked obvious disease-causing PTPN11 mutations, and that PTPN11 mutations are not a common cause of Ollier disease or Maffucci syndrome.

Dissection of the complex phenotype in cuticular mutants of Arabidopsis reveals a role of SERRATE as a mediator.

Mutations in LACERATA (LCR), FIDDLEHEAD (FDH), and BODYGUARD (BDG) cause a complex developmental syndrome that is consistent with an important role for these Arabidopsis genes in cuticle biogenesis. The genesis of their pleiotropic phenotypes is, however, poorly understood. We provide evidence that neither distorted depositions of cutin, nor deficiencies in the chemical composition of cuticular lipids, account for these features, instead suggesting that the mutants alleviate the functional disorder of the cuticle by reinforcing their defenses. To better understand how plants adapt to these mutations, we performed a genome-wide gene expression analysis. We found that apparent compensatory transcriptional responses in these mutants involve the induction of wax, cutin, cell wall, and defense genes. To gain greater insight into the mechanism by which cuticular mutations trigger this response in the plants, we performed an overlap meta-analysis, which is termed MASTA (MicroArray overlap Search Tool and Analysis), of differentially expressed genes. This suggested that different cell integrity pathways are recruited in cesA cellulose synthase and cuticular mutants. Using MASTA for an in silico suppressor/enhancer screen, we identified SERRATE (SE), which encodes a protein of RNA-processing multi-protein complexes, as a likely enhancer. In confirmation of this notion, the se lcr and se bdg double mutants eradicate severe leaf deformations as well as the organ fusions that are typical of lcr and bdg and other cuticular mutants. Also, lcr does not confer resistance to Botrytis cinerea in a se mutant background. We propose that there is a role for SERRATE-mediated RNA signaling in the cuticle integrity pathway.

Phosphorylation regulates FOXC2-mediated transcription in lymphatic endothelial cells.

One of the key mechanisms linking cell signaling and control of gene expression is reversible phosphorylation of transcription factors. FOXC2 is a forkhead transcription factor that is mutated in the human vascular disease lymphedema-distichiasis and plays an essential role in lymphatic vascular development. However, the mechanisms regulating FOXC2 transcriptional activity are not well understood. We report here that FOXC2 is phosphorylated on eight evolutionarily conserved proline-directed serine/threonine residues. Loss of phosphorylation at these sites triggers substantial changes in the FOXC2 transcriptional program. Through genome-wide location analysis in lymphatic endothelial cells, we demonstrate that the changes are due to selective inhibition of FOXC2 recruitment to chromatin. The extent of the inhibition varied between individual binding sites, suggesting a novel rheostat-like mechanism by which expression of specific genes can be differentially regulated by FOXC2 phosphorylation. Furthermore, unlike the wild-type protein, the phosphorylation-deficient mutant of FOXC2 failed to induce vascular remodeling in vivo. Collectively, our results point to the pivotal role of phosphorylation in the regulation of FOXC2-mediated transcription in lymphatic endothelial cells and underscore the importance of FOXC2 phosphorylation in vascular development.

Asymmetric and differential gene introgression at a contact zone between two highly divergent lineages of field voles (Microtus agrestis).

Secondary contact zones have the potential to shed light on the mode and rate at which reproductive isolation accumulates during allopatric speciation. We investigated the population genetics of a contact zone between two highly divergent lineages of field voles (Microtus agrestis) in the Swiss Jura mountains. To shed light on the processes underlying introgression, we used maternally, paternally, and bi-parentally inherited markers. Though the two lineages maintained a strong genetic structure, we found some hybrids and evidence of gene flow. The extent of introgression varied with the mode of inheritance, being highest for mtDNA and absent for the Y chromosome. In addition, introgression was asymmetric, occurring only from the Northern to the Southern lineage. Both patterns seem parsimoniously explained by neutral processes linked to differences in effective sizes and sex-biased dispersal rates. The lineage with lower effective population size was also the more introgressed, and the mode-of-inheritance effect correlated with the male-biased dispersal rate of microtine rodents. We cannot exclude, however, that Haldane's effect contributed to the latter, as we found a marginally significant deficit in males (the heterogametic sex) among hybrids. We propose a possible demographic scenario to account for the patterns documented, and empirical extensions to further investigate this contact zone.

The peroxisome proliferator-activated receptor alpha regulates amino acid metabolism.

The peroxisome proliferator-activated receptor alpha is a ligand-activated transcription factor that plays an important role in the regulation of lipid homeostasis. PPARalpha mediates the effects of fibrates, which are potent hypolipidemic drugs, on gene expression. To better understand the biological effects of fibrates and PPARalpha, we searched for genes regulated by PPARalpha using oligonucleotide microarray and subtractive hybridization. By comparing liver RNA from wild-type and PPARalpha null mice, it was found that PPARalpha decreases the mRNA expression of enzymes involved in the metabolism of amino acids. Further analysis by Northern blot revealed that PPARalpha influences the expression of several genes involved in trans- and deamination of amino acids, and urea synthesis. Direct activation of PPARalpha using the synthetic PPARalpha ligand WY14643 decreased mRNA levels of these genes, suggesting that PPARalpha is directly implicated in the regulation of their expression. Consistent with these data, plasma urea concentrations are modulated by PPARalpha in vivo. It is concluded that in addition to oxidation of fatty acids, PPARalpha also regulates metabolism of amino acids in liver, indicating that PPARalpha is a key controller of intermediary metabolism during fasting.

Type I signal peptidase from Leishmania is a target of the immune response in human cutaneous and visceral leishmaniasis.

The gene encoding type I signal peptidase (Lmjsp) has been cloned from Leishmania major. Lmjsp encodes a protein of 180 amino residues with a predicted molecular mass of 20.5 kDa. Comparison of the protein sequence with those of known type I signal peptidases indicates homology in five conserved domains A-E which are known to be important, or essential, for catalytic activity. Southern blot hybridisation analysis indicates that there is a single copy of the Lmjsp gene. A recombinant SPase protein and a synthetic peptide of the L. major signal peptidase were used to examine the presence of specific antibodies in sera from either recovered or active individuals of both cutaneous and visceral leishmaniasis. This evaluation demonstrated that sera from cutaneous and visceral forms of leishmaniasis are highly reactive to both the recombinant and synthetic signal peptidase antigens. Therefore, the Leishmania signal peptidase, albeit localised intracellularly, is a significant target of the Leishmania specific immune response and highlights its potential use for serodiagnosis of cutaneous and visceral leishmaniasis.

Secreted subtilisin gene family in Trichophyton rubrum.

Secreted proteases constitute potential virulence factors of dermatophytes. A total of seven genes encoding putative serine proteases of the subtilisin family (SUB) were isolated in Trichophyton rubrum. Based on sequence data and intron-exon structure, a phylogenetic analysis of subtilisins from T. rubrum and other fungi revealed a presumed ancestral lineage comprising T. rubrum SUB2 and Aspergillus SUBs. All other SUBs (SUB1, SUB3-7) are dermatophyte-specific and have apparently emerged more recently, through successive gene duplication events. We showed that two subtilisins, Sub3 and Sub4, were detected in culture supernatants of T. rubrum grown in a medium containing soy protein as a sole nitrogen source. Both recombinant enzymes produced in Pichia pastoris are highly active on keratin azure suggesting that these proteases play an important role in invasion of keratinised tissues by the fungus. The set of deduced amino acid sequences of T. rubrum SUB ORFs allowed the identification of orthologous Subs secreted by other dermatophyte species using proteolysis and mass spectrometry.

A species delimitation approach in the Trochulus sericeus/hispidus complex reveals two cryptic species within a sharp contact zone.

BACKGROUND: Mitochondrial DNA sequencing increasingly results in the recognition of genetically divergent, but morphologically cryptic lineages. Species delimitation approaches that rely on multiple lines of evidence in areas of co-occurrence are particularly powerful to infer their specific status. We investigated the species boundaries of two cryptic lineages of the land snail genus Trochulus in a contact zone, using mitochondrial and nuclear DNA marker as well as shell morphometrics. RESULTS: Both mitochondrial lineages have a distinct geographical distribution with a small zone of co-occurrence. In the same area, we detected two nuclear genotype clusters, each being highly significantly associated to one mitochondrial lineage. This association however had exceptions: a small number of individuals in the contact zone showed intermediate genotypes (4%) or cytonuclear disequilibrium (12%). Both mitochondrial lineage and nuclear cluster were statistically significant predictors for the shell shape indicating morphological divergence. Nevertheless, the lineage morphospaces largely overlapped (low posterior classification success rate of 69% and 78%, respectively): the two lineages are truly cryptic. CONCLUSION: The integrative approach using multiple lines of evidence supported the hypothesis that the investigated Trochulus lineages are reproductively isolated species. In the small contact area, however, the lineages hybridise to a limited extent. This detection of a hybrid zone adds an instance to the rare reported cases of hybridisation in land snails.

The sigma factor AlgU (AlgT) controls exopolysaccharide production and tolerance towards desiccation and osmotic stress in the biocontrol agent Pseudomonas fluorescens CHA0.

A variety of stress situations may affect the activity and survival of plant-beneficial pseudomonads added to soil to control root diseases. This study focused on the roles of the sigma factor AlgU (synonyms, AlgT, RpoE, and sigma(22)) and the anti-sigma factor MucA in stress adaptation of the biocontrol agent Pseudomonas fluorescens CHA0. The algU-mucA-mucB gene cluster of strain CHA0 was similar to that of the pathogens Pseudomonas aeruginosa and Pseudomonas syringae. Strain CHA0 is naturally nonmucoid, whereas a mucA deletion mutant or algU-overexpressing strains were highly mucoid due to exopolysaccharide overproduction. Mucoidy strictly depended on the global regulator GacA. An algU deletion mutant was significantly more sensitive to osmotic stress than the wild-type CHA0 strain and the mucA mutant were. Expression of an algU'-'lacZ reporter fusion was induced severalfold in the wild type and in the mucA mutant upon exposure to osmotic stress, whereas a lower, noninducible level of expression was observed in the algU mutant. Overexpression of algU did not enhance tolerance towards osmotic stress. AlgU was found to be essential for tolerance of P. fluorescens towards desiccation stress in a sterile vermiculite-sand mixture and in a natural sandy loam soil. The size of the population of the algU mutant declined much more rapidly than the size of the wild-type population at soil water contents below 5%. In contrast to its role in pathogenic pseudomonads, AlgU did not contribute to tolerance of P. fluorescens towards oxidative and heat stress. In conclusion, AlgU is a crucial determinant in the adaptation of P. fluorescens to dry conditions and hyperosmolarity, two major stress factors that limit bacterial survival in the environment.