Genome-wide analysis of alternative splicing in cow: implications in bovine as a model for human diseases

Functional analysis of the GPAT4 gene mutation predicted to affect splicing

The GPAT4 gene is considered as a potential functional candidate for single nucleotide polymorphism (SNP) studies in dairy cattle breeding due to its association with dairy performance in cattle by encoding an enzyme responsible for the presence of diacylglycerols and triacylglycerols in milk. Using the example of the GPAT4 gene, we applied the minigene splicing assay to analyze the functional consequences of its variant that was predicted to affect normal splicing. The results of functional analysis revealed the sequence variations (rs442541537), transfection experiments in a wild type and mutant cell line model system demonstrated that the investigated mutation in the second intron of the GPAT4 gene was responsible for the presence of a second exon in mature messenger RNA (mRNA). The cases of its absence in the spliced mature mRNA transcript resulted in a truncated dysfunctional protein due to the appearance of a stop codon. Thus, the discovered SNP led to alternative splicing in pre-mRNA by the 'cassette exon' ('exon skipping') mechanism. The studied mutation can potentially be a molecular genetic marker for alternative splicing for the GPAT4 gene and, therefore contributes to economic benefits in cattle breeding programs.

Alternative splicing of pre-mRNA modulates the immune response in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis

Little is known about the role of alternative splicing (AS) in regulating gene expression in Mycobacteria-infected individuals in distinct stages of infection. Pre-mRNA AS consists of the removal of introns and the assembly of exons contained in eukaryotic genes. AS events can influence transcript stability or structure with important physiological consequences. Using RNA-Seq data from peripheral blood (PB) and ileocecal valve (ICV) samples collected from Holstein cattle with focal and diffuse paratuberculosis (PTB)-associated histopathological lesions in gut tissues and without lesions (controls), we detected differential AS profiles between the infected and control groups. Four of the identified AS events were experimentally validated by reverse transcription-digital droplet PCR (RT-ddPCR). AS events in several genes correlated with changes in gene expression. In the ICV of animals with diffuse lesions, for instance, alternatively spliced genes correlated with changes in the expression of genes involved in endocytosis, antigen processing and presentation, complement activation, and several inflammatory and autoimmune diseases in humans. Taken together, our results identified common mechanisms of AS involvement in the pathogenesis of PTB and human diseases and shed light on novel diagnostic and therapeutic interventions to control these diseases.

MAEL gene contributes to bovine testicular development through the m5C-mediated splicing

Knowledge of RNA molecules regulating testicular development and spermatogenesis in bulls is essential for elite bull selection and an ideal breeding program. Herein, we performed direct RNA sequencing (DRS) to explore the functional characterization of RNA molecules produced in the testicles of 9 healthy Simmental bulls at three testicular development stages (prepuberty, puberty, and postpuberty). We identified 5,043 differentially expressed genes associated with testicular weight. These genes exhibited more alternative splicing at sexual maturity, particularly alternative 3' (A3) and 5' (A5) splice sites usage and exon skipping (SE). The expression of hub genes in testicular developmental stages was also affected by both m6A and m5C RNA modifications. We found m5C-mediated splicing events significantly (p < 0.05) increased MAEL gene expression at the isoform level, likely promoting spermatogenesis. Our findings highlight the complexity of RNA processing and expression as well as the regulation of transcripts involved in testicular development and spermatogenesis.

Characterization of peripheral white blood cells transcriptome to unravel the regulatory signatures of bovine subclinical mastitis resistance

