Search

Background Bovine Digital Dermatitis (BDD) is a prevalent infectious disease, causing painful foot skin lesions and lameness in cattle. We describe herein the bovine foot skin microbiota and its associations with BDD using 16S rRNA gene amplicon and shotgun metagenomic sequencing on samples from 259 dairy cows from three UK dairy farms. Results We show evidence of dysbiosis, and differences in taxonomy and functional profiles in the bovine foot skin microbiome of clinically healthy animals that subsequently develop BDD lesions, compared to those that do not. Our results suggest that taxonomical and functional differences together with alterations in ecological interactions between bacteria in the normal foot skin microbiome may predispose an animal to develop BDD lesions. Using genome-wide association and regional heritability mapping approaches, we provide first evidence for interactions between host genotype and certain members of the foot skin microbiota. We show the existence of significant genetic variation in the relative abundance of Treponema spp. and Peptoclostridium spp. and identify regions in the bovine genome that explain a significant proportion of this variation. Conclusions Collectively this work shows early changes in taxonomic and functional profiles of the bovine foot-skin microbiota in clinically healthy animals which are associated with subsequent development of BDD and could be relevant to prevention of disease. The description of host genetic control of members of the foot skin microbiota, combined with the association of the latter with BDD development offer new insights into a complex relationship that can be exploited in selective breeding programmes.

Energy storing tendons such as the human Achilles and equine superficial digital flexor tendon (SDFT) are prone to injury, with incidence increasing with ageing. The interfascicular matrix (IFM), which binds tendon fascicles, plays a key role in energy storing tendon mechanics, and ageing alterations to the IFM negatively impact tendon function. While the mechanical role of the IFM in tendon function is well-established, the biological role of IFM-resident cell populations remains to be elucidated. Therefore, the aim of this study was to identify IFM-resident cell populations and establish how these populations are affected by ageing. Cells from young and old SDFTs were subjected to single cell RNA-sequencing, and immunolabelling for markers of each resulting population used to localise cell clusters. Eleven cell clusters were identified, including tenocytes, endothelial cells, mural cells and immune cells. One tenocyte cluster localised to the fascicular matrix, whereas nine clusters localised to the IFM. Interfascicular tenocytes and mural cells were preferentially affected by ageing, with differential expression of genes related to senescence, dysregulated proteostasis and inflammation. This is the first study to uncover the importance of the IFM niche for a diverse range of cell populations, and to identify age-related alterations specific to IFM-localised cells.

Background Contagious Ovine Digital Dermatitis (CODD) is an emerging and common infectious foot disease of sheep which causes severe welfare and economic problems for the sheep industry. The aetiology of the disease is not fully understood and control of the disease is problematic. The aim of this study was to investigate the polybacterial aetiopathogenesis of CODD and the effects of antibiotic treatment, in a longitudinal study of an experimentally induced disease outbreak using a 16S rRNA gene amplicon sequencing approach. Results CODD was induced in 15/30 experimental sheep. During the development of CODD three distinct phenotypic lesion stages were observed. These were an initial interdigital dermatitis (ID) lesion, followed by a footrot (FR) lesion, then finally a CODD lesion. Distinct microbiota were observed for each lesion in terms of microbial diversity, clustering and composition. Porphyromonadaceae, Family XI, Veillonellaceae and Fusobacteriaceae were significantly associated with the diseased feet. Veillonellaceae and Fusobacteriaceae were most associated with the earlier stages of ID and footrot rather than CODD. Following antibiotic treatment of the sheep, the foot microbiota showed a strong tendency to return to the composition of the healthy state. The microbiota composition of CODD lesions collected by swab and biopsy methods were different. In particular, the Spirochaetaceae family were more abundant in samples collected by the biopsy method, suggesting that these bacteria are present in deeper tissues of the diseased foot. Conclusion In this study, CODD presented as part of a spectrum of poly-bacterial foot disease strongly associated with bacterial families Porphyromonadaceae, Family XI (a family in Clostridiales also known as Clostridium cluster XI), Veillonellaceae and Fusobacteriaceae which are predominately Gram-negative anaerobes. Following antibiotic treatment, the microbiome showed a strong tendency to return to the composition of the healthy state. The composition of the healthy foot microbiome does not influence susceptibility to CODD. Based on the data presented here and that CODD appears to be the severest end stage of sheep infectious foot disease lesions, better control of the initial ID and FR lesions would enable better control of CODD and enable better animal welfare.

BACKGROUND:Clostridium difficile is capable of causing severe enterocolitis in adults. The significance of toxin-producing C. difficile in children with diarrhea is unclear and practice differs on whether to institute treatment. We aimed to characterize the microbiome in relation to the presence of C. difficile and co-infection with other pathogens and to describe host response to infection. METHODS:Participants were children with acute diarrhea, 0-16 years of age, from whom stool samples had been submitted to the hospital laboratory for routine microbiology/virology. Convenience sampling was used for 50 prospective and 150 retrospective samples. No participants were treated for C. difficile. Rates of culture positivity for C. difficile, presence of toxin and PCR-ribotype were compared between age groups. Presence of other potential pathogens, comorbidities and complications were recorded. Microbiotal diversity was measured by 16S profiling. RESULTS:Nineteen of 77 (25%) children 2 years of age were C. difficile positive, of whom 10 (53%) and 9 (69%), respectively, carried toxigenic strains. Increased Shannon diversity was seen in children carrying C. difficile, with altered milieu. Presence of C. difficile was not associated with adverse clinical outcomes. In stools containing both Norovirus and C. difficile, there was increased relative abundance of verrucomicrobia. CONCLUSIONS:Children with diarrhea regularly carried toxigenic and non-toxigenic strains of C. difficile, demonstrating enhanced microbiotal diversity, and change in milieu, without apparent morbidity. This unexpected finding is contrary to that seen in adults with C. difficile disease.

