Phenotypic and genetic parameters compared during repeated measures of longissimus muscle area and subcutaneous fat thickness in Nelore cattle

Real-time ultrasound is currently used for in vivo carcass evaluation of beef cattle. We estimated heritability and repeatability coefficients for ultrasound-measured carcass traits in Nelore cattle. We measured longissimus muscle area, backfat thickness, and rump fat thickness in male and female animals ranging in age from 10 to 26 months. The variance components were estimated by single-trait analysis using the derivative-free restricted maximum likelihood method, under an animal model. The complete data file included 3077 records of 1172 animals born between 2003 and 2008. Two analyses were performed using a repeatability model: a) all records of longissimus muscle area, backfat thickness, and rump fat thickness of animals aged 10 to 26 months, and b) considering two age periods (11 to 17 and 18 to 24 months). The heritability estimates ranged from 0.33 to 0.46 for longissimus muscle area, from 0.20 to 0.26 for backfat, and from 0.26 to 0.29 for rump fat, demonstrating that direct selection for these traits can provide genetic gain. The repeatability estimates were moderate, ranging from 0.42 to 0.73. The highest repeatability estimates were obtained for longissimus muscle area in both the first (0.70) and second (0.73) analyses. The repeatability estimates suggest that ultrasound measures are precise; higher precision was obtained for longissimus muscle area than for subcutaneous fat thickness. The latter could be obtained during periods of greater feed availability to increase precision.     

Genome-wide association study for growth traits in Nelore cattle

The objective of this study was to investigate the association of single nucleotide polymorphisms (SNPs) with birth weight, weight gain from birth to weaning and from weaning to yearling, yearling height and cow weight in Nelore cattle. Data from 5064 animals participating in the DeltaGen and PAINT breeding programs were used. The animals were genotyped with a panel of 777 962 SNPs (Illumina BovineHD BeadChip) and 412 993 SNPs remained after quality control analysis of the genomic data. A genome-wide association study was performed using a single-step methodology. The analyses were processed with the BLUPF90 family of programs. When applied to a genome-wide association studies, the single-step GBLUP methodology is an iterative process that estimates weights for the SNPs. The weights of SNPs were included in all analyses by iteratively applying the single-step GBLUP methodology and repeated twice so that the effect of the SNP and the effect of the animal were recalculated in order to increase the weight of SNPs with large effects and to reduce the weight of those with small effects. The genome-wide association results are reported based on the proportion of variance explained by windows of 50 adjacent SNPs. Considering the two iterations, only windows with an additive genetic variance >1.5% were presented in the results. Associations were observed with birth weight on BTA 14, with weight gain from birth to weaning on BTA 5 and 29, with weight gain from weaning to yearling on BTA 11, and with yearling height on BTA 8, showing the genes TMEM68 (transmembrane protein 8B) associated with birth weight and yearling height, XKR4 (XK, Kell blood group complex subunit-related family, member 4) associated with birth weight, NPR2 (natriuretic peptide receptor B) associated with yearling height, and REG3G (regenerating islet-derived 3-gamma) associated with weight gain from weaning to yearling. These genes play an important role in feed intake, weight gain and the regulation of skeletal growth.     

Accuracy of genotype imputation in Nelore cattle

Background

Genotype imputation from low-density (LD) to high-density single nucleotide polymorphism (SNP) chips is an important step before applying genomic selection, since denser chips tend to provide more reliable genomic predictions. Imputation methods rely partially on linkage disequilibrium between markers to infer unobserved genotypes. Bos indicus cattle (e.g. Nelore breed) are characterized, in general, by lower levels of linkage disequilibrium between genetic markers at short distances, compared to taurine breeds. Thus, it is important to evaluate the accuracy of imputation to better define which imputation method and chip are most appropriate for genomic applications in indicine breeds.

Methods

Accuracy of genotype imputation in Nelore cattle was evaluated using different LD chips, imputation software and sets of animals. Twelve commercial and customized LD chips with densities ranging from 7 K to 75 K were tested. Customized LD chips were virtually designed taking into account minor allele frequency, linkage disequilibrium and distance between markers. Software programs FImpute and BEAGLE were applied to impute genotypes. From 995 bulls and 1247 cows that were genotyped with the Illumina® BovineHD chip (HD), 793 sires composed the reference set, and the remaining 202 younger sires and all the cows composed two separate validation sets for which genotypes were masked except for the SNPs of the LD chip that were to be tested.

