Genomic selection for survival under naturally occurring Saprolegnia oomycete infection in farmed European whitefish Coregonus lavaretus.

Saprolegnia oomycete infection causes serious economic losses and reduces fish health in aquaculture. Genomic selection based on thousands of DNA markers is a powerful tool to improve fish traits in selective breeding programs. Our goal was to develop a single nucleotide polymorphism (SNP) marker panel and to test its use in genomic selection for improved survival against Saprolegnia infection in European whitefish Coregonus lavaretus, the second most important farmed fish species in Finland. We used a double digest restriction site associated DNA (ddRAD) genotyping by sequencing method to produce a SNP panel, and we tested it analyzing data from a cohort of 1,335 fish, which were measured at different times for mortality to Saprolegnia oomycete infection and weight traits. We calculated the genetic relationship matrix (GRM) from the genome-wide genetic data, integrating it in multivariate mixed models used for the estimation of variance components and genomic breeding values (GEBVs), and to carry out Genome-Wide Association Studies for the presence of quantitative trait loci (QTL) affecting the phenotypes in analysis. We identified one major QTL on chromosome 6 affecting mortality to Saprolegnia infection, explaining 7.7% to 51.3% of genetic variance, and a QTL for weight on chromosome 4, explaining 1.8% to 5.4% of genetic variance. Heritability for mortality was 0.20 to 0.43 on the liability scale, and heritability for weight was 0.44 to 0.53. The QTL for mortality showed an additive allelic effect. We tested whether integrating the QTL for mortality as a fixed factor, together with a new GRM calculated excluding the QTL from the genetic data, would improve the accuracy estimation of GEBVs. This test was done through a cross-validation approach, which indicated that the inclusion of the QTL increased the mean accuracy of the GEBVs by 0.28 points, from 0.33 to 0.61, relative to the use of full GRM only. The area under the curve of the receiver-operator curve for mortality increased from 0.58 to 0.67 when the QTL was included in the model. The inclusion of the QTL as a fixed effect in the model increased the correlation between the GEBVs of early mortality with the late mortality, compared to a model that did not include the QTL. These results validate the usability of the produced SNP panel for genomic selection in European whitefish and highlight the opportunity for modeling QTLs in genomic evaluation of mortality due to Saprolegnia infection.

Saprolegnia infection causes serious economic losses and reduces fish health in aquaculture. We created a novel set of genetic markers to use in the selective breeding of European whitefish to reduce mortality due to the fungus. Using genetic markers, we estimated how much different fish traits are determined by genetic variation, and thus what potential traits have to be selected. We observed that resistance to infection was controlled by both a genetic variant with a major effect on mortality and by many other variants with a small effect distributed across the genome. We tested whether we could increase the precision of genomic breeding values used in the selective breeding by explicitly adding the major genetic variant to the analysis, and we observed an increase in precision in our results. We conclude that directly including information about the major genetic variant increases the precision of our predictions, rather than assuming that all genetic variants each explain a small amount of the genetic variation.

Rumen Microbiota Predicts Feed Efficiency of Primiparous Nordic Red Dairy Cows.

Efficient feed utilization in dairy cows is crucial for economic and environmental reasons. The rumen microbiota plays a significant role in feed efficiency, but studies utilizing microbial data to predict host phenotype are limited. In this study, 87 primiparous Nordic Red dairy cows were ranked for feed efficiency during their early lactation based on residual energy intake, and the rumen liquid microbial ecosystem was subsequently evaluated using 16S rRNA amplicon and metagenome sequencing. The study used amplicon data to build an extreme gradient boosting model, demonstrating that taxonomic microbial variation can predict efficiency (rtest = 0.55). Prediction interpreters and microbial network revealed that predictions were based on microbial consortia and the efficient animals had more of the highly interacting microbes and consortia. Rumen metagenome data was used to evaluate carbohydrate-active enzymes and metabolic pathway differences between efficiency phenotypes. The study showed that an efficient rumen had a higher abundance of glycoside hydrolases, while an inefficient rumen had more glycosyl transferases. Enrichment of metabolic pathways was observed in the inefficient group, while efficient animals emphasized bacterial environmental sensing and motility over microbial growth. The results suggest that inter-kingdom interactions should be further analyzed to understand their association with the feed efficiency of animals.

Long-term effects of early-life rumen microbiota modulation on dairy cow production performance and methane emissions.

