Review: Diving into the cow hologenome to reduce methane emissions and increase sustainability.

Interest on methane emissions from livestock has increased in later years as it is an anthropogenic greenhouse gas with an important warming potential. The rumen microbiota has a large influence on the production of enteric methane. Animals harbour a second genome consisting of microbes, collectively referred to as the "microbiome". The rumen microbial community plays an important role in feed digestion, feed efficiency, methane emission and health status. This review recaps the current knowledge on the genetic control that the cow exerts on the rumen microbiota composition. Heritability estimates for the rumen microbiota composition range between 0.05 and 0.40 in the literature, depending on the taxonomical group or microbial gene function. Variables depicting microbial diversity or aggregating microbial information are also heritable within the same range. This study includes a genome-wide association analysis on the microbiota composition, considering the relative abundance of some microbial taxa previously associated to enteric methane in dairy cattle (Archaea, Dialister, Entodinium, Eukaryota, Lentisphaerae, Methanobrevibacter, Neocallimastix, Prevotella and Stentor). Host genomic regions associated with the relative abundance of these microbial taxa were identified after Benjamini-Hoschberg correction (Padj < 0.05). An in-silico functional analysis using FUMA and DAVID online tools revealed that these gene sets were enriched in tissues like brain cortex, brain amigdala, pituitary, salivary glands and other parts of the digestive system, and are related to appetite, satiety and digestion. These results allow us to have greater knowledge about the composition and function of the rumen microbiome in cattle. The state-of-the art strategies to include methane traits in the selection indices in dairy cattle populations is reviewed. Several strategies to include methane traits in the selection indices have been studied worldwide, using bioeconomical models or economic functions under theoretical frameworks. However, their incorporation in the breeding programmes is still scarce. Some potential strategies to include methane traits in the selection indices of dairy cattle population are presented. Future selection indices will need to increase the weight of traits related to methane emissions and sustainability. This review will serve as a compendium of the current state of the art in genetic strategies to reduce methane emissions in dairy cattle.

Maintenance of genetic diversity in subdivided populations using genomic coancestry matrices.

For both undivided and subdivided populations, the consensus method to maintain genetic diversity is the Optimal Contribution (OC) method. For subdivided populations, this method determines the optimal contribution of each candidate to each subpopulation to maximize global genetic diversity (which implicitly optimizes migration between subpopulations) while balancing the relative levels of coancestry between and within subpopulations. Inbreeding can be controlled by increasing the weight given to within-subpopulation coancestry (λ). Here we extend the original OC method for subdivided populations that used pedigree-based coancestry matrices, to the use of more accurate genomic matrices. Global levels of genetic diversity, measured as expected heterozygosity and allelic diversity, their distributions within and between subpopulations, and the migration pattern between subpopulations, were evaluated via stochastic simulations. The temporal trajectory of allele frequencies was also investigated. The genomic matrices investigated were (i) the matrix based on deviations of the observed number of alleles shared by two individuals from the expected number under Hardy-Weinberg equilibrium; and (ii) a matrix based on a genomic relationship matrix. The matrix based on deviations led to higher global and within-subpopulation expected heterozygosities, lower inbreeding and similar allelic diversity than the second genomic and pedigree-based matrices when a relatively high weight was given to the within-subpopulation coancestries (λ ≥ 5). Under this scenario, allele frequencies moved only slightly away from the initial frequencies. Therefore, the recommended strategy is to use the former matrix in the OC methodology giving a high weight to the within-subpopulation coancestry.

A comparison of marker-based estimators of inbreeding and inbreeding depression.

