Longitudinal blood transcriptomic analysis to identify molecular regulatory patterns of bovine respiratory disease in beef cattle.

Bovine respiratory disease (BRD) is the most common disease in beef cattle and leads to considerable economic losses in both beef and dairy cattle. It is important to uncover the molecular mechanisms underlying BRD and to identify biomarkers for early identification of BRD cattle in order to address its impact on production and welfare. In this study, a longitudinal transcriptomic analysis was conducted using blood samples collected from 24 beef cattle at three production stages in the feedlot: 1) arrival (Entry group); 2) when identified as sick (diagnosed as BRD) and separated for treatment (Pulled); 3) prior to marketing (Close-out, representing healthy animals). Expressed genes were significantly different in the same animal among Entry, Pulled and Close-out stages (false discovery rate (FDR) < 0.01 & |Fold Change| > 2). Beef steers at both Entry and Pulled stages presented obvious difference in GO terms (FDR < 0.05) and affected biological functions (FDR < 0.05 & |Z-score| > 2) when compared with animals at Close-out. However, no significant functional difference was observed between Entry and Pulled animals. The interferon signaling pathway showed the most significant difference between animals at Entry/Pulled and Close-out stages (P < .001 & |Z-score| > 2), suggesting the animals initiated antiviral responses at an early stage of infection. Six key genes including IFI6, IFIT3, ISG15, MX1, and OAS2 were identified as biomarkers to predict and recognize sick cattle at Entry. A gene module with 169 co-expressed genes obtained from WGCNA analysis was most positively correlated (R = 0.59, P = 6E-08) with sickness, which was regulated by 11 transcription factors. Our findings provide an initial understanding of the BRD infection process in the field and suggests a subset of novel marker genes for identifying BRD in cattle at an early stage of infection.

Creation of miniature pig model of human Waardenburg syndrome type 2A by ENU mutagenesis.

Human Waardenburg syndrome 2A (WS2A) is a dominant hearing loss (HL) syndrome caused by mutations in the microphthalmia-associated transcription factor (MITF) gene. In mouse models with MITF mutations, WS2A is transmitted in a recessive pattern, which limits the study of hearing loss (HL) pathology. In the current study, we performed ENU (ethylnitrosourea) mutagenesis that resulted in substituting a conserved lysine with a serine (p. L247S) in the DNA-binding domain of the MITF gene to generate a novel miniature pig model of WS2A. The heterozygous mutant pig (MITF +/L247S) exhibits a dominant form of profound HL and hypopigmentation in skin, hair, and iris, accompanied by degeneration of stria vascularis (SV), fused hair cells, and the absence of endocochlear potential, which indicate the pathology of human WS2A. Besides hypopigmentation and bilateral HL, the homozygous mutant pig (MITF L247S/L247S) and CRISPR/Cas9-mediated MITF bi-allelic knockout pigs both exhibited anophthalmia. Three WS2 patients carrying MITF mutations adjacent to the corresponding region were also identified. The pig models resemble the clinical symptom and molecular pathology of human WS2A patients perfectly, which will provide new clues for better understanding the etiology and development of novel treatment strategies for human HL.

Pilot study of large-scale production of mutant pigs by ENU mutagenesis.

N-ethyl-N-nitrosourea (ENU) mutagenesis is a powerful tool to generate mutants on a large scale efficiently, and to discover genes with novel functions at the whole-genome level in Caenorhabditis elegans, flies, zebrafish and mice, but it has never been tried in large model animals. We describe a successful systematic three-generation ENU mutagenesis screening in pigs with the establishment of the Chinese Swine Mutagenesis Consortium. A total of 6,770 G1 and 6,800 G3 pigs were screened, 36 dominant and 91 recessive novel pig families with various phenotypes were established. The causative mutations in 10 mutant families were further mapped. As examples, the mutation of SOX10 (R109W) in pig causes inner ear malfunctions and mimics human Mondini dysplasia, and upregulated expression of FBXO32 is associated with congenital splay legs. This study demonstrates the feasibility of artificial random mutagenesis in pigs and opens an avenue for generating a reservoir of mutants for agricultural production and biomedical research.

Transcriptome microRNA profiling of bovine mammary epithelial cells challenged with Escherichia coli or Staphylococcus aureus bacteria reveals pathogen directed microRNA expression profiles.