Rumen microbiota modulation during the pre-weaning period has been suggested as means to affect animal performance later in life. In this follow-up study, we examined the post-weaning rumen microbiota development differences in monozygotic twin-heifers that were inoculated (T-group) or not inoculated (C-group) (n = 4 each) with fresh adult rumen liquid during their pre-weaning period. We also assessed the treatment effect on production parameters and methane emissions of cows during their 1st lactation period. The rumen microbiota was determined by the 16S rRNA gene, 18S rRNA gene, and ITS1 amplicon sequencing. Animal weight gain and rumen fermentation parameters were monitored from 2 to 12 months of age. The weight gain was not affected by treatment, but butyrate proportion was higher in T-group in month 3 (p = 0.04). Apart from archaea (p = 0.084), the richness of bacteria (p < 0.0001) and ciliate protozoa increased until month 7 (p = 0.004) and anaerobic fungi until month 11 (p = 0.005). The microbiota structure, measured as Bray-Curtis distances, continued to develop until months 3, 6, 7, and 10, in archaea, ciliate protozoa, bacteria, and anaerobic fungi, respectively (for all: p = 0.001). Treatment or age × treatment interaction had a significant (p < 0.05) effect on 18 bacterial, 2 archaeal, and 6 ciliate protozoan taxonomic groups, with differences occurring mostly before month 4 in bacteria, and month 3 in archaea and ciliate protozoa. Treatment stimulated earlier maturation of prokaryote community in T-group before month 4 and earlier maturation of ciliate protozoa at month 2 (Random Forest: 0.75 month for bacteria and 1.5 month for protozoa). No treatment effect on the maturity of anaerobic fungi was observed. The milk production and quality, feed efficiency, and methane emissions were monitored during cow's 1st lactation. The T-group had lower variation in energy-corrected milk yield (p < 0.001), tended to differ in pattern of residual energy intake over time (p = 0.069), and had numerically lower somatic cell count throughout their 1st lactation period (p = 0.081), but no differences between the groups in methane emissions (g/d, g/kg DMI, or g/kg milk) were observed. Our results demonstrated that the orally administered microbial inoculant induced transient changes in early rumen microbiome maturation. In addition, the treatment may influence the later production performance, although the mechanisms that mediate these effects need to be further explored.

SNPs in Sheep: Characterization of Lithuanian Sheep Populations.

In Lithuania, there are two recognised native sheep breeds: old native Lithuanian Coarsewooled and Lithuanian Blackface. In addition, in 2005, primitive Heidschnucke-type Skudde sheep were imported to Lithuania and were argued to possibly represent a lost Lithuanian sheep type. The aim of the study was to investigate the genetic variation in the two Lithuanian native sheep breeds, compare them with the imported Skudde sheep and establish the historical patterns of admixture and the genetic relatedness of Lithuanian sheep to British, Central European and Nordic sheep breeds included in the SheepHapMap study. In total, 72 individuals, representing two Lithuanian native and imported Skudde sheep breeds, were genotyped using a Neogen 12K Illumina Infinium chip. The population analysis was carried out by model-based clustering, principal component analysis and neighbour net analysis, and showed similar patterns for the Lithuanian sheep populations. Lithuanian Coarsewooled and Skudde in Lithuania have unique divergence and possibly some shared ancestry, while the Lithuanian Blackface conforms to a modern synthetic breed. The study clearly showed that the Coarsewooled and the Skudde breeds are distinct from each other. Historical data strongly suggest that the Coarsewooled breed represents a local breed, while the Skudde origin is less directly linked to the geographical area of modern-day Lithuania. Within the modern-day Lithuanian context, the Lithuanian Coarsewooled sheep is very important historical sheep type for conservation.

The plant-based by-product diets for the mass-rearing of Acheta domesticus and Gryllus bimaculatus.

Edible insect rearing could provide one alternative for protein production by having a smaller environmental impact than traditional livestock farming due to insects' ability to convert organic side streams. Currently, the insect rearing industry utilizes soybeans as a major source of protein in the feeds. Protein-rich by-products of food industry could be used to replace them in insect feeds, but it is not known if they also meet the insects' nutritional requirements. Our study evaluated the growth performance of two widely used edible cricket species, Acheta domesticus and Gryllus bimaculatus (Orthoptera: Gryllidae), on 18 experimental diets. The experimental diets included commercial chicken feeds and cricket diets, where soybean was partly and completely replaced with by-products from food industry: potato protein, barley mash, barley feed, compressed leftover of turnip rape and mix of broad bean and pea on three levels of protein. We found that the high- and medium-protein turnip rape and barley mash diets produced the highest yield and an increase in all performance variables. Overall, the high- and medium-protein diets produced the highest yield, growth and fastest development. Our results showed that by-products of food industry could be utilized as a part of the cricket feeds and thus advance the goals of circular economy.

