Genomic clues of the evolutionary history of Bos indicus cattle.

Together with their sister subspecies Bos taurus, zebu cattle (Bos indicus) have contributed to important socioeconomic changes that have shaped modern civilizations. Zebu cattle were domesticated in the Indus Valley 8000 years before present (YBP). From the domestication site, they expanded to Africa, East Asia, southwestern Asia and Europe between 4000 and 1300 YBP, intercrossing with B. taurus to form clinal variations of zebu ancestry across the landmass of Afro-Eurasia. In the past 150 years, zebu cattle reached the Americas and Oceania, where they have contributed to the prosperity of emerging economies. The zebu genome is characterized by two mitochondrial haplogroups (I1 and I2), one Y chromosome haplogroup (Y3) and three major autosomal ancestral groups (Indian-Pakistani, African and Chinese). Phenotypically, zebu animals are recognized by their hump, large ears and excess skin. They are rustic, resilient to parasites and capable of bearing the hot and humid climates of the tropics. Many resources are available to study the zebu genome, including commercial arrays of SNP, reference assemblies and publicly available genotypes and whole-genome sequences. Nevertheless, many of these resources were initially developed to support research and subsidize industrial applications in B. taurus, and therefore they can produce bias in data analysis. The combination of genomics with precision agriculture holds great promise for the identification of genetic variants affecting economically important traits such as tick resistance and heat tolerance, which were naturally selected for millennia and played a major role in the evolution of B. indicus cattle.

Pre- and post-puberty expression of genes and proteins in the uterus of Bos indicus heifers: the luteal phase effect post-puberty.

Progesterone signaling and uterine function are crucial in terms of pregnancy establishment. To investigate how the uterine tissue and its secretion changes in relation to puberty, we sampled tissue and uterine fluid from six pre- and six post-pubertal Brahman heifers. Post-pubertal heifers were sampled in the luteal phase. Gene expression of the uterine tissue was investigated with RNA-sequencing, whereas the uterine fluid was used for protein profiling with mass spectrometry. A total of 4034 genes were differentially expressed (DE) at a nominal P-value of 0.05, and 26 genes were significantly DE after Bonferroni correction (P < 3.1 × 10-6 ). We also identified 79 proteins (out of 230 proteins) that were DE (P < 1 × 10-5 ) in the uterine fluid. When we compared proteomics and transcriptome results, four DE proteins were identified as being encoded by DE genes: OVGP1, GRP, CAP1 and HBA. Except for CAP1, the other three had lower expression post-puberty. The function of these four genes hypothetically related to preparation of the uterus for a potential pregnancy is discussed in the context of puberty. All DE genes and proteins were also used in pathway and ontology enrichment analyses to investigate overall function. The DE genes were enriched for terms related to ribosomal activity. Transcription factors that were deemed key regulators of DE genes are also reported. Transcription factors ZNF567, ZNF775, RELA, PIAS2, LHX4, SOX2, MEF2C, ZNF354C, HMG20A, TCF7L2, ZNF420, HIC1, GTF3A and two novel genes had the highest regulatory impact factor scores. These data can help to understand how puberty influences uterine function.

Genome-wide detection of signatures of selection in Korean Hanwoo cattle.

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26-0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise FST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome-wide smoothed FST , 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty-four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.

Performance of different SNP panels for parentage testing in two East Asian cattle breeds.

The International Society for Animal Genetics (ISAG) proposed a panel of single nucleotide polymorphisms (SNPs) for parentage testing in cattle (a core panel of 100 SNPs and an additional list of 100 SNPs). However, markers specific to East Asian taurine cattle breeds were not included, and no information is available as to whether the ISAG panel performs adequately for these breeds. We tested ISAG's core (100 SNP) and full (200 SNP) panels on two East Asian taurine breeds: the Korean Hanwoo and the Japanese Wagyu, the latter from the Australian herd. Even though the power of exclusion was high at 0.99 for both ISAG panels, the core panel performed poorly with 3.01% false-positive assignments in the Hanwoo population and 3.57% in the Wagyu. The full ISAG panel identified all sire-offspring relations correctly in both populations with 0.02% of relations wrongly excluded in the Hanwoo population. Based on these results, we created and tested two population-specific marker panels: one for the Wagyu population, which showed no false-positive assignments with either 100 or 200 SNPs, and a second panel for the Hanwoo, which still had some false-positive assignments with 100 SNPs but no false positives using 200 SNPs. In conclusion, for parentage assignment in East Asian cattle breeds, only the full ISAG panel is adequate for parentage testing. If fewer markers should be used, it is advisable to use population-specific markers rather than the ISAG panel.

