Differentially expressed genes in the femur cartilage transcriptome clarify the understanding of femoral head separation in chickens.

Locomotor problems are among one of the main concerns in the current poultry industry, causing major economic losses and affecting animal welfare. The most common bone anomalies in the femur are dyschondroplasia, femoral head separation (FHS), and bacterial chondronecrosis with osteomyelitis (BCO), also known as femoral head necrosis (FHN). The present study aimed to identify differentially expressed (DE) genes in the articular cartilage (AC) of normal and FHS-affected broilers by RNA-Seq analysis. In the transcriptome analysis, 12,169 genes were expressed in the femur AC. Of those, 107 genes were DE (FDR < 0.05) between normal and affected chickens, of which 9 were downregulated and 98 were upregulated in the affected broilers. In the gene-set enrichment analysis using the DE genes, 79 biological processes (BP) were identified and were grouped into 12 superclusters. The main BP found were involved in the response to biotic stimulus, gas transport, cellular activation, carbohydrate-derived catabolism, multi-organism regulation, immune system, muscle contraction, multi-organism process, cytolysis, leukocytes and cell adhesion. In this study, the first transcriptome analysis of the broilers femur articular cartilage was performed, and a set of candidate genes (AvBD1, AvBD2, ANK1, EPX, ADA, RHAG) that could trigger changes in the broiler´s femoral growth plate was identified. Moreover, these results could be helpful to better understand FHN in chickens and possibly in humans.

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.

Gene Expression in the Salivary Gland of Rhipicephalus (Boophilus) microplus Fed on Tick-Susceptible and Tick-Resistant Hosts.

The success of cattle tick fixation largely depends on the secretion of substances that alter the immune response of the host. The majority of these substances are expressed by the parasite salivary gland and secreted in tick saliva. It is known that hosts can mount immune responses against ticks and bovine European breeds, and bovine industrial crossbreeds are more susceptible to infestations than are Bos indicus cattle. To identify candidates for the development of novel control strategies for the cattle tick Rhipicephalus (Boophilus) microplus, a salivary gland transcriptome analysis of engorged females fed on susceptible or resistant hosts was performed. Using RNA-Seq, transcriptomes were de novo assembled and produced a total of 235,451 contigs with 93.3% transcriptome completeness. Differential expression analysis identified 137 sequences as differentially expressed genes (DEGs) between ticks raised on tick-susceptible or tick-resistant cattle. DEGs predicted to be secreted proteins include innexins, which are transmembrane proteins that form gap junction channels; the transporters Na+/dicarboxylate, Na+/tricarboxylate, and phosphate transporter and a putative monocarboxylate transporter; a phosphoinositol 4-phosphate adaptor protein; a cysteine-rich protein containing a trypsin inhibitor-like (TIL) domain; a putative defense protein 3 containing a reeler domain; and an F-actin-uncapping protein LRRC16A with a CARMIL_C domain; these genes were upregulated in ticks fed on tick-susceptible cattle. DEGs predicted to be non-secreted proteins included a small heat shock protein and the negative elongation factor B-like, both acting in a coordinated manner to increase HSP transcript levels in the salivary glands of the ticks fed on tick-susceptible cattle; the 26S protease regulatory subunit 6B and another chaperone with similarity to calnexin, also upregulated in ticks fed on tick-susceptible cattle; an EF-hand calcium binding protein and a serine carboxypeptidase (SCP), both involved in the blood coagulation cascade and upregulated in ticks fed on tick-susceptible cattle; and two ribosomal proteins, the 60S acidic ribosomal protein P2 and the 60S ribosomal protein L19. These results help to characterize cattle tick salivary gland gene expression in tick-susceptible and tick-resistant hosts and suggest new putative targets for the control of tick infestations, as those genes involved in the mechanism of stress response during blood feeding.

Predominance of syntrophic bacteria, Methanosaeta and Methanoculleus in a two-stage up-flow anaerobic sludge blanket reactor treating coffee processing wastewater at high organic loading rate.

The effect of the organic loading rate (OLR) on the performance and microbial composition of a two-stage UASB system treating coffee processing wastewater was assessed. The system was operated with OLR up to 18.2 g COD (L d)-1 and effluent recirculation. Methane production and effluent characteristics were monitored. The microbial composition was examined through next-generation sequencing and qPCR from the anaerobic sludge of the first reactor (R1) operated at low and high OLR. The system showed operational stability, obtaining a maximum methane production of 2.2 L CH4 (L d)-1, with a removal efficiency of COD and phenolic compounds of 84 and 73%, respectively. The performance of R1 at high OLR in steady conditions was associated with an appropriate proportion of nutrients (particularly Fe) and a marked increase of the syntrophic bacteria Syntrophus and Candidatus Cloacimonas, and acetoclastic and hydrogenotrophic methanogens, mainly Methanosaeta, Methanoculleus, Methanobacterium and Methanomassiliicoccus.

