Effective population size (Ne) is a key population genetic parameter that describes the amount of genetic drift in a population. Estimating Ne has been subject to much research over the last 80 years. Methods to estimate Ne from linkage disequilibrium (LD) were developed ~40 years ago but depend on the availability of large amounts of genetic marker data that only the most recent advances in DNA technology have made available. Here we introduce SNeP, a multithreaded tool to perform the estimate of Ne using LD using the standard PLINK input file format (.ped and.map files) or by using LD values calculated using other software. Through SNeP the user can apply several corrections to take account of sample size, mutation, phasing, and recombination rate. Each variable involved in the computation such as the binning parameters or the chromosomes to include in the analysis can be modified. When applied to published datasets, SNeP produced results closely comparable with those obtained in the original studies. The use of SNeP to estimate Ne trends can improve understanding of population demography in the recent past, provided a sufficient number of SNPs and their physical position in the genome are available. Binaries for the most common operating systems are available at https://sourceforge.net/projects/snepnetrends/.

Three distinct mitochondrial maternal lineages (haplotype Groups A, B, and C) have been found in the domestic sheep. Group B has been observed primarily in European domestic sheep. The European mouflon carries this haplotype group. This could suggest that European mouflon was independently domesticated in Europe, although archaeological evidence supports sheep domestication in the central part of the Fertile Crescent. To investigate this question, we sequenced a highly variable segment of mitochondrial DNA (mtDNA) in 406 unrelated animals from 48 breeds or local varieties. They originated from a wide area spanning northern Europe and the Balkans to the Altay Mountains in south Siberia. The sample included a representative cross-section of sheep breeds from areas close to the postulated Near Eastern domestication center and breeds from more distant northern areas. Four (A, B, C, and D) highly diverged sheep lineages were observed in Caucasus, 3 (A, B and C) in Central Asia, and 2 (A and B) in the eastern fringe of Europe, which included the area north and west of the Black Sea and the Ural Mountains. Only one example of Group D was detected. The other haplotype groups demonstrated signs of population expansion. Sequence variation within the lineages implied Group A to have expanded first. This group was the most frequent type only in Caucasian and Central Asian breeds. Expansion of Group C appeared most recently. The expansion of Group B involving Caucasian sheep took place at nearly the same time as the expansion of Group A. Group B expansion for the eastern European area started approximately 3,000 years after the earliest inferred expansion. An independent European domestication of sheep is unlikely. The distribution of Group A variation as well as other results are compatible with the Near East being the domestication site. Groups C and D may have been introgressed later into a domestic stock, but larger samples are needed to infer their geographical origin. The results suggest that some mitochondrial lineages arrived in northern Europe from the Near East across Russia.