Marine sponges are hosts to large, diverse communities of microorganisms. These microbiomes are distinct among sponge species and from seawater bacterial communities, indicating a key role of host identity in shaping its resident microbial community. However, the factors governing intraspecific microbiome variability are underexplored and may shed light on the evolutionary and ecological relationships between host and microbiome.Here, we examined the influence of genetic variation and geographic location on the composition of the Ircinia campana microbiome.We developed new microsatellite markers to genotype I. campana from two locations in the Florida Keys, USA, and characterized their microbiomes using V4 16S rRNA amplicon sequencing.We show that microbial community composition and diversity is influenced by host genotype, with more genetically similar sponges hosting more similar microbial communities. We also found that although I. campana was not genetically differentiated between sites, microbiome composition differed by location.Our results demonstrate that both host genetics and geography influence the composition of the sponge microbiome. Host genotypic influence on microbiome composition may be due to stable vertical transmission of the microbial community from parent to offspring, making microbiomes more similar by descent. Alternatively, sponge genotypic variation may reflect variation in functional traits that influence the acquisition of environmental microbes. This study reveals drivers of microbiome variation within and among locations, and shows the importance of intraspecific variability in mediating eco‐evolutionary dynamics of host‐associated microbiomes., For the first time, the authors show that intraspecific genetic variation affects microbiome composition in a marine sponge (Ircinia campana), with positive correlations observed between genetic and microbiome similarity. This has significant implications for our understanding of the ecological and evolutionary relationships between host and microbiome in this important model system.

The etiology of Crohn’s disease (CD) is multifactorial. Bacterial and fungal microbiota are involved in the onset and/or progression of the disease. A bacterial dysbiosis in CD patients is accepted; however, less is known about the mycobiome and the relationships between the two communities. We investigated the interkingdom relationships, their metabolic consequences, and the changes in the fungal community during relapse and remission in CD. Two cohorts were evaluated: a British cohort (n = 63) comprising CD and ulcerative colitis patients, and controls. The fungal and bacterial communities of biopsy and fecal samples were analyzed, with the fecal volatiles; datasets were also integrated; and a Dutch cohort (n = 41) comprising CD patients and healthy controls was analyzed for stability of the gut mycobiome. A dysbiosis of the bacterial community was observed in biopsies and stool. Results suggest Bacteroides is likely key in CD and may modulate Candida colonization. A dysbiosis of the fungal community was observed only in the Dutch cohort; Malassezia and Candida were increased in patients taking immunosuppressants. Longitudinal analysis showed an increase in Cyberlindnera in relapse. Saccharomyces was dominant in all fecal samples, but not in biopsies, some of which did not yield fungal reads; amino acid degradation was the main metabolic change associated with CD and both bacteria and fungi might be implicated. We have shown that Bacteroides and yeasts may play a role in CD; understanding their role and relationship in the disease would shed new light on the development and treatment of CD.

CRISPR-Cas immune systems are widespread in bacteria and archaea, but not ubiquitous. Previous work has demonstrated that CRISPR immunity is associated with an infection-induced fitness cost, which may help explain the patchy distribution observed. However, the mechanistic basis of this cost has remained unclear. Using Pseudomonas aeruginosa PA14 and its phage DMS3vir as a model, we perform a 30-day evolution experiment under phage mediated selection. We demonstrate that although CRISPR is initially selected for, bacteria carrying mutations in the phage receptor rapidly invade the population following subsequent reinfections. We then test three potential mechanisms for the observed cost of CRISPR: (1) autoimmunity from the acquisition of self-targeting spacers, (2) immunopathology or energetic costs from increased cas gene expression and (3) toxicity caused by phage gene expression prior to CRISPR-mediated cleavage. We find that phages can express genes before the immune system clears the infection and that expression of these genes can have a negative effect on host fitness. While infection does not lead to increased expression of cas genes, it does cause differential expression of multiple other host processes that may further contribute to the cost of CRISPR immunity. In contrast, we found little support for infection-induced autoimmunological and immunopathological effects. Phage gene expression prior to cleavage of the genome by the CRISPR-Cas immune system is therefore the most parsimonious explanation for the observed phage-induced fitness cost.

Background The dam is considered an important source of microbes for the calf; consequently, the development of calf microbiota may vary with farming system due to differences between the contact the calf has with the dam. The objective of this study was to characterise the early changes in the composition of oral and faecal microbiota in beef and dairy calves (N = 10) using high-throughput sequencing of the 16S rRNA gene. The microbiota of calves was compared to selected anatomical niches on their dams which were likely to contribute to the vertical transfer of microbes. Results A total of 14,125 amplicon sequence variants (ASVs) were identified and taxonomically assigned. The oral microbiota of calves and their dams were composed of more similar microbes after the first 4 weeks of life than immediately after calving. The faecal microbiota of four-week old calves was composed of microbes which were more similar to those found in the oral microbiota of calves and adult cows than the faecal microbiota of adult cows. Specific ASVs were identified in the oral microbiota of four-week old calves that were also present in cow niches at calving, whereas very few ASVs were present in the calf faecal microbiota at four-weeks of age were present in any adult cow niche at calving. These results were observed in both beef and dairy calves. Conclusions We did not observe any marked differences in the maturation of the oral and faecal microbiota between beef or dairy calves, despite dairy calves having very limited contact with their dam. This suggests the development of gastrointestinal microbiota in calves may not be affected by continued vertical transmission of microbes from the dam. Although the calf faecal microbiota changed over the first four-weeks of life, it was composed of microbes which were phylogenetically closer to those in the oral microbiota of calves and adult cows than the faeces of adult cows. There was little evidence of persistent microbial seeding of the calf faeces from anatomical niches on the cow at calving in either beef or dairy animals.

Phages play a major role in shaping the composition, evolution and function of bacterial communities. While bacteria and phages coexist in many natural environments, their coexistence is often short-lived in the lab due to the evolution of phage resistance. However, fitness costs associated with resistance and mutational loss of resistance alleles may limit the durability of acquired resistances, potentially allowing phages to re-invade the population. Here, we explore this idea in the context of bacteria that evolve CRISPR-based immunity against their phages. Consistent with previous studies, we found that the bacterium Pseudomonas aeruginosa PA14 evolved high levels of CRISPR-based immunity and low levels of surface-based resistance following infection with phage DMS3vir, which led to rapid phage extinction. However, when these pre-immunized bacterial populations were subsequently challenged with the same phage, they failed to clear the infection and instead stably coexisted with the phage. Analysis of bacterial genotypes and phenotypes over time explained why CRISPR-Cas immunity provides only a transient advantage: in the absence of phage (i.e. following the initial phage extinction) formerly CRISPR-immune bacteria regain sensitivity due to evolutionary loss of spacers, whereas in the presence of phage (i.e. upon reinfection) selection favours surface-based resistance over CRISPR immunity. The latter results from an infection-induced fitness cost of CRISPR-immunity that is due to phage gene expression prior to target DNA cleavage by the immune system. Together, these results show that CRISPR-Cas immune systems provide only a transient benefit to bacteria upon phage infection and help to explain why bacteria and phages can coexist in natural environments even when bacteria carry CRISPR-Cas adaptive immune systems that allow for rapid acquisition of immunity against phages.