Advances in tick vaccinology in Brazil: from gene expression to immunoprotection.

Rhipicephalus (Boophilus) microplus has substantial economic impact on the cattle breeding industry and, chemical control and tick resistance development are the major concern. There is a worldwide search for new options, and control using vaccines has been the main focus nowadays. Studies performed in Brazil found that Bm86-based immunization of bovines reduced the infestation of R. (B.) microplus of vaccinated bovines by 45% to 60%. Native Boophilusmicroplus tripsin inhibitors (BmTIs) with trypsin-, kallikrein-, and elastase-inhibiting activities have been used as immunogens in bovines reaching 72.8.% of efficacy. The reverse vaccinology approach has also been used for antigen search using transcriptome analysis to identify and characterize potential antigens. Study has generated more than 600 million sequences using RNA-seq of larvae, nymphs, salivary glands, intestines, and ovaries of the tick R. (B) microplus. Based on the set of transcripts obtained using this strategy, a total of 20,326 protein sequences have been identified. A pipeline analysis built in house identified the protein sequences that were most likely to be immunogenic based on the overall structural characteristic analysis.

Complete mitochondrial genome from South American catfish Pseudoplatystoma reticulatum (Eigenmann & Eigenmann) and its impact in Siluriformes phylogenetic tree.

The cachara (Pseudoplatystoma reticulatum) is a Neotropical freshwater catfish from family Pimelodidae (Siluriformes) native to Brazil. The species is of relative economic importance for local aquaculture production and basic biological information is under development to help boost efforts to domesticate and raise the species in commercial systems. The complete cachara mitochondrial genome was obtained by assembling Illumina RNA-seq data from pooled samples. The full mitogenome was found to be 16,576 bp in length, showing the same basic structure, order, and genetic organization observed in other Pimelodidae, with 13 protein-coding genes, 2 rNA genes, 22 trNAs, and a control region. Observed base composition was 24.63% T, 28.47% C, 31.45% A, and 15.44% G. With the exception of NAD6 and eight tRNAs, all of the observed mitochondrial genes were found to be coded on the H strand. A total of 107 SNPs were identified in P. reticulatum mtDNA, 67 of which were located in coding regions. Of these SNPs, 10 result in amino acid changes. Analysis of the obtained sequence with 94 publicly available full Siluriformes mitogenomes resulted in a phylogenetic tree that generally agreed with available phylogenetic proposals for the order. The first report of the complete Pseudoplatystoma reticulatum mitochondrial genome sequence revealed general gene organization, structure, content, and order similar to most vertebrates. Specific sequence and content features were observed and may have functional attributes which are now available for further investigation.

New insights from molecular characterization of the tick Rhipicephalus (Boophilus) microplus in Brazil.

The Rhipicephalus (Boophilus) microplus complex currently consists of five taxa, namely R. australis, R. annulatus, R. (B.) microplus clade A sensu, R. microplus clade B sensu, and R. (B.) microplus clade C sensu. Mitochondrial DNA-based methods help taxonomists when they are facing the morpho-taxonomic problem of distinguishing members of the R. (B.) microplus complex. The purpose of this study was to perform molecular characterization of ticks in all five regions of Brazil and infer their phylogenetic relationships. Molecular analysis characterized 10 haplotypes of the COX-1 gene. Molecular network analysis revealed that haplotype H-2 was the most dispersed of the studied populations (n = 11). Haplotype H-3 (n = 2) had the greatest genetic differentiation when compared to other Brazilian populations. A Bayesian phylogenetic tree of the COX-1 gene obtained strong support. In addition, it was observed that the population of R. (B.) microplus haplotype H-3 exhibited diverging branches among the other Brazilian populations in the study. The study concludes that the different regions of Brazil have R. (B.) microplus tick populations with distinct haplotypes.

New insights from molecular characterization of the tick Rhipicephalus (Boophilus) microplus in Brazil / Novos pontos de vista sobre a caracterização molecular em carrapatos Rhipicephalus (Boophilus) microplus, no Brasil

