Differential expression of miRNAs associated with pectoral myopathies in young broilers: insights from a comparative transcriptome analysis.

INTRODUCTION: White Striping (WS) and Wooden Breast (WB) pectoral myopathies are relevant disorders for contemporary broiler production worldwide. Several studies aimed to elucidate the genetic components associated with the occurrence of these myopathies. However, epigenetic factors that trigger or differentiate these two conditions are still unclear. The aim of this study was to identify miRNAs differentially expressed (DE) between normal and WS and WB-affected broilers, and to verify the possible role of these miRNAs in metabolic pathways related to the manifestation of these pectoral myopathies in 28-day-old broilers. RESULTS: Five miRNAs were DE in the WS vs control (gga-miR-375, gga-miR-200b-3p, gga-miR-429-3p, gga-miR-1769-5p, gga-miR-200a-3p), 82 between WB vs control and 62 between WB vs WS. Several known miRNAs were associated with WB, such as gga-miR-155, gga-miR-146b, gga-miR-222, gga-miR-146-5p, gga-miR- 29, gga-miR-21-5p, gga-miR-133a-3p and gga-miR-133b. Most of them had not previously been associated with the development of this myopathy in broilers. We also have predicted 17 new miRNAs expressed in the broilers pectoral muscle. DE miRNA target gene ontology analysis enriched 6 common pathways for WS and WB compared to control: autophagy, insulin signaling, FoxO signaling, endocytosis, and metabolic pathways. The WS vs control contrast had two unique pathways, ERBB signaling and the mTOR signaling, while WB vs control had 14 unique pathways, with ubiquitin-mediated proteolysis and endoplasmic reticulum protein processing being the most significant. CONCLUSIONS: We found miRNAs DE between normal broilers and those affected with breast myopathies at 28 days of age. Our results also provide novel evidence of the miRNAs role on the regulation of WS and in the differentiation of both WS and WB myopathies. Overall, our study provides insights into miRNA-mediated and pathways involved in the occurrence of WS and WB helping to better understand these chicken growth disorders in an early age. These findings can help developing new approaches to reduce these complex issues in poultry production possibly by adjustments in nutrition and management conditions. Moreover, the miRNAs and target genes associated with the initial stages of WS and WB development could be potential biomarkers to be used in selection to reduce the occurrence of these myopathies in broiler production.

The downregulation of genes encoding muscle proteins have a potential role in the development of scrotal hernia in pigs.

BACKGROUND: Testicular descent is a physiological process regulated by many factors. Eventually, disturbances in the embryological/fetal development path facilitate the occurrence of scrotal hernia, a congenital malformation characterized by the presence of intestinal portions within the scrotal sac due to the abnormal expansion of the inguinal ring. In pigs, some genes have been related to this anomaly, but the genetic mechanisms involved remain unclear. This study aimed to investigate the expression profile of a set of genes potentially involved with the manifestation of scrotal hernia in the inguinal ring tissue. METHODS AND RESULTS: Tissue samples from the inguinal ring/canal of normal and scrotal hernia-affected male pigs with approximately 30 days of age were used. Relative expression analysis was performed using qPCR to confirm the expression profile of 17 candidate genes previously identified in an RNA-Seq study. Among them, the Myosin heavy chain 1 (MYH1), Desmin (DES), and Troponin 1 (TNNI1) genes were differentially expressed between groups and had reduced levels of expression in the affected animals. These genes encode proteins involved in the formation of muscle tissue, which seems to be important for increasing the resistance of the inguinal ring to the abdominal pressure, which is essential to avoid the occurrence of scrotal hernia. CONCLUSIONS: The downregulation of muscular candidate genes in the inguinal tissue clarifies the genetic mechanisms involved with this anomaly in its primary site, providing useful information for developing strategies to control this malformation in pigs and other mammals.

A joint analysis using exome and transcriptome data identifiescandidate polymorphisms and genes involved with umbilical hernia in pigs.