Microbial ecology studies are often performed through extraction of metagenomic DNA followed by amplification and sequencing of a marker. It is known that each step may bias the results. These biases have been explored for the study of bacterial communities, but rarely for fungi. Our aim was therefore to evaluate methods for the study of the gut mycobiome. We first evaluated DNA extraction methods in fungal cultures relevant to the gut. Afterwards, to assess how these methods would behave with an actual sample, stool from a donor was spiked with cells from the same cultures. We found that different extraction kits favour some species and bias against others. In terms of amplicon sequencing, we evaluated five primer sets, two for ITS2 and one for ITS1, 18S and 28S rRNA. Results showed that 18S rRNA outperformed the other markers: it was able to amplify all the species in the mock community and to discriminate among them. ITS primers showed both amplification and sequencing biases, the latter related to the variable length of the product. We identified several biases in the characterisation of the gut mycobiome and showed how crucial it is to be aware of these before drawing conclusions from the results of these studies.

Antigenic variation is employed by many pathogens to evade the host immune response, andTrypanosoma bruceihas evolved a complex system to achieve this phenotype, involving sequential use of variant surface glycoprotein (VSG) genes encoded from a large repertoire of ∼2,000 alleles.T. bruceiexpress multiple, sometimes closely related, VSGs in a population at any one time, and the ability to resolve and analyse this diversity has been limited. We applied long read sequencing (PacBio) to VSG amplicons generated from blood extracted from batches of mice sacrificed at time points (days 3, 6, 10 and 12) post-infection with T. brucei TREU927. The data showed that long read sequencing is reliable for resolving allelic differences between VSGs, and demonstrated that there is significant expressed diversity (449 VSGs detected across 20 mice) and across the timeframe of study there was a clear semi-reproducible pattern of expressed diversity (median of 27 VSGs per sample at day 3 post infection (p.i.), 82 VSGs at day 6 p.i., 187 VSGs at day 10 p.i. and 132 VSGs by day 12 p.i.). There was also consistent detection of one VSG dominating expression across replicates at days 3 and 6, and emergence of a second dominant VSG across replicates by day 12. The innovative application of ecological diversity analysis to VSG reads enabled characterisation of hierarchical VSG expression in the dataset, and resulted in a novel method for analysing such patterns of variation. Additionally, the long read approach allowed detection of mosaic VSG expression from very few reads – this was observed as early as day 3, the earliest that such events have been detected. Therefore, our results indicate that long read analysis is a reliable tool for resolving diverse allele expression profiles, and provides novel insights into the complexity and nature of VSG expression in trypanosomes, revealing significantly higher diversity than previously shown and identifying mosaic gene formation unprecedentedly early during the infection process.

The authors would like to acknowledge funding from the Natural Environmental Research Council (NERC) grant NE/E006795/1 and the Pilot Competition grant number NBAF908 NBAF-L. PCR amplification, sequencing and taxonomic analysis was carried out at the NERC Biomolecular Analysis Facility (NBAF), Liverpool, UK. This work received funding from the MASTS pooling initiative (The Marine Alliance for Science and Technology for Scotland) and their support is gratefully acknowledged. MASTS is funded by the Scottish Funding Council (grant reference HR09011) and contributing institutions. Marine sediments are important sites for global biogeochemical cycling, mediated by macrofauna and microalgae. However, it is the microorganisms that drive these key processes. There is strong evidence that coastal benthic habitats will be affected by changing environmental variables (rising temperature, elevated CO2), and research has generally focused on the impact on macrofaunal biodiversity and ecosystem services. Despite their importance, there is less understanding of how microbial community assemblages will respond to environmental changes. In this study, a manipulative mesocosm experiment was employed, using next-generation sequencing to assess changes in microbial communities under future environmental change scenarios. Illumina sequencing generated over 11 million 16S rRNA gene sequences sequences (using a primer set biased towards bacteria) and revealed Bacteroidetes and Proteobacteria dominated the total bacterial community of sediment samples. In this study, the sequencing coverage and depth revealed clear changes in species abundance within some phyla. Bacterial community composition was correlated with simulated environmental conditions, and species level community composition was significantly influenced by the mean temperature of the environmental regime (p = 0.002), but not by variation in CO2 or diurnal temperature variation. Species level changes with increasing mean temperature corresponded with changes in NH4 concentration, suggesting there is no functional redundancy in microbial communities for nitrogen cycling. Marine coastal biogeochemical cycling under future environmental conditions is likely to be driven by changes in nutrient availability as a direct result of microbial activity. Publisher PDF