Subclinical bovine mastitis is a pathogenic infection of the breast characterized by a marked decrease in milk production and quality. As it has no obvious clinical symptoms, diagnosis and treatment are challenging. Therefore, searching for biomarkers in cows’ peripheral white blood cells is valuable for preventing and treating subclinical mastitis. Thus, in this study, the transcriptome of peripheral blood from healthy and subclinical mastitis cows was characterized to find the regulatory signatures of bovine subclinical mastitis using RNA-seq. A total of 287 differentially expressed genes (DEGs) and 70 differentially expressed lncRNAs (DELs) were detected, and 37 DELs were documented near known Quantitative Trait Loci (QTL) associated with the mastitis of cows. Bioinformatic analysis indicated that lncRNAs MSTRG25101.2, MSTRG.56327.1, and MSTRG.18968.1, which are adjacent to the SCS QTL and SCC QTL, may be candidate lncRNAs that influence the pathogenesis of mastitis in cows by up-regulating the expression of genes TLR4, NOD2, CXCL8, and OAS2. Moreover, the alternative splicing (AS) pattern of transcriptional sequence differences between healthy cows and subclinical mastitis cows suggested a molecular mechanism of mastitis resistance and susceptibility. A total of 2,212 differential alternative splicing (DAS) events, corresponding to 1,621 unique DAS genes, were identified in both groups and significantly enriched in immune and inflammatory pathways. Of these, 29 DAS genes were subject to regulation by 32 alternative splicing SNPs, showing diverse and specific splicing patterns and events. It is hypothesized that the PIK3C2B and PPRPF8 splice variants associated with AS SNPs (rs42705933 and rs133847062) may be risk factors for susceptibility to bovine subclinical mastitis. Altogether, these key blood markers associated with resistance to subclinical mastitis and SNPs associated with alternative splicing of genes provide the basis for genetic breeding for resistance to subclinical mastitis in cows.

Gene expression adjustment of inflammatory mechanisms in dairy cow mammary gland parenchyma during host defense against staphylococci

The aim of the study was to identify differences in the expression of splice variants of the PRMT2, LTF and C4A genes in the mammary glands of healthy dairy cows and those infected with staphylococci. An expression study was conducted on 38 Polish Holstein-Friesian dairy cows who were removed from the herd owing to subclinical and chronic mastitic or reproductive issues. Two days before slaughter, milk samples were taken for microbiological analysis and examined for the presence of bacteria. The mammary gland parenchyma samples with a predominance of secretory tissue were taken; these were divided into three groups according to the health status of the mammary gland: H (without pathogenic bacteria in milk), CoNS (with coagulase-negative staphylococci in milk), and CoPS (with coagulase-positive staphylococci in milk). Two of the investigated genes, LTF and C4A, demonstrated variants unequivocally expressed in infected tissue. Two LTF gene variants were found to be associated with cow health status, and with the type of bacteria causing mastitis (CoPS or CoNS). In addition, the expression of C4A isoforms differed with regard to mastitis etiology groups. The comprehensive evaluation of PRMT2 transcript suggested that the gene may also be involved in course of mastitis: two of four PRMT2 transcripts showed increased expression in the mammary gland of the CoPS group compared to controls. The obtained results are important for the knowledge on the etiology of bovine mastitis. The effects of the identified mastitis-relevant splice variants need to be further explored on the protein level to verify the suitability of splice variants and recognize their contribution towards the disease phenotypes and course.

Application of transposon systems in the transgenesis of bovine somatic and germ cells

Background

Several DNA transposons including PiggyBac (PB), Sleeping Beauty (SB), and Tol2 have been applied as effective means for of transgenesis in many species. Cattle are not typically experimental animals, and relatively little verification has been presented on this species. Thus, the goal here was to determine the applicability of three transposon systems in somatic and embryo cells in cattle, while also investigating which of the three systems is appropriate for each cell type. Green fluorescent protein (GFP)-expressing transposon systems were used for electroporation and microinjection in the somatic cells and embryo stage, respectively. After transfection, the GFP-positive cells or blastocysts were observed through fluorescence, while the transfection efficiency was calculated by FACS.

Results

In bovine somatic cells, the PB (63.97 ± 11.56) showed the highest efficiency of the three systems (SB: 50.74 ± 13.02 and Tol2: 16.55 ± 5.96). Conversely, Tol2 (75.00%) and SB (70.00%) presented a higher tendency in the embryonic cells compared to PB (42.86%).