BACKGROUND: Umbilical Hernia (UH) is characterized by the passage of part of the intestine through the umbilical canal forming the herniary sac. There are several potential causes that can lead to the umbilical hernia such as bacterial infections, management conditions and genetic factors. Since the genetic components involved with UH are poorly understood, this study aimed to identify polymorphisms and genes associated with the manifestation of umbilical hernia in pigs using exome and transcriptome sequencing in a case and control design. RESULTS: In the exome sequencing, 119 variants located in 58 genes were identified differing between normal and UH-affected pigs, and in the umbilical ring transcriptome, 46 variants were identified, located in 27 genes. Comparing the two methodologies, we obtained 34 concordant variants between the exome and transcriptome analyses, which were located in 17 genes, distributed in 64 biological processes (BP). Among the BP involved with UH it is possible to highlight cell adhesion, cell junction regulation, embryonic morphogenesis, ion transport, muscle contraction, within others. CONCLUSIONS: We have generated the first exome sequencing related to normal and umbilical hernia-affected pigs, which allowed us to identify several variants possibly involved with this disorder. Many of those variants present in the DNA were confirmed with the RNA-Seq results. The combination of both exome and transcriptome sequencing approaches allowed us to better understand the complex molecular mechanisms underlying UH in pigs and possibly in other mammals, including humans. Some variants found in genes and other regulatory regions are highlighted as strong candidates to the development of UH in pigs and should be further investigated.

Genome-wide detection of CNVs and their association with performance traits in broilers.

BACKGROUND: Copy number variations (CNVs) are a major type of structural genomic variants that underlie genetic architecture and phenotypic variation of complex traits, not only in humans, but also in livestock animals. We identified CNVs along the chicken genome and analyzed their association with performance traits. Genome-wide CNVs were inferred from Affymetrix® high density SNP-chip data for a broiler population. CNVs were concatenated into segments and association analyses were performed with linear mixed models considering a genomic relationship matrix, for birth weight, body weight at 21, 35, 41 and 42 days, feed intake from 35 to 41 days, feed conversion ratio from 35 to 41 days and, body weight gain from 35 to 41 days of age. RESULTS: We identified 23,214 autosomal CNVs, merged into 5042 distinct CNV regions (CNVRs), covering 12.84% of the chicken autosomal genome. One significant CNV segment was associated with BWG on GGA3 (q-value = 0.00443); one significant CNV segment was associated with BW35 (q-value = 0.00571), BW41 (q-value = 0.00180) and BW42 (q-value = 0.00130) on GGA3, and one significant CNV segment was associated with BW on GGA5 (q-value = 0.00432). All significant CNV segments were verified by qPCR, and a validation rate of 92.59% was observed. These CNV segments are located nearby genes, such as KCNJ11, MyoD1 and SOX6, known to underlie growth and development. Moreover, gene-set analyses revealed terms linked with muscle physiology, cellular processes regulation and potassium channels. CONCLUSIONS: Overall, this CNV-based GWAS study unravels potential candidate genes that may regulate performance traits in chickens. Our findings provide a foundation for future functional studies on the role of specific genes in regulating performance in chickens.

RNA-seq reveals downregulated osteochondral genes potentially related to tibia bacterial chondronecrosis with osteomyelitis in broilers.

BACKGROUND: Bacterial chondronecrosis with osteomyelitis (BCO) develops in the growth plate (GP) of the proximal femur and tibia and is initiated by damage to the less mineralized chondrocytes followed by colonization of opportunistic bacteria. This condition affects approximately 1% of all birds housed, being considered one of the major causes of lameness in fast growing broilers. Although several studies have been previously performed aiming to understand its pathogenesis, the molecular mechanisms involved with BCO remains to be elucidated. Therefore, this study aimed to generate a profile of global differential gene expression involved with BCO in the tibia of commercial broilers, through RNA sequencing analysis to identity genes and molecular pathways involved with BCO in chickens. RESULTS: Our data showed 192 differentially expressed (DE) genes: 63 upregulated and 129 downregulated in the GP of the tibia proximal epiphysis of BCO-affected broilers. Using all DE genes, six Biological Processes (BP) were associated with bone development (connective tissue development, cartilage development, skeletal system development, organ morphogenesis, system development and skeletal system morphogenesis). The analyses of the upregulated genes did not indicate any significant BP (FDR < 0.05). However, with the downregulated genes, the same BP were identified when using all DE genes in the analysis, with a total of 26 coding genes explaining BCO in the tibia: ACAN, ALDH1A2, CDH7, CHAD, CHADL, COL11A1, COMP, CSGALNACT1, CYR61, FRZB, GAL3ST1, HAPLN1, IHH, KIF26B, LECT1, LPPR1, PDE6B, RBP4A, SERINC5, SFRP1, SOX8, SOX9, TENM2, THBS1, UCHL1 and WFIKKN2. In addition, seven transcription factors were also associated to BCO: NFATC2, MAFB, HIF1A-ARNT, EWSR1-FLI1, NFIC, TCF3 and NF-KAPPAB. CONCLUSIONS: Our data show that osteochondral downregulated genes are potential molecular causes of BCO in broilers, and the bacterial process seems to be, in fact, a secondary condition. Sixteen genes responsible for bone and cartilage formation were downregulated in BCO-affected broilers being strong candidate genes to trigger this disorder.

