Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology

BACKGROUND

The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues.

RESULTS

A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network.

CONCLUSIONS

These validated results show significant improvement over current bovine genome annotations.

Alternative Transcripts Diversify Genome Function for Phenome Relevance to Health and Diseases

Manipulation using alternative exon splicing (AES), alternative transcription start (ATS), and alternative polyadenylation (APA) sites are key to transcript diversity underlying health and disease. All three are pervasive in organisms, present in at least 50% of human protein-coding genes. In fact, ATS and APA site use has the highest impact on protein identity, with their ability to alter which first and last exons are utilized as well as impacting stability and translation efficiency. These RNA variants have been shown to be highly specific, both in tissue type and stage, with demonstrated importance to cell proliferation, differentiation and the transition from fetal to adult cells. While alternative exon splicing has a limited effect on protein identity, its ubiquity highlights the importance of these minor alterations, which can alter other features such as localization. The three processes are also highly interwoven, with overlapping, complementary, and competing factors, RNA polymerase II and its CTD (C-terminal domain) chief among them. Their role in development means dysregulation leads to a wide variety of disorders and cancers, with some forms of disease disproportionately affected by specific mechanisms (AES, ATS, or APA). Challenges associated with the genome-wide profiling of RNA variants and their potential solutions are also discussed in this review.

Imprinted Genes: Genomic Conservation, Transcriptomic Dynamics and Phenomic Significance in Health and Diseases

Since its discovery in 1991, genomic imprinting has been the subject of numerous studies into its mechanisms of establishment and regulation, evolution and function, and presence in multiple genomes. Disturbance of imprinting has been implicated in a range of diseases, ranging from debilitating syndromes to cancers to fetal deficiencies. Despite this, studies done on the prevalence and relevance of imprinting on genes have been limited in scope, tissue types available, and focus, by both availability and resources. This has left a gap in comparative studies. To address this, we assembled a collection of imprinted genes available in current literature covering five species. Here we sought to identify trends and motifs in the imprinted gene set (IGS) in three distinct arenas: evolutionary conservation, across-tissue expression, and health phenomics. Overall, we found that imprinted genes displayed less conservation and higher proportions of non-coding RNA while maintaining synteny. Maternally expressed genes (MEGs) and paternally expressed genes (PEGs) occupied distinct roles in tissue expression and biological pathway use, while imprinted genes collectively showed a broader tissue range, notable preference for tissue specific expression and limited gene pathways than comparable sex differentiation genes. Both human and murine imprinted genes showed the same clear phenotypic trends, that were distinct from those displayed by sex differentiation genes which were less involved in mental and nervous system disease. While both sets had representation across the genome, the IGS showed clearer clustering as expected, with PEGs significantly more represented than MEGs.

Integrated time-series transcriptomic and metabolomic analyses reveal different inflammatory and adaptive immune responses contributing to host resistance to PRRSV

Porcine reproductive and respiratory syndrome virus (PRRSV) is a highly contagious disease that affects the global pig industry. To understand mechanisms of susceptibility/resistance to PRRSV, this study profiled the time-serial white blood cells transcriptomic and serum metabolomic responses to PRRSV in piglets from a crossbred population of PRRSV-resistant Tongcheng pigs and PRRSV-susceptible Large White pigs. Gene set enrichment analysis (GSEA) illustrated that PRRSV infection up-regulated the expression levels of marker genes of dendritic cells, monocytes and neutrophils and inflammatory response, but down-regulated T cells, B cells and NK cells markers. CIBERSORT analysis confirmed the higher T cells proportion in resistant pigs during PRRSV infection. Resistant pigs showed a significantly higher level of T cell activation and lower expression levels of monocyte surface signatures post infection than susceptible pigs, corresponding to more severe suppression of T cell immunity and inflammatory response in susceptible pigs. Differentially expressed genes between resistant/susceptible pigs during the course of infection were significantly enriched in oxidative stress, innate immunity and humoral immunity, cell cycle, biotic stimulated cellular response, wounding response and behavior related pathways. Fourteen of these genes were distributed in 5 different QTL regions associated with PRRSV-related traits. Chemokine CXCL10 levels post PRRSV infection were differentially expressed between resistant pigs and susceptible pigs and can be a promising marker for susceptibility/resistance to PRRSV. Furthermore, the metabolomics dataset indicated differences in amino acid pathways and lipid metabolism between pre-infection/post-infection and resistant/susceptible pigs. The majority of metabolites levels were also down-regulated after PRRSV infection and were significantly positively correlated to the expression levels of marker genes in adaptive immune response. The integration of transcriptome and metabolome revealed concerted molecular events triggered by the infection, notably involving inflammatory response, adaptive immunity and G protein-coupled receptor downstream signaling. This study has increased our knowledge of the immune response differences induced by PRRSV infection and susceptibility differences at the transcriptomic and metabolomic levels, providing the basis for the PRRSV resistance mechanism and effective PRRS control.

Long-term trends of microplastics in seawater and farmed oysters in the Maowei Sea, China

Microplastic pollution in marine environments and organisms has received a great deal of international attention. However, the long-term field studies of microplastics are rare. Here, we assessed annual variation in microplastic abundance in the Maowei Sea, a classic mariculture bay in southern China, and analyzed the long-term accumulation in oyster tissues. U-shaped time trends of microplastics in water were observed from January to December in 2018 in the estuarine region, inner bay, and mouth bay sites, representing an inverse relationship with the local rainfall patterns. The common microplastic particles in Maowei Sea are PET/PE fibers, and polystyrene foams, which are mainly related to textile pollution and fishery activities. After one year of continuous monitoring, we did not find accumulation of microplastics in the whole soft tissues of oyster after 10% KOH digestion. No significant correlation of microplastic abundances between water and oysters was observed. The microplastic abundance in oyster was correlated with some environmental variables (i.e. salinity, pH, nutrients and total organic carbon) of the surrounding water following Spearman correlation analysis. The microplastic levels in oysters could probably be influenced by the environmental variables.

Genome-to-phenome research in rats: progress and perspectives

Because of their relatively short lifespan (<4 years), rats have become the second most used model organism to study health and diseases in humans who may live for up to 120 years. First-, second- and third-generation sequencing technologies and platforms have produced increasingly greater sequencing depth and accurate reads, leading to significant advancements in the rat genome assembly during the last 20 years. In fact, whole genome sequencing (WGS) of 47 strains have been completed. This has led to the discovery of genome variants in rats, which have been widely used to detect quantitative trait loci underlying complex phenotypes based on gene, haplotype, and sweep association analyses. DNA variants can also reveal strain, chromosome and gene functional evolutions. In parallel, phenome programs have advanced significantly in rats during the last 15 years and more than 10 databases host genome and/or phenome information. In order to discover the bridges between genome and phenome, systems genetics and integrative genomics approaches have been developed. On the other hand, multiple level information transfers from genome to phenome are executed by differential usage of alternative transcriptional start (ATS) and polyadenylation (APA) sites per gene. We used our own experiments to demonstrate how alternative transcriptome analysis can lead to enrichment of phenome-related causal pathways in rats. Development of advanced genome-to-phenome assays will certainly enhance rats as models for human biomedical research.

Genome-Wide Characterization and Comparative Analyses of Simple Sequence Repeats among Four Miniature Pig Breeds

Simple Summary

