RNA sequencing differential gene expression analysis of isolated perfused bovine udders experimentally inoculated with Streptococcus agalactiae.

The aim of this study was to elucidate the differential gene expression in the RNA sequencing transcriptome of isolated perfused udders collected from 4 slaughtered Holstein × Zebu crossbred dairy cows experimentally inoculated with Streptococcus agalactiae. We studied 3 different statistical tools (edgeR, baySeq, and Cuffdiff 2). In summary, 2 quarters of each udder were experimentally inoculated with Strep. agalactiae and the other 2 were used as a control. Mammary tissue biopsies were collected at times 0 and 3 h after infection. The total RNA was extracted and sequenced on an Illumina HiSeq 2000 (Illumina Inc., San Diego, CA). Transcripts were assembled from the reads aligned to the bovine UMD 3.1 reference genome, and the statistical analyses were performed using the previously mentioned tools (edgeR, baySeq, and Cuffdiff 2). Finally, the identified genes were submitted to pathway enrichment analysis. A total of 1,756, 1,161, and 3,389 genes with differential gene expression were identified when using edgeR, baySeq, and Cuffdiff 2, respectively. A total of 122 genes were identified by the overlapping of the 3 methods; however, only the platelet activation presented a significantly enriched pathway. From the results, we suggest the FCER1G, GNAI2, ORAI1, and VASP genes shared among the 3 methods in this pathway for posterior biological validation.

Relationship of runs of homozygosity with adaptive and production traits in a paternal broiler line.

Genomic regions under high selective pressure present specific runs of homozygosity (ROH), which provide valuable information on the genetic mechanisms underlying the adaptation to environment imposed challenges. In broiler chickens, the adaptation to conventional production systems in tropical environments lead the animals with favorable genotypes to be naturally selected, increasing the frequency of these alleles in the next generations. In this study, ~1400 chickens from a paternal broiler line were genotyped with the 600 K Affymetrix® Axiom® high-density (HD) genotyping array for estimation of linkage disequilibrium (LD), effective population size (N e ), inbreeding and ROH. The average LD between adjacent single nucleotide polymorphisms (SNPs) in all autosomes was 0.37, and the LD decay was higher in microchromosomes followed by intermediate and macrochromosomes. The N e of the ancestral population was high and declined over time maintaining a sufficient number of animals to keep the inbreeding coefficient of this population at low levels. The ROH analysis revealed genomic regions that harbor genes associated with homeostasis maintenance and immune system mechanisms, which may have been selected in response to heat stress. Our results give a comprehensive insight into the relationship between shared ROH regions and putative regions related to survival and production traits in a paternal broiler line selected for over 20 years. These findings contribute to the understanding of the effects of environmental and artificial selection in shaping the distribution of functional variants in the chicken genome.