BACKGROUND: MicroRNAs (miRNAs) can post-transcriptionally regulate gene expression and have been shown to be critical regulators to the fine-tuning of epithelial immune responses. However, the role of miRNAs in bovine responses to E. coli and S. aureus, two mastitis causing pathogens, is not well understood. RESULTS: The global expression of miRNAs in bovine mammary epithelial cells (MAC-T cells) challenged with and without heat-inactivated Staphylococcus aureus (S. aureus) or Escherichia coli (E. coli) bacteria at 0, 6, 12, 24, and 48 hr was profiled using RNA-Seq. A total of 231 known bovine miRNAs were identified with more than 10 counts per million in at least one of 13 libraries and 5 miRNAs including bta-miR-21-5p, miR-27b, miR-22-3p, miR-184 and let-7f represented more than 50% of the abundance. One hundred and thirteen novel miRNAs were also identified and more than one third of them belong to the bta-miR-2284 family. Seventeen miRNAs were significantly (P < 0.05) differentially regulated by the presence of pathogens. E. coli initiated an earlier regulation of miRNAs (6 miRNAs differentially regulated within the first 6 hrs post challenge as compared to 1 miRNA for S. aureus) while S. aureus presented a delayed response. Five differentially expressed miRNAs (bta-miR-184, miR-24-3p, miR-148, miR-486 and let-7a-5p) were unique to E. coli while four (bta-miR-2339, miR-499, miR-23a and miR-99b) were unique to S. aureus. In addition, our study revealed a temporal differential regulation of five miRNAs (bta-miR-193a-3p, miR-423-5p, miR-30b-5p, miR-29c and miR-un116) in unchallenged cells. Target gene predictions of pathogen differentially expressed miRNAs indicate a significant enrichment in gene ontology functional categories in development/cellular processes, biological regulation as well as cell growth and death. Furthermore, target genes were significantly enriched in several KEGG pathways including immune system, signal transduction, cellular process, nervous system, development and human diseases. CONCLUSION: Using next-generation sequencing, our study identified a pathogen directed differential regulation of miRNAs in MAC-T cells with roles in immunity and development. Our study provides a further confirmation of the involvement of mammary epithelia cells in contributing to the immune response to infecting pathogens and suggests the potential of miRNAs to serve as biomarkers for diagnosis and development of control measures.

MicroRNAs in bovine adipogenesis: genomic context, expression and function.

BACKGROUND: MicroRNAs (miRNAs) are small non-coding RNAs found to regulate several biological processes including adipogenesis. Understanding adipose tissue regulation is critical for beef cattle as fat is an important determinant of beef quality and nutrient value. This study analyzed the association between genomic context characteristics of miRNAs with their expression and function in bovine adipose tissue. Twenty-four subcutaneous adipose tissue biopsies were obtained from eight British-continental crossbred steers at 3 different time points. Total RNA was extracted and miRNAs were profiled using a miRNA microarray with expression further validated by qRT-PCR. RESULTS: A total of 224 miRNAs were detected of which 155 were expressed in all steers (n = 8), and defined as the core miRNAs of bovine subcutaneous adipose tissue. Core adipose miRNAs varied in terms of genomic location (59.5% intergenic, 38.7% intronic, 1.2% exonic, and 0.6% mirtron), organization (55.5% non-clustered and 44.5% clustered), and conservation (49% highly conserved, 14% conserved and 37% poorly conserved). Clustered miRNAs and highly conserved miRNAs were more highly expressed (p < 0.05) and had more predicted targets than non-clustered or less conserved miRNAs (p < 0.001). A total of 34 miRNAs were coordinately expressed, being part of six identified relevant networks. Two intronic miRNAs (miR-33a and miR-1281) were confirmed to have coordinated expression with their host genes, transcriptional factor SREBF2 and EP300 (a transcriptional co-activator of transcriptional factor C/EBPα), respectively which are involved in lipid metabolism, suggesting these miRNAs may also play a role in regulation of bovine lipid metabolism/adipogenesis. Furthermore, a total of 17 bovine specific miRNAs were predicted to be involved in the regulation of energy balance in adipose tissue. CONCLUSIONS: These findings improve our understanding on the behavior of miRNAs in the regulation of bovine adipogenesis and fat metabolism as it reveals that miRNA expression patterns and functions are associated with miRNA genomic location, organization and conservation.

Elucidation of molecular mechanisms of physiological variations between bovine subcutaneous and visceral fat depots under different nutritional regimes.