Conclusions

These results demonstrate that these three transposon systems can be used in bovine somatic cells and embryos as gene engineering experimental methods. Moreover, they demonstrate which type of transposon system to apply depending on the cell type.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03252-1.

Differences in gene expression and variable splicing events of ovaries between large and small litter size in Chinese Xiang pigs

BACKGROUND

Although lots of quantitative trait loci (QTLs) and genes present roles in litter size of some breeds, the information might not make it clear for the huge diversity of reproductive capability in pig breeds. To elucidate the inherent mechanisms of heterogeneity of reproductive capability in litter size of Xiang pig, we performed transcriptome analysis for the expression profile in ovaries using RNA-seq method.

RESULTS

We identified 1,419 up-regulated and 1,376 down-regulated genes in Xiang pigs with large litter size. Among them, 1,010 differentially expressed genes (DEGs) were differently spliced between two groups with large or small litter sizes. Based on GO and KEGG analysis, numerous members of genes were gathered in ovarian steroidogenesis, steroid biosynthesis, oocyte maturation and reproduction processes.

CONCLUSIONS

Combined with gene biological function, twelve genes were found out that might be related with the reproductive capability of Xiang pig, of which, eleven genes were recognized as hub genes. These genes may play a role in promoting litter size by elevating steroid and peptide hormones supply through the ovary and facilitating the processes of ovulation and in vivo fertilization.

Evolution of tissue and developmental specificity of transcription start sites in Bos taurus indicus

To further the understanding of the evolution of transcriptional regulation, we profiled genome-wide transcriptional start sites (TSSs) in two sub-species, Bos taurus taurus and Bos taurus indicus, that diverged approximately 500,000 years ago. Evolutionary and developmental-stage differences in TSSs were detected across the sub-species, including translocation of dominant TSS and changes in TSS distribution. The 16% of all SNPs located in significant differentially used TSS clusters across sub-species had significant shifts in allele frequency (472 SNPs), indicating they may have been subject to selection. In spleen and muscle, a higher relative TSS expression was observed in Bos indicus than Bos taurus for all heat shock protein genes, which may be responsible for the tropical adaptation of Bos indicus.

Forutan et al. measure RNA expression and map transcription start sites (TSSs) between two sub-species of cattle in adult and fetal tissues and demonstrate translocation and changes in TSS distribution between the sub-species. This study provides insight into cattle adaptation and provides support of rapidly evolving TSS between two recently diverged subspecies.

Characterization of Alternative Splicing (AS) Events during Chicken (Gallus gallus) Male Germ-Line Stem Cell Differentiation with Single-Cell RNA-seq

Simple Summary

Studies have shown that alternative splicing (AS) has been utilized in a wide variety of life processes. However, there are very few studies on AS during germ cell development. In this study, we preliminarily investigated the variation of variable shear events during the formation of chicken germ cells through the RNA-seq data analysis of embryonic stem cells (ESCs), gonad PGCs (gPGCs), and spermatogonia stem cells (SSCs), and the critical AS mode for several crucial stage-specific genes, which were identified during germ cell development. The results of this study lay a theoretical foundation for further analysis of the regulation mechanism of key genes involved in germ cell formation.

Abstract

Alternative splicing (AS) is a ubiquitous, co-transcriptional, and post-transcriptional regulation mechanism during certain developmental processes, such as germ cell differentiation. A thorough understanding of germ cell differentiation will help us to open new avenues for avian reproduction, stem cell biology, and advances in medicines for human consumption. Here, based on single-cell RNA-seq, we characterized genome-wide AS events in manifold chicken male germ cells: embryonic stem cells (ESCs), gonad primordial germ cells (gPGCs), and spermatogonia stem cells (SSCs). A total of 38,494 AS events from 15,338 genes were detected in ESCs, with a total of 48,955 events from 14,783 genes and 49,900 events from 15,089 genes observed in gPGCs and SSCs, respectively. Moreover, this distribution of AS events suggests the diverse splicing feature of ESCs, gPGCs, and SSCs. Finally, several crucial stage-specific genes, such as NANOG, POU5F3, LIN28B, BMP4, STRA8, and LHX9, were identified in AS events that were transmitted in ESCs, gPGCs, and SSCs. The gene expression results of the RNA-seq data were validated by qRT-PCR. In summary, we provided a comprehensive atlas of the genome-wide scale of the AS event landscape in male chicken germ-line cells and presented its distribution for the first time. This research may someday improve treatment options for men suffering from male infertility.

