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.

Bone morphogenetic protein receptor 2 in patients with idiopathic portal hypertension

In idiopathic portal hypertension (IPH) typical vascular lesions are present in the branches of the portal vein or in the perisinusoidal area of the liver. Similar histological alterations have been reported in the pulmonary vasculature of patients with idiopathic pulmonary artery hypertension (IPAH). As IPAH is associated with mutations of the bone morphogenetic protein receptor 2 (BMPR2) gene, the aim of this study was to investigate whether this association might also be found in patients with IPH. Twenty-three samples belonging to 21 unrelated caucasian patients with IPH followed in the hepatic haemodynamic laboratory of the Hospital Clinic in Barcelona were included in the study. All patients were studied for the entire open reading frame and splice site of the BMPR2 gene by direct sequencing and multiple ligation probe amplification (MLPA) in order to detect large deletions/duplications. None of the 23 patients had pulmonary artery hypertension. Four patients presented one single nucleotide polymorphism (SNP) in intron 5, four patients had a SNP in exon 12 and a SNP in exon 1 was found in two cases. Two patients had both intron 5 and exon 12 polymorphisms. All SNPs were previously described. Except for these three SNPs, neither mutations nor rearrangements have been identified in the BMPR2 gene in this population. We did not detect mutations or rearrangements in the coding region of the BMPR2 gene in our patients with IPH. These findings suggest that, in contrast to IPAH, mutations in BMPR2 are not involved in the pathogenesis of IPH.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Detection of prion protein gene polymorphisms in Baltic breeds of sheep

A total of 167 sheep belonging to the Estonian whiteheaded mutton, Estonian blackheaded mutton, Lithuanian coarsewool native, Lithuanian blackface and Latvian darkheaded mutton breeds, and a population of sheep kept isolated on the Estonian island of Ruhnu, were sequence-analysed for polymorphisms in the prion protein (PrP) gene, to determine their genotype and the allele frequencies of polymorphisms in PrP known to confer resistance to scrapie. A 939 base pair fragment of exon 3 from the PrP gene was amplified by pcr and analysed by direct sequencing. For animals showing polymorphism at two nucleotide positions, both haplotypes of these double-heterozygous genotypes were further verified by pcr cloning and sequence analysis. Known polymorphisms were observed at codons 136, 154 and 171, and six different haplotypes (arr, ahq, arh, ahr, arq and vrq) were determined. On the basis of these polymorphisms, the six populations of sheep possessed the resistant arr haplotype at different frequencies. The high-risk arq haplotype occurred in high frequencies in all six populations, but vrq, the haplotype carrying the highest risk, occurred at low frequencies and in only three of the populations.

Stejnulxin, a novel snake C-type lectin-like protein from Trimeresurus stejnegeri venom is a potent platelet agonist acting specifically via GPVI

Stejnulxin, a novel snake C-type lectin-like protein with potent platelet activating activity, was purified and characterized from Trimeresurus stejnegeri venom. Under non-reducing conditions, it migrated on a SDS-polyacrylamide gel with an apparent molecular mass of 120 kDa. On reduction, it separated into three polypeptide subunits with apparent molecular masses of 16 kDa (alpha), 20 kDa (beta1) and 22 kDa (beta2), respectively. The complete amino acid sequences of its subunits were deduced from cloned cDNAs. The N-terminal sequencing and cDNA cloning indicated that beta1 and beta2 subunits of stejnulxin have identical amino acid sequences and each contains two N-glycosylation sites. Accordingly, the molecular mass difference between beta1 and beta2 is caused by glycosylation heterogenity. The subunit amino acid sequences of stejnulxin are similar to those of convulxin, with sequence identities of 52.6% and 66.4% for the alpha and beta, respectively. Stejnulxin induced human platelet aggregation in a dose-dependent manner. Antibodies against alphaIIbbeta3 inhibited the aggregation response to stejnulxin, indicating that activation of alphaIIbbeta3 and binding of fibrinogen are involved in stejnulxin-induced platelet aggregation. Antibodies against GPIbalpha or alpha2beta1 as well as echicetin or rhodocetin had no significant effect on stejnulxin-induced platelet aggregation. However, platelet activation induced by stejnulxin was blocked by anti-GPVI antibodies. In addition, stejnulxin induced a tyrosine phosphorylation profile in platelets that resembled that produced by convulxin. Biotinylated stejnulxin bound specifically to platelet membrane GPVI.