Stable reference genes for expression studies in breast muscle of normal and white striping-affected chickens.

The normalization with proper reference genes is a crucial step to obtain accurate mRNA expression levels in quantitative PCR (qPCR) studies. Therefore, in this study, 10 reference candidate genes were evaluated to determine their stability in normal pectoralis major muscle of broilers and those counterparts affected with White Striping (WS) myopathy at 42 days age. Four different tools were used for ranking the most stable genes: GeNorm, NormFinder, BestKeeper and Comparative Ct (ΔCt), and a general ranking was performed using the RankAggreg tool to select the best reference genes among all tools. From the 10 genes evaluated in the breast muscle of broilers, 8 were amplified. Most of the algorithms/tools indicated the same two genes, RPL30 and RPL5, as the most stable in the broilers breast muscle. In addition, there was agreement among the tools for the least stable genes: MRPS27, GAPDH and RPLP1 in the broilers breast muscle. Therefore, it is interesting to note that even with different tools for evaluating gene expression, there was consensus on the most and least stable genes. These results indicate that the Ribosomal protein L30 (RPL30) and Ribosomal protein L5 (RPL5) can be recommended for accurate normalization in qPCR studies with chicken pectoralis major muscle affected with White Striping and other myopathies.

Identification of selection signatures involved in performance traits in a paternal broiler line.

BACKGROUND: Natural and artificial selection leads to changes in certain regions of the genome resulting in selection signatures that can reveal genes associated with the selected traits. Selection signatures may be identified using different methodologies, of which some are based on detecting contiguous sequences of homozygous identical-by-descent haplotypes, called runs of homozygosity (ROH), or estimating fixation index (FST) of genomic windows that indicates genetic differentiation. This study aimed to identify selection signatures in a paternal broiler TT line at generations 7th and 16th of selection and to investigate the genes annotated in these regions as well as the biological pathways involved. For such purpose, ROH and FST-based analysis were performed using whole genome sequence of twenty-eight chickens from two different generations. RESULTS: ROH analysis identified homozygous regions of short and moderate size. Analysis of ROH patterns revealed regions commonly shared among animals and changes in ROH abundance and size between the two generations. Results also suggest that whole genome sequencing (WGS) outperforms SNPchip data avoiding overestimation of ROH size and underestimation of ROH number; however, sequencing costs can limited the number of animals analyzed. FST-based analysis revealed genetic differentiation in several genomic windows. Annotation of the consensus regions of ROH and FST windows revealed new and previously identified genes associated with traits of economic interest, such as APOB, IGF1, IGFBP2, POMC, PPARG, and ZNF423. Over-representation analysis of the genes resulted in biological terms of skeletal muscle, matrilin proteins, adipose tissue, hyperglycemia, diabetes, Salmonella infections and tyrosine. CONCLUSIONS: Identification of ROH and FST-based analyses revealed selection signatures in TT line and genes that have important role in traits of economic interest. Changes in the genome of the chickens were observed between the 7th and 16th generations showing that ancient and recent selection in TT line may have acted over genomic regions affecting diseases and performance traits.

Proximal femoral head transcriptome reveals novel candidate genes related to epiphysiolysis in broiler chickens.