Several conditions associated with reduced gastric acid secretion confer an altered risk of developing a gastric malignancy. Helicobacter pylori-induced atrophic gastritis predisposes to gastric adenocarcinoma, autoimmune atrophic gastritis is a precursor of type I gastric neuroendocrine tumours, whereas proton pump inhibitor (PPI) use does not affect stomach cancer risk. We hypothesised that each of these conditions was associated with specific alterations in the gastric microbiota and that this influenced subsequent tumour risk. 95 patients (in groups representing normal stomach, PPI treated, H. pylori gastritis, H. pylori-induced atrophic gastritis and autoimmune atrophic gastritis) were selected from a cohort of 1400. RNA extracted from gastric corpus biopsies was analysed using 16S rRNA sequencing (MiSeq). Samples from normal stomachs and patients treated with PPIs demonstrated similarly high microbial diversity. Patients with autoimmune atrophic gastritis also exhibited relatively high microbial diversity, but with samples dominated by Streptococcus. H. pylori colonisation was associated with decreased microbial diversity and reduced complexity of co-occurrence networks. H. pylori-induced atrophic gastritis resulted in lower bacterial abundances and diversity, whereas autoimmune atrophic gastritis resulted in greater bacterial abundance and equally high diversity compared to normal stomachs. Pathway analysis suggested that glucose-6-phospahte1-dehydrogenase and D-lactate dehydrogenase were over represented in H. pylori-induced atrophic gastritis versus autoimmune atrophic gastritis, and that both these groups showed increases in fumarate reductase. Autoimmune and H. pylori-induced atrophic gastritis were associated with different gastric microbial profiles. PPI treated patients showed relatively few alterations in the gastric microbiota compared to healthy subjects., Author summary Different conditions such as autoimmune atrophic gastritis and Helicobacter pylori associated atrophic gastritis are associated with different types of gastric cancer, specifically neuroendocrine tumours and adenocarcinoma. Both conditions result in reduced gastric acid secretion, potentially allowing non-H. pylori bacteria to colonise the stomach. However patients receiving proton pump inhibitors (PPI) experience similar levels of acid secretion, but do not develop gastric cancer. The aims of this study were to investigate the contribution of non-H. pylori microbiota to gastric tumour development in the presence of reduced gastric acid secretion. 16S rRNA sequencing identified relatively few alterations in the gastric microbiota in patients receiving PPI therapy, despite reduced acid secretion, but more substantial alterations in those patents who had atrophic gastritis. Significant differences were also found between the patients who had atrophic gastritis of autoimmune and H. pylori associated types. Differences in biochemical pathways that potentially contribute to gastric tumorigenesis were also predicted. This work increases understanding of the mechanisms involved in gastric tumour development, and demonstrates how non-H. pylori bacteria may be important. This work may eventually lead to the development of novel chemopreventive therapies for stomach cancer that are based on altering the composition of the gastric microbiota.

Primary ciliary dyskinesia (PCD) is a genetic disorder, usually autosomal recessive, causing early respiratory disease and later subfertility. Whole exome sequencing may enable efficient analysis for locus heterogeneous disorders such as PCD. We whole-exome-sequenced one consanguineous Saudi Arabian with clinically diagnosed PCD and normal laterality, to attempt ab initio molecular diagnosis. We reviewed 13 known PCD genes and potentially autozygous regions (extended homozygosity) for homozygous exon deletions, non-dbSNP codon, splice-site base variants or small indels. Homozygous non-dbSNP changes were also reviewed exome-wide. One single molecular read representing RSPH9 p.Lys268del was observed, with no wild-type reads, and a notable deficiency of mapped reads at this location. Among all observations, RSPH9 was the strongest candidate for causality. Searching unmapped reads revealed seven more mutant reads. Direct assay for p.Lys268del (MboII digest) confirmed homozygosity in the affected individual, then confirmed homozygosity in three siblings with bronchiectasis. Our finding in southwest Saudi Arabia indicates that p.Lys268del, previously observed in two Bedouin families (Israel, UAE), is geographically widespread in the Arabian Peninsula. Analogous with cystic fibrosis CFTR p.Phe508del, screening for RSPH9 p.Lys268del (which lacks sentinel dextrocardia) in those at risk would help in early diagnosis, tailored clinical management, genetic counselling and primary prevention.

The gut microbiota profoundly affects the biology of its host. The composition of the microbiota is dynamic and is affected by both host genetic and many environmental effects. The gut microbiota of laboratory mice has been studied extensively, which has uncovered many of the effects that the microbiota can have. This work has also shown that the environments of different research institutions can affect the mouse microbiota. There has been relatively limited study of the microbiota of wild mice, but this has shown that it typically differs from that of laboratory mice (and that maintaining wild caught mice in the laboratory can quite quickly alter the microbiota). There is also inter-individual variation in the microbiota of wild mice, with this principally explained by geographical location. In this study we have characterised the gut (both the caecum and rectum) microbiota of wild caught Mus musculus domesticus at three UK sites and have investigated how the microbiota varies depending on host location and host characteristics. We find that the microbiota of these mice are generally consistent with those described from other wild mice. The rectal and caecal microbiotas of individual mice are generally more similar to each other, than they are to the microbiota of other individuals. We found significant differences in the diversity of the microbiotas among mice from different sample sites. There were significant correlations of microbiota diversity and body weight, a measure of age, body-mass index, serum concentration of leptin, and virus, nematode and mite infection.

Nature Communications 5: Article number: 3979 (2014); Published 9 June 2014; Updated 9 April 2015. The original version of this Article contained a typographical error in the spelling of the author Nicholas E. Ilott, which was incorrectly given as Nicholas E. IIott. This has now been corrected in both the PDF and HTML versions of the Article.

Down syndrome critical region gene 1-like 2 (DSCR1L2) belongs to the human DSCR1-like gene family, which also includes DSCR1 and DSCR1L1. Both DSCR1 and DSCR1L1 proteins interact with calcineurin, a calcium/calmodulin-dependent phosphatase. To date, no interactor has been described for DSCR1L2. The aim of this work was to perform a first functional study of DSCR1L2 using yeast two-hybrid analysis conducted on a human heart cDNA library. Here, we report the interaction between DSCR1L2 and the human cardiac troponin I (TNNI3), the heart-specific inhibitory subunit of the troponin complex, a central component of the contractile apparatus. This interaction was confirmed by both yeast cotransformation and GST (glutathione-sepharose transferase) fusion protein assay. Moreover, a new DSCR1L2 mRNA isoform, generated by alternative splicing, was identified and cloned in different tissues: it lacks two central exons, encoding the most conserved domains among the DSCR1-like protein family. A quantitative relative reverse transcription-polymerase chain reaction (RT-PCR) assay showed that in heart tissue the normalized expression level ratio for DSCR1L2 and DSCR1L2-E2E5 mRNA isoforms is 3.5 : 1, respectively. The yeast cotransformation and GST fusion protein assay demonstrated the interaction between this new DSCR1L2 variant and the human cardiac troponin I and the prominent role of DSCR1L2 exon 2 in determining binding between both DSCR1L2 isoforms and TNNI3. These data indicate an entirely new role for a DSCR1-like family gene, suggesting a possible involvement of DSCR1L2 in cardiac contraction.