BACKGROUND: The availability of genome-wide marker data allows estimation of inbreeding coefficients (F, the probability of identity-by-descent, IBD) and, in turn, estimation of the rate of inbreeding depression (ΔID). We investigated, by computer simulations, the accuracy of the most popular estimators of inbreeding based on molecular markers when computing F and ΔID in populations under random mating, equalization of parental contributions, and artificially selected populations. We assessed estimators described by Li and Horvitz (FLH1 and FLH2), VanRaden (FVR1 and FVR2), Yang and colleagues (FYA1 and FYA2), marker homozygosity (FHOM), runs of homozygosity (FROH) and estimates based on pedigree (FPED) in comparison with estimates obtained from IBD measures (FIBD). RESULTS: If the allele frequencies of a base population taken as a reference for the computation of inbreeding are known, all estimators based on marker allele frequencies are highly correlated with FIBD and provide accurate estimates of the mean ΔID. If base population allele frequencies are unknown and current frequencies are used in the estimations, the largest correlation with FIBD is generally obtained by FLH1 and the best estimator of ΔID is FYA2. The estimators FVR2 and FLH2 have the poorest performance in most scenarios. The assumption that base population allele frequencies are equal to 0.5 results in very biased estimates of the average inbreeding coefficient but they are highly correlated with FIBD and give relatively good estimates of ΔID. Estimates obtained directly from marker homozygosity (FHOM) substantially overestimated ΔID. Estimates based on runs of homozygosity (FROH) provide accurate estimates of inbreeding and ΔID. Finally, estimates based on pedigree (FPED) show a lower correlation with FIBD than molecular estimators but provide rather accurate estimates of ΔID. An analysis of data from a pig population supports the main findings of the simulations. CONCLUSIONS: When base population allele frequencies are known, all marker-allele frequency-based estimators of inbreeding coefficients generally show a high correlation with FIBD and provide good estimates of ΔID. When base population allele frequencies are unknown, FLH1 is the marker frequency-based estimator that is most correlated with FIBD, and FYA2 provides the most accurate estimates of ΔID. Estimates from FROH are also very precise in most scenarios. The estimators FVR2 and FLH2 have the poorest performances.

Purging of deleterious burden in the endangered Iberian lynx.

SignificanceThe dynamics of deleterious variation under contrasting demographic scenarios remain poorly understood in spite of their relevance in evolutionary and conservation terms. Here we apply a genomic approach to study differences in the burden of deleterious alleles between the endangered Iberian lynx (Lynx pardinus) and the widespread Eurasian lynx (Lynx lynx). Our analysis unveils a significantly lower deleterious burden in the former species that should be ascribed to genetic purging, that is, to the increased opportunities of selection against recessive homozygotes due to the inbreeding caused by its smaller population size, as illustrated by our analytical predictions. This research provides theoretical and empirical evidence on the evolutionary relevance of genetic purging under certain demographic conditions.

Estimates of recent and historical effective population size in turbot, seabream, seabass and carp selective breeding programmes.

BACKGROUND: The high fecundity of fish species allows intense selection to be practised and therefore leads to fast genetic gains. Based on this, numerous selective breeding programmes have been started in Europe in the last decades, but in general, little is known about how the base populations of breeders have been built. Such knowledge is important because base populations can be created from very few individuals, which can lead to small effective population sizes and associated reductions in genetic variability. In this study, we used genomic information that was recently made available for turbot (Scophthalmus maximus), gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and common carp (Cyprinus carpio) to obtain accurate estimates of the effective size for commercial populations. METHODS: Restriction-site associated DNA sequencing data were used to estimate current and historical effective population sizes. We used a novel method that considers the linkage disequilibrium spectrum for the whole range of genetic distances between all pairs of single nucleotide polymorphisms (SNPs), and thus accounts for potential fluctuations in population size over time. RESULTS: Our results show that the current effective population size for these populations is small (equal to or less than 50 fish), potentially putting the sustainability of the breeding programmes at risk. We have also detected important drops in effective population size about five to nine generations ago, most likely as a result of domestication and the start of selective breeding programmes for these species in Europe. CONCLUSIONS: Our findings highlight the need to broaden the genetic composition of the base populations from which selection programmes start, and suggest that measures designed to increase effective population size within all farmed populations analysed here should be implemented in order to manage genetic variability and ensure the sustainability of the breeding programmes.