Profiling Novel Alternative Splicing within Multiple Tissues Provides Useful Insights into Porcine Genome Annotation

Alternative splicing (AS) is a process during gene expression that results in a single gene coding for different protein variants. AS contributes to transcriptome and proteome diversity. In order to characterize AS in pigs, genome-wide transcripts and AS events were detected using RNA sequencing of 34 different tissues in Duroc pigs. In total, 138,403 AS events and 29,270 expressed genes were identified. An alternative donor site was the most common AS form and accounted for 44% of the total AS events. The percentage of the other three AS forms (exon skipping, alternative acceptor site, and intron retention) was approximately 19%. The results showed that the most common AS events involving alternative donor sites could produce different transcripts or proteins that affect the biological processes. The expression of genes with tissue-specific AS events showed that gene functions were consistent with tissue functions. AS increased proteome diversity and resulted in novel proteins that gained or lost important functional domains. In summary, these findings extend porcine genome annotation and highlight roles that AS could play in determining tissue identity.

Gene co-expression and alternative splicing analysis of key metabolic tissues to unravel the regulatory signatures of fatty acid composition in cattle

Increasing the healthy/unhealthy fatty acid (FA) ratio in meat is one of the urgent tasks required to address consumer concerns. However, the regulatory mechanisms ultimately resulting in FA profiles vary among animals and remain largely unknown. In this study, using ~1.2 Tb high-quality RNA-Seq-based transcriptomic data of 188 samples from four key metabolic tissues (rumen, liver, muscle, and backfat) together with the contents of 49 FAs in backfat, the molecular regulatory mechanisms of these tissues contributing to FA formation in cattle were explored. Using this large dataset, the alternative splicing (AS) events, one of the transcriptional regulatory mechanisms in four tissues were identified. The highly conserved and absent AS events were detected in rumen tissue, which may contribute to its functional differences compared with the other three tissues. In addition, the healthy/unhealthy FA ratio related AS events, differential expressed (DE) genes, co-expressed genes, and their functions in four tissues were analysed. Eight key genes were identified from the integrated analysis of DE, co-expressed, and AS genes between animals with high and low healthy/unhealthy FA ratios. This study provides an applicable pipeline for AS events based on comprehensive RNA-Seq analysis and improves our understanding of the regulatory mechanism of FAs in beef cattle.

Circular RNA TTN Acts As a miR-432 Sponge to Facilitate Proliferation and Differentiation of Myoblasts via the IGF2/PI3K/AKT Signaling Pathway

Circular RNAs (circRNAs) are ubiquitous endogenous RNA found in various organisms that can regulate gene expression in eukaryotes. However, little is known about potential roles for circRNAs in muscle development. We analyzed circRNA sequencing data of bovine skeletal muscle tissue and found differential expression of circTitin (circTTN) in fetal and adult bovine muscle tissue. We then further studied the role of circTTN in bovine muscle development. Overexpression and inhibition of circTTN together elicited its promoting roles in proliferation and differentiation of bovine primary myoblasts. Mechanistically, circTTN showed interaction with miR-432 by luciferase screening and RNA immunoprecipitation (RIP) assays. Additionally, miR-432 is a regulator of insulin-like growth factor 2 (IGF2), as indicated by luciferase activity, quantitative real-time PCR, and western blotting assays. Increased miR-432 expression inhibited the expression of IGF2, but this effect was remitted by circTTN. Conclusively, our results showed that circTTN promoted proliferation and differentiation of bovine primary myoblasts via competitively combining with miR-432 to activate the IGF2/phosphatidylinositol 3-kinase (PI3K)/AKT signaling pathway.

