Population Structure and Diversity in European Honey Bees (Apis mellifera L.)—An Empirical Comparison of Pool and Individual Whole-Genome Sequencing

Background: Whole-genome sequencing has become routine for population genetic studies. Sequencing of individuals provides maximal data but is rather expensive and fewer samples can be studied. In contrast, sequencing a pool of samples (pool-seq) can provide sufficient data, while presenting less of an economic challenge. Few studies have compared the two approaches to infer population genetic structure and diversity in real datasets. Here, we apply individual sequencing (ind-seq) and pool-seq to the study of Western honey bees (Apis mellifera). Methods: We collected honey bee workers that belonged to 14 populations, including 13 subspecies, totaling 1347 colonies, who were individually (139 individuals) and pool-sequenced (14 pools). We compared allele frequencies, genetic diversity estimates, and population structure as inferred by the two approaches. Results: Pool-seq and ind-seq revealed near identical population structure and genetic diversities, albeit at different costs. While pool-seq provides genome-wide polymorphism data at considerably lower costs, ind-seq can provide additional information, including the identification of population substructures, hybridization, or individual outliers. Conclusions: If costs are not the limiting factor, we recommend using ind-seq, as population genetic structure can be inferred similarly well, with the advantage gained from individual genetic information. Not least, it also significantly reduces the effort required for the collection of numerous samples and their further processing in the laboratory.

Transcriptomic dataset for Sardina pilchardus: Assembly, annotation, and expression of nine tissues

European sardine or pilchard is a planktonic small pelagic fish present from the North Sea in Europe to the coast of Senegal in the North of Africa, and across the Mediterranean sea to the Black Sea. Ecologically, sardines are an intermediary link in the trophic network, preying on plankton and being predated by larger fishes, marine mammals, and seabirds. This species is of great nutritional and economic value as a cheap but rich source of protein and fat. It is either consumed directly by humans or fed as fishmeal for aquaculture and farm animals. Despite its importance in the food basket, little is known about the molecular mechanisms involved in protein and lipid synthesis in this species. We collected nine tissues of Sardina pilchardus and reconstructed the transcriptome. In all, 198,597 transcripts were obtained, from which 68,031 are protein-coding. Quality assessment of the transcriptome was performed by back-mapping reads to the transcriptome and by searching for Single Copy Orthologs. Additionally, Gene Ontology and KEGG annotations were retrieved for most of the protein-coding genes. Finally, each library was quantified in terms of Transcripts per Million to disclose their expression patterns.

Authoritative subspecies diagnosis tool for European honey bees based on ancestry informative SNPs

BACKGROUND

With numerous endemic subspecies representing four of its five evolutionary lineages, Europe holds a large fraction of Apis mellifera genetic diversity. This diversity and the natural distribution range have been altered by anthropogenic factors. The conservation of this natural heritage relies on the availability of accurate tools for subspecies diagnosis. Based on pool-sequence data from 2145 worker bees representing 22 populations sampled across Europe, we employed two highly discriminative approaches (PCA and F_ST) to select the most informative SNPs for ancestry inference.

RESULTS

Using a supervised machine learning (ML) approach and a set of 3896 genotyped individuals, we could show that the 4094 selected single nucleotide polymorphisms (SNPs) provide an accurate prediction of ancestry inference in European honey bees. The best ML model was Linear Support Vector Classifier (Linear SVC) which correctly assigned most individuals to one of the 14 subspecies or different genetic origins with a mean accuracy of 96.2% ± 0.8 SD. A total of 3.8% of test individuals were misclassified, most probably due to limited differentiation between the subspecies caused by close geographical proximity, or human interference of genetic integrity of reference subspecies, or a combination thereof.

CONCLUSIONS

The diagnostic tool presented here will contribute to a sustainable conservation and support breeding activities in order to preserve the genetic heritage of European honey bees.

EXFI: Exon and splice graph prediction without a reference genome

For population genetic studies in nonmodel organisms, it is important to use every single source of genomic information. This paper presents EXFI, a Python pipeline that predicts the splice graph and exon sequences using an assembled transcriptome and raw whole-genome sequencing reads. The main algorithm uses Bloom filters to remove reads that are not part of the transcriptome, to predict the intron-exon boundaries, to then proceed to call exons from the assembly, and to generate the underlying splice graph. The results are returned in GFA1 format, which encodes both the predicted exon sequences and how they are connected to form transcripts. EXFI is written in Python, tested on Linux platforms, and the source code is available under the MIT License at https://github.com/jlanga/exfi.

