Genome-wide association study identifies genetic variants underlying footrot in Portuguese Merino sheep.

BACKGROUND: Ovine footrot caused by Dichelobacter nodosus (D. nodosus) is a contagious disease with serious economic and welfare impacts in sheep production systems worldwide. A better understanding of the host genetic architecture regarding footrot resistance/susceptibility is crucial to develop disease control strategies that efficiently reduce infection and its severity. A genome-wide association study was performed using a customized SNP array (47,779 SNPs in total) to identify genetic variants associated to footrot resistance/susceptibility in two Portuguese native breeds, i.e. Merino Branco and Merino Preto, and a population of crossbred animals. A cohort of 1375 sheep sampled across 17 flocks, located in the Alentejo region (southern Portugal), was included in the analyses. RESULTS: Phenotypes were scored from 0 (healthy) to 5 (severe footrot) based on visual inspection of feet lesions, following the Modified Egerton System. Using a linear mixed model approach, three SNPs located on chromosome 24 reached genome-wide significance after a Bonferroni correction (p < 0.05). Additionally, six genome-wide suggestive SNPs were identified each on chromosomes 2, 4, 7, 8, 9 and 15. The annotation and KEGG pathway analyses showed that these SNPs are located within regions of candidate genes such as the nonsense mediated mRNA decay associated PI3K related kinase (SMG1) (chromosome 24) and the RALY RNA binding protein like (RALYL) (chromosome 9), both involved in immunity, and the heparan sulfate proteoglycan 2 (HSPG2) (chromosome 2) and the Thrombospodin 1 (THBS1) (chromosome 7) implicated in tissue repair and wound healing processes. CONCLUSION: This is the first attempt to identify molecular markers associated with footrot in Portuguese Merino sheep. These findings provide relevant information on a likely genetic association underlying footrot resistance/susceptibility and the potential candidate genes affecting this trait. Genetic selection strategies assisted on the information obtained from this study could enhance Merino sheep-breeding programs, in combination with farm management strategies, for a more effective and sustainable long-term solution for footrot control.

A metagenomics approach to characterize the footrot microbiome in Merino sheep.

In the Portuguese Alentejo region, Merino sheep breed is the most common breed, reared for the production of meat, dairy, and wool. Footrot is responsible for lameness, decreased animal welfare, and higher production losses, generating a negative economic impact. The disease is caused by Dichelobacter nodosus that interacts with the sheep foot microbiome, to date largely uncharacterized. In fact, Dichelobacter nodosus is not able to induce footrot by itself being required the presence of a second pathogen known as Fusobacterium necrophorum. To understand and characterize the footrot microbiome dynamics of different footrot lesion scores, a whole metagenome sequencing (WMGS) approach was used. Foot tissue samples were collected from 212 animals with different degrees of footrot lesion scores, ranging from 0 to 5. Distinct bacterial communities were associated with feet with different footrot scores identifying a total of 63 phyla and 504 families. As the severity of footrot infection increases the microorganisms' diversity decreases triggering a shift in the composition of the microbiome from a dominant gram-positive in mild stages to a dominant gram-negative in the severe stages. Several species previously associated with footrot and other polymicrobial diseases affecting the epidermis and provoking inflammatory responses such as Treponema spp., Staphylococcus spp., Streptococcus spp. and Campylobacter spp. were identified proliferating along with the lesions' severity. Although these bacteria are not able to initiate footrot, several evidences have been described supporting their association with the severity and incidence increase of footrot lesions caused by Dichelobacter nodosus and Fusobacterium necrophorum. Further investigation is required to establish the roles of particular taxa and identify which of them play a role in the disease process and which are opportunistic pathogens.

Genome-wide assessment of the population structure and genetic diversity of four Portuguese native sheep breeds.

As the effects of global warming become increasingly complex and difficult to manage, the conservation and sustainable use of locally adapted sheep breeds are gaining ground. Portuguese native sheep breeds are important reservoirs of genetic diversity, highly adapted to harsh environments and reared in low input production systems. Genomic data that would describe the breeds in detail and accelerate the selection of more resilient animals to be able to cope with climatic challenges are still lacking. Here, we sequenced the genomes of 37 animals from four Portuguese native sheep breeds (Campaniça, Bordaleira Serra da Estrela, Merino Branco and Merino Preto) and 19 crossbred sheep to make inferences on their genomic diversity and population structure. Mean genomic diversities were very similar across these breeds (.30 ≤ Ho ≤ .34; .30 ≤ He ≤ .35; 1.7 × 10-3 ≤ π ≤ 3.1 × 10-3) and the levels of inbreeding were negligible (.005 ≤ FIS ≤ .038). The Principal Components, Bayesian clustering and Treemix analyses split the Portuguese breeds in two main groups which are consistent with historical records: one comprising Campaniça and Serra da Estrela together with other European and transboundary dairy breeds; and another of the well-differentiated multi-purpose Merino and Merino-related breeds. Runs of homozygosity analyses yielded 1,690 ROH segments covering an average of 2.27 Gb across the genome in all individuals. The overall genome covered by ROH segments varied from 27,75 Mb in Serra da Estrela to 61,29 Mb in Campaniça. The phylogenetic analysis of sheep mitogenomes grouped the Portuguese native breeds within sub-haplogroup B1a along with two animals of the Akkaraman breed from Turkey. This result provides additional support to a direct influence of Southwest Asian sheep in local breeds from the Iberian Peninsula. Our study is a first step pertaining to the genomic characterization of Portuguese sheep breeds and the results emphasize the potential of genomic data as a valid tool to guide conservation efforts in locally adapted sheep breeds. In addition, the genomic data we generated can be used to identify markers for breed assignment and traceability of certified breed-products.