Simple sequence repeats (SSRs) are present at high densities in regulatory elements, suggesting that they may affect gene function and phenotypic traits. Therefore, SSRs can be exploited in marker-assisted selection. In addition, they can be widely used as molecular markers to study genetic diversity, population structure, and evolution. While SSRs have been widely studied in many mammalian species, very little research has focused on genome-wide SSRs of miniature pigs, a small but special group of pigs that express the dwarf phenotype. Based on the SSR-enriched library building and sequencing, about 30,000 novel polymorphic SSRs for four miniature pig breeds were mapped to the Duroc pig reference genome. The four miniature pig breeds had different numbers and types of SSRs and distributions of repeat units. There were 2518 polymorphic SSRs in the intron or exon regions that were common to all four breeds and functional analyses revealed 17 genes that were associated with body size and other genes that were associated with growth and development. In conclusion, the SSRs detected in the miniature pigs in this study may provide useful genetic markers for the selection of farm animals and the polymorphic SSRs provide valuable insights into the determination of mature body size, as well as the immunity, growth and development of animals.

Abstract

Simple sequence repeats (SSRs) are commonly used as molecular markers in research on genetic diversity and discrimination among taxa or breeds because polymorphisms in these regions contribute to gene function and phenotypically important traits. In this study, we investigated genome-wide characteristics, repeat units, and polymorphisms of SSRs using sequencing data from SSR-enriched libraries created from Wuzhishan (WZS), Bama (BM), inbred Luchuan (LC) and Zangxiang (ZX) miniature pig breeds. The numbers and types of SSRs, distributions of repeat units and polymorphic SSRs varied among the four breeds. Compared to the Duroc pig reference genome, 2518 polymorphic SSRs were unique and common to all four breeds and functional annotation revealed that they may affect the coding and regulatory regions of genes. Several examples, such as FGF23, MYF6, IGF1R, and LEPROT, are associated with growth and development in pigs. Three of the polymorphic SSRs were selected to confirm the polymorphism and the corresponding alleles through fluorescence polymerase chain reaction (PCR) and capillary electrophoresis. Together, this study provides useful insights into the discovery, characteristics and distribution of SSRs in four pig breeds. The polymorphic SSRs, especially those common and unique to all four pig breeds, might affect associated genes and play important roles in growth and development.

A novel quantitative real-time PCR method for the detection of mammalian and poultry species based on a shared single-copy nuclear DNA sequence

Accurate quantification of species fractions is critical to determine meat adulteration. This study aimed to develop a novel quantitative real-time PCR (qRT-PCR) method for detection of mammalian and poultry DNA. A shared single-copy nuclear DNA sequence derived from the first exon of the LcoR gene was identified as a multi-species universal reference for a qRT-PCR assay. The conservation and copy number of the LcoR gene were evaluated among different species. The limit of detection was 0.01 ng DNA or 0.01% meat ingredient, and the limit of quantification was 0.01 ng DNA or 0.05% meat ingredient. Both the relative error (R.E.) and relative standard deviation (R.S.D.) were ≤ 25%. Moreover, modified coefficient k was introduced into this quantitative system to improve the accuracy and reliability of results, with maximum R.E. improved from 19.43% to 16.16%. The quantitative method would contribute to fighting against meat adulteration and maintaining a fair market.

341 Up-regulation of wound healing pathway may trigger adipogenic potentials of intramuscular progenitor cells in Wagyu as compared to Angus cattle

Marbling and muscle growth are significantly different between Angus and Wagyu cattle, but the underlying mechanisms remain unclear. The objective of the present study was to unravel genetic information flows specific to each phenotype. Biopsies of the Biceps femoris muscle were collected from three animals per breed. Intramuscular progenitor cells were isolated from the biopsies and alternative polyadenylation (APA) events were profiled using the WTTS-seq (whole transcriptome termini site sequencing) method. In adipogenic progenitor cells, 271 and 341 differentially expressed (DE) APA sites were identified as upregulated in Angus and Wagyu, respectively (adjusted P &lt; 0.1). Pathway enrichment revealed that cellular wound healing events are unique to Wagyu. For these seven pairs of gene family members (in a broad sense), AKT1, CD47 and CD99, IGFBP3, KIF3A, RIOK1, RUNX1 and THBS1 were upregulated in Angus, while AKT3, CD59 and CD200, IGFBP2 and IGFBP4, KIF3B, RIOK3, RUNX2 and THBS2 were upregulated in Wagyu. There were 313 and 171 upregulated DE-APA sites in satellite cells derived from Angus and Wagyu, respectively (adjusted P &lt; 0.1). Muscle system related pathways were enriched specifically for Angus. The CAPZA1, FAM227B, LDB2 and TTLL4 genes contained at least one DE-APA site that was specifically abundant in Angus or Wagyu. For these six gene family member pairs, ANKRD1, CD82, FBXO17, IGFBP2, PDLIM3 and THBS1 were abundantly expressed in Angus, but ANKRD28, CD9, CD14 and CD164, FBXO33, IGFBP5, PDLIM1 and THBS2 were upregulated in Wagyu. In brief, intramuscular progenitor cells and satellite cells may be characteristically different between Angus and Wagyu. Based on DE-APA events, the myogenesis pathways are enriched in cells from Angus, but the adipogenesis pathways are dominant in cells derived from Wagyu. This research was supported by the USDA-NIFA under Award Numbers 2016-67015-24470/2018-67051-27500 to ZJ and 2015-67015-23219/2016-68006-24634 to MD.

PSXI-4 Composting management to reduce greenhouse gas and ammonia emissions from feedlot manure

Management of livestock manure may recycle nutrients and decrease greenhouse gas (GHG) and ammonia (NH3) emissions. The objectives were to ascertain effects of environmental conditions and turning on methane (CH4), nitrous oxide (N2O), and NH3 emissions and if treatment with 8.5 g of dicyandiamide (DCD), a denitrification agent, altered GHG emissions. Manure and bedding were collected from feedlot pens and used to construct 3 piles (~1.9 m3 volume) each in winter (WI) and spring (SP). WI piles were turned once, and SP piles were turned twice. Methane, N2O, and NH3 emissions were collected. Methane and N2O flux measurements were collected from SP piles using a static chamber (3.7m L x 2.2m W x 0.9m H). Initial dry matter and nitrogen contents were 33.2 and 30.0% and 20.1 and 17.7 g/kg in WI and SP piles, respectively. Average ambient temperatures and wind speeds were 0.3oC and 10.7oC and 1.76 m/s and 1.97 m/s during WI and SP, respectively. Internal temperatures reached 51±3.9oC on d 4–11 and gradually decreased. Normalized CH4 averaged 2.19 mg٠s٠m-4 and N2O emissions averaged 0.84 mg٠s٠m-4, and were not different between the WI and SP piles. Turning did not affect CH4 emissions from WI piles, but were 55% greater (P &lt; 0.05) when SP piles were turned a second time. Emissions of N2O increased 51% when WI and SP piles were turned (P &lt; 0.05). Ammonia emissions were 83.5% greater from WI piles due to their higher initial concentrations of NH4+-N (2.21 vs. 1.11 g/kg; P &lt; 0.05). Turning did not influence CH4 and N2O fluxes. Addition of DCD at pile formation appears to decrease N2O emissions and fluxes 3 and 10 d later. Turning management and season impacted overall CH4, N2O, and NH3 emissions. Fine-tuning manure handling and management during different seasons may effectively reduce GHG and NH3 emissions.

Genome-wide association study shows that microtia in Altay sheep is caused by a 76 bp duplication of HMX1

Microtia is a congenital malformation of the external ear that can be observed in many species including sheep. However, the genetic basis of microtia still remains unclear. Here, a GWAS was conducted to investigate the genetic basis underlying microtia. A total of 55 samples from 26 microtia and 29 normal animals were genotyped with Illumina OvineHD BeadChip. The strongest significant SNP was identified on OAR6, approximating the evolutionarily conserved region of the HMX1 gene, which is related to congenital malformations of the external ear in other species such as cattle and rats. Sequencing an evolutionarily conserved region surrounding HMX1 revealed a duplication of 76 bp, which is concordant with microtia, suggesting a dominant inheritance mode. Identification of this causal mutation in the HMX1 gene indicates the role of this particular gene in the development of the external ear and provides a genetic marker for selection against microtia.

