cDNA Library Generation for the Analysis of Small RNAs by High-Throughput Sequencing

The RNome of a cell is highly diverse and consists besides messenger RNAs (mRNAs), transfer RNAs (tRNAs), and ribosomal RNAs (rRNAs) also of other small and long transcript entities without apparent coding potential. This class of molecules, commonly referred to as non-protein-coding RNAs (ncRNAs), is involved in regulating numerous biological processes and thought to contribute to cellular complexity. Therefore, much effort is put into their identification and further functional characterization. Here we provide a cost-effective and reliable method for cDNA library construction of small RNAs in the size range of 20-500 residues. The effectiveness of the described method is demonstrated by the analysis of ribosome-associated small RNAs in the eukaryotic model organism Trypanosoma brucei.

Whole-Genome Sequencing of a Canine Family Trio Reveals a FAM83G Variant Associated with Hereditary Footpad Hyperkeratosis

Over 250 Mendelian traits and disorders, caused by rare alleles have been mapped in the canine genome. Although each disease is rare in the dog as a species, they are collectively common and have major impact on canine health. With SNP-based genotyping arrays, genome-wide association studies (GWAS) have proven to be a powerful method to map the genomic region of interest when 10-20 cases and 10-20 controls are available. However, to identify the genetic variant in associated regions, fine-mapping and targeted re-sequencing is required. Here we present a new approach using whole-genome sequencing (WGS) of a family trio without prior GWAS. As a proof-of-concept, we chose an autosomal recessive disease known as hereditary footpad hyperkeratosis (HFH) in Kromfohrl änder dogs. To our knowledge, this is the first time this family trio WGS-approach, has successfully been used to identify a genetic variant that perfectly segregates with a canine disorder. The sequencing of three Kromfohrl änder dogs from a family trio (an affected offspring and both its healthy parents) resulted in an average genome coverage of 9.2X per individual. After applying stringent filtering criteria for candidate causative coding variants, 527 single nucleotide variants (SNVs) and 15 indels were found to be homozygous in the affected offspring and heterozygous in the parents. Using the computer software packages ANNOVAR and SIFT to functionally annotate coding sequence differences and to predict their functional effect, resulted in seven candidate variants located in six different genes. Of these, only FAM83G:c155G>C (p.R52P) was found to be concordant in eight additional cases and 16 healthy Kromfohrl änder dogs.

Standard genotyping overestimates transmission of Mycobacterium tuberculosis among immigrants in a low incidence country.

Immigrants from high tuberculosis (TB) incidence regions are a risk group for TB in low-incidence countries such as Switzerland. In a previous analysis of a nationwide collection of 520 Mycobacterium tuberculosis isolates from 2000-2008, we identified 35 clusters comprising 90 patients based on standard genotyping (24-loci MIRU-VNTR and spoligotyping). Here, we used whole genome sequencing (WGS) to revisit these transmission clusters. Genome-based transmission clusters were defined as isolate pairs separated by ≤12 single nucleotide polymorphisms (SNPs). WGS confirmed 17/35 (49%) MIRU-VNTR clusters; the other 18 clusters contained pairs separated by >12 SNPs. Most transmission clusters (3/4) of Swiss-born patients were confirmed by WGS, as opposed to 25% (4/16) of clusters involving only foreign-born patients. The overall clustering proportion using standard genotyping was 17% (90 patients, 95% confidence interval [CI]: 14-21%), but only 8% (43 patients, 95% CI: 6-11%) using WGS. The clustering proportion was 17% (67/401, 95% CI: 13-21%) using standard genotyping and 7% (26/401, 95% CI: 4-9%) using WGS among foreign-born patients, and 19% (23/119, 95% CI: 13-28%) and 14% (17/119, 95% CI: 9-22%), respectively, among Swiss-born patients. Using weighted logistic regression, we found weak evidence for an association between birth origin and transmission (aOR 2.2, 95% CI: 0.9-5.5, comparing Swiss-born patients to others). In conclusion, standard genotyping overestimated recent TB transmission in Switzerland when compared to WGS, particularly among immigrants from high TB incidence regions, where genetically closely related strains often predominate. We recommend the use of WGS to identify transmission clusters in low TB incidence settings.