Changes in Allele Frequencies When Different Genomic Coancestry Matrices Are Used for Maintaining Genetic Diversity.

A main objective in conservation programs is to maintain genetic variability. This can be achieved using the Optimal Contributions (OC) method that optimizes the contributions of candidates to the next generation by minimizing the global coancestry. However, it has been argued that maintaining allele frequencies is also important. Different genomic coancestry matrices can be used on OC and the choice of the matrix will have an impact not only on the genetic variability maintained, but also on the change in allele frequencies. The objective of this study was to evaluate, through stochastic simulations, the genetic variability maintained and the trajectory of allele frequencies when using two different genomic coancestry matrices in OC to minimize the loss of diversity: (i) the matrix based on deviations of the observed number of alleles shared between two individuals from the expected numbers under Hardy-Weinberg equilibrium (θLH); and (ii) the matrix based on VanRaden's genomic relationship matrix (θVR). The results indicate that the use of θLH resulted in a higher genetic variability than the use of θVR. However, the use of θVR maintained allele frequencies closer to those in the base population than the use of θLH.

The value of genomic relationship matrices to estimate levels of inbreeding.

BACKGROUND: Genomic relationship matrices are used to obtain genomic inbreeding coefficients. However, there are several methodologies to compute these matrices and there is still an unresolved debate on which one provides the best estimate of inbreeding. In this study, we investigated measures of inbreeding obtained from five genomic matrices, including the Nejati-Javaremi allelic relationship matrix (FNEJ), the Li and Horvitz matrix based on excess of homozygosity (FL&H), and the VanRaden (methods 1, FVR1, and 2, FVR2) and Yang (FYAN) genomic relationship matrices. We derived expectations for each inbreeding coefficient, assuming a single locus model, and used these expectations to explain the patterns of the coefficients that were computed from thousands of single nucleotide polymorphism genotypes in a population of Iberian pigs. RESULTS: Except for FNEJ, the evaluated measures of inbreeding do not match with the original definitions of inbreeding coefficient of Wright (correlation) or Malécot (probability). When inbreeding coefficients are interpreted as indicators of variability (heterozygosity) that was gained or lost relative to a base population, both FNEJ and FL&H led to sensible results but this was not the case for FVR1, FVR2 and FYAN. When variability has increased relative to the base, FVR1, FVR2 and FYAN can indicate that it decreased. In fact, based on FYAN, variability is not expected to increase. When variability has decreased, FVR1 and FVR2 can indicate that it has increased. Finally, these three coefficients can indicate that more variability than that present in the base population can be lost, which is also unreasonable. The patterns for these coefficients observed in the pig population were very different, following the derived expectations. As a consequence, the rate of inbreeding depression estimated based on these inbreeding coefficients differed not only in magnitude but also in sign. CONCLUSIONS: Genomic inbreeding coefficients obtained from the diagonal elements of genomic matrices can lead to inconsistent results in terms of gain and loss of genetic variability and inbreeding depression estimates, and thus to misleading interpretations. Although these matrices have proven to be very efficient in increasing the accuracy of genomic predictions, they do not always provide a useful measure of inbreeding.

On the estimation of inbreeding depression using different measures of inbreeding from molecular markers.