Results

Imputation accuracy increased with the SNP density of the LD chip. However, the gain in accuracy with LD chips with more than 15 K SNPs was relatively small because accuracy was already high at this density. Commercial and customized LD chips with equivalent densities presented similar results. FImpute outperformed BEAGLE for all LD chips and validation sets. Regardless of the imputation software used, accuracy tended to increase as the relatedness between imputed and reference animals increased, especially for the 7 K chip.

Conclusions

If the Illumina® BovineHD is considered as the target chip for genomic applications in the Nelore breed, cost-effectiveness can be improved by genotyping part of the animals with a chip containing around 15 K useful SNPs and imputing their high-density missing genotypes with FImpute.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0069-1) contains supplementary material, which is available to authorized users.     

Transcriptome analysis of cattle muscle identifies potential markers for skeletal muscle growth rate and major cell types

Background

This study aimed to identify markers for muscle growth rate and the different cellular contributors to cattle muscle and to link the muscle growth rate markers to specific cell types.

Results

The expression of two groups of genes in the longissimus muscle (LM) of 48 Brahman steers of similar age, significantly enriched for “cell cycle” and “ECM (extracellular matrix) organization” Gene Ontology (GO) terms was correlated with average daily gain/kg liveweight (ADG/kg) of the animals. However, expression of the same genes was only partly related to growth rate across a time course of postnatal LM development in two cattle genotypes, Piedmontese x Hereford (high muscling) and Wagyu x Hereford (high marbling). The deposition of intramuscular fat (IMF) altered the relationship between the expression of these genes and growth rate. K-means clustering across the development time course with a large set of genes (5,596) with similar expression profiles to the ECM genes was undertaken. The locations in the clusters of published markers of different cell types in muscle were identified and used to link clusters of genes to the cell type most likely to be expressing them. Overall correspondence between published cell type expression of markers and predicted major cell types of expression in cattle LM was high. However, some exceptions were identified: expression of SOX8 previously attributed to muscle satellite cells was correlated with angiogenesis. Analysis of the clusters and cell types suggested that the “cell cycle” and “ECM” signals were from the fibro/adipogenic lineage. Significant contributions to these signals from the muscle satellite cells, angiogenic cells and adipocytes themselves were not as strongly supported. Based on the clusters and cell type markers, sets of five genes predicted to be representative of fibro/adipogenic precursors (FAPs) and endothelial cells, and/or ECM remodelling and angiogenesis were identified.

Conclusions

Gene sets and gene markers for the analysis of many of the major processes/cell populations contributing to muscle composition and growth have been proposed, enabling a consistent interpretation of gene expression datasets from cattle LM. The same gene sets are likely to be applicable in other cattle muscles and in other species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1403-x) contains supplementary material, which is available to authorized users.     

Transcriptional profiling of skeletal muscle tissue from two breeds of cattle

We used a 9.6 K cattle muscle/fat cDNA microarray to study gene expression differences between the longuissimus dorsi (LD) muscle of Japanese Black (JB) and Holstein (HOL) cattle. JB cattle exhibit an unusual ability to accumulate intramuscular adipose tissue with fat melting points lower than that in other breeds. The LD biopsies from three JB (Tajima strain) and three HOL animals were used in this breed comparison. Seventeen genes were identified as preferentially expressed in LD samples from JB and seven genes were found to be expressed more highly in HOL. The expression of six selected differentially expressed genes was confirmed by quantitative real-time PCR. The genes more highly expressed in JB are associated with unsaturated fatty acid synthesis, fat deposition, and the thyroid hormone pathway. These results are consistent with the increased amounts and proportions of monounsaturated fatty acids observed in the muscle of JB animals. By discovering as yet uncharacterized genes that are differentially regulated in this comparison, the work may lead us to a better understanding of the regulatory pathways involved in the development of intramuscular adipose tissue.     

Transcriptome profiling in crustaceans as a tool for ecotoxicogenomics