A novel human gene has been identified by in-depth bioinformatics analysis of chromosome 21 segment 40/105 (21q21.1), with no coding region predicted in any previous analysis. Brain-derived DNA complementary to RNA (cDNA) sequencing predicts a 154-amino acid product with no similarity to any known protein. The gene has been named cysteine and tyrosine-rich protein 1 gene (symbol cysteine and tyrosine-rich 1, CYYR1). The CYYR1 messenger RNA was found by Northern blot analysis in a broad range of tissues (two transcripts of 3.4 and 2.2 kb). The gene consists of four exons and spans about 107 kb, including a very large intron of 85.8 kb. Analysis of expressed sequence tags shows high CYYR1 expression in cells belonging to the amine precursor uptake and decarboxylation system. We also cloned the cDNA of the murine ortholog Cyyr1, which was mapped by a radiation hybrid panel on chromosome 16 within the region corresponding to that containing the respective human homolog on chromosome 21. Sequence and phylogenetic analysis led to identification of several genes encoding CYYR1 homologous proteins. The most prominent feature identified in the protein family is a central, unique cysteine and tyrosine-rich domain, which is strongly conserved from lower vertebrates (fishes) to humans but is absent in bacteria and invertebrates.

A recently recognized gene family, conserved from yeast to humans, includes Down syndrome candidate region 1 gene (DSCR1), Adapt78 (recognized as the hamster ortholog of the DSCR1 isoform 4), ZAKI-4 (renamed DSCR1-like 1, DSCR1L1) and DSCR1L2 (a novel gene on human chromosome 1), along with yeast and C. elegans single members (Strippoli P., Lenzi L., Petrini M., Carinci P., Zannotti M., 2000. A new gene family including DSCR1 (Down Syndrome Candidate Region 1) and ZAKI-4: characterization from yeast to human and identification of DSCR1-like 2, a novel human member. Genomics 64, 252-263). The proposed family labels were a putative single-strand nucleic acid binding domain similar to the RNA recognition motif, and a unique, highly-conserved serine-proline motif. We have used a bioinformatics-driven molecular biology approach to characterize the murine members of DSCR1-like gene family. Systematic expressed-sequence-tags (EST) database search and reverse-transcription polymerase chain rection (RT-PCR) product sequencing allowed identification of the murine DSCR1, DSCR1L1 and DSCR1L2. The sequences of the respective protein products are of 198, 197 and 241 amino acids, respectively, and are very similar to the corresponding human proteins. The very broad expression pattern of the murine DSCR1 genes is similar to that of the human genes. Using a radiation hybrid panel, we mapped the murine DSCR1-like family members. The murine DSCR1 ortholog is located on the chromosome 16, in a region corresponding to that on human chromosome 21 just upstream of the Down syndrome candidate region. DSCR1L1 and DSCR1L2 murine genes are also located in chromosomal segments of chromosome 17 and 4, respectively, exactly corresponding to those containing the respective human homologs on chromosomes 6 and 1. Description of the mouse orthologs for DSCR1-like genes will allow knockout mice to be obtained for specific family members.

Summary: Extracting the desired data from a database entry for later analysis is a constant need in the biological sequence analysis community; GeneRecords 1.0 is a solution for GenBank biological flat file parsing, as it implements a structured representation of each feature and feature qualifier in GenBank following import in a common database managing system usable in a personal computer (Macintosh and Windows environments). This collection of related databases enables the local management of GenBank records, allowing indexing, retrieval and analysis of both information and sequences on a personal computer. Availability: the current release, including the FileMaker Pro runtime application (built for Windows and Macintosh environments), is freely available at http://apollo11.isto.unibo.it/software/

Housekeeping (HK) genes are constitutively expressed in order to maintain cellular function. They produce the minimal essential transcripts necessary for normal cellular physiology. Wide range expression, stable expression level and high expression level are independent features of a single gene expression and are all desirable for the definition of an “ideal” HK. Recent studies have questioned the possible existence of “ideal” HK mRNAs, mainly because of the wide expression conditions variability. This would imply that for each investigated organism the suitability of a putative HK should be verified. We perform a systematic analysis to identify “optimal” HK genes in Danio rerio (zebrafish), to be used in expression analyses conducted on embryos/larvae at different developmental stages, as well as on differentiated adult tissues from single donors. The expression pattern of candidate genes, selected on the basis of the literature available and of ad hoc bioinformatics analysis, was assessed by quantitative relative RT-PCR in an RNA panel, including six different embryo/larvae developmental stages and six adult tissues. Statistical analysis was performed to identify genes with the lowest expression standard deviation in the studied panel. Our results showed that beta-actin 2 (bactin2) is the mRNA with the lowest variability of expression.

Human RCAN3 (Regulator of calcineurin 3; previously known as DSCR1L2, Down syndrome critical region gene 1-like 2) is a five-exon gene mapped on chromosome 1 and belongs to the human RCAN gene family which also includes RCAN1 and RCAN2. The novel denomination RCAN for genes and proteins, instead of DSCR1L (Down syndrome critical region gene 1-like) has recently been widely discussed. The aim of the present work was to perform a multiple approach analysis of five RCAN3 mRNA and encoded protein isoforms, two of which have been identified for the first time in this research. The two new RCAN3 mRNA isoforms, RCAN3-2,4,5, which lacks exon 3, and RCAN3-2,3,5, which lacks exon 4, were identified during RCAN3 RT-PCR (reverse transcription-polymerase chain reaction) cloning, the product of which unexpectedly revealed the presence of five isoforms as opposed to the three previously known. In order to analyze the expression pattern of the five RCAN3 mRNA isoforms in seven different human tissues, a quantitative relative RT-PCR was performed: interestingly, all isoforms are present in all tissues investigated, with a statistically significant constant prevalence of RCAN3 isoform (the most complete, "reference" isoform). The RCAN3 locus expression level was comparable in all seven tissues analyzed, considering all isoforms, which indicates a ubiquitous expression of this human RCAN family member. To date two possible interactors have been described for this protein: human cardiac troponin I (TNNI3) and calcineurin. Here we report the interaction between the new RCAN3 variants and TNNI3, demonstrated by both yeast cotransformation and by the GST (glutathione-sepharose transferase) fusion protein assay, as was to be expected from the presence of exon 2 whose product has been seen to be sufficient for binding to TNNI3.