The inbreeding coefficient (F) of individuals can be estimated from molecular marker data, such as SNPs, using measures of homozygosity of individual markers or runs of homozygosity (ROH) across the genome. These different measures of F can then be used to estimate the rate of inbreeding depression (ID) for quantitative traits. Some recent simulation studies have investigated the accuracy of this estimation with contradictory results. Whereas some studies suggest that estimates of inbreeding from ROH account more accurately for ID, others suggest that inbreeding measures from SNP-by-SNP homozygosity giving a large weight to rare alleles are more accurate. Here, we try to give more light on this issue by carrying out a set of computer simulations considering a range of population genetic parameters and population sizes. Our results show that the previous studies are indeed not contradictory. In populations with low effective size, where relationships are more tight and selection is relatively less intense, F measures based on ROH provide very accurate estimates of ID whereas SNP-by-SNP-based F measures with high weight to rare alleles can show substantial upwardly biased estimates of ID. However, in populations of large effective size, with more intense selection and trait allele frequencies expected to be low if they are deleterious for fitness because of purifying selection, average estimates of ID from SNP-by-SNP-based F values become unbiased or slightly downwardly biased and those from ROH-based F values become slightly downwardly biased. The noise attached to all these estimates, nevertheless, can be very high in large-sized populations. We also investigate the relationship between the different F measures and the homozygous mutation load, which has been suggested as a proxy of inbreeding depression.

Genetic parameters for uniformity of harvest weight in Pacific white shrimp (Litopenaeus vannamei).

BACKGROUND: Uniformity of body weight is a trait of great economic importance in the production of white shrimp (Litopenaeus vannamei). A necessary condition to improve this trait through selective breeding is the existence of genetic variability for the environmental variance of body weight. Although several studies have reported such variability in other aquaculture species, to our knowledge, no estimates are available for shrimp. Our aim in this study was to estimate the genetic variance for weight uniformity in a farmed population of shrimp to determine the potential of including this trait in the selection program. We also estimated the genetic correlation of weight uniformity between two environments (selection nucleus and commercial population). METHODS: The database contained phenotypic records for body weight on 51,346 individuals from the selection nucleus and 38,297 individuals from the commercial population. A double hierarchical generalized linear model was used to analyse weight uniformity in the two environments. Fixed effects included sex and year for the nucleus data and sex and year-pond combination for the commercial data. Environmental and additive genetic effects were included as random effects. RESULTS: The estimated genetic variance for weight uniformity was greater than 0 (0.06 ± 0.01) in both the nucleus and commercial populations and the genetic coefficient of variation for the residual variance was 0.25 ± 0.01. The genetic correlation between weight and weight uniformity was close to zero in both environments. The estimate of the genetic correlation of weight uniformity between the two environments (selection nucleus and commercial population) was 0.64 ± 0.06. CONCLUSIONS: The existence of genetic variance for weight uniformity suggests that genetic improvement of this trait is possible. Selection for weight uniformity should not decrease weight, given the near zero genetic correlation between these two traits. The strong genetic correlation of weight uniformity between the two environments indicates that response to selection for uniformity in the nucleus will be at least partially transmitted to the commercial population if this trait is included in the breeding goal.

Optimum mating designs for exploiting dominance in genomic selection schemes for aquaculture species.

BACKGROUND: In commercial fish, dominance effects could be exploited by predicting production abilities of the offspring that would be generated by different mating pairs and choosing those pairs that maximise the average offspring phenotype. Consequently, matings would be performed to reduce inbreeding depression. This can be achieved by applying mate selection (MS) that combines selection and mating decisions in a single step. An alternative strategy to MS would be to apply minimum coancestry mating (MCM) after selection based on estimated breeding values. The objective of this study was to evaluate, by computer simulations, the potential benefits that can be obtained by implementing MS or MCM based on genomic data for exploiting dominance effects when creating commercial fish populations that are derived from a breeding nucleus. METHODS: The selected trait was determined by a variable number of loci with additive and dominance effects. The population consisted of 50 full-sib families with 30 offspring each. Males and females with the highest estimated genomic breeding values were selected in the nucleus and paired using the MCM strategy. Both MCM and MS were used to create the commercial population. RESULTS: For a moderate number of SNPs, equal or even higher mean phenotypic values are obtained by selecting on genomic breeding values and then applying MCM than by using MS when the trait exhibited substantial inbreeding depression. This could be because MCM leads to high levels of heterozygosity across the whole genome, even for loci affecting the trait that are in linkage equilibrium with the SNPs. In contrast, MS specifically promotes heterozygosity for SNPs for which a dominance effect has been detected. CONCLUSIONS: In most scenarios, for the management of aquaculture breeding programs it seems advisable to follow the MCM strategy when creating the commercial population, especially for traits with large inbreeding depression. Moreover, MCM has the appealing property of reducing inbreeding levels, with a corresponding reduction in inbreeding depression for traits beyond those included in the selection objective.