Variation in the XKR4 gene was significantly associated with subcutaneous rump fat thickness in indicine and composite cattle.

Variation in the XK, Kell blood group complex subunit-related family, member 4 (XKR4) gene on BTA14 was associated with rump fat thickness in a recent genome-wide association study. This region is also of interest because it is known to show evidence of a signature of population genetic selection. In this study, additional variation in this gene was genotyped in a sample of a total of 1283 animals of the Belmont Red (BEL) and Santa Gertrudis (SGT) breeds. The SNP rs41724387 was significantly (P < 0.001) associated with rump fat thickness and explained 5.9% of the genetic variance for the trait in this sample. Using the 4466 genotypes for the SNP rs42646708 from several data sets to estimate effects in seven breeds, this relatively large quantitative trait locus effect appears to be a result of the variation in indicine and taurine-indicine composite cattle. However, the only DNA variant found in Brahman cattle that altered the predicted amino acid sequence of XKR4 was not associated with rump fat thickness. This suggests that causative mutations lie outside the coding sequence of this gene.

DNA variation in the gene ELTD1 is associated with tick burden in cattle.

Ticks and tick-born diseases are major constraints on cattle production in tropical and subtropical regions in the world. Previously, we identified single nucleotide polymorphisms (SNPs) associated with tick resistance on bovine chromosome 3 at approximately 70 Mb. In this study, we genotyped a dairy (n = 1133) and a beef (n = 774) sample to confirm the association of the intronic SNP rs29019303 and its gene (ELTD1) with tick burden. We genotyped 18 additional SNPs in a region of 181 kb and found that rs29019303 was significantly (P < 0.05) associated with tick burden in both samples with the same favourable allele. A second SNP in this same genomic region was also significantly associated with tick burden in each sample. The associations using haplotypes were stronger than for single markers, including a haplotype of nine tag SNPs that was highly significantly (P = 0.0008) associated with tick counts in the dairy animals. This haplotype and two others were significant after Bonferroni correction for multiple testing. The estimated size of the effects was close to 0.9% of the residual variance in both samples tested.

Myostatin (GDF8) single nucleotide polymorphisms in Nellore cattle.