BACKGROUND: The proximal femoral head separation (FHS) or epiphysiolysis is a prevalent disorder affecting the chicken femur epiphysis, being considered a risk factor to infection which can cause bacterial chondronecrosis with osteomyelitis in broilers. To identify the genetic mechanisms involved in epiphysiolysis, differentially expressed (DE) genes in the femur of normal and FHS-affected broilers were identified using RNA-Seq technology. Femoral growth plate (GP) samples from 35-day-old commercial male broilers were collected from 4 healthy and 4 FHS-affected broilers. Sequencing was performed using an Illumina paired-end protocol. Differentially expressed genes were obtained using the edgeR package based on the False Discovery Rate (FDR < 0.05). RESULTS: Approximately 16 million reads/sample were generated with 2 × 100 bp paired-end reads. After data quality control, approximately 12 million reads/sample were mapped to the reference chicken genome (Galgal5). A total of 12,645 genes were expressed in the femur GP. Out of those, 314 were DE between groups, being 154 upregulated and 160 downregulated in FHS-affected broilers. In the functional analyses, several biological processes (BP) were overrepresented. Among them, those related to cell adhesion, extracellular matrix (ECM), bone development, blood circulation and lipid metabolism, which are more related to chicken growth, are possibly involved with the onset of FHS. On the other hand, BP associated to apoptosis or cell death and immune response, which were also found in our study, could be related to the consequence of the FHS. CONCLUSIONS: Genes with potential role in the epiphysiolysis were identified through the femur head transcriptome analysis, providing a better understanding of the mechanisms that regulate bone development in fast-growing chickens. In this study, we highlighted the importance of cell adhesion and extracellular matrix related genes in triggering FHS. Furthermore, we have shown new insights on the involvement of lipidemia and immune response/inflammation with FHS in broilers. Understanding the changes in the GP transcriptome might support breeding strategies to address poultry robustness and to obtain more resilient broilers.

Heat tolerance responses in a Bos taurus cattle herd raised in a Brazilian climate.

Knowledge of thermoregulatory responses in taurine cattle contribute to identification of animals most adapted to heat and productive when raised under Brazilian climate. The objectives were to verify the morphological and physiological responses related to adaptation to heat of taurine breeds raised under in Brazilian meteorological conditions in different seasons of the year and day periods, and to detect differences within and between breeds to know breed is most adapted. Measurements were made of 74 young bulls (n = 31 Angus; n = 43 Simmental) for the morphological traits: hair length (HL), number of hairs (NH), and coat thickness (CT); and for the physiological traits: respiratory rate (RR) and hair coat surface temperature (ST). The temperature-humidity index (THI) was calculated. The data were subjected to analyses of variance, cluster analysis, and principal component analysis (PCA). The THI (<74) indicates thermal comfort. In the winter, the HL and CT higher than in the spring season (P < 0.0001) in both breeds. Angus exhibited higher HL and CT (P < 0.0001). Within each breed, the animals differed from one another for HL (P < 0.0005). In the spring, CT was similar between the breeds, differing only in the winter season. Angus had higher values (P < 0.0005) of RR and lower values (P < 0.0001) of ST. Both breeds had higher (P < 0.0001) RR and ST in the afternoon. PCA showed that NH and HL better explained variation in adaptation. In general, the breeds have similar morphological responses in the hottest months, but have different physiological responses; Simmental proves to be more physiologically resistant. The afternoon was more stressful than the morning, even though the animals were in a thermal comfort zone. Measuring traits related to hair coat is sufficient for effective evaluation of adaptation, and the season affects the morphological and physiological traits of taurine cattle raised.

A genome-wide association study reveals novel genomic regions and positional candidate genes for fat deposition in broiler chickens.