Specific processing of tenascin-C by the metalloprotease meprinbeta neutralizes its inhibition of cell spreading

The metalloprotease meprin has been implicated in tissue remodelling due to its capability to degrade extracellular matrix components. Here, we investigated the susceptibility of tenascin-C to cleavage by meprinbeta and the functional properties of its proteolytic fragments. A set of monoclonal antibodies against chicken and human tenascin-C allowed the mapping of proteolytic fragments generated by meprinbeta. In chicken tenascin-C, meprinbeta processed all three major splicing variants by removal of 10kDa N-terminal and 38kDa C-terminal peptides, leaving a large central part of subunits intact. A similar cleavage pattern was found for large human tenascin-C variant where two N-terminal peptides (10 or 15kDa) and two C-terminal fragments (40 and 55kDa) were removed from the intact subunit. N-terminal sequencing revealed the exact amino acid positions of cleavage sites. In both chicken and human tenascin-C N-terminal cleavages occurred just before and/or after the heptad repeats involved in subunit oligomerization. In the human protein, an additional cleavage site was identified in the alternative fibronectin type III repeat D. Whereas all these sites are known to be attacked by several other proteases, a unique cleavage by meprinbeta was located to the 7th constant fibronectin type III repeat in both chicken and human tenascin-C, thereby removing the C-terminal domain involved in its anti-adhesive activity. In cell adhesion assays meprinbeta-digested human tenascin-C was not able to interfere with fibronectin-mediated cell spreading, confirming cleavage in the anti-adhesive domain. Whereas the expression of meprinbeta and tenascin-C does not overlap in normal colon tissue, inflamed lesions of the mucosa from patients with Crohn's disease exhibited many meprinbeta-positive leukocytes in regions where tenascin-C was strongly induced. Our data indicate that, at least under pathological conditions, meprinbeta might attack specific functional sites in tenascin-C that are important for its oligomerization and anti-adhesive activity.

Molecular epidemiology of the nasal colonization by methicillin-susceptible Staphylococcus aureus in Swiss children

Nasal carriage of Staphylococcus aureus contributes to an increased risk of developing an infection with the same bacterial strain. Genetic regulatory elements and toxin-expressing genes are virulence factors associated with the pathogenic potential of S. aureus. We undertook an extensive molecular characterization of methicillin-susceptible S. aureus (MSSA) carried by children. MSSA were recovered from the nostrils of children. The presence of Panton-Valentine leukocidin (PVL), exfoliatins A and B (exfoA and exfoB), and the toxic-shock staphylococcal toxin (TSST-1) and agr group typing were determined by quantitative PCR. A multiple-locus variable-number of tandem repeat analysis (MLVA) assay was also performed for genotyping. Five hundred and seventy-two strains of MSSA were analysed. Overall, 30% were positive for toxin-expressing genes: 29% contained one toxin and 1.6% two toxins. The most commonly detected toxin gene was tst, which was present in 145 (25%) strains. The TSST-1 gene was significantly associated with the agr group 3 (OR 56.8, 95% CI 32.0-100.8). MLVA analysis revealed a large diversity of genetic content and no clonal relationship was demonstrated among the analysed MSSA strains. Multilocus sequence typing confirmed this observation of diversity and identified ST45 as a frequent colonizer. This broad diversity in MSSA carriage strains suggests a limited selection pressure in our geographical area.

Deep-sequencing identification of the genomic targets of the cytidine deaminase AID and its cofactor RPA in B lymphocytes