Abstract The Rhipicephalus (Boophilus) microplus complex currently consists of five taxa, namely R. australis, R. annulatus, R. (B.) microplus clade A sensu, R. microplus clade B sensu, and R. (B.) microplus clade C sensu. Mitochondrial DNA-based methods help taxonomists when they are facing the morpho-taxonomic problem of distinguishing members of the R. (B.) microplus complex. The purpose of this study was to perform molecular characterization of ticks in all five regions of Brazil and infer their phylogenetic relationships. Molecular analysis characterized 10 haplotypes of the COX-1 gene. Molecular network analysis revealed that haplotype H-2 was the most dispersed of the studied populations (n = 11). Haplotype H-3 (n = 2) had the greatest genetic differentiation when compared to other Brazilian populations. A Bayesian phylogenetic tree of the COX-1 gene obtained strong support. In addition, it was observed that the population of R. (B.) microplus haplotype H-3 exhibited diverging branches among the other Brazilian populations in the study. The study concludes that the different regions of Brazil have R. (B.) microplus tick populations with distinct haplotypes.

Resumo Carrapatos do complexo R. (B.) microplus se distribuem em cinco taxa: R. australis, R. annulatus, R. (B.) microplus clado A sensu R. microplus clado B sensue e R. (B.) microplus clado C sensu. Métodos baseados no DNA mitocondrial podem auxiliar taxonomistas quando há dificuldades em estabelecer diferenças morfológicas para distinguir membros do complexo R. (B.) microplus. O objetivo deste estudo foi a caracterização molecular e a inferência de relações filogenéticas em carrapatos de todas as cinco regiões geográficas do Brasil. Para o gene COX-1, a análise molecular caracterizou 10 haplótipos. Na análise molecular em rede foi observado que o haplótipo H-2 é o mais disperso entre as populações (n=11). O haplótipo H-3 (n=2) foi o que obteve maior diferenciação genética ao ser comparado com outras populações brasileiras. A árvore filogenética Bayesiana de gene COX-1 gerou suporte robusto e foi observado que a população de R. (B.) microplus haplótipo H-3 apresentou ramificação com divergência entre as outras populações brasileiras apresentadas neste estudo. Conclui-se que as populações brasileiras possuem diversidade haplotípica com divergência entre as diversas populações de R. (B.) microplus no Brasil.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits.

BACKGROUND: Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance. RESULTS: Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively. CONCLUSIONS: Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

Genome-Wide Detection of CNVs and Their Association with Meat Tenderness in Nelore Cattle.

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina® High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing ≈6.5% of the genome. Comparing our results with previous studies revealed an overlap in ≈1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.

Genomic Variants Revealed by Invariably Missing Genotypes in Nelore Cattle.

High density genotyping panels have been used in a wide range of applications. From population genetics to genome-wide association studies, this technology still offers the lowest cost and the most consistent solution for generating SNP data. However, in spite of the application, part of the generated data is always discarded from final datasets based on quality control criteria used to remove unreliable markers. Some discarded data consists of markers that failed to generate genotypes, labeled as missing genotypes. A subset of missing genotypes that occur in the whole population under study may be caused by technical issues but can also be explained by the presence of genomic variations that are in the vicinity of the assayed SNP and that prevent genotyping probes from annealing. The latter case may contain relevant information because these missing genotypes might be used to identify population-specific genomic variants. In order to assess which case is more prevalent, we used Illumina HD Bovine chip genotypes from 1,709 Nelore (Bos indicus) samples. We found 3,200 missing genotypes among the whole population. NGS re-sequencing data from 8 sires were used to verify the presence of genomic variations within their flanking regions in 81.56% of these missing genotypes. Furthermore, we discovered 3,300 novel SNPs/Indels, 31% of which are located in genes that may affect traits of importance for the genetic improvement of cattle production.

Hepatic mRNA expression and plasma levels of insulin-like growth factor-I (IGF-I) in broiler chickens selected for different growth rates

The hepatic expression and plasma concentrations of IGF-I were investigated in three broiler chicken strains selected for different growth rates (HP-Hubbard-Pettersen, a fast growing strain; NN-Naked-neck, a strain with an intermediate growth rate and a heterozygous genotype, and C-Caipira, a slow growing crossbred strain). The chickens were studied at 1, 21 and 42 days of age and had free access to food throughout the study. Hepatic IGF-I mRNA expression was assessed by dot blot analysis using a randomly labeled chicken IGF-I cDNA as the probe and plasma IGF-I concentrations were assayed by radioimmunoassay. The hepatic levels of IGF-I mRNA increased from 1 to 21 days of age in all strains, with NN chickens showing a higher (p < 0.05) IGF-I expression than the other strains. Plasma IGF-I concentrations increased (p < 0.05) with broiler chicken age, but there were no significant differences among the strains. These results indicate that despite differences in the growth rates among the strains, the changes in the expression of IGF-I mRNA in liver and in the plasma levels of IGF-I were independent of broiler chicken strain, but varied with chicken age.