Adipose tissue plays a critical role in energy homeostasis and metabolism. There is sparse understanding of the molecular regulation at the protein level of bovine adipose tissues, especially within different fat depots under different nutritional regimes. The objective of this study was to analyze the differences in protein expression between bovine subcutaneous and visceral fat depots in steers fed different diets and to identify the potential regulatory molecular mechanisms of protein expression. Subcutaneous and visceral fat tissues were collected from 16 British-continental steers (15.5 month old) fed a high-fat diet (7.1% fat, n=8) or a control diet (2.7% fat, n=8). Protein expression was profiled using label free quantification LC-MS/MS and expression of selected transcripts was evaluated using qRT-PCR. A total of 682 proteins were characterized and quantified with fat depot having more impact on protein expression, altering the level of 51.0% of the detected proteins, whereas diet affected only 5.3%. Functional analysis revealed that energy production and lipid metabolism were among the main functions associated with differentially expressed proteins between fat depots, with visceral fat being more metabolically active than subcutaneous fat as proteins associated with lipid and energy metabolism were upregulated. The expression of several proteins was significantly correlated to subcutaneous fat thickness and adipocyte size, indicating their potential as adiposity markers. A poor correlation (r=0.245) was observed between mRNA and protein levels for 9 genes, indicating that many proteins may be subjected to post-transcriptional regulation. A total of 8 miRNAs were predicted to regulate more than 20% of lipid metabolism proteins differentially expressed between fat depots, suggesting that miRNAs play a role in adipose tissue regulation. Our results show that proteomic changes support the distinct metabolic and physiological characteristics observed between subcutaneous and visceral adipose tissue depots in cattle.

Altered microRNA expression in bovine subcutaneous and visceral adipose tissues from cattle under different diet.

BACKGROUND: MicroRNAs (miRNAs) are a class of molecular regulators found to participate in numerous biological processes, including adipogenesis in mammals. This study aimed to evaluate the differences of miRNA expression between bovine subcutaneous (backfat) and visceral fat depots (perirenal fat) and the dietary effect on miRNA expression in these fat tissues. METHODOLOGY/PRINCIPAL FINDINGS: Fat tissues were collected from 16 Hereford×Aberdeen Angus cross bred steers (15.5 month old) fed a high-fat diet (5.85% fat, n = 8) or control diet (1.95% fat, n = 8). Total RNA from each animal was subjected to miRNA microarray analysis using a customized Agilent miRNA microarray containing 672 bovine miRNA probes. Expression of miRNAs was not equal between fat depots as well as diets: 207 miRNAs were detected in both fat depots, while 37 of these were found to be tissue specific; and 169 miRNAs were commonly expressed under two diets while 75 were diet specific. The number of miRNAs detected per animal fed the high fat diet was higher than those fed control diet (p = 0.037 in subcutaneous fat and p = 0.002 visceral fat). Further qRT-PCR analysis confirmed that the expression of some miRNAs was highly influenced by diet (miR-19a, -92a, -92b, -101, -103, -106, -142-5p, and 296) or fat depot (miR-196a and -2454). CONCLUSIONS/SIGNIFICANCE: Our results revealed that the miRNA may differ among adipose depots and level of fat in the diet, suggesting that miRNAs may play a role in the regulation of bovine adipogenesis.

MicroRNA regulation in mammalian adipogenesis.

Adipogenesis, the complex development from preadipocytes or mesenchymal stem cells to mature adipocytes, is essential for fat formation and metabolism of adipose tissues in mammals. It has been reported to be regulated by hormones and various adipogenic transcription factors which are expressed as a transcriptional cascade promoting adipocyte differentiation, leading to the mature adipocyte phenotype. Recent findings indicate that microRNAs (miRNAs), a family of small RNA molecules of approximately 22 nucleotides in length, are involved in the regulatory network of many biological processes, including cell differentiation, through post-transcriptional regulation of transcription factors and/or other genes. In this review, we focus on the recent understanding of the roles of miRNAs in adipogenesis, including the most recent and relevant findings that support the role of several miRNAs as pro- or antiadipogenic factors regulating adipogenesis in mice, human and cattle to propose the future role of miRNA in adipogenesis of farm animal models.

Characterization of microRNA expression in bovine adipose tissues: a potential regulatory mechanism of subcutaneous adipose tissue development.

BACKGROUND: MicroRNAs (miRNAs), a family of small non-coding RNA molecules, appear to regulate animal lipid metabolism and preadipocyte conversion to form lipid-assimilating adipocytes (i.e. adipogenesis). However, no miRNA to date has been reported to modulate adipogenesis and lipid deposition in beef cattle. RESULTS: The expression patterns of 89 miRNAs including four bovine specific miRNAs in subcutaneous adipose tissues from three groups of crossbred steers differing in backfat thickness were compared using qRT-PCR analysis. Eighty-six miRNAs were detectable in all samples, with 42 miRNAs differing among crossbreds (P < 0.05) and 15 miRNAs differentially expressed between tissues with high and low backfat thickness (P < 0.05). The expression levels of 18 miRNAs were correlated with backfat thickness (P < 0.05). The miRNA most differentially expressed and the most strongly associated with backfat thickness was miR-378, with a 1.99-fold increase in high backfat thickness tissues (r = 0.72). CONCLUSIONS: MiRNA expression patterns differed significantly in response to host genetic components. Approximately 20% of the miRNAs in this study were identified as being correlated with backfat thickness. This result suggests that miRNAs may play a regulatory role in white adipose tissue development in beef animals.

