The objective of this study was to perform a genome-wide association study (GWAS) to detect chromosome regions associated with indicator traits of sexual precocity in Nellore cattle. Data from Nellore animals belonging to farms which participate in the DeltaGen® and Paint® animal breeding programs, were used. The traits used in this study were the occurrence of early pregnancy (EP) and scrotal circumference (SC). Data from 72,675 females and 83,911 males with phenotypes were used; of these, 1,770 females and 1,680 males were genotyped. The SNP effects were estimated with a single-step procedure (WssGBLUP) and the observed phenotypes were used as dependent variables. All animals with available genotypes and phenotypes, in addition to those with only phenotypic information, were used. A single-trait animal model was applied to predict breeding values and the solutions of SNP effects were obtained from these breeding values. The results of GWAS are reported as the proportion of variance explained by windows with 150 adjacent SNPs. The 10 windows that explained the highest proportion of variance were identified. The results of this study indicate the polygenic nature of EP and SC, demonstrating that the indicator traits of sexual precocity studied here are probably controlled by many genes, including some of moderate effect. The 10 windows with large effects obtained for EP are located on chromosomes 5, 6, 7, 14, 18, 21 and 27, and together explained 7.91% of the total genetic variance. For SC, these windows are located on chromosomes 4, 8, 11, 13, 14, 19, 22 and 23, explaining 6.78% of total variance. GWAS permitted to identify chromosome regions associated with EP and SC. The identification of these regions contributes to a better understanding and evaluation of these traits, and permits to indicate candidate genes for future investigation of causal mutations.

BackgroundAn important goal of Zebu breeding programs is to improve reproductive performance. A major problem faced with the genetic improvement of reproductive traits is that recording the time for an animal to reach sexual maturity is costly. Another issue is that accurate estimates of breeding values are obtained only a long time after the young bulls have gone through selection. An alternative to overcome these problems is to use traits that are indicators of the reproductive efficiency of the herd and are easier to measure, such as age at first calving. Another problem is that heifers that have conceived once may fail to conceive in the next breeding season, which increases production costs. Thus, increasing heifer’s rebreeding rates should improve the economic efficiency of the herd. Response to selection for these traits tends to be slow, since they have a low heritability and phenotypic information is provided only later in the life of the animal. Genome-wide association studies (GWAS) are useful to investigate the genetic mechanisms that underlie these traits by identifying the genes and metabolic pathways involved.ResultsData from 1853 females belonging to the Agricultural Jacarezinho LTDA were used. Genotyping was performed using the BovineHD BeadChip (777 962 single nucleotide polymorphisms (SNPs)) according to the protocol of Illumina - Infinium Assay II ® Multi-Sample HiScan with the unit SQ ™ System. After quality control, 305 348 SNPs were used for GWAS. Forty-two and 19 SNPs had a Bayes factor greater than 150 for heifer rebreeding and age at first calving, respectively. All significant SNPs for age at first calving were significant for heifer rebreeding. These 42 SNPs were next or within 35 genes that were distributed over 18 chromosomes and comprised 27 protein-encoding genes, six pseudogenes and two miscellaneous noncoding RNAs.ConclusionsThe use of Bayes factor to determine the significance of SNPs allowed us to identify two sets of 42 and 19 significant SNPs for heifer rebreeding and age at first calving, respectively, which explain 11.35 % and 6.42 % of their phenotypic variance, respectively. These SNPs provide relevant information to help elucidate which genes affect these traits.

The liver plays a central role in metabolism and produces important hormones. Hepatic estrogen receptors and the release of insulin-like growth factor 1 (IGF1) are critical links between liver function and the reproductive system. However, the role of liver in pubertal development is not fully understood. To explore this question, we applied transcriptomic analyses to liver samples of pre- and post-pubertal Brahman heifers and identified differentially expressed (DE) genes and genes encoding transcription factors (TFs). Differential expression of genes suggests potential biological mechanisms and pathways linking liver function to puberty. The analyses identified 452 DE genes and 82 TF with significant contribution to differential gene expression by using a regulatory impact factor metric. Brain-derived neurotrophic factor was observed as the most down-regulated gene (P = 0.003) in post-pubertal heifers and we propose this gene influences pubertal development in Brahman heifers. Additionally, co-expression network analysis provided evidence for three TF as key regulators of liver function during pubertal development: the signal transducer and activator of transcription 6, PBX homeobox 2, and polybromo 1. Pathway enrichment analysis identified transforming growth factor-beta and Wnt signaling pathways as significant annotation terms for the list of DE genes and TF in the co-expression network. Molecular information regarding genes and pathways described in this work are important to further our understanding of puberty onset in Brahman heifers.

ABSTRACT. The objective of this study was to identify associations between known polymorphisms in genes related to adipose tissue and sexual precocity in Nellore cattle. A total of 1689 precocious and nonprecocious heifers belonging to farms participating in Conexão Delta G breeding program were studied. SNPs from the Illumina High-Density Bovine SNP BeadChip were used. This chip contains 777,000 SNPs located within the region of the candidate genes at a distance of up to 5 kb, considering that linkage disequilibrium (LD) exists at this distance. Linear models were used for statistical analysis. The fastPHASE and GenomeStudio programs were used for haplotype reconstruction and LD analysis based on r 2 statistics. Fifty-seven candidate genes and 443 SNPs were analyzed: among the latter, 370 SNPs formed 83 haplotypes, while the remaining SNPs were studied separately. Statistical analysis showed that only three haplotypes, one haplotype consisting of two SNPs located in the FABP4 gene and two haplotypes consisting of four and Lipid metabolism-related genes as candidates of precocity two SNPs located in the PPP3CA gene, had a significant effect on sexual precocity at P < 0.05. It can be concluded that the FABP4 and PPP3CA genes influence sexual precocity and may therefore be used in selection programs designed to improve sexual precocity in Nellore cattle.