BACKGROUND: Excess fat content in chickens has a negative impact on poultry production. The discovery of QTL associated with fat deposition in the carcass allows the identification of positional candidate genes (PCGs) that might regulate fat deposition and be useful for selection against excess fat content in chicken's carcass. This study aimed to estimate genomic heritability coefficients and to identify QTLs and PCGs for abdominal fat (ABF) and skin (SKIN) traits in a broiler chicken population, originated from the White Plymouth Rock and White Cornish breeds. RESULTS: ABF and SKIN are moderately heritable traits in our broiler population with estimates ranging from 0.23 to 0.33. Using a high density SNP panel (355,027 informative SNPs), we detected nine unique QTLs that were associated with these fat traits. Among these, four QTL were novel, while five have been previously reported in the literature. Thirteen PCGs were identified that might regulate fat deposition in these QTL regions: JDP2, PLCG1, HNF4A, FITM2, ADIPOR1, PTPN11, MVK, APOA1, APOA4, APOA5, ENSGALG00000000477, ENSGALG00000000483, and ENSGALG00000005043. We used sequence information from founder animals to detect 4843 SNPs in the 13 PCGs. Among those, two were classified as potentially deleterious and two as high impact SNPs. CONCLUSIONS: This study generated novel results that can contribute to a better understanding of fat deposition in chickens. The use of high density array of SNPs increases genome coverage and improves QTL resolution than would have been achieved with low density. The identified PCGs were involved in many biological processes that regulate lipid storage. The SNPs identified in the PCGs, especially those predicted as potentially deleterious and high impact, may affect fat deposition. Validation should be undertaken before using these SNPs for selection against carcass fat accumulation and to improve feed efficiency in broiler chicken production.

Detection and quantification of Babesia bovis and Babesia bigemina using different target genes.

Molecular assays have been widely used for the detection and quantification of bovine babesiosis due to their high sensitivity and specificity. However, variations in the sensitivity of pathogen detection may occur depending on the selected target gene. Thus, this study aimed to compare the detection sensitivity (DS) of Babesia bovis and B. bigemina infection levels in artificially and naturally infected cattle using quantitative PCR (qPCR) and six target genes. For B. bovis, the merozoite surface antigen genes 2b and 2c (msa-2b and msa-2c), and the mitochondrial cytochrome b gene (cybmt) were used. For B. bigemina, the genes encoding the proteins associated with rhoptry 1c (rap-1c), rap-1a, and cybmt were used. Six bovines, free of babesiosis, were artificially infected with 1 × 10-8 red blood cells infected (iRBC) with B. bovis (n = 3) or 1 × 10-6B. bigemina iRBC (n = 3). The animals were evaluated daily until parasitemia was confirmed (≥ 2.0%). The quantity of iRBC present in each animal was determined by examining blood smears. Blood samples were then subjected to DNA extraction, serial dilution, and qPCR analysis to determine the DS of each target gene. In addition, 30 calves naturally infected by Babesia spp. were also evaluated using the same six target genes. Regarding the artificial infection, B. bovis cybmt showed 25-fold higher sensitivity than the msa-2b and msa-2c genes, while the B. bigemina cybmt exhibited 5-fold and 25-fold higher sensitivity than the rap-1a and rap-1c genes, respectively. The rap-1a gene was found to be 5 times more sensitive than the rap-1c gene, while the B. bovis msa-2b and msa-2c genes exhibited similar DS. The positive frequencies of naturally infected calves for the target cybmt, msa-2b, and msa-2c genes (B. bovis) were: 100%, 33.3% and 50%, while cybmt, rap-1a, and rap-1c genes (B. bigemina) were 90%, 83.3%, and 63.3%, respectively. This study may contribute to the selection of suitable genes for molecular monitoring of bovine babesiosis. Mitochondrial genes could be considered as an alternative to improve the sensitivity of B. bovis and B. bigemina detection using qPCR.

Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep.

BACKGROUND: The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control. METHODS: Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content. RESULTS: The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep. CONCLUSIONS: The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.

A small proportion of Zebu genetic background in crossbred calves may not be enough to improve resistance against natural bovine Babesia spp. infections.

