Genome-wide association studies for diarrhoea outcomes identified genomic regions affecting resistance to a severe enteropathy in suckling rabbits.

Selection and breeding strategies to improve resistance to enteropathies are essential to reaching the sustainability of the rabbit production systems. However, disease heterogeneity (having only as major visible symptom diarrhoea) and low disease heritability are two barriers for the implementation of these strategies. Diarrhoea condition can affect rabbits at different life stages, starting from the suckling period, with large negative economic impacts. In this study, from a commercial population of suckling rabbits (derived from 133 litters) that experienced an outbreak of enteropathy, we first selected a few animals that died with severe symptoms of diarrhoea and characterized their microbiota, using 16S rRNA gene sequencing data. Clostridium genus was consistently present in all affected specimens. In addition, with the aim to identify genetic markers in the rabbit genome that could be used as selection tools, we performed genome-wide association studies for symptoms of diarrhoea in the same commercial rabbit population. These studies were also complemented with FST analyses between the same groups of rabbits. A total of 332 suckling rabbits (151 with severe symptoms of diarrhoea, 42 with mild symptoms and 129 without any symptoms till the weaning period), derived from 45 different litters (a subset of the 133 litters) were genotyped with the Affymetrix Axiom OrcunSNP Array. In both genomic approaches, rabbits within litters were paired to constitute two groups (susceptible and resistant, including the mildly affected in one or the other group) and run case and control genome-wide association analyses. Genomic heritability estimated in the designed experimental structure integrated in a commercial breeding scheme was 0.19-0.21 (s.e. 0.09-0.10). A total of eight genomic regions on rabbit chromosome 2 (OCU2), OCU3, OCU7, OCU12, OCU13, OCU16 and in an unassembled scaffold had significant single nucleotide polymorphisms (SNPs) and/or markers that trespassed the FST percentile distribution. Among these regions, three main peaks of SNPs were identified on OCU12, OCU13 and OCU16. The QTL region on OCU13 encompasses several genes that encode members of a family of immunoglobulin Fc receptors (FCER1G, FCRLA, FCRLB and FCGR2A) involved in the immune innate system, which might be important candidate genes for this pathogenic condition. The results obtained in this study demonstrated that resistance to an enteropathy occurring in suckling rabbits is in part genetically determined and can be dissected at the genomic level, providing DNA markers that could be used in breeding programmes to increase resistance to enteropathies in meat rabbits.

Whole genome sequencing identifies candidate genes and mutations that can explain diluted and other colour varieties of domestic canaries (Serinus canaria).

The domestic canary (Serinus canaria) is one of the most common pet birds and has been extensively selected and bred over the last few centuries to constitute many different varieties. Plumage pigmentation is one of the main phenotypic traits that distinguish canary breeds and lines. Feather colours in these birds, similarly to other avian species, are mainly depended on the presence of two major types of pigments: carotenoids and melanins. In this study, we exploited whole genome sequencing (WGS) datasets produced from five canary lines or populations (Black Frosted Yellow, Opal, Onyx, Opal × Onyx and Mogno, some of which carrying different putative dilute alleles), complemented with other WGS datasets retrieved from previous studies, to identify candidate genes that might explain pigmentation variability across canary breeds and varieties. Sequencing data were obtained using a DNA pool-seq approach and genomic data were compared using window-based FST analyses. We identified signatures of selection in genomic regions harbouring genes involved in carotenoid-derived pigmentation variants (CYP2J19, EDC, BCO2 and SCARB1), confirming the results reported by previous works, and identified several other signatures of selection in the correspondence of melanogenesis-related genes (AGRP, ASIP, DCT, EDNRB, KITLG, MITF, MLPH, SLC45A2, TYRP1 and ZEB2). Two putative causative mutations were identified in the MLPH gene that may explain the Opal and Onyx dilute mutant alleles. Other signatures of selection were also identified that might explain additional phenotypic differences between the investigated canary populations.

Analysis of honey environmental DNA indicates that the honey bee (Apis mellifera L.) trypanosome parasite Lotmaria passim is widespread in the apiaries of the North of Italy.