Chemicals released into the environment have the potential to affect various species and it is important to evaluate such chemical effect on ecosystems, including aquatic organisms. Among aquatic organisms, Daphnia magna has been used extensively for acute toxicity or reproductive toxicity tests. Although these types of tests can provide information on hazardous concentrations of chemicals, they provide no information on their mode of action. Recent advances in toxicogenomics, the integration of genomics with toxicology, have the potential to afford a better understanding of the responses of aquatic organisms to pollutants. In a previous study, we developed an oligonucleotide-based DNA microarray with high reproducibility using a Daphnia expressed sequence tag (EST) database. In this study, we increased the number of genes on the array and used it for a careful ecotoxicogenomic assessment of Daphnia magna. The DNA microarray was used to evaluate gene expression profiles of neonate daphnids exposed to beta-naphthoflavone (bNF). Exposure to this chemical resulted in a characteristic gene expression pattern. As the number of the genes on an array was increased, the number of genes that were found to respond to the chemicals was also increased, which made the classification of the toxic chemicals easier and more accurate. This newly developed DNA microarray can be useful for a obtaining a better mechanistic understanding of chemical toxicity effects on a common freshwater organism.     

Signatures of selection in Nelore cattle revealed by whole-genome sequencing data

Nelore cattle breed was farmed worldwide due to its economic importance in the beef market and adaptation to the tropics. In Brazil, purebred Nelore animals (PO) receive a certificate from the breeders' association based on the animal's genealogy and morphological characterization. The top 20 to 30% of the superior animals are eligible to receive the Special Certificate of Identification and Production (CEIP), meaning animals from this category were selected and evaluated in a breeding program to improve economically important traits. We used whole-genome sequencing and approaches based on haplotype differentiation and allelic differentiation to detect regions of selection signatures in Nelore cattle by comparing animals from PO and CEIP categories. From a total of 150 animals, a hierarchical clustering analysis was performed to choose the more unrelated animals from each category (16 PO and 40 CEIP). The hapFLK statistic was performed, and extensions of hapFLK values were investigated considering continuous regions with significant q-values. The Weir and Cockerham's Fst estimator (wcFst) was computed using the GPAT++ software library. The total of 82,326 SNPs with hapFLK values passed the FDR control (q-value<0.05), and 718 segments were target as signatures of selection. A total of 1713 highly differentiated genomic regions were identified based on the segmentFst approach. The signatures of selection were spread across the genome. Annotation of overlapping selection signature regions between the two methods revealed 118 genes in common. A variant located within the 3′ region of the BOLA-DRB3 gene was found as a promising candidate polymorphism. Within genomic regions that deserves attention, we found genes previously associated with adaptation to tropical environments (HELB), growth and navel size (HMGA2), fat deposition and domestication (IRAK3), and feed efficiency and postmortem carcass traits (GABRG3). The genes BOLA-DQA2, BOLA-DQB, BOLA-DQA5, BOLA-DQA1, BOLA-DRB3, ENSBTAG00000038397 on chromosome 23 are part of the Bovine Major Histocompatibility Complex (MHC) Class II gene family, representing good candidates for immune response and adaptation to tropical conditions. The BoLA family genes and the interaction of ROBO1 with SLIT genes appeared in the enrichment results. Genomic regions located in intronic regions were also identified and might play a regulatory role in traits under selection in PO and CEIP subpopulations. The regions here identified contribute to our knowledge regarding genes and variants that have an important role in complex traits selected in this breed     

Ovarian follicular dynamics in Nelore breed (Bos indicus) cattle

The most common beef cattle raised in Brazil is the Nelore breed (Bos indicus). Information obtained by ultrasonography on follicular growth in Bos taurus cattle has been accumulating rapidly. However, there are few publications to date on follicular development in Bos indicus breeds. The follicular dynamics in Nelore heifers and cows during natural or prostaglandin (PG)-induced estrous cycle were studied. From the detection of estrus onward, all animals were examined daily by ultrasonography for one (n = 35) or two (n = 10) consecutive estrous cycles. The follicular dynamic in Nelore cattle was characterized by the predominance of 2 follicular waves in the cows (83.3%, n = 18, P < 0.05) and 3 waves in the heifers (64.7%, n = 16, P < 0.05). Most of the cattle observed over 2 consecutive estrous cycles presented the same pattern of follicular waves in the first and second cycle, and only 30% showed variation in the number of waves from one cycle to the other. Most of the follicular parameters analyzed were not affected by PG treatment or age but were altered by follicular waves. Consequently, data on cows and heifers were combined according to the number of follicular waves. The ovulatory follicle was larger than the other dominant follicles (P < 0.05), and the ovulatory wave was shorter than the preceding waves (P < 0.05). The interovulatory interval was longer in animals showing 3 waves than those exhibiting 2 waves (P < 0.05). Maximum diameter of the dominant follicle (around 11 mm) and of the corpus luteum (CL, approximately 17 mm) were smaller than those reported for European breeds. In conclusion, the results demonstrate that although the dominant follicle and corpus luteum are smaller than in European breeds, the follicular dynamics in Nelore cattle were similar to those observed in European breeds and were characterized by 2 or 3 follicular waves for cows and heifers, respectively, during the natural or prostaglandin-induced estrous cycle.     