Alternative polyadenylation coordinates embryonic development, sexual dimorphism and longitudinal growth in Xenopus tropicalis

RNA alternative polyadenylation contributes to the complexity of information transfer from genome to phenome, thus amplifying gene function. Here, we report the first X. tropicalis resource with 127,914 alternative polyadenylation (APA) sites derived from embryos and adults. Overall, APA networks play central roles in coordinating the maternal-zygotic transition (MZT) in embryos, sexual dimorphism in adults and longitudinal growth from embryos to adults. APA sites coordinate reprogramming in embryos before the MZT, but developmental events after the MZT due to zygotic genome activation. The APA transcriptomes of young adults are more variable than growing adults and male frog APA transcriptomes are more divergent than females. The APA profiles of young females were similar to embryos before the MZT. Enriched pathways in developing embryos were distinct across the MZT and noticeably segregated from adults. Briefly, our results suggest that the minimal functional units in genomes are alternative transcripts as opposed to genes.

Alternative polyadenylation analysis in animals and plants: newly developed strategies for profiling, processing and validation

Alternative polyadenylation is an essential RNA processing event that contributes significantly to regulation of transcriptome diversity and functional dynamics in both animals and plants. Here we review newly developed next generation sequencing methods for genome-wide profiling of alternative polyadenylation (APA) sites, bioinformatics pipelines for data processing and both wet and dry laboratory approaches for APA validation. The library construction methods LITE-Seq (Low-Input 3'-Terminal sequencing) and PAC-seq (PolyA Click sequencing) tag polyA+ cDNA, while BAT-seq (BArcoded, three-prime specific sequencing) and PAPERCLIP (Poly(A) binding Protein-mediated mRNA 3'End Retrieval by CrossLinking ImmunoPrecipitation) enrich polyA+ RNA. Interestingly, only WTTS-seq (Whole Transcriptome Termini Site sequencing) targets both polyA+ RNA and polyA+ cDNA. Varieties of bioinformatics pipelines are well established to pursue read quality control, mapping, clustering, characterization and pathway analysis. The RHAPA (RNase H alternative polyadenylation assay) and 3'RACE-seq (3' rapid amplification of cDNA end sequencing) methods directly validate APA sites, while WTSS-seq (whole transcriptome start site sequencing), RNA-seq (RNA sequencing) and public APA databases can serve as indirect validation methods. We hope that these tools, pipelines and resources trigger huge waves of interest in the research community to investigate APA events underlying physiological, pathological and psychological changes and thus understand the information transfer events from genome to phenome relevant to economically important traits in both animals and plants.

The role of pparγ in embryonic development of Xenopus tropicalis under triphenyltin-induced teratogenicity

Evidence has shown that triphenyltin (TPT) triggers severe malformations in Xenopus tropicalis embryos, partly due to activation of PPARγ (peroxisome proliferator activated receptor γ) protein. In the present study, we investigated how abundance of pparγ and TPT exposure interact and affect X. tropicalis embryonic development. We observed pparγ expression signals appeared in the neural crest and neural fold, as well as in the brain, eyes and spinal cord organs. Both pparγ overexpression and its Morpholino (MO) knockdown inhibited pax6 (paired box 6) expression, a marker of eye development, and significantly up- and down-regulated lipid and glucose homeostasis related genes, such as lpl (lipoprotein lipase), slc2a4 (solute carrier family 2 (facilitated glucose transporter), member 4) and pck1 (phosphoenolpyruvate carboxykinase 1, cytosolic), thus inducing eye phenotypes. Overexpression of pparγ induced small eye phenotype, while pparγ MO induced small eye plus turbid eye lens microencephaly and enlarged trunk. In contrast, 5-20μgSn/L (stannum/L) TPT exposure reversed some impacts induced by pparγ overexpression, i.e., no small eye, up-regulation of pax6 and down-regulation of pparγ, lpl, slc2a4 and pck1. Meanwhile, microinjection of pparγ MO combined with exposure to 20μgSn/L TPT caused 85% mortality. In brief, our work clearly indicates that pparγ is essential to eye development and inhibition of its expression combined with TPT exposure can be extremely harmful to X. tropicalis embryo.

Genome Wide Sampling Sequencing for SNP Genotyping: Methods, Challenges and Future Development

Genetic polymorphisms, particularly single nucleotide polymorphisms (SNPs), have been widely used to advance quantitative, functional and evolutionary genomics. Ideally, all genetic variants among individuals should be discovered when next generation sequencing (NGS) technologies and platforms are used for whole genome sequencing or resequencing. In order to improve the cost-effectiveness of the process, however, the research community has mainly focused on developing genome-wide sampling sequencing (GWSS) methods, a collection of reduced genome complexity sequencing, reduced genome representation sequencing and selective genome target sequencing. Here we review the major steps involved in library preparation, the types of adapters used for ligation and the primers designed for amplification of ligated products for sequencing. Unfortunately, currently available GWSS methods have their drawbacks, such as inconsistency in the number of reads per sample library, the number of sites/targets per individual, and the number of reads per site/target, all of which result in missing data. Suggestions are proposed here to improve library construction, genotype calling accuracy, genome-wide marker density and read mapping rate. In brief, optimized GWSS library preparation should generate a unique set of target sites with dense distribution along chromosomes and even coverage per site across all individuals.

Genome Wide Distributions and Functional Characterization of Copy Number Variations between Chinese and Western Pigs

Copy number variations (CNVs) refer to large insertions, deletions and duplications in the genomic structure ranging from one thousand to several million bases in size. Since the development of next generation sequencing technology, several methods have been well built for detection of copy number variations with high credibility and accuracy. Evidence has shown that CNV occurring in gene region could lead to phenotypic changes due to the alteration in gene structure and dosage. However, it still remains unexplored whether CNVs underlie the phenotypic differences between Chinese and Western domestic pigs. Based on the read-depth methods, we investigated copy number variations using 49 individuals derived from both Chinese and Western pig breeds. A total of 3,131 copy number variation regions (CNVRs) were identified with an average size of 13.4 Kb in all individuals during domestication, harboring 1,363 genes. Among them, 129 and 147 CNVRs were Chinese and Western pig specific, respectively. Gene functional enrichments revealed that these CNVRs contribute to strong disease resistance and high prolificacy in Chinese domestic pigs, but strong muscle tissue development in Western domestic pigs. This finding is strongly consistent with the morphologic characteristics of Chinese and Western pigs, indicating that these group-specific CNVRs might have been preserved by artificial selection for the favored phenotypes during independent domestication of Chinese and Western pigs. In this study, we built high-resolution CNV maps in several domestic pig breeds and discovered the group specific CNVs by comparing Chinese and Western pigs, which could provide new insight into genomic variations during pigs’ independent domestication, and facilitate further functional studies of CNV-associated genes.

Molecular characterization of the porcine S100A6 gene and analysis of its expression in pigs infected with highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV)

Our previous microarray study revealed that S100A6 was significantly upregulated in porcine alveolar macrophages (PAMs) infected with highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV). In the present study, we cloned both cDNA and genomic DNA sequences of the gene. Transient transfection indicated that the porcine S100A6 protein was located in the nucleus and cytoplasm. Reverse transcription polymerase chain reaction (RT-PCR) revealed that the porcine S100A6 gene was highly expressed in the kidney and subcutaneous fat. Polyinosinic-polycytidylic acid [poly (I:C)] induced porcine S100A6 gene expression in PK-15 cells. Quantitative real-time PCR (Q-PCR) analysis further showed that the porcine S100A6 gene was upregulated in different cells and tissues of Tongcheng pigs infected with HP-PRRSV. Chromosome walking obtained the porcine S100A6 promoter region and then luciferase reporter assays confirmed its regulatory activities. We observed a putative NF-κB binding site in the core promoter region, which may explain the upregulation of porcine S100A6 in response to PRRSV. Transfection of NF-κB (p65 subunit) intensely induced the promoter activity of the porcine S100A6 gene, while an NF-κB inhibitor, pyrrolidine dithiocarbamate (PDTC), inhibited this activity. Furthermore, compared to its wild type, the promoter activity was significantly reduced when it contained a mutant NF-κB binding site. All these results provide a solid foundation to further investigate how S100A6 is involved in PRRSV infection.