ChIP-Seq reveals that QsMYB1 directly targets genes involved in lignin and suberin biosynthesis pathways in cork oak (Quercus suber).

BACKGROUND: Gene activity is largely controlled by transcriptional regulation through the action of transcription factors and other regulators. QsMYB1 is a member of the R2R3-MYB transcription factor family related to secondary growth, and in particular, with the cork development process. In order to identify the putative gene targets of QsMYB1 across the cork oak genome we developed a ChIP-Seq strategy. RESULTS: Results provide direct evidence that QsMY1B targets genes encoding for enzymes involved in the lignin and suberin pathways as well as gene encoding for ABCG transporters and LTPs implicated in the transport of monomeric suberin units across the cellular membrane. These results highlight the role of QsMYB1 as a regulator of lignin and suberin biosynthesis, transport and assembly. CONCLUSION: To our knowledge, this work constitutes the first ChIP-Seq experiment performed in cork oak, a non-model plant species with a long-life cycle, and these results will contribute to deepen the knowledge about the molecular mechanisms of cork formation and differentiation.

Characterization of the cork formation and production transcriptome in Quercus cerris × suber hybrids.

Cork oak is the main cork-producing species worldwide, and plays a significant economic, ecological and social role in the Mediterranean countries, in particular in Portugal and Spain. The ability to produce cork is limited to a few species, hence it must involve specific regulation mechanisms that are unique to these species. However, to date, these mechanisms remain largely understudied, especially with approaches involving the use of high-throughput sequencing technology. In this study, the transcriptome of cork-producing and non-cork-producing Quercus cerris × suber hybrids was analyzed in order to elucidate the differences between the two groups of trees displaying contrasting phenotypes for cork production. The results revealed the presence of a significant number of genes exclusively associated with cork production, in the trees that developed cork. Moreover, several gene ontology subcategories, such as cell wall biogenesis, lipid metabolic processes, metal ion binding and apoplast/cell wall, were only detected in the trees with cork production. These results indicate the existence, at the transcriptome level, of mechanisms that seem to be unique and necessary for cork production, which is an advancement in our knowledge regarding the genetic regulation behind cork formation and production.

The draft genome sequence of cork oak.

Cork oak (Quercus suber) is native to southwest Europe and northwest Africa where it plays a crucial environmental and economical role. To tackle the cork oak production and industrial challenges, advanced research is imperative but dependent on the availability of a sequenced genome. To address this, we produced the first draft version of the cork oak genome. We followed a de novo assembly strategy based on high-throughput sequence data, which generated a draft genome comprising 23,347 scaffolds and 953.3 Mb in size. A total of 79,752 genes and 83,814 transcripts were predicted, including 33,658 high-confidence genes. An InterPro signature assignment was detected for 69,218 transcripts, which represented 82.6% of the total. Validation studies demonstrated the genome assembly and annotation completeness and highlighted the usefulness of the draft genome for read mapping of high-throughput sequence data generated using different protocols. All data generated is available through the public databases where it was deposited, being therefore ready to use by the academic and industry communities working on cork oak and/or related species.

Draft Genome Sequence of a Rare Pigmented Mycobacterium avium subsp. paratuberculosis Type C Strain.

Mycobacterium avium subsp. paratuberculosis is the causative agent of paratuberculosis. We report here the draft genome sequence of a rare pigmented M. avium subsp. paratuberculosis type C strain, comprising 58 contigs and having a genome size of 4,851,414 bp. The genome will assist in the execution of pigmentation and virulence studies on this mycobacterium.

Genome-wide linkage analysis of global gene expression in loin muscle tissue identifies candidate genes in pigs.

BACKGROUND: Nearly 6,000 QTL have been reported for 588 different traits in pigs, more than in any other livestock species. However, this effort has translated into only a few confirmed causative variants. A powerful strategy for revealing candidate genes involves expression QTL (eQTL) mapping, where the mRNA abundance of a set of transcripts is used as the response variable for a QTL scan. METHODOLOGY/PRINCIPAL FINDINGS: We utilized a whole genome expression microarray and an F(2) pig resource population to conduct a global eQTL analysis in loin muscle tissue, and compared results to previously inferred phenotypic QTL (pQTL) from the same experimental cross. We found 62 unique eQTL (FDR <10%) and identified 3 gene networks enriched with genes subject to genetic control involved in lipid metabolism, DNA replication, and cell cycle regulation. We observed strong evidence of local regulation (40 out of 59 eQTL with known genomic position) and compared these eQTL to pQTL to help identify potential candidate genes. Among the interesting associations, we found aldo-keto reductase 7A2 (AKR7A2) and thioredoxin domain containing 12 (TXNDC12) eQTL that are part of a network associated with lipid metabolism and in turn overlap with pQTL regions for marbling, % intramuscular fat (% fat) and loin muscle area on Sus scrofa (SSC) chromosome 6. Additionally, we report 13 genomic regions with overlapping eQTL and pQTL involving 14 local eQTL. CONCLUSIONS/SIGNIFICANCE: Results of this analysis provide novel candidate genes for important complex pig phenotypes.