The main objective of cattle breeders in tropical and subtropical regions is to acquire animals with taurine-productive traits adapted to the broad weather range of these regions. However, one of the main challenges on using taurine genetics in these areas is the high susceptibility of these animals to tick-borne diseases. Consequently, the present study evaluated from 10 November 2021-19 April 2022, the over 13 assessments, the Babesia bovis and Babesia bigemina DNA loads and the IgG anti-B. bovis and anti-B. bigemina levels in Angus (n = 17, 100% Taurine) and Ultrablack (n = 14, ∼82% taurine and 18% Zebu) calves. Data were analyzed using a multivariate mixed model with repeated measures of the same animal including the fixed effects of evaluation, genetic group, sex, Babesia spp., and their interactions. The repeatability values were estimated from the (co)variances matrix and expressed for each species. The correlations between the DNA loads (CNlog) and IgG titers (S/P) values for the two species were also estimated using the same model. Regarding the specific IgG antibody titers for both Babesia spp., no significant differences were observed between the two genetic groups. However, for B. bovis and B. bigemina DNA loads, Ultrablack calves presented significantly higher values than Angus calves. Under the conditions evaluated in this study, our findings suggest that the low percentage of Zebu genetic in the Ultrablack breed was insufficient to improve resistance against babesiosis. Further studies must demonstrate if the low percentages of Zebu genetics in Taurine breeds can modify the susceptibility to babesiosis infections.

CYP24A1 and TRPC3 Gene Expression in Kidneys and Their Involvement in Calcium and Phosphate Metabolism in Laying Hens.

Ca and P homeostasis across the egg-laying cycle is a complex process involving absorption in the small intestine, reabsorption/excretion in the kidneys, and eggshell gland secretion. Diets with inadequate calcium and phosphorus can interfere with their absorption and digestibility, resulting in eggshell quality losses and reduced productive life, affecting egg production and welfare. A better understanding of gene expression profiles in the kidneys of laying hens during the late egg-laying period could clarify the renal role in mineral metabolism at this late stage. Therefore, the performance, egg quality and bone integrity-related traits, and expression profiles of kidney candidate genes were evaluated in 73-week-old laying hens receiving different Ca and P ratios in their diet: a high Ca/P ratio (HR, 22.43), a low ratio (LR, 6.71), and a medium ratio (MR, 11.43). The laying hens receiving the HR diet had improved egg production and eggshell quality traits compared to the other two groups. Humerus length was shorter in the HR than in the other groups. The CYP24A1 and TRPC3 genes were differentially expressed (p.adj ≤ 0.05) among the groups. Therefore, their expression profiles could be involved in calcium and phosphate transcellular transport in 73-week-old laying hens as a way to keep mineral absorption at adequate levels.

Bioinformatics and expression analyses of the Ixodes scapularis tick cystatin family.

The cystatins are inhibitors of papain- and legumain-like cysteine proteinases, classified in MEROPS subfamilies I25A-I25C. This study shows that 84 % (42/50) of tick cystatins are putatively extracellular in subfamily I25B and the rest are putatively intracellular in subfamily I25A. On the neighbor joining phylogeny guide tree, subfamily I25A members cluster together, while subfamily I25B cystatins segregate among prostriata or metastriata ticks. Two Ixodes scapularis cystatins, AAY66864 and ISCW011771 that show 50-71 % amino acid identity to metastriata tick cystatins may be linked to pathways that are common to all ticks, while ISCW000447 100 % conserved in I. ricinus is important among prostriata ticks. Likewise metastriata tick cystatins, Dermacentor variabilis-ACF35512, Rhipicephalus microplus-ACX53850, A. americanum-AEO36092, R. sanguineus-ACX53922, D. variabilis-ACF35514, R. sanguineus-ACX54033 and A. maculatum-AEO35155 that show 73-86 % amino acid identity may be essential to metastriata tick physiology. RT-PCR expression analyses revealed that I. scapularis cystatins were constitutively expressed in the salivary glands, midguts and other tissues of unfed ticks and ticks that were fed for 24-120 h, except for ISCW017861 that are restricted to the 24 h feeding time point. On the basis of mRNA expression patterns, I. scapularis cystatins, ISCW017861, ISCW011771, ISCW002215 and ISCW0024528 that are highly expressed at 24 h are likely involved in regulating early stage tick feeding events such as tick attachment onto host skin and creation of the feeding lesion. Similarly, ISCW018602, ISCW018603 and ISCW000447 that show 2-3 fold transcript increase by 120 h of feeding are likely associated with blood meal up take, while those that maintain steady state expression levels (ISCW018600, ISCW018601 and ISCW018604) during feeding may not be associated with tick feeding regulation. We discuss our findings in the context of advancing our knowledge of tick molecular biology.