Relationships between residual feed intake and hepatic mitochondrial function in growing beef cattle

The objective of this study was to evaluate the relationship between hepatic mitochondrial function and residual feed intake (RFI) in growing beef cattle. In Trial 1, RFI was measured in 29 Angus heifers (initial BW = 258.0 ± 24.9 kg) from divergent IGF-I selection lines created at the Eastern Agricultural Research Station (The Ohio State University) fed a grain-based diet (calculated ME = 2.85 Mcal/kg DM). In Trial 2, RFI was measured in 119 Santa Gertrudis steers (initial BW = 308.4 ± 28.1 kg) fed a roughage-based diet (calculated ME = 2.21 Mcal/kg DM). At the end of the RFI measurement period, cattle in Trial 1 (n = 7 low RFI and n = 7 high RFI) and in Trial 2 (n = 6 low RFI and n = 8 high RFI) with measures of RFI exceeding 0.5 (Trial 1) or 1.0 (Trial 2) SD from the mean RFI were selected to measure mitochondrial function. Overall ADG, DMI, and RFI were 1.19 ± 0.15, 9.31 ± 1.12, and 0.00 ± 0.63 kg/d and 0.83 ± 0.16, 9.48 ± 1.00, and 0.00 ± 0.86 kg/d in Trial 1 and 2, respectively. Cattle with low RFI consumed 13 and 24% less (P < 0.05) DM and had 14 and 56% greater (P < 0.05) G:F than cattle with high RFI in Trial 1 and 2, respectively, even though ADG and BW were similar (P > 0.10). In Trial 1, cattle with low RFI tended (P = 0.06) to have greater state 3 respiration rates than cattle with high RFI, but state 3 respiration rates were similar (P > 0.10) between cattle with low and high RFI in Trial 2. In both trials, cattle with low RFI had greater (P < 0.05) acceptor control ratios than their high RFI counterparts. The respiratory control ratio tended (P = 0.09) to be greater for cattle with low RFI compared with high RFI cattle in Trial 1, but no difference (P > 0.10) was observed in Trial 2. Proton-leak kinetics were similar (P > 0.05) between cattle with low and high RFI in both trials. These data suggest that ADP has greater control of oxidative phosphorylation in liver mitochondrial of cattle with low RFI compared to their high RFI counterparts.

Genome wide screening of candidate genes for improving piglet birth weight using high and low estimated breeding value populations

Birth weight is an economically important trait in pig production because it directly impacts piglet growth and survival rate. In the present study, we performed a genome wide survey of candidate genes and pathways associated with individual birth weight (IBW) using the Illumina PorcineSNP60 BeadChip on 24 high (HEBV) and 24 low estimated breeding value (LEBV) animals. These animals were selected from a reference population of 522 individuals produced by three sires and six dam lines, which were crossbreds with multiple breeds. After quality-control, 43,257 SNPs (single nucleotide polymorphisms), including 42,243 autosomal SNPs and 1,014 SNPs on chromosome X, were used in the data analysis. A total of 27 differentially selected regions (DSRs), including 1 on Sus scrofa chromosome 1 (SSC1), 1 on SSC4, 2 on SSC5, 4 on SSC6, 2 on SSC7, 5 on SSC8, 3 on SSC9, 1 on SSC14, 3 on SSC18, and 5 on SSCX, were identified to show the genome wide separations between the HEBV and LEBV groups for IBW in piglets. A DSR with the most number of significant SNPs (including 7 top 0.1% and 31 top 5% SNPs) was located on SSC6, while another DSR with the largest genetic differences in F_ST was found on SSC18. These regions harbor known functionally important genes involved in growth and development, such as TNFRSF9 (tumor necrosis factor receptor superfamily member 9), CA6 (carbonic anhydrase VI) and MDFIC (MyoD family inhibitor domain containing). A DSR rich in imprinting genes appeared on SSC9, which included PEG10 (paternally expressed 10), SGCE (sarcoglycan, epsilon), PPP1R9A (protein phosphatase 1, regulatory subunit 9A) and ASB4 (ankyrin repeat and SOCS box containing 4). More importantly, our present study provided evidence to support six quantitative trait loci (QTL) regions for pig birth weight, six QTL regions for average birth weight (ABW) and three QTL regions for litter birth weight (LBW) reported previously by other groups. Furthermore, gene ontology analysis with 183 genes harbored in these 27 DSRs suggested that protein, metal, ion and ATP binding, viral process and innate immune response present important pathways for deciphering their roles in fetal growth or development. Overall, our study provides useful information on candidate genes and pathways for regulating birth weight in piglets, thus improving our understanding of the genetic mechanisms involved in porcine embryonic or fetal development.

Molecular heterogeneities of adipose depots - potential effects on adipose-muscle cross-talk in humans, mice and farm animals

Adipose tissue is considered as a major endocrine organ that secretes numerous proteins called adipokines. The heterogeneous nature of adipose tissue in different parts of the body suggests respective heterogeneity of proteomes and secretomes. This review consolidates knowledge from recent studies targeting the diversity of different adipose depots affecting the pattern of secreted adipokines and discusses potential consequences for the cross-talk between adipose and skeletal muscle in humans, rodent models and farm animals. Special attention is paid to muscle-associated fat depots like inter- and intramuscular fat that become focus of attention in the context of the rather new notion of skeletal muscle as a major endocrine organ. Understanding the complexity of communication between adipocytes and skeletal muscle cells will allow developing strategies for improvement of human health and for sustainable production of high quality meat.

Obesity gene atlas in mammals

Obesity in humans has increased at an alarming rate over the past two decades and has become one of the leading public health problems worldwide. Studies have revealed a large number of genes/markers that are associated with obesity and/or obesity-related phenotypes, indicating an urgent need to develop a central database for helping the community understand the genetic complexity of obesity. In the present study, we collected a total of 1,736 obesity associated loci and created a freely available obesity database, including 1,515 protein-coding genes and 221 microRNAs (miRNAs) collected from four mammalian species: human, cattle, rat, and mouse. These loci were integrated as orthologs on comparative genomic views in human, cattle, and mouse. The database and genomic views are freely available online at: http://www.integratomics-time.com/fat_deposition. Bioinformatics analyses of the collected data revealed some potential novel obesity related molecular markers which represent focal points for testing more targeted hypotheses and designing experiments for further studies. We believe that this centralized database on obesity and adipogenesis will facilitate development of comparative systems biology approaches to address this important health issue in human and their potential applications in animals.

Genome-wide genetic diversity and differentially selected regions among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep

Sheep are among the major economically important livestock species worldwide because the animals produce milk, wool, skin, and meat. In the present study, the Illumina OvineSNP50 BeadChip was used to investigate genetic diversity and genome selection among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds from the United States. After quality-control filtering of SNPs (single nucleotide polymorphisms), we used 48,026 SNPs, including 46,850 SNPs on autosomes that were in Hardy-Weinberg equilibrium and 1,176 SNPs on chromosome × for analysis. Phylogenetic analysis based on all 46,850 SNPs clearly separated Suffolk from Rambouillet, Columbia, Polypay, and Targhee, which was not surprising as Rambouillet contributed to the synthesis of the later three breeds. Based on pair-wise estimates of F_ST, significant genetic differentiation appeared between Suffolk and Rambouillet (F_ST = 0.1621), while Rambouillet and Targhee had the closest relationship (F_ST = 0.0681). A scan of the genome revealed 45 and 41 differentially selected regions (DSRs) between Suffolk and Rambouillet and among Rambouillet-related breed populations, respectively. Our data indicated that regions 13 and 24 between Suffolk and Rambouillet might be good candidates for evaluating breed differences. Furthermore, ovine genome v3.1 assembly was used as reference to link functionally known homologous genes to economically important traits covered by these differentially selected regions. In brief, our present study provides a comprehensive genome-wide view on within- and between-breed genetic differentiation, biodiversity, and evolution among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds. These results may provide new guidance for the synthesis of new breeds with different breeding objectives.