The concept of effective population size loses its meaning in the context of optimal management of diversity using molecular markers.

Effective population size is a key parameter in conservation genetics. In the management of conservation programs using pedigree information, there is a consensus that the optimal method for maximizing effective population size is to calculate the contribution of each potential parent (the number of offspring that each individual leaves to the next generation) by minimizing the global pedigree-based coancestry between potential parents weighted by their contributions. When using molecular data, the optimal method for managing genetic diversity will remain the same but now the molecular coancestry calculated from markers will replace the pedigree-based coancestry. However, in this situation, the concept of effective population size loses its meaning because with optimal molecular management, genetic diversity increases in early generations and therefore effective population size takes negative values. Furthermore, in the long term, the molecular effective population size does not attain an asymptotic value but it shows an unpredictable behaviour.

Identification of genomic regions regulating sex determination in Atlantic salmon using high density SNP data.

BACKGROUND: A complete understanding of the genetic basis for sexual determination and differentiation is necessary in order to implement efficient breeding schemes at early stages of development. Atlantic salmon belongs to the family Salmonidae of fishes and represents a species of great commercial value. Although the species is assumed to be male heterogametic with XY sex determination, the precise genetic basis of sexual development remains unclear. The complexity is likely associated to the relatively recent salmonid specific whole genome duplication that may be responsible for certain genome instability. This instability together with the capacity of the sex-determining gene to move across the genome as reported by previous studies, may explain that sexual development genes are not circumscribed to the same chromosomes in all members of the species. In this study, we have used a 220 K SNP panel developed for Atlantic salmon to identify the chromosomes explaining the highest proportion of the genetic variance for sex as well as candidate regions and genes associated to sexual development in this species. RESULTS: Results from regional heritability analysis showed that the chromosomes explaining the highest proportion of variance in these populations were Ssa02 (heritability = 0.42, SE = 0.12) and Ssa21 (heritability = 0.26, SE = 0.11). After pruning by linkage disequilibrium, genome-wide association analyses revealed 114 SNPs that were significantly associated with sex, being Ssa02 the chromosome containing a greatest number of regions. Close examination of the candidate regions evidenced important genes related to sex in other species of Class Actinopterygii, including SDY, genes from family SOX, RSPO1, ESR1, U2AF2A, LMO7, GNRH-R, DND and FIGLA. CONCLUSIONS: The combined results from regional heritability analysis and genome-wide association have provided new advances in the knowledge of the genetic regulation of sex determination in Atlantic salmon, supporting that Ssa02 is the candidate chromosome for sex in this species and suggesting an alternative population lineage in Spanish wild populations according to the results from Ssa21.

Disentangling Genetic Variation for Resistance and Endurance to Scuticociliatosis in Turbot Using Pedigree and Genomic Information.