Legacies of domestication, trade and herder mobility shape extant male zebu cattle diversity in South Asia and Africa.

All tropically adapted humped cattle (Bos indicus or "zebu"), descend from a domestication process that took place >8,000 years ago in South Asia. Here we present an intercontinental survey of Y-chromosome diversity and a comprehensive reconstruction of male-lineage zebu cattle history and diversity patterns. Phylogenetic analysis revealed that all the zebu Y-chromosome haplotypes in our dataset group within three different lineages: Y3A, the most predominant and cosmopolitan lineage; Y3B, only observed in West Africa; and Y3C, predominant in South and Northeast India. The divergence times estimated for these three Zebu-specific lineages predate domestication. Coalescent demographic models support either de novo domestication of genetically divergent paternal lineages or more complex process including gene flow between wild and domestic animals. Our data suggest export of varied zebu lineages from domestication centres through time. The almost exclusive presence of Y3A haplotypes in East Africa is consistent with recent cattle restocking in this area. The cryptic presence of Y3B haplotypes in West Africa, found nowhere else, suggests that these haplotypes might represent the oldest zebu lineage introduced to Africa ca. 3,000 B.P. and subsequently replaced in most of the world. The informative ability of Interspersed Multilocus Microsatellites and Y-specific microsatellites to identify genetic structuring in cattle populations is confirmed.

Taxon abundance, diversity, co-occurrence and network analysis of the ruminal microbiota in response to dietary changes in dairy cows.

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.

Co-infections determine patterns of mortality in a population exposed to parasite infection.

Many individual hosts are infected with multiple parasite species, and this may increase or decrease the pathogenicity of the infections. This phenomenon is termed heterologous reactivity and is potentially an important determinant of both patterns of morbidity and mortality and of the impact of disease control measures at the population level. Using infections with Theileria parva (a tick-borne protozoan, related to Plasmodium) in indigenous African cattle [where it causes East Coast fever (ECF)] as a model system, we obtain the first quantitative estimate of the effects of heterologous reactivity for any parasitic disease. In individual calves, concurrent co-infection with less pathogenic species of Theileria resulted in an 89% reduction in mortality associated with T. parva infection. Across our study population, this corresponds to a net reduction in mortality due to ECF of greater than 40%. Using a mathematical model, we demonstrate that this degree of heterologous protection provides a unifying explanation for apparently disparate epidemiological patterns: variable disease-induced mortality rates, age-mortality profiles, weak correlations between the incidence of infection and disease (known as endemic stability), and poor efficacy of interventions that reduce exposure to multiple parasite species. These findings can be generalized to many other infectious diseases, including human malaria, and illustrate how co-infections can play a key role in determining population-level patterns of morbidity and mortality due to parasite infections.

SNeP: a tool to estimate trends in recent effective population size trajectories using genome-wide SNP data.

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/.

Genetic susceptibility to infectious disease in East African Shorthorn Zebu: a genome-wide analysis of the effect of heterozygosity and exotic introgression.

BACKGROUND: Positive multi-locus heterozygosity-fitness correlations have been observed in a number of natural populations. They have been explained by the correlation between heterozygosity and inbreeding, and the negative effect of inbreeding on fitness (inbreeding depression). Exotic introgression in a locally adapted population has also been found to reduce fitness (outbreeding depression) through the breaking-up of co-adapted genes, or the introduction of non-locally adapted gene variants. In this study we examined the inter-relationships between genome-wide heterozygosity, introgression, and death or illness as a result of infectious disease in a sample of calves from an indigenous population of East African Shorthorn Zebu (crossbred Bos taurus x Bos indicus) in western Kenya. These calves were observed from birth to one year of age as part of the Infectious Disease in East African Livestock (IDEAL) project. Some of the calves were found to be genetic hybrids, resulting from the recent introgression of European cattle breed(s) into the indigenous population. European cattle are known to be less well adapted to the infectious diseases present in East Africa. If death and illness as a result of infectious disease have a genetic basis within the population, we would expect both a negative association of these outcomes with introgression and a positive association with heterozygosity. RESULTS: In this indigenous livestock population we observed negative associations between heterozygosity and both death and illness as a result of infectious disease and a positive association between European taurine introgression and episodes of clinical illness. CONCLUSION: We observe the effects of both inbreeding and outbreeding depression in the East African Shorthorn Zebu, and therefore find evidence of a genetic component to vulnerability to infectious disease. These results indicate that the significant burden of infectious disease in this population could, in principle, be reduced by altered breeding practices.

Genetic and expression analysis of cattle identifies candidate genes in pathways responding to Trypanosoma congolense infection.

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.