The extent of linkage disequilibrium (LD) and effective population size in Finnish Landrace and Finnish Yorkshire pig populations were studied using a whole genome SNP panel (Illumina PorcineSNP60 BeadChip) and pedigree data. Genotypic data included 86 Finnish Landrace and 32 Finnish Yorkshire boars. Pedigree data included 608,138 Finnish Landrace 554,237 and Finnish Yorkshire pigs, and on average 15 ancestral generations were known for the reference animals, born in 2005 to 2009. The breeding animals of the 2 populations have been kept separate in the breeding programs. Based on the pedigree data, the current effective population size for Finnish Landrace is 91 and for Finnish Yorkshire 61. Linkage disequilibrium measures (D' and r(2)) were estimated for over 1.5 million pairs of SNP. Average r(2) for SNP 30 kb apart was 0.47 and 0.49 and for SNP 5 Mb apart 0.09 and 0.12 for Finnish Landrace and Finnish Yorkshire, respectively. Average LD (r(2)) between adjacent SNP in the Illumina PorcineSNP60 BeadChip was 0.43 (57% of the adjacent SNP pairs had r(2) > 0.2) for Finnish Landrace and 0.46 (60% of the adjacent SNP pairs had r(2) > 0.2) for Finnish Yorkshire, and average r(2) > 0.2 extended to 1.0 and 1.5 Mb for Finnish Landrace and Finnish Yorkshire, respectively. Effective population size estimates based on the decay of r(2) with distance were similar to those based on the pedigree data: 80 and 55 for Finnish Landrace and Finnish Yorkshire, respectively. Thus, the results indicate that the effective population size of Finnish Yorkshire is smaller than of Finnish Landrace and has a clear effect on the extent of LD. The current effective population size of both breeds is above the recommended minimum of 50 but may get smaller than that in the near future, if no action is taken to balance the inbreeding rate and selection response. Because a moderate level of LD extends over a long distance, selection based on whole genome SNP markers (genomic selection) is expected to be efficient for both breeds.

African bovine trypanosomiasis caused by Trypanosoma sp., is a major constraint on cattle productivity in sub-Saharan Africa. Some African Bos taurus breeds are highly tolerant of infection, but the potentially more productive Bos indicus zebu breeds are much more susceptible. Zebu cattle are well adapted for plowing and haulage, and increasing their tolerance of trypanosomiasis could have a major impact on crop cultivation as well as dairy and beef production. We used three strategies to obtain short lists of candidate genes within QTL that were previously shown to regulate response to infection. We analyzed the transcriptomes of trypanotolerant N'Dama and susceptible Boran cattle after infection with Trypanosoma congolense. We sequenced EST libraries from these two breeds to identify polymorphisms that might underlie previously identified quantitative trait loci (QTL), and we assessed QTL regions and candidate loci for evidence of selective sweeps. The scan of the EST sequences identified a previously undescribed polymorphism in ARHGAP15 in the Bta2 trypanotolerance QTL. The polymorphism affects gene function in vitro and could contribute to the observed differences in expression of the MAPK pathway in vivo. The expression data showed that TLR and MAPK pathways responded to infection, and the former contained TICAM1, which is within a QTL on Bta7. Genetic analyses showed that selective sweeps had occurred at TICAM1 and ARHGAP15 loci in African taurine cattle, making them strong candidates for the genes underlying the QTL. Candidate QTL genes were identified in other QTL by their expression profile and the pathways in which they participate. nagana | positional cloning | sustainable agriculture | trypanosomosis

The ruminal microbiome, comprising large numbers of bacteria, ciliate protozoa, archaea and fungi, responds to diet and dietary additives in a complex way. The aim of this study was to investigate the benefits of increasing the depth of the community analysis in describing and explaining responses to dietary changes. Quantitative PCR, ssu rRNA amplicon based taxa composition, diversity and co-occurrence network analyses were applied to ruminal digesta samples obtained from four multiparous Nordic Red dairy cows fitted with rumen cannulae. The cows received diets with forage:concentrate ratio either 35:65 (diet H) or 65:35 (L), supplemented or not with sunflower oil (SO) (0 or 50 g/kg diet dry matter), supplied in a 4 × 4 Latin square design with a 2 × 2 factorial arrangement of treatments and four 35-day periods. Digesta samples were collected on days 22 and 24 and combined. QPCR provided a broad picture in which a large fall in the abundance of fungi was seen with SO in the H but not the L diet. Amplicon sequencing showed higher community diversity indices in L as compared to H diets and revealed diet specific taxa abundance changes, highlighting large differences in protozoal and fungal composition. Methanobrevibacter ruminantium and Mbb. gottschalkii dominated archaeal communities, and their abundance correlated negatively with each other. Co-occurrence network analysis provided evidence that no microbial domain played a more central role in network formation, that some minor-abundance taxa were at nodes of highest centrality, and that microbial interactions were diet specific. Networks added new dimensions to our understanding of the diet effect on rumen microbial community interactions.

Sequence variation present within the mitochondrial genome was used to investigate genetic diversity within sheep breeds from Asia and Europe. Comparison of 2027 bp of sequence from 121 animals revealed 44 phylogenetically informative nucleotide positions and a single insertion/deletion. A total of 57 haplotypes were observed which formed two distinct clades. Type A haplotypes were found in breeds from Asia (India, Indonesia, Mongolia, and Tibet), while type B haplotypes were observed at the highest frequency in breeds sourced from Europe (nine breeds from Austria, Aland, Finland, Spain, and northwestern Russia). The distribution of haplotypes indicates sheep appear to have the weakest population structure and the highest rate of intercontinental dispersal of any domestic animal reported to date. Only 2.7% of the sequence variation observed was partitioned between continents, which is lower than both goat (approximately 10%) and cattle (approximately 50%). Diagnostic restriction fragment length polymorphism polymerase chain reaction (RFLP-PCR) tests which distinguish type A and B haplotypes were used to test an additional 223 animals from 17 breeds of European and Asian origin. A mixture of the two lineages was found in every breed except Suffolk and the Indian Garole, indicating introgression has played a major part during breed development and subsequent selection.

Population contribution to genetic diversity can be estimated using neutral variation. However, population expansion or hybridization of diverged ancestries may weaken correlation between neutral and non-neutral variation. Microsatellite variation was studied at 25 loci in 20 native and 12 modern or imported northern European sheep breeds. Breed contributions to total gene diversity, allelic richness and mean allele-sharing distance between individuals were measured. Indications of changes in population size and admixtures of divergent ancestries were investigated and the extent of inbreeding was estimated. The northern European sheep demonstrated signs of reduction in effective population size. Many old, small populations made a substantial positive contribution to total molecular variation, but populations with several divergent major ancestries did not contribute substantially to molecular variation, with the exception of the Norwegian Rygja sheep. However, several diverged major ancestries may cause it to contribute less to non-neutral variation than expected from the microsatellite data. Breed uniqueness and within-breed variability generally had opposite effects on breed contributions to molecular diversity. The degree of inbreeding did not reflect the breed contribution to total gene diversity or allelic richness, but inbred populations increased the mean allele-sharing distance between individuals. Our study indicates breed conservation to be especially important in maintaining allelic variation in northern European sheep and supports the evolutionary importance of peripheral populations.

Effective population size is an important parameter for the assessment of genetic diversity within a livestock population and its development over time. If pedigree information is not available, linkage disequilibrium (LD) analysis might offer an alternative perspective for the estimation of effective population size. In this study, 128 individuals of the Swiss Eringer breed were genotyped using the Illumina BovineSNP50 beadchip. We set bin size at 50 kb for LD analysis, assuming that LD for proximal single nucleotide polymorphism (SNP)-pairs reflects distant breeding history while LD from distal SNP-pairs would reflect near history. Recombination rates varied among different regions of the genome. The use of physical distances as an approximation of genetic distances (e.g. setting 1 Mb = 0.01 Morgan) led to an upward bias in LD-based estimates of effective population size for generations beyond 50, while estimates for recent history were unaffected. Correction for restricted sample size did not substantially affect these results. LD-based actual effective population size was estimated in the range of 87-149, whereas pedigree-based effective population size resulted in 321 individuals. For conservation purposes, requiring knowledge of recent history (<50 generations), approximation assuming constant recombination rate seemed adequate.