Selective breeding for improving host responses to infectious pathogens is a promising option for disease control. In fact, disease resilience, the ability of a host to survive or cope with infectious challenge, has become a highly desirable breeding goal. However, resilience is a complex trait composed of two different host defence mechanisms, namely resistance (the ability of a host to avoid becoming infected or diseased) and endurance (the ability of an infected host to survive the infection). While both could be targeted for genetic improvement, it is currently unknown how they contribute to survival, as reliable estimates of genetic parameters for both traits obtained simultaneously are scarce. A difficulty lies in obtaining endurance phenotypes for genetic analyses. In this study, we present the results from an innovative challenge test carried out in turbot whose design allowed disentangling the genetic basis of resistance and endurance to Philasterides dicentrarchi, a parasite causing scuticociliatosis that leads to substantial economic losses in the aquaculture industry. A noticeable characteristic of the parasite is that it causes visual signs that can be used for disentangling resistance and endurance. Our results showed the existence of genetic variation for both traits (heritability = 0.26 and 0.12 for resistance and endurance, respectively) and for the composite trait resilience (heritability = 0.15). The genetic correlation between resistance and resilience was very high (0.90) indicating that both are at a large extent the same trait, but no significant genetic correlation was found between resistance and endurance. A total of 18,125 SNPs obtained from 2b-RAD sequencing enabled genome-wide association analyses for detecting QTLs controlling the three traits. A candidate QTL region on linkage group 19 that explains 33% of the additive genetic variance was identified for resilience. The region contains relevant genes related to immune response and defence mechanisms. Although no significant associations were found for resistance, the pattern of association was the same as for resilience. For endurance, one significant association was found on linkage group 2. The accuracy of genomic breeding values was also explored for resilience, showing that it increased by 12% when compared with the accuracy of pedigree-based breeding values. To our knowledge, this is the first study in turbot disentangling the genetic basis of resistance and endurance to scuticociliatosis.

Genetic differences in host infectivity affect disease spread and survival in epidemics.