Production of an aberrant splice variant of CCL5 is not caused by genetic mutation in the mammary glands of mastitis‑infected Holstein cows

Genetic mutations, including single nucleotide polymorphisms (SNPs), result in aberrant alternatively splicing of gene and involves in susceptibility of inflammatory diseases, including bovine mastitis. C‑C motif chemokine ligand 5 (CCL5) is an immune‑associated gene, but its alternative splicing (AS) mechanism of gene expression has not yet been understood. The present study identified the splice variant of CCL5 and the compared differential expression of various transcripts between healthy and mastitic mammary gland tissue from cows. A novel transcript lacking exon 2 with a deletion of 112 bp (referred to as CCL5‑AS) was identified in the mammary gland. The expression of CCL5‑AS was lower compared with CCL5‑reference in the healthy and mastitic mammary tissues. A total of two novel SNPs (g.1647 C>T and g.1804 G>A) were identified in exon 2 of CCL5. Using the splicing mini‑gene reporter assay in bovine mammary epithelial cells (MAC‑T) and 293T cells, it was confirmed that the production of CCL5‑AS was not caused by the two SNPs. The present findings suggested that alternative splicing is one of the mechanisms of CCL5 expression regulation and is involved in mastitis infection, but that genetic mutation was not responsible for the generation of the abnormal transcript of CCL5 in cows.

A survey of transcriptome complexity in Sus scrofa using single-molecule long-read sequencing

Alternative splicing (AS) and fusion transcripts produce a vast expansion of transcriptomes and proteomes diversity. However, the reliability of these events and the extend of epigenetic mechanisms have not been adequately addressed due to its limitation of uncertainties about the complete structure of mRNA. Here we combined single-molecule real-time sequencing, Illumina RNA-seq and DNA methylation data to characterize the landscapes of DNA methylation on AS, fusion isoforms formation and lncRNA feature and further to unveil the transcriptome complexity of pig. Our analysis identified an unprecedented scale of high-quality full-length isoforms with over 28,127 novel isoforms from 26,881 novel genes. More than 92,000 novel AS events were detected and intron retention predominated in AS model, followed by exon skipping. Interestingly, we found that DNA methylation played an important role in generating various AS isoforms by regulating splicing sites, promoter regions and first exons. Furthermore, we identified a large of fusion transcripts and novel lncRNAs, and found that DNA methylation of the promoter and gene body could regulate lncRNA expression. Our results significantly improved existed gene models of pig and unveiled that pig AS and epigenetic modify were more complex than previously thought.

Transcriptome analyses reveal genes of alternative splicing associated with muscle development in chickens

Alternative splicing (AS) of pre-mRNA is a central mode of genetic regulation in higher eukaryotes. High-throughput experimental verification of alternative splice forms, functional characterization, and regulation of alternative splicing are key directions for research. However, little information is available on the transcriptome-wide changes during different ages in different chicken breeds. In this study, the sequencing reads of chicken muscle tissues collected from White feather broiler (day 42) and Luning Chicken (day 70, 120, 150) were mapped to the chicken genome. Results showed that a total of 16,958 genes were annotated, with 2230 differentially expressed genes (DEGs) when comparing White feather broiler and Luning Chicken, and an average of 700 DEGs when comparing different ages in Luning Chicken. Functional classification by Gene Ontology (GO) and pathways analysis for selecting the genes showed most DEGs were related to muscle development and immune response. Of the 16,958 genes, a total of 6249 genes (36.85%) underwent AS events, and over 40% were specifically expressed in each library. Additionally, 6 DEGs (SRPK3, ENSGALG00000022884, CCL4, GATM, SESN1, PTTG1IP) involved in muscle development and immunity response were found to be alternatively spliced among all the four muscle tissues. In conclusion, we present a complete dataset involving the spatial and temporal transcriptome of chicken muscle tissue using RNA -seq. These data will facilitate the understanding of the effects of breed and age on the development of muscle and uncover that AS events of candidate genes may have important functional roles in regulation of muscle development in chicken.