The cytidine deaminase AID hypermutates immunoglobulin genes but can also target oncogenes, leading to tumorigenesis. The extent of AID's promiscuity and its predilection for immunoglobulin genes are unknown. We report here that AID interacted broadly with promoter-proximal sequences associated with stalled polymerases and chromatin-activating marks. In contrast, genomic occupancy of replication protein A (RPA), an AID cofactor, was restricted to immunoglobulin genes. The recruitment of RPA to the immunoglobulin loci was facilitated by phosphorylation of AID at Ser38 and Thr140. We propose that stalled polymerases recruit AID, thereby resulting in low frequencies of hypermutation across the B cell genome. Efficient hypermutation and switch recombination required AID phosphorylation and correlated with recruitment of RPA. Our findings provide a rationale for the oncogenic role of AID in B cell malignancy.

PTCH promoter methylation at low level in sporadic basal cell carcinoma analysed by three different approaches

Please cite this paper as: PTCH promoter methylation at low level in sporadic basal cell carcinoma analysed by three different approaches. Experimental Dermatology 2010. Abstract: Basal cell carcinoma (BCC) is the most common form of skin cancer. Mutations of the PTCH hallmark gene are detected in about 50-60% of BCCs, which raises the question whether other mechanisms such as promoter methylation can inactivate PTCH. Therefore, we performed methylation analysis of the PTCH promoter in a total of 56 BCCs. The sensitivity of three different methods, including direct bisulphite sequencing PCR, MethyLight and high-resolution melting (HRM), was applied and compared. We found that HRM analysis of DNA from fresh tissue [rather than formalin-fixed and paraffin-embedded tissue (FFPE)] was the most sensitive method to detect methylation. Low-level methylation of the PTCH promoter was detected in five out of 16 analysed BCCs (31%) on DNA from fresh tissue but only in two (13%) samples on short-time stored FFPE DNA from the very same tumors. In contrast, we were unable to detect methylation by HRM on long-time stored DNA in any of the remaining 40 BCC samples. Our data suggest that (i) HRM on DNA extracted from fresh tissue is the most sensitive method to detect methylation and (ii) methylation of the PTCH promoter may only play a minor role in BCC carcinogenesis.

Multilocus sequence analysis and rpoB sequencing of Mycobacterium abscessus (sensu lato) strains

Mycobacterium abscessus, Mycobacterium bolletii, and Mycobacterium massiliense (Mycobacterium abscessus sensu lato) are closely related species that currently are identified by the sequencing of the rpoB gene. However, recent studies show that rpoB sequencing alone is insufficient to discriminate between these species, and some authors have questioned their current taxonomic classification. We studied here a large collection of M. abscessus (sensu lato) strains by partial rpoB sequencing (752 bp) and multilocus sequence analysis (MLSA). The final MLSA scheme developed was based on the partial sequences of eight housekeeping genes: argH, cya, glpK, gnd, murC, pgm, pta, and purH. The strains studied included the three type strains (M. abscessus CIP 104536(T), M. massiliense CIP 108297(T), and M. bolletii CIP 108541(T)) and 120 isolates recovered between 1997 and 2007 in France, Germany, Switzerland, and Brazil. The rpoB phylogenetic tree confirmed the existence of three main clusters, each comprising the type strain of one species. However, divergence values between the M. massiliense and M. bolletii clusters all were below 3% and between the M. abscessus and M. massiliense clusters were from 2.66 to 3.59%. The tree produced using the concatenated MLSA gene sequences (4,071 bp) also showed three main clusters, each comprising the type strain of one species. The M. abscessus cluster had a bootstrap value of 100% and was mostly compact. Bootstrap values for the M. massiliense and M. bolletii branches were much lower (71 and 61%, respectively), with the M. massiliense cluster having a fuzzy aspect. Mean (range) divergence values were 2.17% (1.13 to 2.58%) between the M. abscessus and M. massiliense clusters, 2.37% (1.5 to 2.85%) between the M. abscessus and M. bolletii clusters, and 2.28% (0.86 to 2.68%) between the M. massiliense and M. bolletii clusters. Adding the rpoB sequence to the MLSA-concatenated sequence (total sequence, 4,823 bp) had little effect on the clustering of strains. We found 10/120 (8.3%) isolates for which the concatenated MLSA gene sequence and rpoB sequence were discordant (e.g., M. massiliense MLSA sequence and M. abscessus rpoB sequence), suggesting the intergroup lateral transfers of rpoB. In conclusion, our study strongly supports the recent proposal that M. abscessus, M. massiliense, and M. bolletii should constitute a single species. Our findings also indicate that there has been a horizontal transfer of rpoB sequences between these subgroups, precluding the use of rpoB sequencing alone for the accurate identification of the two proposed M. abscessus subspecies.