Lotmaria passim is a trypanosomatid that infects honey bees. In this study, we established an axenic culture of L. passim from Italian isolates and then used its DNA as a control in subsequent analyses that investigated environmental DNA (eDNA) to detect this trypasonosomatid. The source of eDNA was honey, which has been already demonstrated to be useful to detect honey bee parasites. DNA from a total of 164 honey samples collected in the North of Italy was amplified with three L. passim specific PCR primers and 78% of the analysed samples gave positive results. These results indicated a high prevalence rate of this trypanosomatid in the North of Italy, where it might be considered another threat to honey bee health.

Signatures of Admixture and Genetic Uniqueness in the Autochthonous Greek Black Pig Breed Deduced from Gene Polymorphisms Affecting Domestication-Derived Traits.

The Greek Black Pig (or Greek Pig) is the only recognized autochthonous pig breed raised in Greece, usually in extensive or semi-extensive production systems. According to its name, the characteristic breed coat color is solid black. In this study, with the aim to start a systematic genetic characterization of the Greek Black Pig breed, we investigated polymorphisms in major genes well known to affect exterior and production traits (MC1R, KIT, NR6A1, VRTN and IGF2) and compared these data with population genetic information available in other Mediterranean and Western Balkan pig breeds and wild boars. None of the investigated gene markers were fixed for one allele, suggesting that, in the past, this breed experienced introgression from wild boars and admixture from cosmopolitan pig breeds, enriching the breed genetic pool that should be further investigated to design appropriate conservation genetic strategies. We identified a new MC1R allele, containing two missense mutations already reported in two other independent alleles, but here present in the same haplotype. This allele might be useful to disclose biological information that can lead to better understanding the cascade transmission of signals to produce melanin pigments. This study demonstrated that autochthonous genetic resources can be an interesting reservoir of unexpected genetic variants.

Honey Environmental DNA Can Be Used to Detect and Monitor Honey Bee Pests: Development of Methods Useful to Identify Aethina tumida and Galleria mellonella Infestations.

Environmental DNA (eDNA) contained in honey derives from the organisms that directly and indirectly have been involved in the production process of this matrix and that have played a role in the hive ecosystems where the honey has been produced. In this study we set up PCR-based assays to detect the presence of DNA traces left in the honey by two damaging honey bee pests: the small hive beetle (Aethina tumida) and the greater wax moth (Galleria mellonella). DNA was extracted from 82 honey samples produced in Italy and amplified using two specific primer pairs that target the mitochondrial gene cytochrome oxidase I (COI) of A. tumida and two specific primer pairs that target the same gene in G. mellonella. The limit of detection was tested using sequential dilutions of the pest DNA. Only one honey sample produced in Calabria was positive for A. tumida whereas about 66% of all samples were positively amplified for G. mellonella. The use of honey eDNA could be important to establish early and effective measures to contain at the local (e.g., apiary) or regional scales these two damaging pests and, particularly for the small hive beetle, to prevent its widespread diffusion.

A genotyping by sequencing approach can disclose Apis mellifera population genomic information contained in honey environmental DNA.

Awareness has been raised over the last years on the genetic integrity of autochthonous honey bee subspecies. Genomic tools available in Apis mellifera can make it possible to measure this information by targeting individual honey bee DNA. Honey contains DNA traces from all organisms that contributed or were involved in its production steps, including the honey bees of the colony. In this study, we designed and tested a genotyping by sequencing (GBS) assay to analyse single nucleotide polymorphisms (SNPs) of A. mellifera nuclear genome using environmental DNA extracted from honey. A total of 121 SNPs (97 SNPs informative for honey bee subspecies identification and 24 SNPs associated with relevant traits of the colonies) were used in the assay to genotype honey DNA, which derives from thousands of honey bees. Results were integrated with information derived from previous studies and whole genome resequencing datasets. This GBS method is highly reliable in estimating honey bee SNP allele frequencies of the whole colony from which the honey derived. This assay can be used to identify the honey bee subspecies of the colony that produced the honey and, in turn, to authenticate the entomological origin of the honey.

Application of Next Generation Semiconductor-Based Sequencing for the Identification of Apis mellifera Complementary Sex Determiner (csd) Alleles from Honey DNA.