Analysis of alternative splicing events by RNA sequencing in the ovaries of Xiang pig at estrous and diestrous

Alternative splicing (AS) is an important mechanism for controlling gene expression, which regulates multiple biological processes in higher organisms. Chinese indigenous Xiang pig has distinctive biological characteristics, such as small size, early sexual maturity, lower litter size and not very clear exhibition of estrous behaviors. To further understand how AS responds to estrous cycles in Xiang pig, the genome-wide analysis of AS events was performed by RNA-seq method in Xiang pig ovaries at diestrous and estrous. Using ASprofile program, we analyzed twelve basic AS events in Xiang pig ovaries and identified 68,775 AS events in 15,142 genes from diestrous ovaries and 69,493 AS events in 15,291 genes from the estrous ovaries with average 4.54 splicing events. 94.4-95.5% of expressed genes underwent alternative splicing in this tissues. The frequencies of AS events were similar to each other at diestrous and estrous. Transcription start site (TSS) was the predominant type of AS events, followed by transcription terminal site (TTS), and skipped exon (SKIP). The remaining type of AS events, e.g., intron retention (IR) and alternative exon ends (AE), showed the lower frequencies. Further comparison analysis of gene expression indicated that 4,433 genes had at least one splice variant differentially expressed during estrous, whereas only 2,382 of them were differentially expressed at gene level. Numerous genes involved in gonad development and hormone metabolism were differentially regulated through AS. Twelve genes with different types of alternatively splicing were validated by using RT-PCR method. The GO annotation and KEGG pathway analysis clearly revealed that a lot of DEGs (differentially expressed genes) and DSGs (differentially spliced genes) were involved in follicular development and ovarian steroid biosynthesis. A large number of DSGs, although not differentially expressed, were enriched in circadian rhythm and several signaling pathways. These pathways potentially regulated the female animal reproductive function in gene and/or AS level. Our results suggested that alternative splicing play an essential role in regulation of gene expression in female pigs during estrous. Numerous genes involved in gonad development, hormone metabolism, circadian rhythm were differentially regulated through alternative splicing.

Differential expression of mRNA isoforms in the skeletal muscle of pigs with distinct growth and fatness profiles

Background

The identification of genes differentially expressed in the skeletal muscle of pigs displaying distinct growth and fatness profiles might contribute to identify the genetic factors that influence the phenotypic variation of such traits. So far, the majority of porcine transcriptomic studies have investigated differences in gene expression at a global scale rather than at the mRNA isoform level. In the current work, we have investigated the differential expression of mRNA isoforms in the gluteus medius (GM) muscle of 52 Duroc HIGH (increased backfat thickness, intramuscular fat and saturated and monounsaturated fatty acids contents) and LOW pigs (opposite phenotype, with an increased polyunsaturated fatty acids content).

Results

Our analysis revealed that 10.9% of genes expressed in the GM muscle generate alternative mRNA isoforms, with an average of 2.9 transcripts per gene. By using two different pipelines, one based on the CLC Genomics Workbench and another one on the STAR, RSEM and DESeq2 softwares, we have identified 10 mRNA isoforms that both pipelines categorize as differentially expressed in HIGH vs LOW pigs (P-value < 0.01 and ±0.6 log₂fold-change). Only five mRNA isoforms, produced by the ITGA5, SEMA4D, LITAF, TIMP1 and ANXA2 genes, remain significant after correction for multiple testing (q-value < 0.05 and ±0.6 log₂fold-change), being upregulated in HIGH pigs.

Conclusions

