Whole-genome sequence analysis reveals differences in population management and selection of European low-input pig breeds

Background

A major concern in conservation genetics is to maintain the genetic diversity of populations. Genetic variation in livestock species is threatened by the progressive marginalisation of local breeds in benefit of high-output pigs worldwide. We used high-density SNP and re-sequencing data to assess genetic diversity of local pig breeds from Europe. In addition, we re-sequenced pigs from commercial breeds to identify potential candidate mutations responsible for phenotypic divergence among these groups of breeds.

Results

Our results point out some local breeds with low genetic diversity, whose genome shows a high proportion of regions of homozygosis (>50%) and that harbour a large number of potentially damaging mutations. We also observed a high correlation between genetic diversity estimates using high-density SNP data and Next Generation Sequencing data (r = 0.96 at individual level). The study of non-synonymous SNPs that were fixed in commercial breeds and also in any local breed, but with different allele, revealed 99 non-synonymous SNPs affecting 65 genes. Candidate mutations that may underlie differences in the adaptation to the environment were exemplified by the genes AZGP1 and TAS2R40. We also observed that highly productive breeds may have lost advantageous genotypes within genes involve in immune response – e.g. IL12RB2 and STAB1–, probably as a result of strong artificial in the intensive production systems in pig.

Conclusions

The high correlation between genetic diversity computed with the 60K SNP and whole genome re-sequence data indicates that the Porcine 60K SNP Beadchip provides reliable estimates of genomic diversity in European pig populations despite the expected bias. Moreover, this analysis gave insights for strategies to the genetic characterization of local breeds. The comparison between re-sequenced local pigs and re-sequenced commercial pigs made it possible to report candidate mutations to be responsible for phenotypic divergence among those groups of breeds. This study highlights the importance of low input breeds as a valuable genetic reservoir for the pig production industry. However, the high levels of ROHs, inbreeding and potentially damaging mutations emphasize the importance of the genetic characterization of local breeds to preserve their genomic variability.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-601) contains supplementary material, which is available to authorized users.     

Frequent Seizures Are Associated with a Network of Gray Matter Atrophy in Temporal Lobe Epilepsy with or without Hippocampal Sclerosis

Objective

Patients with temporal lobe epilepsy (TLE) with hippocampal sclerosis (HS) have diffuse subtle gray matter (GM) atrophy detectable by MRI quantification analyses. However, it is not clear whether the etiology and seizure frequency are associated with this atrophy. We aimed to evaluate the occurrence of GM atrophy and the influence of seizure frequency in patients with TLE and either normal MRI (TLE-NL) or MRI signs of HS (TLE-HS).

Methods

We evaluated a group of 172 consecutive patients with unilateral TLE-HS or TLE-NL as defined by hippocampal volumetry and signal quantification (122 TLE-HS and 50 TLE-NL) plus a group of 82 healthy individuals. Voxel-based morphometry was performed with VBM8/SPM8 in 3T MRIs. Patients with up to three complex partial seizures and no generalized tonic-clonic seizures in the previous year were considered to have infrequent seizures. Those who did not fulfill these criteria were considered to have frequent seizures.

Results

Patients with TLE-HS had more pronounced GM atrophy, including the ipsilateral mesial temporal structures, temporal lobe, bilateral thalami and pre/post-central gyri. Patients with TLE-NL had more subtle GM atrophy, including the ipsilateral orbitofrontal cortex, bilateral thalami and pre/post-central gyri. Both TLE-HS and TLE-NL showed increased GM volume in the contralateral pons. TLE-HS patients with frequent seizures had more pronounced GM atrophy in extra-temporal regions than TLE-HS with infrequent seizures. Patients with TLE-NL and infrequent seizures had no detectable GM atrophy. In both TLE-HS and TLE-NL, the duration of epilepsy correlated with GM atrophy in extra-hippocampal regions.

Conclusion

Although a diffuse network GM atrophy occurs in both TLE-HS and TLE-NL, this is strikingly more evident in TLE-HS and in patients with frequent seizures. These findings suggest that neocortical atrophy in TLE is related to the ongoing seizures and epilepsy duration, while thalamic atrophy is more probably related to the original epileptogenic process.     

Avian Influenza Virus (H11N9) in Migratory Shorebirds Wintering in the Amazon Region,Brazil

Aquatic birds are the natural reservoir for avian influenza viruses (AIV). Habitats in Brazil provide stopover and wintering sites for water birds that migrate between North and South America. The current study was conducted to elucidate the possibility of the transport of influenza A viruses by birds that migrate annually between the Northern and Southern Hemispheres. In total, 556 orotracheal/cloacal swab samples were collected for influenza A virus screening using real-time RT-PCR (rRT-PCR). The influenza A virus-positive samples were subjected to viral isolation. Four samples were positive for the influenza A matrix gene by rRT-PCR. From these samples, three viruses were isolated, sequenced and characterized. All positive samples originated from a single bird species, the ruddy turnstone (Arenaria interpres), that was caught in the Amazon region at Caeté Bay, Northeast Pará, at Ilha de Canelas. To our knowledge, this is the first isolation of H11N9 in the ruddy turnstone in South America.     

Effects of double ligation of Stensen's duct on the rabbit parotid gland