Background CYYR1 is a recently identified gene located on human chromosome 21 whose product has no similarity to any known protein and is of unknown function. Analysis of expressed sequence tags (ESTs) have revealed high human CYYR1 expression in cells belonging to the diffuse neuroendocrine system (DNES). These cells may be the origin of neuroendocrine (NE) tumors. The aim of this study was to conduct an initial analysis of sequence, splicing and expression of the CYYR1 mRNA in human NE tumors. Methods The CYYR1 mRNA coding sequence (CDS) was studied in 32 NE tumors by RT-PCR and sequence analysis. A subtle alternative splicing was identified generating two isoforms of CYYR1 mRNA differing in terms of the absence (CAG- isoform, the first described mRNA for CYYR1 locus) or the presence (CAG+ isoform) of a CAG codon. When present, this specific codon determines the presence of an alanine residue, at the exon 3/exon 4 junction of the CYYR1 mRNA. The two mRNA isoform amounts were determined by quantitative relative RT-PCR in 29 NE tumors, 2 non-neuroendocrine tumors and 10 normal tissues. A bioinformatic analysis was performed to search for the existence of the two CYYR1 isoforms in other species. Results The CYYR1 CDS did not show differences compared to the reference sequence in any of the samples, with the exception of an NE tumor arising in the neck region. Sequence analysis of this tumor identified a change in the CDS 333 position (T instead of C), leading to the amino acid mutation P111S. NE tumor samples showed no significant difference in either CYYR1 CAG- or CAG+ isoform expression compared to control tissues. CYYR1 CAG- isoform was significantly more expressed than CAG+ isoform in NE tumors as well as in control samples investigated. Bioinformatic analysis revealed that only the genomic sequence of Pan troglodytes CYYR1 is consistent with the possible existence of the two described mRNA isoforms. Conclusion A new "subtle" splicing isoform (CAG+) of CYYR1 mRNA, the sequence and the expression of this gene were defined in a large series of NE tumors.