Genotype by environment interaction in Nelore cattle from five Brazilian states

Records from 75,941 Nelore cattle were used to determine the importance of genotype by environment interaction (GEI) in five Brazilian states. (Co)variance components were estimated by single-trait analysis (with yearling weight, W450, considered to be the same trait in all states) and multiple-trait analysis (with the record from each state considered to be a different trait). The direct heritability estimates for yearling weight were 0.51, 0.39, 0.44, 0.37 and 0.41 in the states of Goiás, Mato Grosso, São Paulo, Mato Grosso do Sul and Minas Gerais, respectively. The across-state genetic correlation estimates between Goiás and Mato Grosso, Goiás and Minas Gerais, São Paulo and Minas Gerais, and Mato Grosso do Sul and Minas Gerais ranged from 0.67 to 0.75. These estimates indicate that GEIs are biologically important. No interactions were observed between Goiás and São Paulo, Goiás and Mato Grosso do Sul, Mato Grosso and São Paulo, Mato Grosso and Mato Grosso do Sul, Mato Grosso and Minas Gerais, or São Paulo and Mato Grosso do Sul (0.82 to 0.97). Comparison of single and multiple-trait analyses showed that selection based on the former was less efficient in the presence of GEI, with substantial losses (up to 10%) during selection.     

Genomic prediction ability for beef fatty acid profile in Nelore cattle using different pseudo-phenotypes