Reactomes of porcine alveolar macrophages infected with porcine reproductive and respiratory syndrome virus

Porcine reproductive and respiratory syndrome (PRRS) has devastated pig industries worldwide for many years. It is caused by a small RNA virus (PRRSV), which targets almost exclusively pig monocytes or macrophages. In the present study, five SAGE (serial analysis of gene expression) libraries derived from 0 hour mock-infected and 6, 12, 16 and 24 hours PRRSV-infected porcine alveolar macrophages (PAMs) produced a total 643,255 sequenced tags with 91,807 unique tags. Differentially expressed (DE) tags were then detected using the Bayesian framework followed by gene/mRNA assignment, arbitrary selection and manual annotation, which determined 699 DE genes for reactome analysis. The DAVID, KEGG and REACTOME databases assigned 573 of the DE genes into six biological systems, 60 functional categories and 504 pathways. The six systems are: cellular processes, genetic information processing, environmental information processing, metabolism, organismal systems and human diseases as defined by KEGG with modification. Self-organizing map (SOM) analysis further grouped these 699 DE genes into ten clusters, reflecting their expression trends along these five time points. Based on the number one functional category in each system, cell growth and death, transcription processes, signal transductions, energy metabolism, immune system and infectious diseases formed the major reactomes of PAMs responding to PRRSV infection. Our investigation also focused on dominant pathways that had at least 20 DE genes identified, multi-pathway genes that were involved in 10 or more pathways and exclusively-expressed genes that were included in one system. Overall, our present study reported a large set of DE genes, compiled a comprehensive coverage of pathways, and revealed system-based reactomes of PAMs infected with PRRSV. We believe that our reactome data provides new insight into molecular mechanisms involved in host genetic complexity of antiviral activities against PRRSV and lays a strong foundation for vaccine development to control PRRS incidence in pigs.

Interferon induced IFIT family genes in host antiviral defense

Secretion of interferons (IFNs) from virus-infected cells is a hallmark of host antiviral immunity and in fact, IFNs exert their antiviral activities through the induction of antiviral proteins. The IFN-induced protein with tetratricopeptide repeats (IFITs) family is among hundreds of IFN-stimulated genes. This family contains a cluster of duplicated loci. Most mammals have IFIT1, IFIT2, IFIT3 and IFIT5; however, bird, marsupial, frog and fish have only IFIT5. Regardless of species, IFIT5 is always adjacent to SLC16A12. IFIT family genes are predominantly induced by type I and type III interferons and are regulated by the pattern recognition and the JAK-STAT signaling pathway. IFIT family proteins are involved in many processes in response to viral infection. However, some viruses can escape the antiviral functions of the IFIT family by suppressing IFIT family genes expression or methylation of 5' cap of viral molecules. In addition, the variants of IFIT family genes could significantly influence the outcome of hepatitis C virus (HCV) therapy. We believe that our current review provides a comprehensive picture for the community to understand the structure and function of IFIT family genes in response to pathogens in human, as well as in animals.

Molecular characterization of porcine SARM1 and its role in regulating TLRs signaling during highly pathogenic porcine reproductive and respiratory syndrome virus infection in vivo

Toll-like receptors (TLRs) are important pattern-recognition receptors (PRRs) that trigger innate immune response and mediate acquired immunity. Evidence has shown that SARM1 (sterile-α and TIR motif containing protein 1) is one of five TIR domain-containing adaptor proteins involved in TLRs signaling transduction. In the present study, a full-length cDNA sequence was cloned for the porcine SARM1 gene, which contains nine exons. Using the radiation hybrid mapping approach, we assigned the porcine gene to SSC12 q13. Under the normal condition, porcine SARM1 was highly expressed in brain and spleen. Polyinosinic-polycytidylic acid (poly (I:C)) weakly induced the porcine SARM1 expression in the early stimulation. We found that porcine SARM1 protein is localized in mitochondria and attenuates NF-κB activation induced by stimulation and infection. The quantitative real-time PCR (Q-PCR) analysis showed that the expression of porcine SARM1 significantly decreased in several tissues of Tongcheng pigs infected with highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV). Gene-interaction network analysis for porcine SARM1 in porcine alveolar macrophages (PAMs) showed that down-regulation of SARM1 gene in infected Tongcheng pig may modulate TRIF-depend TLRs signaling and regulate the expression of disease-resistant genes and inflammatory genes. Our findings provide evidence that porcine SARM1 may play an important role in immune regulation with PRRSV infection.

Quantitative genomics of 30 complex phenotypes in Wagyu x Angus F₁ progeny