Background All standard methods for cDNA cloning are affected by a potential inability to effectively clone the 5' region of mRNA. The aim of this work was to estimate mRNA open reading frame (ORF) 5' region sequence completeness in the model organism Danio rerio (zebrafish). Results We implemented a novel automated approach (5'_ORF_Extender) that systematically compares available expressed sequence tags (ESTs) with all the zebrafish experimentally determined mRNA sequences, identifies additional sequence stretches at 5' region and scans for the presence of all conditions needed to define a new, extended putative ORF. Our software was able to identify 285 (3.3%) mRNAs with putatively incomplete ORFs at 5' region and, in three example cases selected (selt1a, unc119.2, nppa), the extended coding region at 5' end was cloned by reverse transcription-polymerase chain reaction (RT-PCR). Conclusion The implemented method, which could also be useful for the analysis of other genomes, allowed us to describe the relevance of the "5' end mRNA artifact" problem for genomic annotation and functional genomic experiment design in zebrafish. Open peer review This article was reviewed by Alexey V. Kochetov (nominated by Mikhail Gelfand), Shamil Sunyaev, and Gáspár Jékely. For the full reviews, please go to the Reviewers' Comments section.

The CF lung facilitates the cohabitation of diverse microbial organisms. These communities diversify as a result of interactions with both the environment and each other. Cystic fibrosis-related diabetes (CFRD) is an increasingly recognised and complex issue in the management of CF patients. There has been little research into the effect that this change in the host has on the associated microbiome. Objectives The main objective of this study was to determine the bacterial community structure in patients with CFRD. Methods and Results Clinically significant dysglycemia was confirmed using continuous glucose monitoring showing interstitial glucose readings of 7.8 mmol/L for more than 4.5% of the time. Amplicon-sequencing using Illumina Miseq of the 16S rRNA gene was performed on DNA extracted from sputum samples taken during outpatient visits. Analysis was performed using uSearch/Qiime pipeline and R package. Sequencing of 19 diabetic and 22 non-diabetic patients revealed a diverse microbial population. In CFRD patients, a significantly higher level of S. maltophilia was detected and confirmed using qPCR. In vitro studies were performed using physiologically relevant levels of glucose (0–10 µM). S. maltophilia displays a significant increase in growth in glucose rich environments suggesting that this bacterium is adept at exploiting this niche. The importance of glucose control in CFRD cannot be underestimated as it has potential to predispose the growth of resistant organisms. Increased understanding of these complex interactions between the microbiome and host may ultimately lead to altered CF patient management. We acknowledge the Leverhulme Trust for funding.

The activation of the insulin-like growth factor 1/IGF1 receptor system (IGF1/IGF1R) is a critical event in the transformation and tumorigenicity processes in a wide variety of human tumors. The IGF1/IGF1R system has been recently studied in carcinoid tumors that often arise in the gastrointestinal tract; these tumors are characterized by hypersecretion of bioamines and neuropeptides, leading to functional tumor disease. Two alternatively spliced IGF1R mRNA transcripts have been described to differ by only three nucleotides (CAG) in the coding sequence, resulting in an amino-acid change from the originally described Thr-Gly to an Arg in the extracellular portion of the receptor beta subunit. In transfected Chinese hamster ovary cells, the form without CAG (CAG-) exhibited an approximate 2-fold increase in IGF1 stimulation of activities required for its mitogenic properties. In this study, we examine the relative expression of the two IGF1R mRNA isoforms by a semiquantitative RT-PCR approach using highly standardized conditions, beta-2 microglobulin (B2M) as a reference gene and gel imaging analysis. We analyzed a large series of human neuroendocrine tumors (32 samples) and 9 normal tissues. A significant higher expression of both isoforms in the tumor samples (approximately 2-fold increase) was found, while a constant CAG+/CAG- IGF1R mRNA isoforms of an approximate 3:1 ratio was observed in all tumoral and normal cell types studied. The phylogenetic study of the IGF1R locus in several species suggests that human IGF1R CAG- mRNA isoform is evolutionarily more recent compared to the IGF1R CAG+ mRNA isoform and it could be used by the splicing apparatus at this intron/exon junction with a lower efficiency. This study highlights the relevance of IGF1R mRNA expression in neuroendocrine tumor cells, and the constant presence of 'subtle' alternative splicing for the IGF1R locus.

Summary: UniGene Tabulator 1.0 provides a solution for full parsing of UniGene flat file format; it implements a structured graphical representation of each data field present in UniGene following import into a common database managing system usable in a personal computer. This database includes related tables for sequence, protein similarity, sequence-tagged site (STS) and transcript map interval (TXMAP) data, plus a summary table where each record represents a UniGene cluster. UniGene Tabulator enables full local management of UniGene data, allowing parsing, querying, indexing, retrieving, exporting and analysis of UniGene data in a relational database form, usable on Macintosh (OS X 10.3.9 or later) and Windows (2000, with service pack 4, XP, with service pack 2 or later) operating systems-based computers. Availability: The current release, including both the FileMaker runtime applications, is freely available at Contact: pierluigi.strippoli@unibo.it Supplementary information: We also distribute a precalculated implementation for current Homo sapiens (build #190, March 2006) and Danio rerio (zebrafish, build #90, March 2006) UniGene data.

Background Formal description of a cell's genetic information should provide the number of DNA molecules in that cell and their complete nucleotide sequences. We pose the formal problem: can the genome sequence forming the genotype of a given living cell be known with absolute certainty so that the cell's behaviour (phenotype) can be correlated to that genetic information? To answer this question, we propose a series of thought experiments. Results We show that the genome sequence of any actual living cell cannot physically be known with absolute certainty, independently of the method used. There is an associated uncertainty, in terms of base pairs, equal to or greater than μs (where μ is the mutation rate of the cell type and s is the cell's genome size). Conclusion This finding establishes an "uncertainty principle" in genetics for the first time, and its analogy with the Heisenberg uncertainty principle in physics is discussed. The genetic information that makes living cells work is thus better represented by a probabilistic model rather than as a completely defined object.

Research has widely supported the efficacy of screening for colorectal cancer in reducing mortality. A blood-based assay potentially represents a more accessible early detection tool for the identification of solid tumor cells originating from a primary tumor site in the body. We demonstrate a relatively easy and highly reproducible technique for the detection of mRNA expression of genes as markers of malignancy in blood samples of patients with colon cancer. The present study aims to identify a set of specific mRNAs expressed in epithelial cells but not in blood cells, which may be useful as markers for early detection of circulating colon cancer cells by a simple, qualitative RT-PCR assay following semi-automated RNA extraction from peripheral blood samples. Our approach includes a systematic search for candidate markers using digital differential display, search on UniGene colon EST libraries and analysis of published data on colon cancer gene expression. A final list included the following genes: bone morphogenetic protein 4 (BMP4), cyclin D (CycD), family with sequence similarity 3, member D (FAM3D), gastrin (GAS), glycoprotein A33 transmembrane (GPA33), glutathione peroxidase 2 gastrointestinal (GPX2), galactoside-binding, soluble, 4 (galectin 4) (LGALS4), non-SMC, structural maintenance of chromosomes, element 1 protein (NSE1), tumor-associated calcium signal transducer 1 (TACSTD1), telomerase reverse transcriptase (hTERT), trefoil factor 3 intestinal (TFF3), transmembrane 4 superfamily member 3 (TM4SF3), UDP glycosyltransferase 1 family, polypeptide A9 (UGT1A9), villin 1 (VIL1), and the novel gene FLJ20127. The mRNA expression of these genes was evaluated in a pool of 16 samples from subjects diagnosed with colon cancer and from 16 normal-controls. We observed expression in 13 of the 15 investigated genes from the blood samples of the vast majority of patients considered, but also in a certain percentage of the controls (from 14.3 to 100%). This finding confirms that the extreme sensitivity of RT-PCR is able to detect minimal amounts of mRNA expressed in a non tissue-specific manner ('illegitimate transcription'). On the contrary, NSE1 and GAS mRNAs were not detected either in patient or in control blood samples; however, they were abundantly expressed in normal and cancerous colon mucosa, encouraging further search for useful markers able to detect epithelial cells in peripheral blood.

The amino acid sequence of gene products is routinely deduced from the nucleotide sequence of the relative cloned cDNA, according to the rules for recognition of start codon (first-AUG rule, optimal sequence context) and the genetic code. From this prediction stem most subsequent types of product analysis, although all standard methods for cDNA cloning are affected by a potential inability to effectively clone the 5′ region of mRNA. Revision by bioinformatics and cloning methods of 109 known genes located on human chromosome 21 (HC 21) shows that 60 mRNAs lack any in-frame stop upstream of the first-AUG, and that in five cases (DSCR1, KIAA0184, KIAA0539, SON, and TFF3) the coding region at the 5′ end was incompletely characterized in the original descriptions. We describe the respective consequences for genomic annotation, domain and ortholog identification, and functional experiments design. We have also analyzed the sequences of 13,124 human mRNAs (RefSeq databank), discovering that in 6448 cases (49%), an in-frame stop codon is present upstream of the initiation codon, while in the other 6676 mRNAs (51%), identification of additional bases at the mRNA 5′ region could well reveal some new upstream in-frame AUG codons in the optimal context. Proportionally to the HC 21 data, about 550 known human genes might thus be affected by this 5′ end mRNA artifact. © 2003 Elsevier B.V. All rights reserved.

Background : The ADARB1 gene is located in 21q22.3 region, previously linked to familial bipolar disorder, and its product has a documented action in the editing of the pre-mRNA of glutamate receptor B subunit. Dysfunction of glutamatergic neurotransmission could play an important role in the patophysiology of bipolar disorder (BD). Glutamate excitatory neurotransmission regulation is a possible mechanism of the initial effect of anticonvulsants in regulating mood. Methods : To investigate the hypothesis of an involvement of ADARB1 gene in the BD, the ADARB1 cDNA has been cloned and sequenced in seven selected bipolar I disorder patients with evidence of familiarity of mood disorders. A detailed investigation of the gene nucleotide sequence in the open reading frame has been performed. Results : No alteration in the sequence of the ADARB1 gene cDNA was found in any patient, except a common neutral polymorphism in three out of seven patients. Conclusions : Mutations in ADARB1 gene are not commonly associated with bipolar I disorder, therefore other genes in the 21q22 region could be associated with bipolar illness in some families, likely in the context of a multifactorial transmission model.

A new gene family has been identified on the basis of in-depth bioinformatics analysis of the Down syndrome candidate region 1 ( DSCR1 ) gene, located on 21q22.1. We have determined the complete coding sequences of similar genes in Saccharomyces cerevisiae and Caenorhabditis elegans, as well as that of a novel human gene, named DSCR1L2 ( DSCR1 -like 2). Peripheral blood leukocyte cDNA sequencing predicts as its product a 241-amino-acid protein highly similar to products of the human genes DSCR1 and ZAKI-4 (HGMW-approved symbol DSCR1L1 ). The highest level of expression of DSCR1L2 mRNA was found by Northern blot analysis in heart and skeletal muscles, liver, kidney, and peripheral blood leukocytes (three transcripts of 3.2, 5.2, and 7.5 kb). The gene consists of four exons and spans about 22 kb on chromosome 1 (1p33–p35.3) (Human Chromosome 1, Sanger Centre). Exon/intron organization is highly conserved between DSCR1 and DSCR1L2. Two alternative DSCR1L2 mRNA splicing forms have been recognized, with one lacking 10 amino acids in the middle of the protein. Analysis of expressed sequence tags (ESTs) shows DSCR1L2 expression in fetal tissues (heart, liver, and spleen) and in adenocarcinomas. ESTs related to the murine DSCR1L2 orthologue are found in the 2-cell stage mouse embryo, in developing brain stem and spinal cord, and in thymus and T cells. The most prominent feature identified in the protein family is a central short, unique serine–proline motif (including an ISPPXSPP box), which is strongly conserved from yeast to human but is absent in bacteria. Moreover, homology with the RNA-binding domain was weakly but consistently detected in a stretch of 80 amino acids at the amino-terminus by fine sequence analysis based on tools utilizing both hidden Markov models and BLAST. The identification of this new gene family should allow a better understanding of the functions of the genes belonging to it.

Background Several tools have been developed to perform global gene expression profile data analysis, to search for specific chromosomal regions whose features meet defined criteria as well as to study neighbouring gene expression. However, most of these tools are tailored for a specific use in a particular context (e.g. they are species-specific, or limited to a particular data format) and they typically accept only gene lists as input. Results TRAM (Transcriptome Mapper) is a new general tool that allows the simple generation and analysis of quantitative transcriptome maps, starting from any source listing gene expression values for a given gene set (e.g. expression microarrays), implemented as a relational database. It includes a parser able to assign univocal and updated gene symbols to gene identifiers from different data sources. Moreover, TRAM is able to perform intra-sample and inter-sample data normalization, including an original variant of quantile normalization (scaled quantile), useful to normalize data from platforms with highly different numbers of investigated genes. When in 'Map' mode, the software generates a quantitative representation of the transcriptome of a sample (or of a pool of samples) and identifies if segments of defined lengths are over/under-expressed compared to the desired threshold. When in 'Cluster' mode, the software searches for a set of over/under-expressed consecutive genes. Statistical significance for all results is calculated with respect to genes localized on the same chromosome or to all genome genes. Transcriptome maps, showing differential expression between two sample groups, relative to two different biological conditions, may be easily generated. We present the results of a biological model test, based on a meta-analysis comparison between a sample pool of human CD34+ hematopoietic progenitor cells and a sample pool of megakaryocytic cells. Biologically relevant chromosomal segments and gene clusters with differential expression during the differentiation toward megakaryocyte were identified. Conclusions TRAM is designed to create, and statistically analyze, quantitative transcriptome maps, based on gene expression data from multiple sources. The release includes FileMaker Pro database management runtime application and it is freely available at http://apollo11.isto.unibo.it/software/, along with preconfigured implementations for mapping of human, mouse and zebrafish transcriptomes.

A novel human gene denominated TruB pseudouridine (psi) synthase homolog 1 (E. coli) (approved symbol, TRUB1) has been identified and characterized. Spanning approximately 40 kb on chromosome 10 and including 8 exons, TRUB1 is the first described human ortholog of bacterial TruB/psi55, a gene involved in tRNA pseudouridinilation. TRUB1 gene encodes a 349-amino acid product, with a VFAVHKPKGPTSA box in positions 71-83 corresponding to motif I of the TruB family (probably involved in conserving protein structure). The TruB domain of TRUB1 lies between W104 and I255, and contains another short motif, GGTLDS AARGVLVV, including the highly conserved D residue that characterizes motif II (involved in uridine recognition and in catalytic function of psi synthases). Northern blot analysis revealed that TRUB1 mRNA is widely expressed in various human tissues (especially heart, skeletal muscle and liver). Phylogenetic analysis of the TruB domain revealed another human gene (approved symbol TRUB2) encoding a conserved TruB domain, located on human chromosome 9. Thus, the human TruB family includes at least three members: i.e. DKC1 (previously identified), TRUB1 and TRUB2. The TRUB1 and TRUB2 products could be the hitherto unidentified human tRNA psi synthases. Although TRUB1 is not highly similar to DKC1/dyskerin (whose mutations cause X-linked dyskeratosis congenita) and putatively affects tRNA rather than rRNA modification, it is the most similar human protein to dyskerin. Study of TRUB1 (and TRUB2) should facilitate understanding of the molecular mechanisms of RNA modification and the involvement of psi synthases in human pathology, including dyskeratosis-like diseases.

Members of the bone morphogenetic protein (BMP) family and their receptors (BMPRs and activin receptors-ActRs) promote the development of bones with a fine regulation of their expression. Mutations in BMPs or BMPRs cause several diseases, as shown in knockout mice, such as skeletal defects, familial primary pulmonary hypertension and neoplasias. Osteosarcoma is the most frequent primary malignant tumor of bone. Due to their importance in bone development, BMPs, BMPRs and ActRs could also play a role in osteosarcoma growth and development. Previous data have shown that the overexpression of the BMPR-II was related to poor prognosis in malignant and metastatic bone tumors. We evaluate by reverse transcription-linked polymerase chain reaction analysis (RT-PCR) the expression pattern of BMPs, BMPRs and ActRs in five different human osteosarcoma cell lines (MG63, G292, HOS, SaOS and U2). Moreover, we performed the mutational screening of the complete BMPR-II mRNA by automated sequencing of the correspondent cDNA to evaluate the presence of point mutations in osteosarcoma cell lines. All the osteosarcoma cell lines studied simultaneously expressed the BMPs, BMPRs and ActRs investigated. No mutations were detected in the BMPR-II cDNA. Our results suggest the presence of a mechanism involving the simultaneous activation of the BMPs and their receptors in osteosarcoma cell lines.