The myostatin gene, also known as GDF8 (growth differentiation factor 8), is located on bovine chromosome 2 (BTA2); it has three exons and two introns. Myostatin is specifically expressed during embryonic development and in adult skeletal muscle, functioning as a negative regulatory protein. Several cattle breeds (Piedmontese, Belgian Blue and Blond'Aquitaine, and others) show polymorphisms in this gene; these polymorphisms are directly related to the double muscling phenotype. We looked for polymorphisms in the Nellore cattle myostatin gene and compared them with those known for taurine breeds. Seven regions, covering the three exons of this gene, were amplified by polymerase chain reaction and sequenced, including the untranslated region. DNA from 30 adult Nellore animals was collected; DNA sequencing revealed three, seven and four polymorphisms in exons 1, 2 and 3, respectively. We found previously reported polymorphisms, as well as several new ones; for instance, 37 polymorphisms were found in the untranslated region segment, and in introns 1 and 2 there were one and three polymorphisms, respectively. The high degree of allelic heterogeneity in the myostatin gene could be related to its high mutation rate; it also could be the result of a long history of artificial selection for meat production, which has probably favored such modifications and maintained them in cattle populations. These polymorphisms identified in Nellore cattle could be useful for breeding programs.

Genomic inbreeding depression for climatic adaptation of tropical beef cattle.

Inbreeding has the potential to negatively impact animal performance. Strategies to monitor and mitigate inbreeding depression require that it can be accurately estimated. Here, we used genomewide SNP data to explore 3 alternative measures of genomic inbreeding: the diagonal elements of the genomic relationship matrix (FGRM), the proportion of homozygous SNP (FHOM), and the proportion of the genome covered by runs of homozygosity (FROH). We used 2,111 Brahman (BR) and 2,550 Tropical Composite (TC) cattle with phenotypes recorded for 10 traits of relevance to tropical adaptation. We further explored 3 marker densities ranging from a high-density chip (729,068 SNP), a medium-density chip (71,726 SNP) specifically designed for cattle, and a low-density chip (18,860 SNP) associated with the measures of inbreeding. Measures of FGRM were highly correlated across the 3 SNP densities and negatively correlated with FHOM and FROH in the BR population. In both populations, there was a strong positive correlation for each measure of inbreeding across the 3 SNP panels. We found significant ( < 0.01) inbreeding depression for various traits, particularly when using the highest-density SNP chip in the BR population, where inbreeding was negatively associated with coat color and coat type such that inbred animals presented shorter, slicker, and lighter coats. Based on FGRM using the medium-density chip, we found that a 1% increase in inbreeding in the BR and TC populations was associated with a decrease of 0.514 and 0.579 kg BW, respectively, in yearlings. In the TC population, a 1% increase in FHOM was associated with a decrease in BCS of -0.636% ( < 0.001). The low-density chip, comprising SNP associated with inbreeding, captured genes, and regions with pleiotropic effects ( < 0.001). However, it did not improve our ability to identify inbreeding depression, relative to the use of higher-density panels. We conclude that where heterogeneous populations are present, such as in tropical environments where composite animals abound, measures of inbreeding that do not depend on allele frequencies, such as FHOM and FROH, are preferable for estimating genomic inbreeding. Finally, the sustainable intensification of livestock systems in tropical regions will rely on genetic safeguards to ensure that productivity is improved while also adapting animals to cope with climate change. The results of this study are a step toward achieving that goal.

Global differential gene expression in the pituitary gland and the ovaries of pre- and postpubertal Brahman heifers.

To understand genes, pathways, and networks related to puberty, we characterized the transcriptome of two tissues: the pituitary gland and ovaries. Samples were harvested from pre- and postpubertal Brahman heifers (same age group). Brahman heifers () are older at puberty compared with , a productivity issue. With RNA sequencing, we identified differentially expressed (DEx) genes and important transcription factors (TF) and predicted coexpression networks. The number of DEx genes detected in the pituitary gland was 284 ( < 0.05), and was the most DEx gene (fold change = 4.12, = 0.01). The gene promotes bone mineralization through transforming growth factor-ß (TGFß) signaling. Further studies of the link between bone mineralization and puberty could target . In ovaries, 3,871 genes were DEx ( < 0.05). Four highly DEx genes were noteworthy for their function: (a γ-aminobutyric acid [GABA] transporter), (), and () and its receptor . These genes had higher ovarian expression in postpubertal heifers. The GABA and its receptors and transporters were expressed in the ovaries of many mammals, suggesting a role for this pathway beyond the brain. The pathway has been known to influence the timing of puberty in rats, via modulation of GnRH. The effects of at the hypothalamus, pituitary gland, and ovaries have been documented. and its receptors are known factors in the release of GnRH, similar to and GABA, although their roles in ovarian tissue are less clear. Pathways previously related to puberty such as TGFß signaling ( = 6.71 × 10), Wnt signaling ( = 4.1 × 10), and peroxisome proliferator-activated receptor (PPAR) signaling ( = 4.84 × 10) were enriched in our data set. Seven genes were identified as key TF in both tissues: , , , , , , and a novel gene. An ovarian subnetwork created with TF and significant ovarian DEx genes revealed five zinc fingers as regulators: , , , , and . Recent work of hypothalamic gene expression also pointed to zinc fingers as TF for bovine puberty. Although some zinc fingers may be ubiquitously expressed, the identification of DEx genes in common across tissues points to key regulators of puberty. The hypothalamus and pituitary gland had eight DEx genes in common. The hypothalamus and ovaries had 89 DEx genes in common. The pituitary gland and ovaries had 48 DEx genes in common. Our study confirmed the complexity of puberty and suggested further investigation on genes that code zinc fingers.

Genomic analyses of tropical beef cattle fertility based on genotyping pools of Brahman cows with unknown pedigree.

We introduce an innovative approach to lowering the overall cost of obtaining genomic EBV (GEBV) and encourage their use in commercial extensive herds of Brahman beef cattle. In our approach, the DNA genotyping of cow herds from 2 independent properties was performed using a high-density bovine SNP chip on DNA from pooled blood samples, grouped according to the result of a pregnancy test following their first and second joining opportunities. For the DNA pooling strategy, 15 to 28 blood samples from the same phenotype and contemporary group were allocated to pools. Across the 2 properties, a total of 183 pools were created representing 4,164 cows. In addition, blood samples from 309 bulls from the same properties were also taken. After genotyping and quality control, 74,584 remaining SNP were used for analyses. Pools and individual DNA samples were related by means of a "hybrid" genomic relationship matrix. The pooled genotyping analysis of 2 large and independent commercial populations of tropical beef cattle was able to recover significant and plausible associations between SNP and pregnancy test outcome. We discuss 24 SNP with significant association ( < 1.0 × 10) and mapped within 40 kb of an annotated gene. We have established a method to estimate the GEBV in young herd bulls for a trait that is currently unable to be predicted at all. In summary, our novel approach allowed us to conduct genomic analyses of fertility in 2 large commercial Brahman herds managed under extensive pastoral conditions.

Transcriptome analyses identify five transcription factors differentially expressed in the hypothalamus of post- versus prepubertal Brahman heifers.

Puberty onset is a developmental process influenced by genetic determinants, environment, and nutrition. Mutations and regulatory gene networks constitute the molecular basis for the genetic determinants of puberty onset. The emerging knowledge of these genetic determinants presents opportunities for innovation in the breeding of early pubertal cattle. This paper presents new data on hypothalamic gene expression related to puberty in (Brahman) in age- and weight-matched heifers. Six postpubertal heifers were compared with 6 prepubertal heifers using whole-genome RNA sequencing methodology for quantification of global gene expression in the hypothalamus. Five transcription factors (TF) with potential regulatory roles in the hypothalamus were identified in this experiment: , , , , and . These TF genes were significantly differentially expressed in the hypothalamus of postpubertal versus prepubertal heifers and were also identified as significant according to the applied regulatory impact factor metric ( < 0.05). Two of these 5 TF, and , were zinc fingers, belonging to a gene family previously reported to have a central regulatory role in mammalian puberty. The gene belongs to the family of homologues of Drosophila sine oculis () genes implicated in transcriptional regulation of gonadotrope gene expression. Tumor-related genes such as and are known to affect basic cellular processes that are relevant in both cancer and developmental processes. Mutations in were associated with puberty in humans. Mutations in these TF, together with other genetic determinants previously discovered, could be used in genomic selection to predict the genetic merit of cattle (i.e., the likelihood of the offspring presenting earlier than average puberty for Brahman). Knowledge of key mutations involved in genetic traits is an advantage for genomic prediction because it can increase its accuracy.

Genome-wide association for the outcome of fixed-time artificial insemination of Brahman heifers in northern Australia.

Fixed-time AI (FTAI) is a powerful tool for genetic improvement of extensively managed beef cattle. A genomewide association study (GWAS) was conducted to investigate genes and genetic markers associated with the outcome (pregnant or not pregnant) of FTAI in 614 commercial Brahman heifers genotyped for 18,895 SNP and imputed to 51,588 SNP. The likelihood of Brahman heifers becoming pregnant after hormonal treatment to synchronize ovulation followed by FTAI was influenced by the content of their genomes, as determined by a principal component analysis. The principal component analysis involved comparisons between the studied heifers and populations of known and ancestry. The heritability of FTAI outcome was = 0.18, which is higher than for most other reproductive outcome traits. The number of SNP associated with FTAI outcome was 101 ( < 0.001, false discovery rate = 0.53). Compared with all SNP tested, associated SNP had a tendency for highly divergent allelic frequencies between and . Associated SNP were located in nearly all chromosomes, a result that shows a complex genetic architecture that is typical of highly complex traits with low heritability. Considering this and previous GWAS that examined Brahman heifer puberty and postpartum anestrus interval, 3 genomic regions emerge as important for overall Brahman heifer fertility, which mapped to chromosomes 1, 7, and 9. Further analyses, including improved genome annotation, are required to elucidate the link between these regions and heifer fertility. Additional studies are needed to confirm SNP and gene associations reported herein and further elucidate the genetics of FTAI outcome. Future GWAS should target other Braham populations and additional cattle breeds with FTAI records, including breeds with higher ancestry.

A marker-derived gene network reveals the regulatory role of PPARGC1A, HNF4G, and FOXP3 in intramuscular fat deposition of beef cattle.

High intramuscular fat (IMF) awards price premiums to beef producers and is associated with meat quality and flavor. Studying gene interactions and pathways that affect IMF might unveil causative physiological mechanisms and inform genomic selection, leading to increased accuracy of predictions of breeding value. To study gene interactions and pathways, a gene network was derived from genetic markers associated with direct measures of IMF, other fat phenotypes, feedlot performance, and a number of meat quality traits relating to body conformation, development, and metabolism that might be plausibly expected to interact with IMF biology. Marker associations were inferred from genomewide association studies (GWAS) based on high density genotypes and 29 traits measured on 10,181 beef cattle animals from 3 breed types. For the network inference, SNP pairs were assessed according to the strength of the correlation between their additive association effects across the 29 traits. The co-association inferred network was formed by 2,434 genes connected by 28,283 edges. Topological network parameters suggested a highly cohesive network, in which the genes are strongly functionally interconnected. Pathway and network analyses pointed towards a trio of transcription factors (TF) as key regulators of carcass IMF: PPARGC1A, HNF4G, and FOXP3. Importantly, none of these genes would have been deemed as significantly associated with IMF from the GWAS. Instead, a total of 313 network genes show significant co-association with the 3 TF. These genes belong to a wide variety of biological functions, canonical pathways, and genetic networks linked to IMF-related phenotypes. In summary, our GWAS and network predictions are supported by the current literature and suggest a cooperative role for the 3 TF and other interacting genes including CAPN6, STC2, MAP2K4, EYA1, COPS5, XKR4, NR2E1, TOX, ATF1, ASPH, TGS1, and TTPA as modulators of carcass and meat quality traits in beef cattle.

Comparative study of in vivo and in vitro phagocytosis including germicidal capacity in Odontaster validus (Koehler, 1906) at 0 degree C.

The phagocytosis and germicidal capacity of Saccharomyces cerevisiae by phagocytic amoebocytes (PA) of the Antarctic starfish Odontaster validus were studied in vivo (after incubation periods of 1, 2, and 4 h) and in vitro (after incubation periods of 1, 2, 4, 8, and 12 h) at 0 degree C. The total number of PA and the phagocytic capacity (PC), phagocytic index (PI), and germicidal capacity (GC) of the PA were calculated. Results showed significant variability of the total PA number in different animals. There was a significant increase in PC and no significant differences in PI and GC for different in vitro incubation times. In vivo, experiments showed no significant difference of PC and PI, but there was a significant increase in GC as incubation periods increased. Comparison between in vitro and in vivo results revealed that PI and PC were significantly higher in vitro and that GC was significantly higher in vivo. The present study shows for the first time the phagocytosis and GC of an Antarctic invertebrate in vivo at low temperature (0 degree C), and the results are comparing with the available literature for echinoderms.