The venom composition of the parasitic wasp Chelonus inanitus resolved by combined expressed sequence tags analysis and proteomic approach

Background Parasitic wasps constitute one of the largest group of venomous animals. Although some physiological effects of their venoms are well documented, relatively little is known at the molecular level on the protein composition of these secretions. To identify the majority of the venom proteins of the endoparasitoid wasp Chelonus inanitus (Hymenoptera: Braconidae), we have randomly sequenced 2111 expressed sequence tags (ESTs) from a cDNA library of venom gland. In parallel, proteins from pure venom were separated by gel electrophoresis and individually submitted to a nano-LC-MS/MS analysis allowing comparison of peptides and ESTs sequences. Results About 60% of sequenced ESTs encoded proteins whose presence in venom was attested by mass spectrometry. Most of the remaining ESTs corresponded to gene products likely involved in the transcriptional and translational machinery of venom gland cells. In addition, a small number of transcripts were found to encode proteins that share sequence similarity with well-known venom constituents of social hymenopteran species, such as hyaluronidase-like proteins and an Allergen-5 protein. An overall number of 29 venom proteins could be identified through the combination of ESTs sequencing and proteomic analyses. The most highly redundant set of ESTs encoded a protein that shared sequence similarity with a venom protein of unknown function potentially specific of the Chelonus lineage. Venom components specific to C. inanitus included a C-type lectin domain containing protein, a chemosensory protein-like protein, a protein related to yellow-e3 and ten new proteins which shared no significant sequence similarity with known sequences. In addition, several venom proteins potentially able to interact with chitin were also identified including a chitinase, an imaginal disc growth factor-like protein and two putative mucin-like peritrophins. Conclusions The use of the combined approaches has allowed to discriminate between cellular and truly venom proteins. The venom of C. inanitus appears as a mixture of conserved venom components and of potentially lineage-specific proteins. These new molecular data enrich our knowledge on parasitoid venoms and more generally, might contribute to a better understanding of the evolution and functional diversity of venom proteins within Hymenoptera.

Epigenetic regulation of somatic angiotensin-converting enzyme by DNA methylation and histone acetylation

Somatic angiotensin-converting enzyme (sACE) is crucial in cardiovascular homeostasis and displays a tissue-specific profile. Epigenetic patterns modulate genes expression and their alterations were implied in pathologies including hypertension. However, the influence of DNA methylation and chromatin condensation state on the expression of sACE is unknown. We examined whether such epigenetic mechanisms could participate in the control of sACE expression in vitro and in vivo. We identified two CpG islands in the human ace-1 gene 3 kb proximal promoter region. Their methylation abolished the luciferase activity of ace-1 promoter/reporter constructs transfected into human liver (HepG2), colon (HT29), microvascular endothelial (HMEC-1) and lung (SUT) cell lines (p < 0.001). Bisulphite sequencing revealed a cell-type specific basal methylation pattern of the ace-1 gene -1,466/+25 region. As assessed by RT-qPCR, inhibition of DNA methylation by 5-aza-2'-deoxycytidine and/or of histone deacetylation by trichostatin A highly stimulated sACE mRNA expression cell-type specifically (p < 0.001 vs. vehicle treated cells). In the rat, in vivo 5-aza-cytidine injections demethylated the ace-1 promoter and increased sACE mRNA expression in the lungs and liver (p = 0.05), but not in the kidney. In conclusion, the expression level of somatic ACE is modulated by CpG-methylation and histone deacetylases inhibition. The basal methylation pattern of the promoter of the ace-1 gene is cell-type specific and correlates to sACE transcription. DNMT inhibition is associated with altered methylation of the ace-1 promoter and a cell-type and tissue-specific increase of sACE mRNA levels. This study indicates a strong influence of epigenetic mechanisms on sACE expression.

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.