The complementary sex determiner (csd) gene plays an essential role in the sex determination of Apis mellifera L. Females develop only if fertilized eggs have functional heterozygous genotypes at this gene whereas males, being haploids, are hemizygous. Two identical csd alleles produce non vital males. In light of the recent decline in honey bee populations, it is therefore important to monitor the allele variability at this gene. In this study, we tested the application of next generation semiconductor-based sequencing technology (Ion Torrent) coupled with environmental honey DNA as a source of honey bee genome information to retrieve massive sequencing data for the analysis of variability at the hypervariable region (HVR) of the csd gene. DNA was extracted from 12 honey samples collected from honeycombs directly retrieved from 12 different colonies. A specifically designed bioinformatic pipeline, applied to analyze a total of about 1.5 million reads, identified a total of 160 different csd alleles, 55% of which were novel. The average number of alleles per sample was compatible with the number of expected patrilines per colony, according to the mating behavior of the queens. Allele diversity at the csd could also provide information useful to reconstruct the history of the honey.

Distribution of the Main Apis mellifera Mitochondrial DNA Lineages in Italy Assessed Using an Environmental DNA Approach.

Growing interest has been emerging on the need to monitor the genetic integrity of the European Apis mellifera subspecies that could be threatened by the human-mediated dispersion of non-native populations and lines. Mitochondrial DNA (mtDNA) lineages can provide useful information for this purpose. In this study, we took advantage of the environmental DNA (eDNA) contained in the honey, which can be analyzed to detect the main groups of mitotypes of the honey bees that produced it. In this study, we applied this eDNA to produce a distribution map all over the Italian peninsula and the two major islands (Sicily and Sardinia) of the following three honey bee mtDNA lineages: A, C and M. A total of 607 georeferenced honey samples, produced in all Italian regions, was analyzed to detect these lineages. The A lineage was widespread in Sicily, as expected, considering that A. m. siciliana carries the African lineage. Surprisingly, this lineage was also reported in about 14% of all other samples produced in almost all continental regions, and in Sardinia. The applied method obtained an updated distribution map of honey bee mtDNA lineages that could be useful to design policies for the conservation of Italian honey bee genetic resources.

Redefinition of the Mora Romagnola Pig Breed Herd Book Standard Based on DNA Markers Useful to Authenticate Its "Mono-Breed" Products: An Example of Sustainable Conservation of a Livestock Genetic Resource.

Mora Romagnola is an autochthonous pig breed, raised in the north of Italy. Mono-breed pork products of this breed are part of important niche value chain that is intrinsically linked to the conservation of this local genetic resources that can only survive due to the premium price that these products can obtain on the market. However, the added value attracts fraudsters that unscrupulously sell mis-labelled Mora Romagnola products, causing consumer distrust that, in turn, undermines the conservation strategy of this breed. To monitor and better characterise this local breed, we phenotyped 826 Mora Romagnola pigs for three breed-specific traits. Then, we genotyped almost all living sows and boars registered to the Herd Book (n. = 357 animals) for polymorphisms in the MC1R and NR6A1 genes (affecting coat colour and vertebral number, respectively). The results were used to re-define the breed descriptors of the Mora Romagnala breed that included information on the allowed genotypes at these two genes. A few pigs that did not carry the allowed genotypes were excluded from its Herd Book. Finally, we evaluated the usefulness of these DNA markers to authenticate Mora Romagnola meat against meat derived from other 11 pig breeds and wild boars. To our knowledge, the Mora Romagnola Herd Book is one of the first examples that established a direct link between a genetic standard of a breed with the possibility to authenticate mono-breed products using DNA markers with the specific purpose to combat frauds and, indirectly, support the conservation of a livestock genetic resource.

Describing variability in pig genes involved in coronavirus infections for a One Health perspective in conservation of animal genetic resources.

Coronaviruses silently circulate in human and animal populations, causing mild to severe diseases. Therefore, livestock are important components of a "One Health" perspective aimed to control these viral infections. However, at present there is no example that considers pig genetic resources in this context. In this study, we investigated the variability of four genes (ACE2, ANPEP and DPP4 encoding for host receptors of the viral spike proteins and TMPRSS2 encoding for a host proteinase) in 23 European (19 autochthonous and three commercial breeds and one wild boar population) and two Asian Sus scrofa populations. A total of 2229 variants were identified in the four candidate genes: 26% of them were not previously described; 29 variants affected the protein sequence and might potentially interact with the infection mechanisms. The results coming from this work are a first step towards a "One Health" perspective that should consider conservation programs of pig genetic resources with twofold objectives: (i) genetic resources could be reservoirs of host gene variability useful to design selection programs to increase resistance to coronaviruses; (ii) the described variability in genes involved in coronavirus infections across many different pig populations might be part of a risk assessment including pig genetic resources.