Whole-Exome Sequencing Indicated New Candidate Genes Associated with Unilateral Cryptorchidism in Pigs.

INTRODUCTION: Cryptorchidism is a hereditary anomaly characterized by the incomplete descent of one or both testicles to the scrotum. One of the challenges of this anomaly is that the retained testicle maintains its endocrine function. As a consequence, cryptorchid animals produce hormone-tainted meat in comparison to castrated animals and are likely to be more aggressive. Cryptorchidism can lead to reduced animal welfare outcomes and cause economic losses. Identifying genetic markers for cryptorchidism is an essential step toward mitigating these negative outcomes and may facilitate genome manipulation to reduce the occurrence of cryptorchidism. Attempts to identify such markers have used genome-wide association studies. Using whole-exome sequencing, we aimed to identify single nucleotide polymorphisms (SNPs) in the coding regions of cryptorchid pigs and to characterize functional pathways concerning these SNPs. METHODS: DNA was extracted and sequenced from 5 healthy and 5 cryptorchid animals from the Landrace breed, using the Illumina HiSeq 2500 platform. Data were pre-processed using the SeqyClean tool and further mapped against the swine reference genome (Sus scrofa 11.1) using BWA software. GATK was used to identify polymorphisms (SNPs and InDels), which were annotated using the VEP tool. Network prediction and gene ontology enrichment analysis were conducted using the Cytoscape platform, and STRING software was used for visualization. RESULTS: A total of 63 SNPs were identified across the genes PIGB, CCPG1, COMMD9, LDLRAD3, TRIM44, MYLPF, SEPTIN, ZNF48, TIA1, FAIM2, KRT18, FBP1, FBP2, CTSL, DAPK1, DHX8, GPR179, DEPDC1B, ENSSSCG00000049573, ENSSSCG00000016384, ENSSSCG00000022657, ENSSSCG00000038825, and ENSSSCG00000001229. Using pathway enrichment analyses and network prospection, we have identified the following significant adjusted p value threshold of 0.001 involved with the biological function pathways of estrogen signaling, cytoskeleton organization, and the pentose phosphate pathway. CONCLUSION: Our data suggest the involvement of new SNPs and genes in developing cryptorchidism in pigs. However, further studies are needed to validate our results in a larger cohort population. Variations in the GPR179 gene, with implications at the protein level, may be associated with the appearance of this anomaly in the swine. Finally, we are showing that the estrogen signaling pathway may be involved in the pathophysiological mechanisms of this congenital anomaly as previously reported in GWAS.

Upregulated genes in articular cartilage may help to counteract femoral head separation in broilers with 21 days of age.

Femoral head necrosis (FHN) is one of the most common conditions in fast growing broilers, being characterized by separation of articular cartilage from epiphysis and classified as femoral head separation (FHS) or FHS with laceration (FHSL) depending on severity. Although molecular mechanisms involved with this disorder have been observed, its etiology is still unclear. Therefore, the expression of 15 candidate genes, chosen based on previous transcriptomic studies, was evaluated in the articular cartilage (AC) of normal and FHS-affected broilers at 21 days of age. Samples were collected based on the absence or presence of FHS for physical-chemical and qPCR analysis. The AvBD2, RHAG, COL28A1, ADA and ANGPTL7 were upregulated in FHS-affected broilers compared to the healthy group. These genes are involved in immune response, defense against pathogens, inflammation, cellular migration and adhesion, indicating different molecular mechanisms to control FHS progression at early age. Our results can contribute to improve the knowledge on FHN etiology in chickens and other species, such as horse and pigs that are severely affected by bone disorders.

Comprehensive Analyses of Bone and Cartilage Transcriptomes Evince Ion Transport, Inflammation and Cartilage Development-Related Genes Involved in Chickens' Femoral Head Separation.