ABSTRACT. Fertility traits, such as heifer pregnancy, are economically important in cattle production systems, and are therefore, used in genetic selection programs. The aim of this study was to identify single nucleotide polymorphisms (SNPs) using RNA-sequencing (RNA-Seq) data from ovary, uterus, endometrium, pituitary gland, hypothalamus, liver, longissimus dorsi muscle, and adipose tissue in 62 candidate genes associated with heifer puberty in cattle. RNA-Seq reads were assembled to the bovine reference genome (UMD 3.1.1) and analyzed in five cattle breeds; Brangus, Brahman, Nellore, Angus, and Holstein. Two approaches used the Brangus data for SNP discovery 1) pooling all samples, and 2) within each individual sample. These approaches revealed 1157 SNPs. These were compared with those identified in the pooled samples of the other breeds. Overall, 172 SNPs within 13 genes (CPNE5, FAM19A4, FOXN4, KLF1, LOC777593, MGC157266, NEBL, NRXN3, PEPT-1, PPP3CA, SCG5, TSG101, and TSHR) were concordant in the five breeds. Using Ensembl's Variant Effector Predictor, we determined that 12% of SNPs were in exons (71% synonymous, 29% nonsynonymous), 1% were in untranslated regions (UTRs), 86% were in introns, and 1% were in intergenic regions. Since these SNPs were discovered in RNA, the variants were predicted to be within exons or UTRs. Overall, 160 novel transcripts in 42 candidate genes and five novel genes overlapping five candidate genes were observed. In conclusion, 1157 SNPs were identified in 62 candidate genes associated with puberty in Brangus cattle, of which, 172 were concordant in the five cattle breeds. Novel transcripts and genes were also identified.

BackgroundQTL mapping through genome-wide association studies (GWAS) is challenging, especially in the case of low heritability complex traits and when few animals possess genotypic and phenotypic information. When most of the phenotypic information is from non-genotyped animals, GWAS can be performed using the weighted single-step GBLUP (WssGBLUP) method, which permits to combine all available information, even that of non-genotyped animals. However, it is not clear to what extent phenotypic information from non-genotyped animals increases the power of QTL detection, and whether factors such as the extent of linkage disequilibrium (LD) in the population and weighting SNPs in WssGBLUP affect the importance of using information from non-genotyped animals in GWAS. These questions were investigated in this study using real and simulated data.ResultsAnalysis of real data showed that the use of phenotypes of non-genotyped animals affected SNP effect estimates and, consequently, QTL mapping. Despite some coincidence, the most important genomic regions identified by the analyses, either using or ignoring phenotypes of non-genotyped animals, were not the same. The simulation results indicated that the inclusion of all available phenotypic information, even that of non-genotyped animals, tends to improve QTL detection for low heritability complex traits. For populations with low levels of LD, this trend of improvement was less pronounced. Stronger shrinkage on SNPs explaining lower variance was not necessarily associated with better QTL mapping.ConclusionsThe use of phenotypic information from non-genotyped animals in GWAS may improve the ability to detect QTL for low heritability complex traits, especially in populations in which the level of LD is high.

ABSTRACT. We evaluated the genetic association of growth traits [weight adjusted to 205 days of age (W205), 365 days of age (W365), and 550 days of age (W550); weight gain between 205 days of age and 365 days of age (WG1) and between 365 days of age and 550 days of age (WG2)] and reproductive traits [age at first calving (AFC); first calving interval (FCI)] with stayability in the herd (STAY), using Bayesian inference in linear and threshold models. We defined STAY as the probability of a cow calving three or more times before the age of 76 months, given that she had calved at least once. We assigned binary codes (0, failure; 1, success) to each female. We used a sire model for analysis and formed different contemporary groups for the investigated traits. We analyzed the results by applying a two-trait sire model that included STAY (threshold trait) and linear traits (W205, W365, W550, WG1, WG2, AFC, and FCI). We used 14957 Genetic association of growth traits with stayability ©FUNPEC-RP www.funpecrp.com.br Genetics and Molecular Research 14 (4): 14956-14966 (2015) Gibbs sampling to estimate variance components and heritabilities. In all the analyses, we found that the mean heritability estimates for STAY were of moderate magnitude (0.20-0.25). The mean heritabilities for W205, W365, W550, WG1, WG2, AFC, and FCI were 0.20, 0.23, 0.39, 0.08, 0.14, 0.12, and 0.11, respectively. We observed wide variation in the posterior distributions of genetic correlations; however, with the exception of those obtained for the reproductive traits, the mean estimates were of low magnitude. Selection for WG2 can results in favorable correlated response in STAY.

Study of pressor response to graded, increased doses of infused norepinephrine in patients with essential hypertension, their normotensive siblings, and normotensive control subjects unrelated to the patients and without a family history of hypertension indicated an increased response in the two former groups. Comparison of the dose-response curves in the three groups showed that the difference in response was due to a reduced threshold to norepinephrine in patients and their siblings and not to differences in the slopes of the dose-response curves. These alterations were not paralleled by differences in heart rate responses. (Hypertension 1990;15(suppl I):I-137-I-139)