The multifunctional FUS protein: Characterization of the biological effects of its depletion

FUS/TLS (fused in sarcoma/translocated in liposarcoma) protein, a ubiquitously expressed RNA-binding protein, has been linked to a variety of cellular processes, such as RNA metabolism, microRNA biogenesis and DNA repair. However, the precise role of FUS protein remains unclear. Recently, FUS has been linked to Amyotrophic Lateral Sclerosis (ALS), a neurodegenerative disorder characterized by the dysfunction and death of motor neurons. Based on the observation that some mutations in the FUS gene induce cytoplasmic accumulation of FUS aggregates, we decided to explore a loss-of-function situation (i.e. inhibition of FUS’ nuclear function) to unravel the role of this protein. To this purpose, we have generated a SH-SY5Y human neuroblastoma cell line which expresses a doxycycline induced shRNA targeting FUS and that specifically depletes the protein. In order to characterize this cell line, we have performed a whole transcriptome analysis by RNA deep sequencing. Preliminary results show that FUS depletion affects both expression and alternative splicing levels of several RNAs. When FUS is depleted we observed 330 downregulated and 81 upregulated genes. We also found that 395 splicing isoforms were downregulated, while 426 were upregulated. Currently, we are focusing our attention on the pathways which are mostly affected by FUS depletion. In addition, to further characterize the FUS-depleted cell line we have performed growth proliferation and survival assays. From these experiments emerge that FUS-depleted cells display growth proliferation alteration. In order to explain this observation, we have tested different hypothesis (e.g. apoptosis, senescence or slow-down growth). We observed that FUS-depleted cells growth slower than controls. Currently, we are looking for putative candidate targets causing this phenotype. Finally, since MEFs and B-lymphocytes derived from FUS knockdown mice display major sensitivity to ionizing radiation and chromosomal aberrations [1,2], we are exploring the effects of DNA damage in FUS-depleted cells by monitoring important components of DNA Damage Response (DDR). Taken together, these studies may contribute to our knowledge of the role of FUS in these cellular processes and will allow us to draw a clearer picture of mechanisms of neurodegenerative diseases.

The FUS protein is required for cell proliferation

FUS/TLS (fused in sarcoma/translocated in liposarcoma) protein, a ubiquitously expressed and highly conserved RNA binding protein, has been linked to a variety of cellular processes from mRNA processing to DNA repair. However, the precise function of FUS is not well understood. Recently, mutations in the FUS gene have been identified in familial and sporadic patients of Amyotrophic Lateral Sclerosis, a fatal neurodegenerative disorder characterized by dysfunction and death of motor neurons. Based on the observation that some mutations in the FUS gene induce cytoplasmic accumulation of FUS aggregates, we decided to explore a loss-of-function situation (i.e. inhibition of FUS’ nuclear function) to unravel the role of this protein. To this purpose, we have generated a SH-SY5Y human neuroblastoma cell line which expresses a doxycycline induced shRNA targeting FUS that efficiently depletes the protein. In order to characterize this cell line, we have characterized the poly(A) fraction by RNA deep sequencing. Preliminary results show that FUS depletion affects both mRNA expression and alternative splicing. Upon FUS depletion 330 genes are downregulated and 81 are upregulated. We also found that 395 splicing isoforms were downregulated, while 426 were upregulated. Currently, we are focusing our attention on the pathways which are mostly affected by FUS depletion. In addition, we are currently characterizing how FUS depletion affects cell proliferation and survival. We find that the lack of FUS impairs cell proliferation but does not induce apoptosis. Finally, since MEFs and B-lymphocytes derived from FUS knockdown mice display major sensitivity to ionizing radiation and chromosomal aberrations [1,2], we are exploring the effects of DNA damage in FUS-depleted cells by monitoring important components of DNA Damage Response (DDR). Taken together, these studies may contribute to our knowledge of the role of FUS in these cellular processes and will allow us to draw a clearer picture of mechanisms of neurodegenerative diseases.