Salivary gland duct ligation is an alternative to gland excision for treating sialorrhea or reducing salivary gland size prior to tumor excision. Duct ligation also is used as an approach to study salivary gland aging, regeneration, radiotherapy, sialolithiasis and sialadenitis. Reports conflict about the contribution of each salivary cell population to gland size reduction after ductal ligation. Certain cell populations, especially acini, reportedly undergo atrophy, apoptosis and proliferation during reduction of gland size. Acini also have been reported to de-differentiate into ducts. These contradictory results have been attributed to different animal or salivary gland models, or to methods of ligation. We report here a bilateral double ligature technique for rabbit parotid glands with histologic observations at 1, 7, 14, 30, 60 days after ligation. A large battery of special stains and immunohistochemical procedures was employed to define the cell populations. Four stages with overlapping features were observed that led to progressive shutdown of gland activities: 1) marked atrophy of the acinar cells occurred by 14 days, 2) response to and removal of the secretory material trapped in the acinar and ductal lumens mainly between 30 and 60 days, 3) reduction in the number of parenchymal (mostly acinar) cells by apoptosis that occurred mainly between 14–30 days, and 4) maintenance of steady-state at 60 days with a low rate of fluid, protein, and glycoprotein secretion, which greatly decreased the number of leukocytes engaged in the removal of the luminal contents. The main post- ligation characteristics were dilation of ductal and acinar lumens, massive transient infiltration of mostly heterophils (rabbit polymorphonuclear leukocytes), acinar atrophy, and apoptosis of both acinar and ductal cells. Proliferation was uncommon except in the larger ducts. By 30 days, the distribution of myoepithelial cells had spread from exclusively investing the intercalated ducts pre-ligation to surrounding a majority of the residual duct-like structures, many of which clearly were atrophic acini. Thus, both atrophy and apoptosis made major contributions to the post-ligation reduction in gland size. Structures also occurred with both ductal and acinar markers that suggested acini differentiating into ducts. Overall, the reaction to duct ligation proceeded at a considerably slower pace in the rabbit parotid glands than has been reported for the salivary glands of the rat.     

Breeding for flavour of fresh market tomato: sources for increasing acid content

Breeding programs in tomato for fresh consumption have concluded in very productive varieties/hybrids with an extraordinary external quality. However, internal quality has not been a priority objective in these breeding programs, so present products have a lack of important internal quality properties. Thus, internal and nutritional improvement including taste and organoleptyc characteristics are important breeding objectives at present. A screening trial to identify sources for internal quality improvement by increasing acid content was carried out. 12 accessions of Lycopersicon esculentum and 8 of L. pimpinellifolium were tested with that purpose. Content in citric, malic, oxalic and fumaric acid by HPLC, soluble solids content (SSC)(1 Brix) by refractometry and total acidity by titration with NaOH were measured. Sources for high citric, malic and fumaric acid content were found to begin those breeding programs. Results could possibly suggest an independent genetic control for every acid.     

Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR–HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR–HLA interactions among all described worldwide populations, and that 83–97% of their KIR–HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR–HLA coevolution and its impact on human health and survival.     

6-Thioguanine-resistant T-lymphocyte autoradiographic assay. Determination of variation frequencies in individuals suspected of radiation exposure

We used the autoradiographic assay to assess human in vivo somatic cell gene mutation at the hypoxanthine guanine phosphoribosyl transferase (HGPRT) locus in T-lymphocytes. Cells able to incorporate tritiated thymidine in vitro in the presence of 6-thioguanine were enumerated in order to determine 6-thioguanine-resistance (TGr) variant frequencies in cryopreserved lymphocytes from 17 normal control individuals, from 3 persons suspected to have been exposed to 60Co in an accident in Cd. Juárez (Mexico), studied 24 months after the accident, and from 4 individuals who were in Kiev during the radiation accident in Chernobyl (U.S.S.R.); 2 of them were studied 1 month after the accident, and again 1 year after the first sampling, the other 2 were studied 13 months after the accident. The data obtained indicate that this assay may be useful in any laboratory of cytogenetics for human population monitoring and that its use following accidental exposure to ionizing radiation should be further evaluated.     

Stable levels of <Emphasis Type="Italic">Coxiella burnetii</Emphasis> prevalence in dairy sheep flocks but changes in genotype distribution after a 10-year period in northern Spain

Bulk tank milk (BTM) samples were collected from 81 sheep flocks in the Basque Country, Spain, in 2015 and were analysed for antibodies against Coxiella burnetii by ELISA and for C. burnetii DNA by real-time PCR. Thirty-two percent of the flocks had BTM antibodies against C. burnetii. Presence of C. burnetii DNA in BTM was detected in 23% of the flocks, suggesting recent C. burnetii infections. Retrospective data of BTM samples obtained from 154 sheep flocks investigated in 2005 in the same geographic area were compiled to assess temporal changes in C. burnetii infection. The overall percentage of infected sheep flocks did not significantly change after the 10-year period. Among the 46 flocks sampled in both periods, 11 flocks that were negative in 2005 were positive in 2015, 18 maintained their initial status (positive or negative), and 17 positive flocks were negative in 2015. These findings indicate that C. burnetii infection is a dynamic process in dairy sheep in northern Spain. Single nucleotide polymorphism (SNP) genotyping of positive samples identified three genotypes, SNP1 being the most prevalent in 2015 and SNP8 in 2005; SNP4 was only detected once in 2005. These results suggest possible changes in the pattern of genotype infection over time.