The increased levels of specific ITGA5, LITAF, TIMP1 and ANXA2 mRNA isoforms in HIGH pigs is consistent with reports indicating that the overexpression of these four genes is associated with obesity and metabolic disorders in humans. A broader knowledge about the functional attributes of these mRNA variants would be fundamental to elucidate the consequences of transcript diversity on the determinism of porcine phenotypes of economic interest.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4515-2) contains supplementary material, which is available to authorized users.

Transcriptome Bioinformatical Analysis of Vertebrate Stages of Schistosoma japonicum Reveals Alternative Splicing Events

Alternative splicing is a molecular process that contributes greatly to the diversification of proteome and to gene functions. Understanding the mechanisms of stage-specific alternative splicing can provide a better understanding of the development of eukaryotes and the functions of different genes. Schistosoma japonicum is an infectious blood-dwelling trematode with a complex lifecycle that causes the tropical disease schistosomiasis. In this study, we analyzed the transcriptome of Schistosoma japonicum to discover alternative splicing events in this parasite, by applying RNA-seq to cDNA library of adults and schistosomula. Results were validated by RT-PCR and sequencing. We found 11,623 alternative splicing events among 7,099 protein encoding genes and average proportion of alternative splicing events per gene was 42.14%. We showed that exon skip is the most common type of alternative splicing events as found in high eukaryotes, whereas intron retention is the least common alternative splicing type. According to intron boundary analysis, the parasite possesses same intron boundaries as other organisms, namely the classic "GT-AG" rule. And in alternative spliced introns or exons, this rule is less strict. And we have attempted to detect alternative splicing events in genes encoding proteins with signal peptides and transmembrane helices, suggesting that alternative splicing could change subcellular locations of specific gene products. Our results indicate that alternative splicing is prevalent in this parasitic worm, and that the worm is close to its hosts. The revealed secretome involved in alternative splicing implies new perspective into understanding interaction between the parasite and its host.

Molecular characterization, alternative splicing and expression analysis of bovine DBC1

Deleted in breast cancer 1 (DBC1, KIAA1967, p30 DBC) is a novel protein that has been recently shown to bind and regulate SIRT1. Loss of function of DBC1 increased SIRT1 deacetylase activity, which promotes "browning" of WAT by deacetylating peroxisome proliferator-activated receptor (PPARγ) on Lys268 and Lys293. In the present study, we have cloned and characterized the bovine DBC1 gene. Two transcript variants of bovine DBC1 were identified, designated DBC1-A and DBC1-B, respectively, which were both located in nucleus. Protein sequence analysis showed that DBC1-A was well conserved across species. Expression analysis of DBC1 in seven different tissues of calves and bulls by RT-PCR indicated that the two transcripts were ubiquitously expressed. However, the relatively level of DBC1-A was higher when compared to DBC1-B in all examined tissues. Surprisingly, the expression of DBC1-A was extraordinary high in calves adipose tissue, which implicated its potential key role in regulating calve adipocyte development. These findings provide new insight into our understanding of the biochemical characteristics and physiological role of bovine DBC1.

Increased expression of a novel splice variant of the complement component 4 (C4A) gene in mastitis-infected dairy cattle

The complement system helps in the direct lysis of invading pathogens and modulates phagocytic, humoral and cellular immune responses. Complement 4 is a critical component in complement activity and protection against many bacterial pathogens because it is essential to classical and lectin activation pathways. We used reverse transcription and PCR to investigate alternative splicing and expression of the complement component 4 (C4A) gene in Chinese Holstein cattle. The PCR products were cloned and sequenced. A novel splice variant involving intron 10 was identified, which we named C4A-AS. To examine how C4A gene activity is affected by bovine mastitis, six Chinese Holstein cattle were divided into healthy (non-mastitic) and Staphylococcus aureus-induced mastitic groups. Real-time quantitative PCR (qRT-PCR) revealed that the C4A-complete and C4A-AS transcripts are expressed at significantly different levels in healthy cows, while there were no significant differences in the mastitic group (P = 0.257). Expression of C4A-AS increased significantly when mastitis developed. We also examined the expression of C4A-complete and C4A-AS in several tissues (liver, heart, spleen, lung, kidney, tongue, and muscle). The two transcripts were expressed in all of these tissues but there were no significant differences in expression between healthy and mastitic cows. We therefore conclude that the C4A-complete transcript is the main transcript under normal physiological conditions, while C4A-AS is augmented when mastitis develops.