Honey as a Source of Environmental DNA for the Detection and Monitoring of Honey Bee Pathogens and Parasites.

Environmental DNA (eDNA) has been proposed as a powerful tool to detect and monitor cryptic, elusive, or invasive organisms. We recently demonstrated that honey constitutes an easily accessible source of eDNA. In this study, we extracted DNA from 102 honey samples (74 from Italy and 28 from 17 other countries of all continents) and tested the presence of DNA of nine honey bee pathogens and parasites (Paenibacillus larvae, Melissococcus plutonius, Nosema apis, Nosema ceranae, Ascosphaera apis,Lotmaria passim, Acarapis woodi, Varroa destructor, and Tropilaelaps spp.) using qualitative PCR assays. All honey samples contained DNA from V. destructor, confirming the widespread diffusion of this mite. None of the samples gave positive amplifications for N. apis, A. woodi, and Tropilaelaps spp. M. plutonius was detected in 87% of the samples, whereas the other pathogens were detected in 43% to 57% of all samples. The frequency of Italian samples positive for P. larvae was significantly lower (49%) than in all other countries (79%). The co-occurrence of positive samples for L. passim and A. apis with N. ceranae was significant. This study demonstrated that honey eDNA can be useful to establish monitoring tools to evaluate the sanitary status of honey bee populations.

The choice of anaesthesia for glioblastoma surgery does not impact the time to recurrence.

Anaesthetics used during cancer surgery may influence tumour cells and immunological response. The aim of this study was to evaluate a potential influence of the anaesthetic method (inhaled anaesthetics versus total-intravenous anaesthesia using propofol) on recurrence-free and overall survival in glioblastoma patients. We retrospectively identified patients undergoing resection of contrast enhancing glioblastoma under general anaesthesia followed by standard adjuvant treatment between January 2010 and February 2017 at two University Hospitals. Matched pairs of patients receiving either balanced with volatile anaesthetics or total intravenous anaesthesia were generated according to the known prognostic factors (extent of resection, methyl-guanine-methyl-transferase (MGMT) promoter methylation, age, Karnofsky performance score). Groups were compared using chi-square and Whitney-Man-U test. Time to recurrence was calculated using Kaplan Meier estimates. Log Rank test was used to assess the influence of the anaesthetic method. One hundred and fifty-eight (79:79) patients were included. Groups showed no significant difference in recurrence-free (volatiles: 8.0 (95% CI 6.5-9.8) vs. propofol: 8.4 (95% CI 7.9-10.1) months; p = 0.54) or overall survival (propofol: 17.4 (95% CI 14.0-20.7) vs. volatiles: 16.9 (95% CI 13.9-20.1) months; p = 0.85). In contrast to potential beneficial effects in some other solid tumours, the choice of anaesthetic method had no impact on survival in patients with glioblastoma in a well-defined cohort.

Microscopic ossicle analyses and the complete mitochondrial genome sequence of Holothuria (Roweothuria) polii (Echinodermata; Holothuroidea) provide new information to support the phylogenetic positioning of this sea cucumber species.

Sea cucumbers (Holothuroidea) are ecologically important organisms for their bioturbation and alkalinization activities of the seabed. These species are extensively fished as they are considered luxury food. Sea cucumbers are also relevant for biomedical studies and the production of bioactive compounds. A few initiatives are recently evaluating sea cucumbers as novel aquaculture species. The aim of this study was to provide morphological and genetic information useful for the identification of Holothuria polii, the white spot sea cucumber (a common species of the Mediterranean Sea). We generated the complete sequence of the mitochondrial DNA (mtDNA) genome of this species and combined it with a detailed ossicle characterization of the sequenced specimen by scanning electron microscopic analysis. Ossicles (known also as sclerites) are anatomical features that can discriminate Holothuroidea species, including the closely related ones of the genus Holothuria. The complete mitochondrial genome was assembled, functionally annotated and then used to evaluate the phylogenetic relationship of H. polii against the other few Holothuroidea species for which the whole mtDNA was available. The 15,907 bp H. polii mtDNA sequence has the same gene order already reported for H. scabra, H. forskali and other species of the same class. Cox1 and 16S gene sequences were informative for species identification across the genus and could be used for the authentication of commercialized Holothuria spp. The mitochondrial genome sequence presented here provides the basis to a future analysis of the variability of H. polii populations in the Mediterranean region.