The aim of the present study was to compare the predictive ability of SNP-BLUP model using different pseudo-phenotypes such as phenotype adjusted for fixed effects, estimated breeding value, and genomic estimated breeding value, using simulated and real data for beef FA profile of Nelore cattle finished in feedlot. A pedigree with phenotypes and genotypes of 10,000 animals were simulated, considering 50% of multiple sires in the pedigree. Regarding to phenotypes, two traits were simulated, one with high heritability (0.58), another with low heritability (0.13). Ten replicates were performed for each trait and results were averaged among replicates. A historical population was created from generation zero to 2020, with a constant size of 2000 animals (from generation zero to 1000) to produce different levels of linkage disequilibrium (LD). Therefore, there was a gradual reduction in the number of animals (from 2000 to 600), producing a “bottleneck effect” and consequently, genetic drift and LD starting in the generation 1001 to 2020. A total of 335,000 markers (with MAF greater or equal to 0.02) and 1000 QTL were randomly selected from the last generation of the historical population to generate genotypic data for the test population. The phenotypes were computed as the sum of the QTL effects and an error term sampled from a normal distribution with zero mean and variance equal to 0.88. For simulated data, 4000 animals of the generations 7, 8, and 9 (with genotype and phenotype) were used as training population, and 1000 animals of the last generation (10) were used as validation population. A total of 937 Nelore bulls with phenotype for fatty acid profiles (Sum of saturated, monounsaturated, omega 3, omega 6, ratio of polyunsaturated and saturated and polyunsaturated fatty acid profile) were genotyped using the Illumina BovineHD BeadChip (Illumina, San Diego, CA) with 777,962 SNP. To compare the accuracy and bias of direct genomic value (DGV) for different pseudo-phenotypes, the correlation between true breeding value (TBV) or DGV with pseudo-phenotypes and linear regression coefficient of the pseudo-phenotypes on TBV for simulated data or DGV for real data, respectively. For simulated data, the correlations between DGV and TBV for high heritability traits were higher than obtained with low heritability traits. For simulated and real data, the prediction ability was higher for GEBV than for Yc and EBV. For simulated data, the regression coefficient estimates (b_(Yc,DGV)), were on average lower than 1 for high and low heritability traits, being inflated. The results were more biased for Yc and EBV than for GEBV. For real data, the GEBV displayed less biased results compared to Yc and EBV for SFA, MUFA, n-3, n-6, and PUFA/SFA. Despite the less biased results for PUFA using the EBV as pseudo-phenotype, the b_(Yi,DGV estimates obtained for the different pseudo-phenotypes (Yc, EBV and GEBV) were very close. Genomic information can assist in improving beef fatty acid profile in Zebu cattle, since the use of genomic information yielded genomic values for fatty acid profile with accuracies ranging from low to moderate. Considering both simulated and real data, the ssGBLUP model is an appropriate alternative to obtain more reliable and less biased GEBVs as pseudo-phenotype in situations of missing pedigree, due to high proportion of multiple sires, being more adequate than EBV and Yc to predict direct genomic value for beef fatty acid profile.     

Mitochondrial abundance and function in skeletal muscle and liver from Simmental beef cattle divergent for residual feed intake

Cellular mitochondrial function has been suggested to contribute to variation in feed efficiency (FE) among animals. The objective of this study was to determine mitochondrial abundance and activities of various mitochondrial respiratory chain complexes (complex I (CI) to complex IV (CIV)) in liver and muscle tissue from beef cattle phenotypically divergent for residual feed intake (RFI), a measure of FE. Individual DM intake (DMI) and growth were measured in purebred Simmental heifers (n = 24) and bulls (n = 28) with an initial mean BW (SD) of 372 kg (39.6) and 387 kg (50.6), respectively. All animals were offered concentrates ad libitum and 3 kg of grass silage daily, and feed intake was recorded for 70 days. Residuals of the regression of DMI on average daily gain (ADG), mid-test BW^0.75 and backfat (BF), using all animals, were used to compute individual RFI coefficients. Animals were ranked within sex, by RFI into high (inefficient; top third of the population), medium (middle third of population) and low (efficient; bottom third of the population) terciles. Statistical analysis was carried out using the MIXED procedure of SAS v 9.3. Overall mean ADG (SD) and daily DMI (SD) for heifers were 1.2 (0.4) and 9.1 (0.5) kg, respectively, and for bulls were 1.8 (0.3) and 9.5 (1.02) kg, respectively. Heifers and bulls ranked as high RFI consumed 10% and 15% more (P < 0.05), respectively, than their low RFI counterparts. There was no effect of RFI on mitochondrial abundance in either liver or muscle (P > 0.05). An RFI × sex interaction was apparent for CI activity in muscle. High RFI animals had an increased activity (P < 0.05) of CIV in liver tissue compared to their low RFI counterparts; however, the relevance of that observation is not clear. Our data provide no clear evidence that cellular mitochondrial function within either skeletal muscle or hepatic tissue has an appreciable contributory role to overall variation in FE among beef cattle.     

Genetic parameter estimates and response to selection for weight and testicular traits in Nelore cattle

We estimated genetic parameters for various phases of body and testicular growth until 550 days of age in Nelore cattle, using Bayesian inference, including correlation values and error estimates. Weight and scrotal records of 54,182 Nelore animals originating from 18 farms participating in the Brazilian Nelore Breeding Program (PMGRN) were included. The following traits were measured: weight at standard ages of 120 (W120), 210 (W210), 365 (W365), 450 (W450), and 550 (W550) days; weight gain between 120/210 (WG1), 210/365 (WG2), 365/450 (WG3), 450/550 (WG4), 120/365 (WG5), 120/450 (WG6), 120/550 (WG7), 210/450 (WG8), 210/550 (WG9), and 365/550 (WG10) days of age; scrotal circumference at 365 (SC365), 450 (SC450) and 550 (SC550) days of age, and testicular growth between 365/450 (TG1), 450/550 (TG2) and 365/550 (TG3) days of age. The model included contemporary group (current farm, year and two-month period of birth, sex, and management group) and age of dam at calving, divided into classes as fixed effects. The model also included random effects for direct additive, maternal additive and maternal permanent environmental, and residual effects. The direct heritability estimates ranged from 0.23 to 0.39, 0.13 to 0.39 and 0.32 to 0.56 for weights at standard ages, weight gains and testicular measures, respectively. The genetic correlations between weights (0.69 to 0.94) and scrotal circumferences (0.91 to 0.97) measured at standard ages were higher than those between weight gain and testicular growth (0.18 to 0.97 and 0.36 to 0.77, respectively). The weights at standard ages responded more effectively to selection, and also gave strong correlations with the other traits.     

Low levels of taurine introgression in the current Brazilian Nelore and Gir indicine cattle populations

Background

Nelore and Gir are the two most important indicine cattle breeds for production of beef and milk in Brazil. Historical records state that these breeds were introduced in Brazil from the Indian subcontinent, crossed to local taurine cattle in order to quickly increase the population size, and then backcrossed to the original breeds to recover indicine adaptive and productive traits. Previous investigations based on sparse DNA markers detected taurine admixture in these breeds. High-density genome-wide analyses can provide high-resolution information on the genetic composition of current Nelore and Gir populations, estimate more precisely the levels and nature of taurine introgression, and shed light on their history and the strategies that were used to expand these breeds.

Results

We used the high-density Illumina BovineHD BeadChip with more than 777 K single nucleotide polymorphisms (SNPs) that were reduced to 697 115 after quality control filtering to investigate the structure of Nelore and Gir populations and seven other worldwide populations for comparison. Multidimensional scaling and model-based ancestry estimation clearly separated the indicine, European taurine and African taurine ancestries. The average level of taurine introgression in the autosomal genome of Nelore and Gir breeds was less than 1% but was 9% for the Brahman breed. Analyses based on the mitochondrial SNPs present in the Illumina BovineHD BeadChip did not clearly differentiate taurine and indicine haplotype groupings.

Conclusions

The low level of taurine ancestry observed for both Nelore and Gir breeds confirms the historical records of crossbreeding and supports a strong directional selection against taurine haplotypes via backcrossing. Random sampling in production herds across the country and subsequent genotyping would be useful for a more complete view of the admixture levels in the commercial Nelore and Gir populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0109-5) contains supplementary material, which is available to authorized users.     