Transcriptome analysis of the model protozoan, Tetrahymena thermophila, using Deep RNA sequencing

BACKGROUND

The ciliated protozoan Tetrahymena thermophila is a well-studied single-celled eukaryote model organism for cellular and molecular biology. However, the lack of extensive T. thermophila cDNA libraries or a large expressed sequence tag (EST) database limited the quality of the original genome annotation.

METHODOLOGY/PRINCIPAL FINDINGS

This RNA-seq study describes the first deep sequencing analysis of the T. thermophila transcriptome during the three major stages of the life cycle: growth, starvation and conjugation. Uniquely mapped reads covered more than 96% of the 24,725 predicted gene models in the somatic genome. More than 1,000 new transcribed regions were identified. The great dynamic range of RNA-seq allowed detection of a nearly six order-of-magnitude range of measurable gene expression orchestrated by this cell. RNA-seq also allowed the first prediction of transcript untranslated regions (UTRs) and an updated (larger) size estimate of the T. thermophila transcriptome: 57 Mb, or about 55% of the somatic genome. Our study identified nearly 1,500 alternative splicing (AS) events distributed over 5.2% of T. thermophila genes. This percentage represents a two order-of-magnitude increase over previous EST-based estimates in Tetrahymena. Evidence of stage-specific regulation of alternative splicing was also obtained. Finally, our study allowed us to completely confirm about 26.8% of the genes originally predicted by the gene finder, to correct coding sequence boundaries and intron-exon junctions for about a third, and to reassign microarray probes and correct earlier microarray data.

CONCLUSIONS/SIGNIFICANCE

RNA-seq data significantly improve the genome annotation and provide a fully comprehensive view of the global transcriptome of T. thermophila. To our knowledge, 5.2% of T. thermophila genes with AS is the highest percentage of genes showing AS reported in a unicellular eukaryote. Tetrahymena thus becomes an excellent unicellular model eukaryote in which to investigate mechanisms of alternative splicing.

Extending Asia Pacific bioinformatics into new realms in the "-omics" era

The 2009 annual conference of the Asia Pacific Bioinformatics Network (APBioNet), Asia's oldest bioinformatics organisation dating back to 1998, was organized as the 8^th International Conference on Bioinformatics (InCoB), Sept. 7-11, 2009 at Biopolis, Singapore. Besides bringing together scientists from the field of bioinformatics in this region, InCoB has actively engaged clinicians and researchers from the area of systems biology, to facilitate greater synergy between these two groups. InCoB2009 followed on from a series of successful annual events in Bangkok (Thailand), Penang (Malaysia), Auckland (New Zealand), Busan (South Korea), New Delhi (India), Hong Kong and Taipei (Taiwan), with InCoB2010 scheduled to be held in Tokyo, Japan, Sept. 26-28, 2010. The Workshop on Education in Bioinformatics and Computational Biology (WEBCB) and symposia on Clinical Bioinformatics (CBAS), the Singapore Symposium on Computational Biology (SYMBIO) and training tutorials were scheduled prior to the scientific meeting, and provided ample opportunity for in-depth learning and special interest meetings for educators, clinicians and students. We provide a brief overview of the peer-reviewed bioinformatics manuscripts accepted for publication in this supplement, grouped into thematic areas. In order to facilitate scientific reproducibility and accountability, we have, for the first time, introduced minimum information criteria for our pubilcations, including compliance to a Minimum Information about a Bioinformatics Investigation (MIABi). As the regional research expertise in bioinformatics matures, we have delineated a minimum set of bioinformatics skills required for addressing the computational challenges of the "-omics" era.