In the present study, a total of 91 genes involved in various pathways were investigated for their associations with six carcass traits and twenty-four fatty acid composition phenotypes in a Wagyu×Angus reference population, including 43 Wagyu bulls and their potential 791 F₁ progeny. Of the 182 SNPs evaluated, 102 SNPs that were in Hardy-Weinberg equilibrium with minor allele frequencies (MAF>0.15) were selected for parentage assignment and association studies with these quantitative traits. The parentage assignment revealed that 40 of 43 Wagyu sires produced over 96.71% of the calves in the population. Linkage disequilibrium analysis identified 75 of 102 SNPs derived from 54 genes as tagged SNPs. After Bonferroni correction, single-marker analysis revealed a total of 113 significant associations between 44 genes and 29 phenotypes (adjusted P<0.05). Multiple-marker analysis confirmed single-gene associations for 10 traits, but revealed two-gene networks for 9 traits and three-gene networks for 8 traits. Particularly, we observed that TNF (tumor necrosis factor) gene is significantly associated with both beef marbling score (P=0.0016) and palmitic acid (C16:0) (P=0.0043), RCAN1 (regulator of calcineurin 1) with rib-eye area (P=0.0103), ASB3 (ankyrin repeat and SOCS box-containing 3) with backfat (P=0.0392), ABCA1 (ATP-binding cassette A1) with both palmitic acid (C16:0) (P=0.0025) and oleic acid (C18:1n9) (P=0.0114), SLC27A1(solute carrier family 27 A1) with oleic acid (C18:1n9) (P=0.0155), CRH (corticotropin releasing hormone) with both linolenic acid (OMEGA-3) (P=0.0200) and OMEGA 6:3 RATIO (P=0.0054), SLC27A2 (solute carrier family 27 A2) with both linoleic acid (OMEGA-6) (P=0.0121) and FAT (P=0.0333), GNG3 (guanine nucleotide binding protein gamma 3 with desaturase 9 (P=0.0115), and EFEMP1 (EGF containing fibulin-like extracellular matrix protein 1), PLTP (phospholipid transfer protein) and DSEL (dermatan sulfate epimerase-like) with conjugated linoleic acid (P=0.0042-0.0044), respectively, in the Wagyu x Angus F₁ population. In addition, we observed an interesting phenomenon that crossbreeding of different breeds might change gene actions to dominant and overdominant modes, thus explaining the origin of heterosis. The present study confirmed that these important families or pathway-based genes are useful targets for improving meat quality traits and healthful beef products in cattle.

Development and initial characterization of a HAPPY panel for mapping the X. tropicalis genome

HAPPY mapping was designed to pursue the analysis of approximately random HAPloid DNA breakage samples using the PolYmerase chain reaction for mapping genomes. In the present study, we improved the method and integrated two other molecular techniques into the process: whole genome amplification and the Sequenom SNP (single nucleotide polymorphism) genotyping assay in order to facilitate whole genome mapping of X. tropicalis. The former technique amplified enough DNA materials to genotype a large number of markers, while the latter allowed for relatively high throughput marker genotyping with multiplex assays on the HAPPY lines. A total of 58 X. tropicalis genes were genotyped on an initial panel of 383 HAPPY lines, which contributed to formation of a working panel of 146 lines. Further genotyping of 29 markers on the working panel led to construction of a HAPPY map for the X. tropicalis genome. We believe that our improved HAPPY method described in the present study has paved the way for the community to map different genomes with a simple, but powerful approach.

The heparan and heparin metabolism pathway is involved in regulation of fatty acid composition

Six genes involved in the heparan sulfate and heparin metabolism pathway, DSEL (dermatan sulfate epimerase-like), EXTL1 (exostoses (multiple)-like 1), HS6ST1 (heparan sulfate 6-O-sulfotransferase 1), HS6ST3 (heparan sulfate 6-O-sulfotransferase 3), NDST3 (N-deacetylase/N-sulfotransferase (heparan glucosaminyl) 3), and SULT1A1 (sulfotransferase family, cytosolic, 1A, phenol-preferring, member 1), were investigated for their associations with muscle lipid composition using cattle as a model organism. Nineteen single nucleotide polymorphisms (SNPs)/multiple nucleotide length polymorphisms (MNLPs) were identified in five of these six genes. Six of these mutations were then genotyped on 246 Wagyu x Limousin F₂ animals, which were measured for 5 carcass, 6 eating quality and 8 fatty acid composition traits. Association analysis revealed that DSEL, EXTL1 and HS6ST1 significantly affected two stearoyl-CoA desaturase activity indices, the amount of conjugated linoleic acid (CLA), and the relative amount of saturated fatty acids (SFA) and monounsaturated fatty acids (MUFA) in skeletal muscle (P<0.05). In particular, HS6ST1 joined our previously reported SCD1 and UQCRC1 genes to form a three gene network for one of the stearoyl-CoA desaturase activity indices. These results provide evidence that genes involved in heparan sulfate and heparin metabolism are also involved in regulation of lipid metabolism in bovine muscle. Whether the SNPs affected heparan sulfate proteoglycan structure is unknown and warrants further investigation.

A global view of 54,001 single nucleotide polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and their transferability to water buffalo

The Illumina BovineSNP50 BeadChip features 54,001 informative single nucleotide polymorphisms (SNPs) that uniformly span the entire bovine genome. Among them, 52,255 SNPs have locations assigned in the current genome assembly (Btau_4.0), including 19,294 (37%) intragenic SNPs (i.e., located within genes) and 32,961 (63%) intergenic SNPs (i.e., located between genes). While the SNPs represented on the Illumina Bovine50K BeadChip are evenly distributed along each bovine chromosome, there are over 14,000 genes that have no SNPs placed on the current BeadChip. Kernel density estimation, a non-parametric method, was used in the present study to identify SNP-poor and SNP-rich regions on each bovine chromosome. With bandwidth = 0.05 Mb, we observed that most regions have SNP densities within 2 standard deviations of the chromosome SNP density mean. The SNP density on chromosome X was the most dynamic, with more than 30 SNP-rich regions and at least 20 regions with no SNPs. Genotyping ten water buffalo using the Illumina BovineSNP50 BeadChip revealed that 41,870 of the 54,001 SNPs are fully scored on all ten water buffalo, but 6,771 SNPs are partially scored on one to nine animals. Both fully scored and partially/no scored SNPs are clearly clustered with various sizes on each chromosome. However, among 43,687 bovine SNPs that were successfully genotyped on nine and ten water buffalo, only 1,159 were polymorphic in the species. These results indicate that the SNPs sites, but not the polymorphisms, are conserved between two species. Overall, our present study provides a solid foundation to further characterize the SNP evolutionary process, thus improving understanding of within- and between-species biodiversity, phylogenetics and adaption to environmental changes.

The reverse cholesterol transport pathway improves understanding of genetic networks for fat deposition and muscle growth in beef cattle

In the present study, thirteen genes involved in the reverse cholesterol transport (RCT) pathway were investigated for their associations with three fat depositions, eight fatty acid compositions and two growth-related phenotypes in a Wagyu x Limousin reference population, including 6 F₁ bulls, 113 F₁ dams, and 246 F₂ progeny. A total of 37 amplicons were used to screen single nucleotide polymorphisms (SNPs) on 6 F₁ bulls. Among 36 SNPs detected in 11 of these 13 genes, 19 were selected for genotyping by the Sequenom assay design on all F₂ progeny. Single-marker analysis revealed seven SNPs in ATP binding cassette A1, apolipoproteins A1, B and E, phospholipid transfer protein and paraoxinase 1 genes significantly associated with nine phenotypes (P<0.05). Previously, we reported genetic networks associated with 19 complex phenotypes based on a total of 138 genetic polymorphisms derived from 71 known functional genes. Therefore, after Bonferroni correction, these significant (adjusted P<0.05) and suggestive (adjusted P<0.10) associations were then used to identify genetic networks related to the RCT pathway. Multiple-marker analysis suggested possible genetic networks involving the RCT pathway for kidney-pelvic-heart fat percentage, rib-eye area, and subcutaneous fat depth phenotypes with markers derived from paraoxinase 1, apolipoproteins A1 and E, respectively. The present study confirmed that genes involved in cholesterol homeostasis are useful targets for investigating obesity in humans as well as for improving meat quality phenotypes in a livestock production.

Water buffalo genome science comes of age

The water buffalo is vital to the lives of small farmers and to the economy of many countries worldwide. Not only are they draught animals, but they are also a source of meat, horns, skin and particularly the rich and precious milk that may be converted to creams, butter, yogurt and many cheeses. Genome analysis of water buffalo has advanced significantly in recent years. This review focuses on currently available genome resources in water buffalo in terms of cytogenetic characterization, whole genome mapping and next generation sequencing. No doubt, these resources indicate that genome science comes of age in the species and will provide knowledge and technologies to help optimize production potential, reproduction efficiency, product quality, nutritional value and resistance to diseases. As water buffalo and domestic cattle, both members of the Bovidae family, are closely related, the vast amount of cattle genetic/genomic resources might serve as shortcuts for the buffalo community to further advance genome science and biotechnologies in the species.

Discovery of novel genetic networks associated with 19 economically important traits in beef cattle

Quantitative or complex traits are determined by the combined effects of many loci, and are affected by genetic networks or molecular pathways. In the present study, we genotyped a total of 138 mutations, mainly single nucleotide polymorphisms derived from 71 functional genes on a Wagyu x Limousin reference population. Two hundred forty six F₂ animals were measured for 5 carcass, 6 eating quality and 8 fatty acid composition traits. A total of 2,280 single marker-trait association runs with 120 tagged mutations selected based on the HAPLOVIEW analysis revealed 144 significant associations (P < 0.05), but 50 of them were removed from the analysis due to the small number of animals (≤ 9) in one genotype group or absence of one genotype among three genotypes. The remaining 94 single-trait associations were then placed into three groups of quantitative trait modes (QTMs) with additive, dominant and overdominant effects. All significant markers and their QTMs associated with each of these 19 traits were involved in a linear regression model analysis, which confirmed single-gene associations for 4 traits, but revealed two-gene networks for 8 traits and three-gene networks for 5 traits. Such genetic networks involving both genotypes and QTMs resulted in high correlations between predicted and actual values of performance, thus providing evidence that the classical Mendelian principles of inheritance can be applied in understanding genetic complexity of complex phenotypes. Our present study also indicated that carcass, eating quality and fatty acid composition traits rarely share genetic networks. Therefore, marker-assisted selection for improvement of one category of these traits would not interfere with improvement of another.

Lipoproteins, cholesterol homeostasis and cardiac health

Cholesterol is an essential substance involved in many functions, such as maintaining cell membranes, manufacturing vitamin D on surface of the skin, producing hormones, and possibly helping cell connections in the brain. When cholesterol levels rise in the blood, they can, however, have dangerous consequences. In particular, cholesterol has generated considerable notoriety for its causative role in atherosclerosis, the leading cause of death in developed countries around the world. Homeostasis of cholesterol is centered on the metabolism of lipoproteins, which mediate transport of the lipid to and from tissues. As a synopsis of the major events and proteins that manage lipoprotein homeostasis, this review contributes to the substantial attention that has recently been directed to this area. Despite intense scrutiny, the majority of phenotypic variation in total cholesterol and related traits eludes explanation by current genetic knowledge. This is somewhat disappointing considering heritability estimates have established these traits as highly genetic. Thus, the continued search for candidate genes, mutations, and mechanisms is vital to our understanding of heart disease at the molecular level. Furthermore, as marker development continues to predict risk of vascular illness, this knowledge has the potential to revolutionize treatment of this leading human disease.

Genome assembly anchored QTL map of bovine chromosome 14

Bovine chromosome 14 (BTA14) has been widely explored for quantitative trait loci (QTL) and genes related to economically important traits in both dairy and beef cattle. We reviewed more than 40 investigations and anchored 126 QTL to the current genome assembly (Btau 4_0). Using this anchored QTL map, we observed that, in dairy cattle, the region spanning 0 – 10 Mb on BTA14 has the highest density QTL map with a total of 56 QTL, mainly for milk production traits. It is very likely that both somatic cell score (SCS) and clinical mastitis share some common QTL in two regions: 61.48 Mb - 73.84 Mb and 7.86 Mb – 39.55 Mb, respectively. As well, both ovulation rate and twinning rate might share a common QTL region from 34.16 Mb to 65.38 Mb. However, there are no common QTL locations in three pregnancy related phenotypes: non-return rate, pregnancy rate and daughter pregnancy rate. In beef cattle, the majority of QTL are located in a broad region of 15 Mb – 45 Mb on the chromosome. Functional genes, such as CRH, CYP11B1, DGAT1, FABP4 and TG, as potential candidates for some of these QTL, were also reviewed. Therefore, our review provides a standardized QTL map anchored within the current genome assembly, which would enhance the process of selecting positional and physiological candidate genes for many important traits in cattle.

Significant associations of stearoyl-CoA desaturase (SCD1) gene with fat deposition and composition in skeletal muscle

Gene expression studies in humans and animals have shown that elevated stearoyl-CoA desaturase (SCD1) activity is associated with increased fat accumulation and monounsaturation of saturated fatty acids in skeletal muscle. However, results of the two reported association studies in humans are inconsistent. In the present study, we annotated the bovine SCD1 gene and identified 3 single nucleotide polymorphisms (SNPs) in its 3'untranslated region (UTR). Genotyping these SNPs on a Wagyu x Limousin reference population revealed that the SCD1 gene was significantly associated with six fat deposition and fatty acid composition traits in skeletal muscle, but not with subcutaneous fat depth and percent kidney-pelvic-heart fat. In particular, we confirmed that the high stearoyl-CoA desaturase activities/alleles were positively correlated with beef marbling score, amount of monounsaturated fatty acids and conjugated linoleic acid content, but negatively with amount of saturated fatty acids. The inconsistent associations between human studies might be caused by using different sets of markers because we observed that most associated markers are located near the end of 3'UTR. We found that the proximity of the polyadenylation signal site is highly conserved among human, cattle and pig, indicating that the region might contain functional elements involved in posttranscriptional control of SCD1 activity. In conclusion, our cross species study provided solid evidence to support SCD1 gene as a critical player in skeletal muscle fat metabolism.

Comparative understanding of UTS2 and UTS2R genes for their involvement in type 2 diabetes mellitus

Several reports have shown that urotensin 2 (UTS2) and its receptor (UTS2R) are involved in glucose metabolism and insulin resistance, which lead to development of type 2 diabetes mellitus (T2DM) in humans. In the present study, we annotated both bovine UTS2 and UTS2R genes and identified 5 single nucleotide polymorphisms (SNPs) for the former gene and 14 mutations for the latter gene. Four mutations were genotyped on a Wagyu x Limousin reference population, including 6 F₁ bulls, 113 F₁ dams and ~250 F₂ progeny. Among 12 phenotypes related to fat deposition and fatty acid composition, we observed that the UTS2 gene was significantly associated with the amount of skeletal saturated fatty acids, while its receptor (UTS2R) gene had significant effects on amounts of saturated and monounsaturated fatty acids, Δ⁹ desaturase activity for converting 16:0 into 16:1, muscle fat (marbling) score and Longissimus Dorsi muscle area. However, in this population, these markers were not associated with subcutaneous fat depth or percent kidney, pelvic and heart fat. We also found that mutations in the promoter regions altered the promoter activities in both genes and coding SNPs might affect the mRNA stability in the UTS2R gene. Overall, our present study provides the first evidence that both UTS2 and UTS2R genes regulate skeletal muscle fat accumulation and fatty acid metabolism, thus indicating their potential pathological functions related to obesity and T2DM in humans.

Functional UQCRC1 polymorphisms affect promoter activity and body lipid accumulation

Obesity and type 2 diabetes constitute leading public health problems worldwide. Studies have shown that insulin resistance affiliated with these conditions is associated with skeletal muscle lipid accumulation, while the latter is associated with mitochondrial dysfunctions. However, the initiation and regulation of mitochondrial biogenesis rely heavily on approximately 1000 nuclear-encoded mitochondrial regulatory proteins. In this study, we targeted the ubiquinol-cytochrome c reductase core protein I gene, a nuclear-encoded component of mitochondrial complex III, for its association with subcutaneous fat depth (SFD) and skeletal muscle lipid accumulation (SMLA) using cattle as a model. Four promoter polymorphisms were identified and genotyped on approximately 250 Wagyu x Limousin F2 progeny. Statistical analysis revealed that two completely linked polymorphic sites, g.13487C>T and g.13709G>C (r2 = 1), were significantly associated with both SFD (p < 0.01) and SMLA (p < 0.0001). The difference between TTCC and CCGG haplotypes was 0.178 cm for SFD and 0.624 scores for SMLA. Interestingly, the former haplotype produced higher promoter activities than the latter by 43% to 49% in three cell lines (p < 0.05). In addition to Rett syndrome and breast/ovarian cancer observed in other studies, we report evidence for the first time, to our knowledge, that overexpression of ubiquinol-cytochrome c reductase core protein I might affect mitochondrial morphology and/or physiology and lead to development of obesity and related conditions.

Corticotropin releasing hormone is a promising candidate gene for marbling and subcutaneous fat depth in beef cattle

The gene corticotropin releasing hormone (CRH) is mapped on bovine chromosome 14 (BTA14), where more than 30 fat-related quantitative trait loci (QTLs) have been reported in dairy and beef cattle. The gene product regulates secretion of adrenocorticotrophin hormone, the hypothalamic–pituitary–adrenal axis, and multiple hypothalamic functions; therefore, we hypothesized that CRH is a promising candidate gene for beef marbling score (BMS) and subcutaneous fat depth (SFD) in a Wagyu × Limousin F2 population. Two pairs of primers were designed and a total of 5 single nucleotide polymorphisms (SNPs) were identified: g.9657C>T, c.10718G>C, c.10841G>A, c.10893A>C, and c.10936G>C (AAFC03076794.1). Among the 4 cSNPs, c.10718G>C, c.10841G>A, and c.10936G>C are missense mutations leading to amino acid changes from arginine to proline, from serine to asparagine, and from aspartic acid to histidine, respectively. These 5 SNPs were genotyped on ~250 F2 progeny, but only 4 were selected as tagging SNPs for association analysis because no historical recombination was observed between c.10718G>C and c.10893A>C. Statistical analysis showed that g.9657C>T, c.10718G>C, and c.10936G>C and their haplotypes had significant effects on SFD, but only c.10936G>C was significantly associated with BMS. The SNP in the promoter (g.9657C>T) led to gain/loss of a CpG site and 4 potential regulatory binding sites. Different haplotypes among the 4 cSNPs significantly affected mRNA secondary structures but were not associated with phenotypes. Overall, our results provide further evidence that CRH is a promising candidate gene for a concordant QTL related to lipid metabolism in mammals.

A novel type of sequence variation: multiple-nucleotide length polymorphisms discovered in the bovine genome

Three types of sequence variations--single-nucleotide polymorphisms (SNPs), insertions and deletions (indels), and short tandem repeats (STRs)--have been extensively reported in mammalian genomes. In this study, we discovered a novel type of sequence variation, i.e., multiple-nucleotide length polymorphisms (MNLPs) in bovine UCN3 (Urocortin 3) and its receptor CRHR2 (corticotropin-releasing hormone receptor 2) genes. Both MNLPs featured involvement of multiple-nucleotide length polymorphisms (5-18 bases), low sequence identity, and 1.7- to 11-fold changes in promoter activity between two alleles. Therefore, this novel genetic complexity would contribute significantly to the evolutionary, functional, and phenotypic complexity of genomes within or among species.

A simplified QTL mapping approach for screening and mapping of novel AFLP markers associated with beef marbling

Genome screening of quantitative trait loci (QTL) for a complex trait is usually costly and highly laborious, as it requires a large number of markers spanning the whole genome. Here we present a simplified approach for screening and mapping of QTL-linked markers for beef marbling using a WagyuxLimousin F(2) reference population. This simplified approach involves integration of the amplified fragment length polymorphism (AFLP) with DNA pooling and selective genotyping and comparative bioinformatics tools. AFLP analysis on two high and two low marbling DNA pools yielded ten visually different markers. Among them, four were confirmed based on individual AFLP validation. Sequencing and in silico characterization assigned two of these AFLP markers to bovine chromosomes 1 (BTA1) and 13 (BTA13), which are orthologous to human chromosomes HSA21q22.2 and HSA10p11.23 with both regions harboring QTL for obesity-related phenotypes. Both AFLP markers showed significantly large additive genetic effects (0.28+/-0.11 on BTA1 and 0.54+/-0.21 on BTA13) on beef-marbling score (BMS) (P<0.05). Overall, this approach is less time consuming, inexpensive and in particular, suitable for screening and mapping QTL-linked markers when targeting one or a few complex traits.

Cross species association examination of UCN3 and CRHR2 as potential pharmacological targets for antiobesity drugs

BACKGROUND

Obesity now constitutes a leading global public health problem. Studies have shown that insulin resistance affiliated with obesity is associated with intramyocellular lipid (IMCL) accumulation. Therefore, identification of genes associated with the phenotype would provide a clear target for pharmaceutical intervention and care for the condition. We hypothesized that urocortin 3 (UCN3) and corticotropin-releasing hormone receptor 2 (CRHR2) are associated with IMCL and subcutaneous fat depth (SFD), because the corticotropin-releasing hormone family of peptides are capable of strong anorectic and thermogenic effects.

METHODOLOGY/PRINCIPAL FINDINGS

We annotated both bovine UCN3 and CRHR2 genes and identified 12 genetic mutations in the former gene and 5 genetic markers in the promoter region of the latter gene. Genotyping of these 17 markers on Wagyu times Limousin F(2) progeny revealed significant associations between promoter polymorphisms and SFD (P = 0.0203-0.0685) and between missense mutations of exon 2 and IMCL (P = 0.0055-0.0369) in the bovine UCN3 gene. The SFD associated promoter SNPs caused a gain/loss of 12 potential transcription regulatory binding sites, while the IMCL associated coding SNPs affected the secondary structure of UCN3 mRNA. However, none of five polymorphisms in CRHR2 gene clearly co-segregated with either trait in the population (P>0.6000).

CONCLUSIONS/SIGNIFICANCE

Because UCN3 is located on human chromosome 10p15.1 where quantitative trait loci for obesity have been reported, our cross species study provides further evidence that it could be proposed as a potential target for developing antiobesity drugs. None of the markers in CRHR2 was associated with obesity-type traits in cattle, which is consistent with findings in human. Therefore, CRHR2 does not lend itself to the development of antiobesity drugs.

The bovine fatty acid binding protein 4 gene is significantly associated with marbling and subcutaneous fat depth in Wagyu x Limousin F2 crosses

Fatty acid binding protein 4 (FABP4), which is expressed in adipose tissue, interacts with peroxisome proliferator-activated receptors and binds to hormone-sensitive lipase and therefore, plays an important role in lipid metabolism and homeostasis in adipocytes. The objective of this study was to investigate associations of the bovine FABP4 gene with fat deposition. Both cDNA and genomic DNA sequences of the bovine gene were retrieved from the public databases and aligned to determine its genomic organization. Primers targeting two regions of the FABP4 gene were designed: from nucleotides 5433-6106 and from nucleotides 7417-7868 (AAFC01136716). Direct sequencing of polymerase chain reaction (PCR) products on two DNA pools from high- and low-marbling animals revealed two single nucleotide polymorphisms (SNPs): AAFC01136716.1:g.7516G>C and g.7713G>C. The former SNP, detected by PCR-restriction fragment length polymorphism using restriction enzyme MspA1I, was genotyped on 246 F2 animals in a Waygu x Limousin F2 reference population. Statistical analysis showed that the FABP4 genotype significantly affected marbling score (P = 0.0398) and subcutaneous fat depth (P = 0.0246). The FABP4 gene falls into a suggestive/significant quantitative trait loci interval for beef marbling that was previously reported on bovine chromosome 14 in three other populations.

A novel nuclear-encoded mitochondrial poly(A) polymerase PAPD1 is a potential candidate gene for the extreme obesity related phenotypes in mammals

People with obesity, especially extreme obesity, are at risk for many health problems. However, the responsible genes remain unknown in >95% of severe obesity cases. Our previous genome-wide scan of Wagyu x Limousin F₂ cattle crosses with extreme phenotypes revealed a molecular marker significantly associated with intramuscular fat deposition. Characterization of this marker showed that it is orthologous to the human gene KIAA1462 located on HSA10p11.23, where a major quantitative trait locus for morbid obesity has been reported. The newly identified mitochondrial poly(A) polymerase associated domain containing 1 (PAPD1) gene, which is located near this marker, is particularly interesting because the polymerase is required for the polyadenylation and stabilization of mammalian mitochondrial mRNAs. In the present study, both cDNA and genomic DNA sequences were annotated for the bovine PAPD1 gene and ten genetic markers were detected in the promoter and exon 1 region. Among seven markers assayed on ~ 250 Wagyu x Limousin F₂ animals, two single nucleotide polymorphisms (SNPs) in the promoter region were significantly associated with intramuscular fat (P<0.05). However, there was a significant interaction (P<0.05) between a third SNP, which causes an amino acid change in coding exon 1, and each of these two promoter SNPs on intramuscular fat deposition. In particular, the differences between double heterozygous animals at two polymorphic sites and the slim genotype animals exceeded 2.3 standard deviations for the trait in both cases. Our study provides evidence for a new mechanism – the involvement of compound heterosis in extreme obesity, which warrants further examination.