Genomic selection signatures in sheep from the Western Pyrenees

BACKGROUND

The current large spectrum of sheep phenotypic diversity results from the combined product of sheep selection for different production traits such as wool, milk and meat, and its natural adaptation to new environments. In this study, we scanned the genome of 25 Sasi Ardi and 75 Latxa sheep from the Western Pyrenees for three types of regions under selection: (1) regions underlying local adaptation of Sasi Ardi semi-feral sheep, (2) regions related to a long traditional dairy selection pressure in Latxa sheep, and (3) regions experiencing the specific effect of the modern genetic improvement program established for the Latxa breed during the last three decades.

RESULTS

Thirty-two selected candidate regions including 147 annotated genes were detected by using three statistical parameters: pooled heterozygosity H, Tajima's D, and Wright's fixation index F_st. For Sasi Ardi sheep, chromosomes Ovis aries (OAR)4, 6, and 22 showed the strongest signals and harbored several candidate genes related to energy metabolism and morphology (BBS9, ELOVL3 and LDB1), immunity (NFKB2), and reproduction (H2AFZ). The major genomic difference between Sasi Ardi and Latxa sheep was on OAR6, which is known to affect milk production, with highly selected regions around the ABCG2, SPP1, LAP3, NCAPG, LCORL, and MEPE genes in Latxa sheep. The effect of the modern genetic improvement program on Latxa sheep was also evident on OAR15, on which several olfactory genes are located. We also detected several genes involved in reproduction such as ESR1 and ZNF366 that were affected by this selection program.

CONCLUSIONS

Natural and artificial selection have shaped the genome of both Sasi Ardi and Latxa sheep. Our results suggest that Sasi Ardi traits related to energy metabolism, morphological, reproductive, and immunological features have been under positive selection to adapt this semi-feral sheep to its particular environment. The highly selected Latxa sheep for dairy production showed clear signatures of selection in genomic regions related to milk production. Furthermore, our data indicate that the selection criteria applied in the modern genetic improvement program affect immunity and reproduction traits.

Training medical assistants for surgery

A successful programme is reported from Mozambique for training middle-level health workers to perform fairly advanced surgical procedures in remote areas where the services of consultants are virtually unobtainable. Manpower and financial constraints obliged Mozambique to train medical assistants to perform surgical work in rural areas, where three broad priorities were identified: pregnancy-related complications, trauma-related complications, and emergency inflammatory conditions. Since 1984, 20 health workers have emerged from three-year courses to become técnicos de cirurgía (assistant medical officers), and it is expected that there will be 46 by 1999. The training comprises two years of lectures and practical sessions in the Maputo Central Hospital, and a practical internship lasting a year at a provincial hospital. Three workshops organized since 1989 suggest that the upgraded personnel are performing well. More detailed evaluation and follow-up are in progress. Throughout 1995 a follow-up was conducted on 14 assistant medical officers. They performed 10,258 surgical operations, some 70% of which were emergency interventions. Low rates of complication occurred and postoperative mortality amounted to 0.4% and 0.1% in emergency and elective interventions respectively.

Type III-B open tibial fractures in Mozambique. A prospective study of 50 cases

We have studied prospectively 50 cases of open tibial fractures Type-III B. Treatment was based on the principles of debridement, wound irrigation, fracture stabilisation (generally by external fixation) and 'early bone coverage' without skin sutures. The cases were divided into two groups: Group 1: Eight fractures with segmental bone defects treated by bone transport using the Ilizarov technique or by intertibial-fibular bone grafting, and group 2: Forty-two fractures with defects of soft tissue coverage without segmental bone defects treated with simple closure or myoplasty. Union was achieved in Group 1 in an average of 16 months (range 6-21.5), and in an average of 5.6 months (range 3-10) in Group 2. The most frequent complication was malunion in angulation which occurred in 11 cases (22%). There was an equinus contracture of the ankle in 8 cases (16%) and chronic osteomyelitis occurred in 5 (10%). Adequate union was achieved in all but one case.

Hip dislocations in chronic osteomyelitis

We present two cases of spontaneous dislocation of the hip joint due to chronic osteomyelitis of the upper femur. Spontaneous hip dislocations occur in certain pathological conditions such as poliomyelitis (Ingram, 1980), cerebral palsy (Howard et al., 1985), osteomyelitis and neurofibromatosis.