Survival during an epidemic is partly determined by host genetics. While quantitative genetic studies typically consider survival as an indicator for disease resistance (an individual's propensity to avoid becoming infected or diseased), mortality rates of populations undergoing an epidemic are also affected by endurance (the propensity of diseased individual to survive the infection) and infectivity (i.e. the propensity of an infected individual to transmit disease). Few studies have demonstrated genetic variation in disease endurance, and no study has demonstrated genetic variation in host infectivity, despite strong evidence for considerable phenotypic variation in this trait. Here we propose an experimental design and statistical models for estimating genetic diversity in all three host traits. Using an infection model in fish we provide, for the first time, direct evidence for genetic variation in host infectivity, in addition to variation in resistance and endurance. We also demonstrate how genetic differences in these three traits contribute to survival. Our results imply that animals can evolve different disease response types affecting epidemic survival rates, with important implications for understanding and controlling epidemics.

Abdominal obesity is a common finding in normal and overweight subjects of Chile and is associated with increased frequency of cardiometabolic risk factors.

BACKGROUND/OBJECTIVES: Abdominal obesity (AO) is associated with elevated risk for cardiovascular diseases; however, this association is less clear for non-obese people. We estimated the association of AO and cardiovascular risk factors (CVRF) and disease in non-obese adult individuals from Chile. SUBJECTS/METHODS: 5248 adults (15 years of age or older) of both sexes from the Chilean National Health Survey (October 2009 -September 2010, response rate 85%.) were included. Information on myocardial infarction and stroke was self-reported. BMI, waist circumference (WC), arterial pressure, plasma glucose, and cholesterol levels were measured. Predictive accuracy of WC was evaluated by area under curve of receiver operating characteristic analysis and cut off points were established by Youden Index. Relationship between AO and CVRF was analyzed by Chi-squared tests. RESULTS: Normal weight/overweight/obesity were present in 34.4%/45.2%/18.1% of men and 33.4%/33.6%/27.5% of women. Predictive accuracy of WC to identify at least one CVRF was 0.70/0.67 and optimal cutoff points for WC in non-obese subjects were 91/83 cm in men/women, respectively. AO was present in 98.2%/99.1% of obese, 70.5%/77.4% of overweight and 12.4%/16.4% of normal weight men/women. AO was associated with increased frequency of CVRF in overweight men (6/8 and stroke) and women (4/8) and higher frequency in normal weight men (8/8 and myocardial infarction/stroke) and women (6/8 and myocardial infarction). CONCLUSIONS: WC cutoff points calculated for non-obese chilean population discriminate more differences in CVRF in normal weight woman. AO significantly increases the frequency of CVRF and diseases in overweight and especially normal weight individuals. WC can be used as a low cost, feasible and reproducible predictor for CVRF in non-obese individuals in most clinical settings.

Novel efficient genome-wide SNP panels for the conservation of the highly endangered Iberian lynx.

BACKGROUND: The Iberian lynx (Lynx pardinus) has been acknowledged as the most endangered felid species in the world. An intense contraction and fragmentation during the twentieth century left less than 100 individuals split in two isolated and genetically eroded populations by 2002. Genetic monitoring and management so far have been based on 36 STRs, but their limited variability and the more complex situation of current populations demand more efficient molecular markers. The recent characterization of the Iberian lynx genome identified more than 1.6 million SNPs, of which 1536 were selected and genotyped in an extended Iberian lynx sample. METHODS: We validated 1492 SNPs and analysed their heterozygosity, Hardy-Weinberg equilibrium, and linkage disequilibrium. We then selected a panel of 343 minimally linked autosomal SNPs from which we extracted subsets optimized for four different typical tasks in conservation applications: individual identification, parentage assignment, relatedness estimation, and admixture classification, and compared their power to currently used STR panels. RESULTS: We ascribed 21 SNPs to chromosome X based on their segregation patterns, and identified one additional marker that showed significant differentiation between sexes. For all applications considered, panels of autosomal SNPs showed higher power than the currently used STR set with only a very modest increase in the number of markers. CONCLUSIONS: These novel panels of highly informative genome-wide SNPs provide more powerful, efficient, and flexible tools for the genetic management and non-invasive monitoring of Iberian lynx populations. This example highlights an important outcome of whole-genome studies in genetically threatened species.

Neuroprotective Effect of Hydroxytyrosol in Experimental Diabetes Mellitus.

The aim of the study was to analyze the possible neuroprotective effect of hydroxytyrosol (HT) in diabetic animals in a model of hypoxia-reoxygenation. Rats (10 animals/group) were distributed in five groups: nondiabetic rats, control diabetic rats (DR), and DR rats treated for 2 months with 1, 5, or 10 mg/kg/day po HT. At the end of follow-up, an experimental model of hypoxia-reoxygenation in brain slices was tested. The DR group showed increased cell death, oxidative and nitrosative stress, and an increase in brain inflammatory mediators. These alterations were significantly greater in DR than in normoglycemic animals. HT significantly reduced oxidative (38.5-52.4% lipid peroxidation) and nitrosative stress (48.0-51.0% nitric oxide and 43.9-75.2% peroxynitrite concentration) and brain inflammatory mediators (18.6-40.6% prostaglandin E2 and 17.0-65.0% interleukin 1ß concentration). Cell death was reduced by 25.9, 37.5, and 41.0% after the administration of 1, 5, or 10 mg/kg/day. The administration of HT in rats with experimental diabetes thus had a neuroprotective effect.

The use of genomic coancestry matrices in the optimisation of contributions to maintain genetic diversity at specific regions of the genome.

BACKGROUND: Optimal contribution methods have proved to be very efficient for controlling the rates at which coancestry and inbreeding increase and therefore, for maintaining genetic diversity. These methods have usually relied on pedigree information for estimating genetic relationships between animals. However, with the large amount of genomic information now available such as high-density single nucleotide polymorphism (SNP) chips that contain thousands of SNPs, it becomes possible to calculate more accurate estimates of relationships and to target specific regions in the genome where there is a particular interest in maximising genetic diversity. The objective of this study was to investigate the effectiveness of using genomic coancestry matrices for: (1) minimising the loss of genetic variability at specific genomic regions while restricting the overall loss in the rest of the genome; or (2) maximising the overall genetic diversity while restricting the loss of diversity at specific genomic regions. RESULTS: Our study shows that the use of genomic coancestry was very successful at minimising the loss of diversity and outperformed the use of pedigree-based coancestry (genetic diversity even increased in some scenarios). The results also show that genomic information allows a targeted optimisation to maintain diversity at specific genomic regions, whether they are linked or not. The level of variability maintained increased when the targeted regions were closely linked. However, such targeted management leads to an important loss of diversity in the rest of the genome and, thus, it is necessary to take further actions to constrain this loss. Optimal contribution methods also proved to be effective at restricting the loss of diversity in the rest of the genome, although the resulting rate of coancestry was higher than the constraint imposed. CONCLUSIONS: The use of genomic matrices when optimising contributions permits the control of genetic diversity and inbreeding at specific regions of the genome through the minimisation of partial genomic coancestry matrices. The formula used to predict coancestry in the next generation produces biased results and therefore it is necessary to refine the theory of genetic contributions when genomic matrices are used to optimise contributions.

Evaluation of the linkage-disequilibrium method for the estimation of effective population size when generations overlap: an empirical case.

BACKGROUND: Within the genetic methods for estimating effective population size (N e ), the method based on linkage disequilibrium (LD) has advantages over other methods, although its accuracy when applied to populations with overlapping generations is a matter of controversy. It is also unclear the best way to account for mutation and sample size when this method is implemented. Here we have addressed the applicability of this method using genome-wide information when generations overlap by profiting from having available a complete and accurate pedigree from an experimental population of Iberian pigs. Precise pedigree-based estimates of N e were considered as a baseline against which to compare LD-based estimates. METHODS: We assumed six different statistical models that varied in the adjustments made for mutation and sample size. The approach allowed us to determine the most suitable statistical model of adjustment when the LD method is used for species with overlapping generations. A novel approach used here was to treat different generations as replicates of the same population in order to assess the error of the LD-based N e estimates. RESULTS: LD-based N e estimates obtained by estimating the mutation parameter from the data and by correcting sample size using the 1/2n term were the closest to pedigree-based estimates. The N e at the time of the foundation of the herd (26 generations ago) was 20.8 ± 3.7 (average and SD across replicates), while the pedigree-based estimate was 21. From that time on, this trend was in good agreement with that followed by pedigree-based N e. CONCLUSIONS: Our results showed that when using genome-wide information, the LD method is accurate and broadly applicable to small populations even when generations overlap. This supports the use of the method for estimating N e when pedigree information is unavailable in order to effectively monitor and manage populations and to early detect population declines. To our knowledge this is the first study using replicates of empirical data to evaluate the applicability of the LD method by comparing results with accurate pedigree-based estimates.

Genome-wide estimates of coancestry, inbreeding and effective population size in the Spanish Holstein population.

Estimates of effective population size in the Holstein cattle breed have usually been low despite the large number of animals that constitute this breed. Effective population size is inversely related to the rates at which coancestry and inbreeding increase and these rates have been high as a consequence of intense and accurate selection. Traditionally, coancestry and inbreeding coefficients have been calculated from pedigree data. However, the development of genome-wide single nucleotide polymorphisms has increased the interest of calculating these coefficients from molecular data in order to improve their accuracy. In this study, genomic estimates of coancestry, inbreeding and effective population size were obtained in the Spanish Holstein population and then compared with pedigree-based estimates. A total of 11,135 animals genotyped with the Illumina BovineSNP50 BeadChip were available for the study. After applying filtering criteria, the final genomic dataset included 36,693 autosomal SNPs and 10,569 animals. Pedigree data from those genotyped animals included 31,203 animals. These individuals represented only the last five generations in order to homogenise the amount of pedigree information across animals. Genomic estimates of coancestry and inbreeding were obtained from identity by descent segments (coancestry) or runs of homozygosity (inbreeding). The results indicate that the percentage of variance of pedigree-based coancestry estimates explained by genomic coancestry estimates was higher than that for inbreeding. Estimates of effective population size obtained from genome-wide and pedigree information were consistent and ranged from about 66 to 79. These low values emphasize the need of controlling the rate of increase of coancestry and inbreeding in Holstein selection programmes.