Bos indicus introgression into (peri-)alpine cattle breeds - evidence from the analysis of bovine whey protein variants.

The major bovine whey proteins, α-lactalbumin (α-LA) and β-lactoglobulin (β-LG), exhibit breed-specific genetic variation. The aim of this study was to identify possible new protein variants and determine the distribution of variants across a variety of 18 taurine and indicine cattle breeds applying a DNA-based sequencing approach. To this end, the open reading frames of the respective genes (LALBA and LGB) were sequenced in 476 animals. Within the LALBA gene, a previously unknown synonymous and a previously undesignated non-synonymous nucleotide exchange were identified. Furthermore, two known α-LA variants (A and B) and four known β-LG variants (A, B, C and W) were determined. The occurrence of typical indicine variants in some taurine cattle breeds, such as Suisse Eringer, German Hinterwälder and Hungarian Grey Steppe, further supports the hypothesis of ancient Bos indicus introgression into (peri-)alpine cattle breeds.

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.

Genome and transcriptome sequencing identifies breeding targets in the orphan crop tef (Eragrostis tef)

Background: Tef (Eragrostis tef), an indigenous cereal critical to food security in the Horn of Africa, is rich in minerals and protein, resistant to many biotic and abiotic stresses and safe for diabetics as well as sufferers of immune reactions to wheat gluten. We present the genome of tef, the first species in the grass subfamily Chloridoideae and the first allotetraploid assembled de novo. We sequenced the tef genome for marker-assisted breeding, to shed light on the molecular mechanisms conferring tef's desirable nutritional and agronomic properties, and to make its genome publicly available as a community resource. Results: The draft genome contains 672 Mbp representing 87% of the genome size estimated from flow cytometry. We also sequenced two transcriptomes, one from a normalized RNA library and another from unnormalized RNASeq data. The normalized RNA library revealed around 38000 transcripts that were then annotated by the SwissProt group. The CoGe comparative genomics platform was used to compare the tef genome to other genomes, notably sorghum. Scaffolds comprising approximately half of the genome size were ordered by syntenic alignment to sorghum producing tef pseudo-chromosomes, which were sorted into A and B genomes as well as compared to the genetic map of tef. The draft genome was used to identify novel SSR markers, investigate target genes for abiotic stress resistance studies, and understand the evolution of the prolamin family of proteins that are responsible for the immune response to gluten. Conclusions: It is highly plausible that breeding targets previously identified in other cereal crops will also be valuable breeding targets in tef. The draft genome and transcriptome will be of great use for identifying these targets for genetic improvement of this orphan crop that is vital for feeding 50 million people in the Horn of Africa.

Whole-exome sequencing reveals a C-terminal germline variant in CEBPA-associated acute myeloid leukemia: 45-year follow-up of a large family.

Familial acute myeloid leukemia is rare and linked to germline mutations in RUNX1, GATA2 or CCAAT/enhancer binding protein-α (CEBPA). We re-evaluated a large family with acute myeloid leukemia originally seen at NIH in 1969. We utilized whole-exome sequencing to study this family, and conducted in silico bioinformatics analysis, protein structural modeling and laboratory experiments to assess the impact of the identified CEBPA Q311P mutation. Unlike most previously identified germline mutations in CEBPA, which were N-terminal frameshift mutations, we identified a novel Q311P variant that was located in the C-terminal bZip domain of C/EBPα. Protein structural modeling suggested that the Q311P mutation alters the ability of the CEBPA dimer to bind DNA. Electrophoretic mobility shift assays showed that the Q311P mutant had attenuated binding to DNA, as predicted by the protein modeling. Consistent with these findings, we found that the Q311P mutation has reduced transactivation, consistent with a loss-of-function mutation. From 45 years of follow-up, we observed incomplete penetrance (46%) of CEBPA Q311P. This study of a large multi-generational pedigree reveals that a germline mutation in the C-terminal bZip domain can alter the ability of C/EBP-α to bind DNA and reduces transactivation, leading to acute myeloid leukemia.