Genome-wide CNV analysis reveals variants associated with growth traits in Bos indicus
- Yang Zhou†1, 2,
- Yuri T. Utsunomiya†3,
- Lingyang Xu†1, 4,
- El Hamidi abdel Hay1,
- Derek M. Bickhart1,
- Pamela Almeida Alexandre5,
- Benjamin D. Rosen1,
- Steven G. Schroeder1,
- Roberto Carvalheiro6,
- Haroldo Henrique de Rezende Neves6,
- Tad S. Sonstegard1, 8,
- Curtis P. Van Tassell1,
- José Bento Sterman Ferraz5,
- Heidge Fukumasu5,
- Jose Fernando Garcia3, 7, 9Email author and
- George E. Liu1Email author
© Zhou et al. 2016
Received: 17 August 2015
Accepted: 11 February 2016
Published: 1 June 2016
Apart from single nucleotide polymorphism (SNP), copy number variation (CNV) is another important type of genetic variation, which may affect growth traits and play key roles for the production of beef cattle. To date, no genome-wide association study (GWAS) for CNV and body traits in beef cattle has been reported, so the present study aimed to investigate this type of association in one of the most important cattle subspecies: Bos indicus (Nellore breed).
We have used intensity data from over 700,000 SNP probes across the bovine genome to detect common CNVs in a sample of 2230 Nellore cattle, and performed GWAS between the detected CNVs and nine growth traits. After filtering for frequency and length, a total of 231 CNVs ranging from 894 bp to 4,855,088 bp were kept and tested as predictors for each growth trait using linear regression analysis with principal components correction. There were 49 significant associations identified among 17 CNVs and seven body traits after false discovery rate correction (P < 0.05). Among the 17 CNVs, three were significant or marginally significant for all the traits. We have compared the locations of associated CNVs with quantitative trait locus and the RefGene database, and found two sets of 9 CNVs overlapping with either known QTLs or genes, respectively. The gene overlapping with CNV100, KCNJ12, is a functional candidate for muscle development and plays critical roles in muscling traits.
This study presents the first CNV-based GWAS of growth traits using high density SNP microarray data in cattle. We detected 17 CNVs significantly associated with seven growth traits and one of them (CNV100) may be involved in growth traits through KCNJ12.
KeywordsCopy number variation (CNV) EBV Association Body size Bos indicus
Growth traits can directly affect production of beef cattle, which is an economically important sector of global agriculture. Both body weight measurements and visual trait scores of conformation, precocity and muscling are used to select calves and steers for breeding and to predict carcass parameters. Various factors can regulate growth traits through the modulation of gene expression levels, and genomic variation is one of the mechanisms, which is not only persistent but also inherited. To date, several genomic variations, such as STR (short tandem repeat, or microsatellite) and SNP (single nucleotide polymorphism) have been employed to identify genetic markers related to cattle growth traits [1, 2]. Copy number variation (CNV) belongs to another family of genetic markers, which could affect gene expression and consequently phenotypes by changes in gene structure and dosage . In human and mouse studies, CNV was found to account for 18 to 30 % of the genetic variation in gene expression and showed critical roles in both normal phenotypic variability and disease susceptibility [4, 5]. However, only few studies have been published using CNVs as candidate markers for genetic selection and disease association at the genome level because the traditional CNV discovery studies tried to detect as many variable regions as possible, instead of emphasizing common CNVs shared by different individuals .
The application of CNV information in animal breeding had long been held back due to the lack of reliable genotyping methods for CNVs. In livestock, most CNV studies had been limited to CNV detection methods [2, 7–10]. Studies performing CNV association with body traits were mainly focused on single rather than genome-wide CNVs [11–13]. Several studies attempted to predict the effects of CNVs on phenotypes based on the function of genes within the CNV regions. For example, an earlier study in Holstein cattle reported CNV loci enriched for solute carriers and ABC transporters that may affect milk production traits . Exploration of genome-wide association study (GWAS) using CNVs and phenotypes has just gotten started recently. Our previous studies have indicated that CNVs could be associated with resistance or susceptibility to gastrointestinal nematodes in Angus cattle  and residual feed intake in Holstein cows . A study on Chinese cattle has shown that one CNV region (CNVR) had significantly negative effects on expression of PLA2G2D gene and two CNVRs were associated with body measurements . Subsequently, we have reported and characterized 34 significantly associated CNVs with milk production traits in Holsteins  and found one deletion polymorphism associated with resistance to gastrointestinal nematodes in Angus cattle . The increase in resolution of the SNP microarray data dramatically improves the accuracy and the number of CNV detected. Up to now, no systematic study of the relationship between CNVs and beef growth traits using SNP microarray data with high resolution has been published.
Few CNV studies including Bos indicus data have been reported. We once performed a high resolution CNV analysis in 674 animals of 27 breeds, which included 113 animals with Bos indicus background . Nellore is one the most important sub-breeds of Bos indicus cattle, which is predominant in Brazil, the world’s second largest producer of cattle. They are farmed for beef that can be seriously affected by growth traits. The present study is the first using common CNVs to detect association with production traits in Bos indicus. Using BovineHD Genotyping BeadChip assay data of 2230 Nellore bulls and cows, we detected common CNVs and performed a systematic GWAS analysis between CNVs and nine growth traits in order to identify candidate CNVs and genes that could be used in cattle genetic improvement programs.
Results and Discussion
CNV segments and genotyping
Traits’ correlation analysis
CNV significant association analysis
Characterization of the association CNVs
P value after FDR correction
Significant CNV characterization
We compared the 17 significantly associated CNVs discovered with those previously described in Hou’s study with the same cattle genome reference . Only 6 CNVs were not supported by Hou’s study. The frequencies of the 17 CNVs ranged from 8.25 % to 95.65 %, with the lowest frequency CNV204 located at chr3:54,541,000-54,963,401 and the highest frequency CNV17 located at chr1:150,048,097-150,057,730. The longest CNV was also CNV204 with 422,402 bp and the shortest CNV175 (chr24: 24,432,820-24,434,451) was 1632 bp. To further analyze the association between discovered CNVs and body traits, we overlapped the 17 significant CNVs with previously reported QTL at the Cattle QTL database (http://www.animalgenome.org/cgi-bin/QTLdb/BT/index). Considering all together, 8 CNVs were found to be overlapped with QTLs (Table 1). There were three CNVs (CNV93, CNV100 and CNV167) overlapped with QTLs that regulate intramuscular fat which is highly related to the growth traits. Iron content is an important factor affecting the quality of beef, including color and flavor. We found CNV107 significantly associated with MW, overlapping with QTLs related to iron content, suggesting that there may be some relations between iron content and MW. It is noted that although we used comparatively strict criteria to examine the overlap of CNV with QTL, some of them may be due to chance.
CNVs can have dramatic phenotypic consequences caused by altering gene dosage, disrupting gene coding sequence or perturbing long-range gene regulation. We compared the locations of the 17 significantly associated CNVs with the bovine refGene database (http://hgdownload.soe.ucsc.edu/goldenPath/bosTau6/database/). There were 9 CNVs overlapping with 13 genes (Fig. 1 and Table 1). Among these 13 genes, five (KCNJ12, SETD4, CBR1, DYNC1LI1 and PRAME) were overlapped by at least one exon on their 5’ terminal end. CNV98 and CNV121 were completely located at the last introns of EPHB3 and COL26A, respectively. Conversely, five genes were completely located within two CNV regions (GBP6 and GBP4 in CNV204, HSFY2 and ZNF280B in CNV96, TSPY in CNV99). CNV99 was also overlapped with the 9th intron of PANCC. Nucleotide sequence variation in the promoter region can also affect the corresponding gene’s expression. We also examined whether any CNVs overlapped with regions up to 2000 bp upstream of gene transcript start sites but none was found.
KCNJ12 as a candidate gene for muscling
Several of the genes found in this study might contribute to body traits, when associated with overlapped CNVs, according to their Gene Ontology and pathway annotation in different species. We selected KCNJ12 (potassium inwardly-rectifying channel, subfamily J, member 12), located at chr19: 35,955,280-35,993,796 as a candidate to do further analysis. Its first exon overlaps with CNV100 (chr19: 35933815-35957534), which was significantly associated with the two highly related traits (muscling and precocity), as well as being associated with the WG and the visual conformation traits. Both muscling and precocity visual traits are designed to predict carcass features in the live animal, as an approximation of the muscle (meat) contained in the carcass itself after slaughtering. Previous studies have shown that KCNJ12 gene family members were expressed in brain and muscle . RNA-seq analysis reported by previous studies revealed that bovine KCNJ12 gene expression was significantly up-regulated from fetal to adult stage in longissimus muscle , pointing out its role in muscle development.
However, to date, there is little knowledge about this gene in cattle. We first compared the amino acids sequence of KCNJ12 in different species including cattle, human, mouse, rat, monkey, pig and dog, showing that KCNJ12 was fairly conserved except for some variable amino acids in the 3’ termination among all the mammals (Additional file 2: Figure S1). This means that the function of KCNJ12 gene may also be conserved in the above species. According to the Gene Ontology Annotation (UniProt-GOA) Database, the cattle KCNJ12 gene has 9 GO terms and it modulates the activity of inward rectifier potassium channel, which regulates the electrical activity of neuronal cells, insulin secretion and epithelial K+ transport . Since the muscle fiber number does not increase after birth , the increase of muscle mass is mainly due to increase in length and cross-sectional area of muscle fibers stimulated by neural, hormonal and nutritional factors after birth. In human, KCNJ12 was found to participate in muscling contraction process . This may be one of the mechanisms by which KCNJ12 regulates the amount of muscle mass through the modulation of muscle contraction in cattle. KCNJ12 was annotated in the GABAB receptor activation pathway, in which GABAB receptors (GABABR) are metabotropic transmembrane receptors for gamma-aminobutyric acid (GABA) which are linked via G-proteins to potassium channels. Previous studies reported that they were required for normal development and that, together with ionotropic GABAA receptors, GABABR could help GABA to stimulate feeding which is another mechanism to promote muscle or even fat accumulation [27, 28]. It’s worth mentioning that the CNV100 was also overlapped with QTLs related to intramuscular fat and calf size. Moreover, we checked the KCNJ12 expression pattern in Bovine Gene Atlas data (http://bovineatlas.arl.arizona.edu/) , hosted by the innatedb (http://www.innatedb.com/) . We found that it was expressed in 45 out of 87 tissues. It was highly expressed in multiple muscle types or muscle related tissues (Additional file 2: Figure S3). We examined the effect of the CNV100’s copy number on the gene expression of KCNJ12 in muscle tissues using qPCR. The Nellore cattle with 3 copies of CNV100 had a significantly higher RNA level of the KCNJ12 gene than the Nellore cattle with 2 copies had (p < 0.05, ANOVA examination, Additional file 2: Figure S4). Therefore, KCNJ12 might also help muscle accumulation by increasing food intake through cooperation with GABAB receptors. In summary, the bovine KCNJ12 gene could be selected as a candidate gene for muscling traits in beef cattle in future studies.
This is the first study applying GWAS approaches between the common CNVs and Bos indicus (Nellore breed) growth traits using bovine high density SNP microarray data. The use of high-density SNP panels enabled the increase in four times the number of CNVs found when compared to Holstein cattle using low density SNP microarray data . We identified 231 CNVs with high confidence and 17 of them were significantly associated with seven moderately or highly correlated growth traits. Three CNVs (CNV93, CNV97 and CNV100) were found significantly associated with all the seven growth traits. Our association study was also supported by the results of overlapping the CNV locations with known QTLs and RefGene databases. We selected KCNJ12, overlapped with CNV100, to functionally assess and concluded that it may be a candidate gene for muscling through the modulation of muscle contraction and food intake.
A total of 952 Nellore bulls and 1278 Nellore cows were genotyped for 777,962 SNPs with the Illumina® BovineHD Genotyping BeadChip assay (Illumina Inc.; 2011 BovineHD Genotyping BeadChip Data Sheet: DNA Analysis. http://www.illumina.com/Documents/products/datasheets/datasheet_bovineHD.pdf), according to the manufacturer’s protocol. This dataset builds on the data reported by previous studies , and comprises part of the genomic selection reference population from a commercial breeding program (DeltaGen) ran by an alliance of Nellore cattle breeders from Brazil. Genotyping data is available upon request from co-corresponding authors for research purposes. All animal protocols were approved by the Institutional Animal Care and Use Committee of Faculty of Food Engineering and Animal Sciences, University of São Paulo (FZEA-USP – protocol number 14.1.636.74.1). Additionally, genomic DNA was extracted either from commercialized semen straws (bulls) or stored hair (cows) samples. To examine the effect of CNV100’s copy number on the gene expression of KCNJ12, DNA and RNA samples were collected from 24 Nellore cattle muscle tissues.
Phenotypic and EBV values
Estimated breeding values (EBVs) were based on Best Linear Unbiased Predictor (BLUP) estimates of single-trait animal models obtained from routine genetic evaluations using performance and pedigree data from the database (available at: http://www.gensys.com.br/home/show_page.php?id=701). Phenotypes used to fit the models comprised records from 542,918 animals born between 1985 and 2011, and raised in 243 grazing-based herds. The evaluated traits included birth weight (BW), post weaning gain (PWG), weaning gain (WG), carcass conformation at weaning (CW), muscling at weaning (MW), carcass finishing precocity at weaning (PW), carcass conformation at yearling (CY), muscling at yearling (MY) and carcass finishing precocity at yearling (PY). Conformation, finishing precocity and muscling traits (CPM) were based on recorded visual scores assigned in a discrete ordered scale, relative to the animals of the same management group. For each trait, only EBVs of animals whose accuracy (i.e., square root of reliability, calculated based on prediction error variance estimates) was > 0.50 were analyzed (Additional file 2: Figure S2 and Additional file 1: Table S3).
CNV segmentation and genotyping
The high intensity data generated by Illumina BovineHD Genotyping BeadChip assay was applied to detect the common CNVs shared by different Nellore cattle. After exporting the DSF file from GenomeStudio Software, we imported its Log R Ratios (LRR) into Golden Helix SNP & Variation Suite (SVS) 8.3.0 (Golden Helix Inc., Bozeman, MT, USA) and successfully mapped 753,182 SNPs (96.81 %) onto the 29 autosomes of Bos taurus genome assembly UMD 3.1. The LRR was then normalized using the default GC correlation file to correct the waviness caused by the GC content. To define the CNV segments, the multivariate method was employed with following thresholds: (1) Max 20 segments per 10 k markers; (2) at least 3 markers per segment; (3) max pairwise segment p-value = 0.005. The three state covariates with a comparatively strict threshold (segment mean ±0.4) was used to define the CNVs as three types (loss, neutral and gain events) across all the samples.
PCA-corrected association testing
Where y i is the EBV of the ith individual, x ij is the CNV genotype of the ith individual, β j is the CNV effect, n is the number of CNVs and e i is the residual. To correct the batch effects/stratification, we selected the Principle Component Analysis (PCA) option. The significance level of (P < 0.05) after FDR correction was chose to define the significant CNVs .
Overlapping with QTL and gene analysis
Gene content of cattle CNV regions was assessed using RefGene annotation file in UCSC database (http://hgdownload.soe.ucsc.edu/goldenPath/bosTau6/database/). QTL database was downloaded from animal QTL database (http://www.animalgenome.org/cgi-bin/QTLdb/index, on UMD3.1 database). Overlaps with genes were detected using R script and defined as at least one bp overlap. We filtered the QTLs with confidence interval (CI) over 30 Mb because of over large CIs for some QTLs and used strict threshold to define the overlap as at least 50 % of the CNV length were covered by QTLs .
CNV validation by qPCR
Primers were designed for qPCR validation using the NCBI Primer-BLAST webtool (http://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi?LINK_LOC=BlastHome). Primer information can be seen in (Additional file 1: Table S4). SYBR green chemistry was used to carry out qPCR in triplicate reactions of 25 μl. All reactions were amplified on a bioRad MyIQ thermocycler. A common Nellore sample was chosen as the reference for all qPCR experiments and the BTF3 and GAPDH were used as the internal DNA and RNA control, respectively. We examined the effect of CNV100’s copy number on the expression of KCNJ12. The CNV100’s copy numbers were determined using qPCR and two samples were found to be 3 copies. We selected twelve samples with 2 copies and the two samples with 3 copies to compare the RNA levels of the KCNJ12 gene in muscle. Each sample was tested three times. The 2-ΔΔCT method was employed to analyze the qPCR result .
We thank Reuben Anderson, and Alexandre Dimtchev for technical assistance. This work was supported in part by Agriculture and Food Research Initiative (AFRI) competitive grant No. 2011-67015-30183 from the USDA National Institute of Food and Agriculture (NIFA) Animal Genome Program. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the US Department of Agriculture.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Curi R, Oliveira H, Silveira A, Lopes C. Effects of polymorphic microsatellites in the regulatory region of IGF1 and GHR on growth and carcass traits in beef cattle. Anim Genet. 2005;36(1):58–62.View ArticlePubMedGoogle Scholar
- Utsunomiya YT, Do Carmo AS, Carvalheiro R, Neves HH, Matos MC, Zavarez LB, O’Brien AMP, Sölkner J, McEwan JC, Cole JB. Genome-wide association study for birth weight in Nellore cattle points to previously described orthologous genes affecting human and bovine height. BMC Genet. 2013;14(1):52.
- Zhang F, Gu W, Hurles ME, Lupski JR. Copy number variation in human health, disease, and evolution. Annu Rev Genom Hum G. 2009;10:451–81.
- Henrichsen CN, Vinckenbosch N, Zöllner S, Chaignat E, Pradervand S, Schütz F, Ruedi M, Kaessmann H, Reymond A. Segmental copy number variation shapes tissue transcriptomes. Nat Genet. 2009;41(4):424–9.
- Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, Redon R, Bird CP, De Grassi A, Lee C. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007;315(5813):848–53.
- Xu L, Cole JB, Bickhart DM, Hou Y, Song J, VanRaden PM, Sonstegard TS, Van Tassell CP, Liu GE. Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins. BMC Genomics. 2014;15(1):683.
- Bae JS, Cheong HS, Kim LH, NamGung S, Park TJ, Chun J-Y, Kim JY, Pasaje CF, Lee JS, Shin HD. Identification of copy number variations and common deletion polymorphisms in cattle. BMC Genomics. 2010;11(1):232.
- Bickhart DM, Hou Y, Schroeder SG, Alkan C, Cardone MF, Matukumalli LK, Song J, Schnabel RD, Ventura M, Taylor JF. Copy number variation of individual cattle genomes using next-generation sequencing. Genome Res. 2012;22(4):778–90.
- Choi J-W, Lee K-T, Liao X, Stothard P, An H-S, Ahn S, et al. Genome-wide copy number variation in Hanwoo, Black Angus, and Holstein cattle. Mamm Genome. 2013;24(3-4):151–63.
- Cicconardi F, Chillemi G, Tramontano A, Marchitelli C, Valentini A, Ajmone-Marsan P, Nardone A. Massive screening of copy number population-scale variation in Bos taurus genome. BMC Genomics. 2013;14(1):124.
- Xu Y, Shi T, Cai H, Zhou Y, Lan X, Zhang C, Lei C, Qi X, Chen H. Associations of myh3 gene copy number variations with transcriptional expression and growth traits in chinese cattle. Gene. 2014;535(2):106–11.
- Glick G, Shirak A, Seroussi E, Zeron Y, Ezra E, Weller JI, Ron M. Fine mapping of a QTL for Fertility on BTA7 and its association with a CNV in the Israeli Holsteins. G3: Genes Genom Genet. 2011;1(1):65–74.
- Kadri NK, Sahana G, Charlier C, Iso-Touru T, Guldbrandtsen B, Karim L, Nielsen US, Panitz F, Aamand GP, Schulman N. A 660-kb deletion with antagonistic effects on fertility and milk production segregates at high frequency in nordic red cattle: additional evidence for the common occurrence of balancing selection in livestock. PLoS Genet. 2014;10(1):e1004049.
- Seroussi E, Glick G, Shirak A, Yakobson E, Weller JI, Ezra E, Zeron Y. Analysis of copy loss and gain variations in Holstein cattle autosomes using BeadChip SNPs. BMC Genomics. 2010;11(1):673.
- Hou Y, Liu GE, Bickhart DM, Matukumalli LK, Li C, Song J, Gasbarre LC, Van Tassell CP, Sonstegard TS. Genomic regions showing copy number variations associate with resistance or susceptibility to gastrointestinal nematodes in Angus cattle. Funct Integr Genomic. 2012;12(1):81–92.
- Hou Y, Bickhart DM, Chung H, Hutchison JL, Norman HD, Connor EE, Liu GE. Analysis of copy number variations in Holstein cows identify potential mechanisms contributing to differences in residual feed intake. Funct Integr Genomic. 2012;12(4):717–23.
- Zhang L, Jia S, Yang M, Xu Y, Li C, Sun J, Huang Y, Lan X, Lei C, Zhou Y. Detection of copy number variations and their effects in Chinese bulls. BMC Genomics. 2014;15(1):480.
- Xu L, Hou Y, Bickhart DM, Song J, Van Tassell CP, Sonstegard TS, Liu GE. A genome-wide survey reveals a deletion polymorphism associated with resistance to gastrointestinal nematodes in Angus cattle. Funct Integr Genomic. 2014;14(2):333–9.
- Hou Y, Bickhart DM, Hvinden ML, Li C, Song J, Boichard DA, Fritz S, Eggen A, DeNise S, Wiggans GR. Fine mapping of copy number variations on two cattle genome assemblies using high density SNP array. BMC Genomics. 2012;13(1):376.
- Liu GE, Hou Y, Zhu B, Cardone MF, Jiang L, Cellamare A, Mitra A, Alexander LJ, Coutinho LL, Dell'Aquila ME. Analysis of copy number variations among diverse cattle breeds. Genome Res. 2010;20(5):693–703.
- Mills RE, Walter K, Stewart C, Handsaker RE, Chen K, Alkan C, Abyzov A, Yoon SC, Ye K, Cheetham RK. Mapping copy number variation by population-scale genome sequencing. Nature. 2011;470(7332):59–65.
- De Boer T, Houtman M, Compier M, Van der Heyden M. The mammalian KIR2. x inward rectifier ion channel family: expression pattern and pathophysiology. Acta Physiol. 2010;199(3):243–56.Google Scholar
- He H, Liu X. Characterization of transcriptional complexity during longissimus muscle development in bovines using high-throughput sequencing. PLoS One. 2013;8(6):e64356.View ArticlePubMedPubMed CentralGoogle Scholar
- Reimann F, Ashcroft FM. Inwardly rectifying potassium channels. Curr Opin Cell Biol. 1999;11(4):503–8.View ArticlePubMedGoogle Scholar
- Allen RE, Merkel RA, Young RB. Cellular aspect of muscle growth: myogenic cell proliferation. J Anim Sci. 1979;49(1):115–27.PubMedGoogle Scholar
- Wible BA, De Biasi M, Majumder K, Taglialatela M, Brown AM. Cloning and functional expression of an inwardly rectifying K+ channel from human atrium. Circ Res. 1995;76(3):343–50.View ArticlePubMedGoogle Scholar
- Dzitoyeva S, Gutnov A, Imbesi M, Dimitrijevic N, Manev H. Developmental role of GABA B (1) receptors in Drosophila. Dev Brain Res. 2005;158(1):111–4.View ArticleGoogle Scholar
- Bäckberg M, Collin M, Ovesjö ML, Meister B. Chemical coding of GABAB receptor‐immunoreactive neurones in hypothalamic regions regulating body weight. J Neuroendocrinol. 2003;15(1):1–14.View ArticlePubMedGoogle Scholar
- Harhay GP, Smith T, Alexander LJ, Haudenschild CD, Keele JW, Matukumalli LK, Schroeder SG, Van Tassell CP, Gresham CR, Bridges SM. An atlas of bovine gene expression reveals novel distinctive tissue characteristics and evidence for improving genome annotation. Genome Biol. 2010;11(10):R102.
- Breuer K, Foroushani AK, Laird MR, Chen C, Sribnaia A, Lo R, Winsor GL, Hancock RE, Brinkman FS, Lynn DJ. InnateDB: systems biology of innate immunity and beyond—recent updates and continuing curation. Nucleic Acids Res. 2013;41(Database issue):D1228-33.