Femoral head separation (FHS) is usually a subclinical condition characterized by the detachment of articular cartilage from the bone. In this study, a comprehensive analysis identifying shared and exclusive expression profiles, biological processes (BP) and variants related to FHS in the femoral articular cartilage and growth plate in chickens was performed through RNA sequencing analysis. Thirty-six differentially expressed (DE) genes were shared between femoral articular cartilage (AC) and growth plate (GP) tissues. Out of those, 23 genes were enriched in BP related to ion transport, translation factors and immune response. Seventy genes were DE exclusively in the AC and 288 in the GP. Among the BP of AC, the response against bacteria can be highlighted, and for the GP tissue, the processes related to chondrocyte differentiation and cartilage development stand out. When the chicken DE genes were compared to other datasets, eight genes (SLC4A1, RHAG, ANK1, MKNK2, SPTB, ADA, C7 and EPB420) were shared between chickens and humans. Furthermore, 89 variants, including missense in the SPATS2L, PRKAB1 and TRIM25 genes, were identified between groups. Therefore, those genes should be more explored to validate them as candidates to FHS/FHN in chickens and humans.

Genome-Wide Association Study for Haemonchus contortus Resistance in Morada Nova Sheep.

Among the gastrointestinal nematodes affecting sheep, Haemonchus contortus is the most prevalent and virulent, resulting in health problems and production losses. Therefore, selecting sheep resistant to H. contortus is a suitable and sustainable strategy for controlling endoparasites in flocks. Here, 287 lambs of the native Brazilian Morada Nova hair sheep breed were subjected to two consecutive artificial infections with H. contortus and assessed for fecal egg count (FEC), packed cell volume (PCV), and live weight (LW). Forty-four animals ranked as having extreme resistance phenotypes were genotyped using the Illumina OvineSNP50v3 chip. A case−control genome-wide association study (GWAS) detected 37 significant (p < 0.001) markers in 12 ovine chromosomes in regions harboring quantitative trait loci (QTL) for FEC, Trichostrongylus spp. adults and larvae, weight, and fat; and candidate genes for immune responses, mucins, hematological parameters, homeostasis, and growth. Four single-nucleotide polymorphisms (SNP; OAR1_rs427671974, OAR2_rs419988472, OAR5_rs424070217, and OAR17_rs401006318) genotyped by qPCR followed by high-resolution melting (HRM) were associated with FEC and LW. Therefore, molecular markers detected by GWAS for H. contortus resistance in Morada Nova sheep may support animal selection programs aimed at controlling gastrointestinal nematode infections in flocks. Furthermore, genotyping of candidate genes using HRM qPCR may provide a rapid and efficient tool for animal identification.

Downregulation of growth plate genes involved with the onset of femoral head separation in young broilers.

Femoral head separation (FHS) is characterized by the detachment of growth plate (GP) and articular cartilage, occurring in tibia and femur. However, the molecular mechanisms involved with this condition are not completely understood. Therefore, genes and biological processes (BP) involved with FHS were identified in 21-day-old broilers through RNA sequencing of the femoral GP. 13,487 genes were expressed in the chicken femoral head transcriptome of normal and FHS-affected broilers. From those, 34 were differentially expressed (DE; FDR ≤0.05) between groups, where all of them were downregulated in FHS-affected broilers. The main BP were enriched in receptor signaling pathways, ossification, bone mineralization and formation, skeletal morphogenesis, and vascularization. RNA-Seq datasets comparison of normal and FHS-affected broilers with 21, 35 and 42 days of age has shown three shared DE genes (FBN2, C1QTNF8, and XYLT1) in GP among ages. Twelve genes were exclusively DE at 21 days, where 10 have already been characterized (SHISA3, FNDC1, ANGPTL7, LEPR, ENSGALG00000049529, OXTR, ENSGALG00000045154, COL16A1, RASD2, BOC, GDF10, and THSD7B). Twelve SNPs were associated with FHS (p < 0.0001). Out of those, 5 were novel and 7 were existing variants located in 7 genes (RARS, TFPI2, TTI1, MAP4K3, LINK54, and AREL1). We have shown that genes related to chondrogenesis and bone differentiation were downregulated in the GP of FHS-affected young broilers. Therefore, these findings evince that candidate genes pointed out in our study are probably related to the onset of FHS in broilers.