Meningitis caused by Mycoplasma mycoides subspecies capri in a goat

A 2-year-old, female goat from Connecticut was submitted for necropsy with a 5-day history of pyrexia and intermittent neurologic signs, including nystagmus, seizures, and circling. Postmortem examination revealed suppurative meningitis. Histologic examination of the brain revealed that the meninges were diffusely infiltrated by moderate numbers of lymphocytes, macrophages, and fibrin, with scattered foci of dense neutrophilic infiltrate. Culture of pus and brainstem yielded typical mycoplasma colonies. DNA sequencing of the 16S ribosomal RNA gene revealed 99% sequence homology with Mycoplasma mycoides subspecies capri and Mycoplasma mycoides subspecies mycoides Large Colony biotype, which are genetically indistinguishable and likely to be combined as a single subspecies labeled M. mycoides subsp. capri. The present case is unusual in that not only are mycoplasma an uncommon cause of meningitis in animals, but additionally, in that all other reported cases of mycoplasma meningitis in goats, systemic lesions were also present. In the present case, meningitis was the only lesion, thus illustrating the need to consider mycoplasma as a differential diagnosis for meningitis in goats.

Genotyping of Mycoplasma hyopneumoniae in wild boar lung samples

Mycoplasma hyopneumoniae is the etiologic agent of enzootic pneumonia (EP), an important cause of disease-associated losses in swine production and a role of wild boar in recurrent infections can be supposed. Genotypes of M. hyopneumoniae from wild boar are unknown but could indicate its role as a potential reservoir. Therefore, 34 lung samples being PCR-positive for M. hyopneumoniae from wild boar from the Geneva region in Switzerland were assayed by genotyping using the p146 and multi-locus sequence typing (MLST) approaches and compared to data from outbreak cases from domestic swine in Switzerland. Successful genotyping was dependent on a sufficiently high concentration of M. hyopneumoniae DNA in the samples as assessed by different real-time PCR assays. The p146 genotyping was more successful with 24 samples (70.5%) being typeable whereas only 6 samples (17.6%) could be genotyped using the MLST approach. Variability of genotypes was high but identical types were found in geographically related animals. Genotypes from wild boar showed phylogenetic relatedness to those from domestic pigs but no matching types could be identified. Results show that direct genotyping from wild boar lung samples is possible and provides a promising approach to investigate future EP outbreak related samples from wild boar. (C) 2011 Elsevier B.V. All rights reserved.

A new DPYD genotyping assay for improving the safety of 5-fluorouracil therapy

Chemotherapeutic use of 5-fluorouracil (5FU) is compromised by 10-20% of patients developing severe toxicity. Recently described genetic variation in dihydropyrimidine dehydrogenase (DPYD) has been shown to be a major predictor of 5FU toxicity. Here, we describe a new genotyping assay for routine clinical use that covers all the major DPYD risk variants.

Bioinformatics tools and databases for the study of human growth hormone

Advances in novel molecular biological diagnostic methods are changing the way of diagnosis and study of metabolic disorders like growth hormone deficiency. Faster sequencing and genotyping methods require strong bioinformatics tools to make sense of the vast amount of data generated by modern laboratories. Advances in genome sequencing and computational power to analyze the whole genome sequences will guide the diagnostics of future. In this chapter, an overview of some basic bioinformatics resources that are needed to study metabolic disorders are reviewed and some examples of bioinformatics analysis of human growth hormone gene, protein and structure are provided.