Breeding value accuracy estimates for growth traits using random regression and multi-trait models in Nelore cattle

We quantified the potential increase in accuracy of expected breeding value for weights of Nelore cattle, from birth to mature age, using multi-trait and random regression models on Legendre polynomials and B-spline functions. A total of 87,712 weight records from 8144 females were used, recorded every three months from birth to mature age from the Nelore Brazil Program. For random regression analyses, all female weight records from birth to eight years of age (data set I) were considered. From this general data set, a subset was created (data set II), which included only nine weight records: at birth, weaning, 365 and 550 days of age, and 2, 3, 4, 5, and 6 years of age. Data set II was analyzed using random regression and multi-trait models. The model of analysis included the contemporary group as fixed effects and age of dam as a linear and quadratic covariable. In the random regression analyses, average growth trends were modeled using a cubic regression on orthogonal polynomials of age. Residual variances were modeled by a step function with five classes. Legendre polynomials of fourth and sixth order were utilized to model the direct genetic and animal permanent environmental effects, respectively, while third-order Legendre polynomials were considered for maternal genetic and maternal permanent environmental effects. Quadratic polynomials were applied to model all random effects in random regression models on B-spline functions. Direct genetic and animal permanent environmental effects were modeled using three segments or five coefficients, and genetic maternal and maternal permanent environmental effects were modeled with one segment or three coefficients in the random regression models on B-spline functions. For both data sets (I and II), animals ranked differently according to expected breeding value obtained by random regression or multi-trait models. With random regression models, the highest gains in accuracy were obtained at ages with a low number of weight records. The results indicate that random regression models provide more accurate expected breeding values than the traditionally finite multi-trait models. Thus, higher genetic responses are expected for beef cattle growth traits by replacing a multi-trait model with random regression models for genetic evaluation. B-spline functions could be applied as an alternative to Legendre polynomials to model covariance functions for weights from birth to mature age.     

EEF1A1 transcription cofactor gene polymorphism is associated with muscle gene expression and residual feed intake in Nelore cattle

Mammalian Genome - Cis-acting effects of noncoding variants on gene expression and regulatory molecules constitute a significant factor for phenotypic variation in complex traits. To provide new...