Characterization of bovine miRNAs by sequencing and bioinformatics analysis.

BACKGROUND: MicroRNAs (miRNAs) are a family of approximately 22 nucleotide small RNA molecules which regulate gene expression by fully or partially binding to their complementary sequences in mRNAs or promoters. A large number of miRNAs and their expression patterns have been reported in human, mouse and rat. However, miRNAs and their expression patterns in live stock species such as beef cattle are not well studied. RESULTS: We constructed and sequenced small-RNA libraries to yield a total of 13,541 small-RNA sequences from 11 bovine tissues including brain, subcutaneous fat, muscle, liver, kidney, spleen and thymus. In total, 228 miRNAs including 29 novel miRNA candidates were identified. Of the 199 miRNAs, 101 have been previously reported as bovine miRNAs and the other 98 are bovine orthologs of known miRNAs that have been identified in at least one other mammalian species. Of the 29 novel miRNA candidates, 17 appeared at this point in time to be bovine specific, while the remaining 12 had evidence of evolutionary conservation in other mammalian species. Five miRNAs (miR-23a, -23b, -99a, -125b and -126-5p) were very abundant across the 11 tissues, accounting for 44.3% of all small RNA sequences. The expression analysis of selected miRNAs using qRT-PCR also showed that miR-26a and -99a were highly expressed in all tissues, while miR-122 and miR-133a were predominantly expressed in liver and muscle, respectively. CONCLUSION: The miRNA expression patterns among 11 tissues from beef cattle revealed that most miRNAs were ubiquitously expressed in all tissues, while only a few miRNAs were tissue specific. Only 60% miRNAs in this study were found to display strand bias, suggesting that there are some key factors for mature miRNA selection other than internal stability. Most bovine miRNAs are highly conserved in other three mammalian species, indicating that these miRNAs may have a role in different species that are potential molecular markers for evolution.

Genome sequence variation analysis of two SARS coronavirus isolates after passage in Vero cell culture.

SARS coronavirus is an RNA virus whose replication is error-prone, which provides possibility for escape of host defenses, and even leads to evolution of new viral strains during the passage or the transmission. Lots of variations have been detected among different SARS-CoV strains. And a study on these variations is helpful for development of efficient vaccine. Moreover, the test of nucleic acid characterization and genetic stability of SARS-CoV is important in the research of inactivated vaccine. The whole genome sequences of two SARS coronavirus strains after passage in Vero cell culture were determined and were compared with those of early passages, respectively. Results showed that both SARS coronavirus strains have high genetic stability, although nearly 10 generations were passed. Four nucleotide variations were observed between the second passage and the 11th passage of Sino1 strain for identification of SARS inactivated vaccine. Moreover, only one nucleotide was different between the third passage and the 10th passage of Sino3 strain for SARS inactivated vaccine. Therefore, this study suggested it was possible to develop inactivated vaccine against SARS-CoV in the future.

Genome sequencing and characterization analysis of a Beijing isolate of chicken coronavirus infectious bronchitis virus.

Avian infectious bronchitis virus (AIBV) is classified as a member of the genus coronavirus in the family coronaviridae. The enveloped virus has a positive-sense, single-stranded RNA genome of approximately 28 kilo-bases, which has a 5' cap structure and 3' polyadenylation tract. The complete genome sequence of infectious bronchitis virus (IBV), Beijing isolate, was determined by cloning sequencing and primer walking. The whole genome is 27733 nucleotides in length, has ten open reading frames: 5' -orf1a-orf1ab-s-3a-3b-e-m-6a-6b-n-3'. Alignments of the genome sequence of IBV Beijing isolate with those of two AIBV strains and one SARS coronavirus were performed respectively. The genome sequence of IBV Beijing isolate compared with that of the IBV strain LX4 (uncompleted, 19440 bp in size) was 91.2% similarity. However, the full-length genome sequence of IBV Beijing isolate was 85.2% identity to that of IBV Strain Beaudette, and was only 50.8% homology to that of SARS coronavirus. The results showed that the genome of IBV has remarkable variation. And IBV Beijing isolate is not closely related to SARS coronavirus. Phylogenetic analyses based on the whole genome sequence, S protein, M protein and N protein, also showed that AIBV Beijing isolate is lone virus in group III and is distant from SARS coronavirus. In conclusion, this study will contribute to the studies of diagnosis and diseases control on IBV in China.