Microsatellite-based genetic diversity and population structure of domestic sheep in northern Eurasia.

BACKGROUND: Identification of global livestock diversity hotspots and their importance in diversity maintenance is essential for making global conservation efforts. We screened 52 sheep breeds from the Eurasian subcontinent with 20 microsatellite markers. By estimating and weighting differently within- and between-breed genetic variation our aims were to identify genetic diversity hotspots and prioritize the importance of each breed for conservation, respectively. In addition we estimated how important within-species diversity hotspots are in livestock conservation. RESULTS: Bayesian clustering analysis revealed three genetic clusters, termed Nordic, Composite and Fat-tailed. Southern breeds from close to the region of sheep domestication were more variable, but less genetically differentiated compared with more northern populations. Decreasing weight for within-breed diversity component led to very high representation of genetic clusters or regions containing more diverged breeds, but did not increase phenotypic diversity among the high ranked breeds. Sampling populations throughout 14 regional groups was suggested for maximized total genetic diversity. CONCLUSIONS: During initial steps of establishing a livestock conservation program populations from the diversity hot-spot area are the most important ones, but for the full design our results suggested that approximately equal population presentation across environments should be considered. Even in this case, higher per population emphasis in areas of high diversity is appropriate. The analysis was based on neutral data, but we have no reason to think the general trend is limited to this type of data. However, a comprehensive valuation of populations should balance production systems, phenotypic traits and available genetic information, and include consideration of probability of success.

Estimation of relatedness among non-pedigreed Yakutian cryo-bank bulls using molecular data: implications for conservation and breed management.

BACKGROUND: Yakutian cattle, the last remaining native cattle breed in Siberia, are well adapted to the extreme sub-arctic conditions. Nowadays only ca. 1200 purebred animals are left in Yakutia. The semen of six Yakutian bulls was stored in a cryo-bank without any pedigree documentation because of the traditional free herding style of the population. METHODS: To clarify the genetic relatedness between these bulls and to provide recommendations to use their semen in future conservation and breed management programs, we have analysed 30 autosomal microsatellites and mitochondrial DNA sequences in 60 individuals including the six for which semen has been stored. Four relatedness estimators were calculated. In addition, we assessed the value of the cryo-bank bulls for the preservation of genetic variation of the contemporary Yakutian cattle by calculating allelic and gene diversity estimates and mean molecular coancestries. RESULTS: On the basis of microsatellite variability, including the Yakutian cryo-bank bulls increases the allelic variation in the contemporary population by 3% and in the male subpopulation by 13%. In terms of the mean molecular coancestries, they are less related to the contemporary cow population than the breeding bulls and therefore could be used to reduce inbreeding in the living population. Although 30 loci are insufficient to resolve definitely their relatedness categories, the data suggest four pairs of cryo-bank bulls as possible half-sibs. CONCLUSIONS: Our results show that even relatively limited cryo-bank storage of semen can carry allelic variation through a bottleneck. We propose a breeding scheme based on the rotation of breeding females and the division of cryo-bank bulls into three groups. Thus, if molecular data (e.g. autosomal microsatellite genotypes) for the contemporary population are available and based on relatively small-scale laboratory analyses, it is possible to avoid serious mistakes in their use for breeding applications. The approach suggested here based on the use of Yakutian cryo-bank semen can be easily extended to cryo-bank materials of other animals in future breeding programs.

Genetic differentiation between the Old and New types of Serbian Tsigai sheep.

Two Tsigai sheep populations exist in Serbia: the Old type, called Cokan, and the New type. It is assumed that the New type results from upgrading Tsigai sheep with exotic genetic material. We investigated genetic diversity and differentiation of these types by analysing 23 autosomal microsatellites. Tests for Hardy-Weinberg proportions, linkage equilibrium between genotypes across loci and the calculation of inbreeding coefficients were performed and the deficiency in the number of alleles within the Tsigai types was examined using a Wilcoxon sign-rank test. The New type displayed a higher level of genetic variability than the Cokan in terms of allele numbers, but the New Tsigai showed a pattern of heterozygosity deficiency. The positive f value for the Cokan suggests the occurrence of inbreeding in this type. The proportion of linkage disequilibrium was below that expected by chance. Exclusion of two loci in Hardy-Weinberg disequilibrium did not alter our conclusions based on the entire data set i.e. the two Tsigai types are clearly differentiated and the New Tsigai type has been influenced by crossbreeding. Therefore, the Cokan Tsigai should be considered as a distinct endangered breed in the FAO classification.

Sheep mitochondrial DNA variation in European, Caucasian, and Central Asian areas.

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.