Biobanking animal genetic resources: critical infrastructure and growth opportunities.

National animal gene banks have acquired substantial quantities of germplasm that protect and preserve a wide range of livestock breeds. New challenges and growth opportunities are emerging. A key challenge will be increased gene bank use, but this requires increased characterization of phenotypes and genotypes for populations and collections.

Whole-Genome Resequencing of Worldwide Wild and Domestic Sheep Elucidates Genetic Diversity, Introgression, and Agronomically Important Loci.

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3'-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

Assessment of FecGE genotypes on reproductive traits in Brazilian Morada Nova and Santa Inês sheep.

The aim of this study was to evaluate the effect of the polymorphic FecGE allele on reproductive traits in Santa Inês and Morada Nova ewes. The traits evaluated were as follows: total progeny weights at birth (PWB) and weaning (PWW) and progeny survival rates at birth (PSRB) and weaning (PSRW). A total of 389 animals, belonging to two Santa Inês herds and one Morada Nova herd, were genotyped. There was a difference between the averages for all the traits studied regarding type of parturition, herd/breed, genotype/herd, and genotype/type of parturition. For each additional progeny, if the female was FecGE/E, the PWB decreased by 1.02 kg and the PWW by 3.16 kg, also with a 0.04% reduction in PSRB and no change in PSRW. If the female was FecGE/+, the reduction in PWB was 0.24 kg, with an increase in PSRW by 0.11%, but no change in PWW and PSRB. In general, these results demonstrate that FecG+/+ females have a better ability to increase their number of progenies without reducing PWB and PWW (also similar to FecGE/+). Thus, it is suggested that further studies on the association between the traits of interest and candidate genes in sheep should be carried out so that the regions which have the greatest effect on the expression of these traits are actually identified. It was not possible to verify the effect of the FecGE allele on the PWB, PWW, PSRB, and PSRW in these Morada Nova and Santa Inês herds.

Selection signatures for heat tolerance in Brazilian horse breeds.

Since domestication, horse breeds have adapted to their environments and differentiated from one another. This paper uses two methods to detect selection signatures in 23 horse breeds, eight of which are Brazilian (610 animals), both cold-blooded and warm-blooded, from temperate and tropical regions. These animals were genotyped using the GGP Equine BeadChip and we analysed the data by Principal Component Analysis (PCA). The samples were separated into groups based on their geographical area of origin and PCA results studied. The genomic regions under selection were detected by hapFLK and PCAdapt methodologies, identifying six regions under selection with at least one Brazilian horse breed. These regions contain genes associated with heat tolerance, skin colour, body size, energy production/metabolism, genes involved in protein degradation/turnover/DNA repair, genes reducing the impact of oxidative stress/cellular repair, and transcriptional regulation. This work confirmed LCORL and NCAPG gene regions in previous studies associated with body size on Equine Chromosome Autosome 3 (ECA3). On the same ECA3, a region implicating genes linked to coat colour was identified, also previously related to heat stress. Regions with genes coding heat shock proteins were found on ECA1 and 2, and many candidate genes for oxidation-reduction which are a natural response to heat stress. However, a larger sample size and whole-genome SNPs are needed to understand better and identify new candidate regions as well as their functional relation with heat tolerance.

Use of logistic models to evaluate the response of superovulation treatment and embryo production in Santa Inês ewes.

The study aimed to verify the influence of the FecGE mutation in superovulated ewes and to evaluate the probability of logistic models to determine the response capacity of these ewes to superovulatory treatment. Santa Inês ewes (n = 29) were genotyped for the FecGE mutation and separated for their genotype group in carriers of the mutant E allele (FecGE/E, FecG+/E) and non-carrier (FecG+/+) alleles. The ewes underwent hormonal treatment for superovulation. Aside from the genotypes, variables included in the statistical model were reproductive status (empty, early lactation, or late lactation), age (> or < 6 years), and number of births (nulliparous, primiparous, multiparous). The carriers of the mutation could be discriminated from the non-carriers based on the number of corpora lutea, rate of frozen embryos, and fecundity. Recovery rate was significantly higher (P < 0.05) in FecGE/E (94.31%) compared to FecG+/E (63.15%) and FecG+/+ (61.90%) (P < 0.05), whereas fecundity rate of FecG+/+ ewes (50.76%) was significantly higher than FecG+/E (18.96%) and FecGE/E (32.53%) (P < 0.05). We determined in this study that the response to superovulation and embryo production can be discriminated between FecGE/E and FecG+/E ewes in relation to the FecG+/+ genotype. Logistic models that included reproductive status and mutation, or reproductive status and age, or reproductive status and number of births were effective in predicting the response to superovulatory treatment.

Fine-scale genetic diversity of the Brazilian Pantaneiro horse breed adapted to flooded regions.

Among the animal species first introduced in Brazil during the country's discovery, horses (Equus caballus) stand out because of their evolutionary history and relationship with humans. Among the Brazilian horse breeds, the Pantaneiro draws attention due to its adaptative traits. Blood samples of 116 Pantaneiro horses were divided into six populations based on their sampling location, aiming to identify the existence of genetic structure and quantify genetic diversity within and between them. Populations were compared to elucidate genetic variability and differentiation better and assess the impact of Pantanal's natural geographic barriers on gene flow between populations. Data from the GGP Equine BeadChip (Geneseek-Neogen, 65.157 SNPs) was used to assess basic diversity parameters, genetic distance (FST), principal component analysis (PCA), and population structure (ADMIXTURE) for the sampled animals. Mantel test was also performed to investigate the correlation between the populations' genetic and geographic distances. Results showed high genetic variability in all populations, with elevated levels of admixture in their structure. High levels of admixture make it challenging to establish a racial pattern and, consequently, populations within the breed, being that only one of the populations differentiated itself from the others. No significant correlations between genetic and geographic distances were observed, indicating that environmental barriers did not hinder gene flow between populations, and neither farmers' selection practices might have change breed genetic composition significantly. Low genetic distance and similar heterozygosity values were observed among populations, suggesting strong genetic proximity and low differentiation. Thereby, the Pantaneiro breed does not exhibit genetic subpopulations and could be considered, for conservation purposes, a single big population in the Panatnal region. This study will support sampling strategies for National genebank.

Integration of georeferenced and genetic data for the management of biodiversity in sheep genetic resources in Brazil.

There are few animal germplasm/gene bank collections in Brazil, and basic studies are needed to attend the future internal and external demands from international partners. The aim of this work was to validate a "proof of concept" that integrates spatial (georeferenced data) and genetic data regarding the local of origin from 3518 DNA samples from 17 different genetic groups or breeds of sheep in the Brazilian Germplasm bank. Spatialisation shows that not all genetic groups have samples in the bank, and collection is concentrated in the conservation nuclei spread nationwide. Only 21% of states with a specific breed have samples in the gene bank. The mean number of animals sampled per collection was 32, while the mean distance travelled to collect samples was 262 km from the conservation nuclei. For example, the Brazilian Somali were only collected in the conservation nucleus in Ceará State. No samples were collected to date for the Cariri breed, which is recognised by the Brazilian Ministry of Agriculture. Only two farms and one breed in the bank are from the northern region. Of the 27 states, there are samples in the gene bank of sheep from 13, so several states have no samples, requiring collection from herds outside the official system of conservation to make sure that studies using this germplasm realised are not biased. Significant genetic differences are seen above 332 km, which should guide future sampling efforts. Suggestions are given for improving the quantity, quality and diversity of samples in the gene bank.

Mixed ancestry from wild and domestic lineages contributes to the rapid expansion of invasive feral swine.

Invasive alien species are a significant threat to both economic and ecological systems. Identifying the processes that give rise to invasive populations is essential for implementing effective control strategies. We conducted an ancestry analysis of invasive feral swine (Sus scrofa, Linnaeus, 1758), a highly destructive ungulate that is widely distributed throughout the contiguous United States, to describe introduction pathways, sources of newly emergent populations and processes contributing to an ongoing invasion. Comparisons of high-density single nucleotide polymorphism genotypes for 6,566 invasive feral swine to a comprehensive reference set of S. scrofa revealed that the vast majority of feral swine were of mixed ancestry, with dominant genetic associations to Western heritage breeds of domestic pig and European populations of wild boar. Further, the rapid expansion of invasive feral swine over the past 30 years was attributable to secondary introductions from established populations of admixed ancestry as opposed to direct introductions of domestic breeds or wild boar. Spatially widespread genetic associations of invasive feral swine to European wild boar deviated strongly from historical S. scrofa introduction pressure, which was largely restricted to domestic pigs with infrequent, localized wild boar releases. The deviation between historical introduction pressure and contemporary genetic ancestry suggests wild boar-hybridization may contribute to differential fitness in the environment and heightened invasive potential for individuals of admixed domestic pig-wild boar ancestry.

Phylogeography of feral Monteiro pig in the Brazilian Pantanal Ecosystem.

The Monteiro is a feral pig found in the Brazilian Pantanal ecosystem. The goal of this research is to generate data and knolewdge related to animal populations wich can be used for management and development of an in vitro conservation program for animal resourses at Pantanal ecosystem. The present study evaluated animals sampled from 10 distinct locations within the region, using 19 microsatellite markers (N = 189) and the control region of mitochondrial DNA (mtDNA) (N = 392). Low genetic differences were found between populations with the microsatellite data. The FST range was between 0.009 and 0.063 (p-value < 0.05). The Mantel test corroborated with previous results, as low correlations between genetic and geographic distances were observed (r2 = 0.2309, p = 0.06). Bayesian analysis for genetic structure identification placed the Monteiro pigs into three main clusters (MOB, Pop 1 and all others Pantanal populations). Most of the Monteiro pigs share a single European haplotype as seen by mtDNA analyses. This haplotype is not exclusive, as it is shared with other swine populations (commercial and other locally adapted breeds). Monteiro populations from different geographic locations within Pantanal are not isolated and can be considered as a large unique population. Since animals roam freely to seek food and water, or even due to seasonal flooding of their habitat, the Monteiro populations presented absence of major genetic structure and evidence of high gene flow. These results can be used to create a management plan and in situ and ex situ conservation program for conservation and use of the Monteiro breed in the Pantanal ecosystem.

Habitat partitioning during character displacement between the sexes.

Ecological differences between the sexes are often interpreted as evidence of within-species ecological character displacement (ECD), a hypothesis with almost no direct tests. Here, we experimentally test two predictions that are direct corollaries of ECD between the sexes, in a salamander. First, we find support for the prediction that each sex has a growth rate advantage in the aquatic microhabitat where it is most commonly found. Second, we test the prediction that selection for ECD in the breeding environment may affect partial migration out of this environment. We found that phenotype-dependent migration resulted in a shift in the phenotypic distribution across treatments, with the highest sexual dimorphism occurring among residents at high founding density, suggesting that migration and ECD can both be driven by competition. Our work illustrates how complex patterns of habitat partitioning evolve during ECD between the sexes and suggest ECD and partial migration can interact to effect both ecological dynamics and evolution of sexual dimorphism.

Complete mitochondrial genome from South American catfish Pseudoplatystoma reticulatum (Eigenmann & Eigenmann) and its impact in Siluriformes phylogenetic tree.

The cachara (Pseudoplatystoma reticulatum) is a Neotropical freshwater catfish from family Pimelodidae (Siluriformes) native to Brazil. The species is of relative economic importance for local aquaculture production and basic biological information is under development to help boost efforts to domesticate and raise the species in commercial systems. The complete cachara mitochondrial genome was obtained by assembling Illumina RNA-seq data from pooled samples. The full mitogenome was found to be 16,576 bp in length, showing the same basic structure, order, and genetic organization observed in other Pimelodidae, with 13 protein-coding genes, 2 rNA genes, 22 trNAs, and a control region. Observed base composition was 24.63% T, 28.47% C, 31.45% A, and 15.44% G. With the exception of NAD6 and eight tRNAs, all of the observed mitochondrial genes were found to be coded on the H strand. A total of 107 SNPs were identified in P. reticulatum mtDNA, 67 of which were located in coding regions. Of these SNPs, 10 result in amino acid changes. Analysis of the obtained sequence with 94 publicly available full Siluriformes mitogenomes resulted in a phylogenetic tree that generally agreed with available phylogenetic proposals for the order. The first report of the complete Pseudoplatystoma reticulatum mitochondrial genome sequence revealed general gene organization, structure, content, and order similar to most vertebrates. Specific sequence and content features were observed and may have functional attributes which are now available for further investigation.

D-loop haplotype diversity in Brazilian horse breeds.

The first horses were brought to Brazil by the colonizers after 1534. Over the centuries, these animals evolved and adapted to local environmental conditions usually unsuitable for exotic breeds, thereby originating locally adapted Brazilian breeds. The present work represents the first description of maternal genetic diversity in these horse breeds based on D-loop sequences. A D-Loop HSV-I fragment of 252 bp, from 141 horses belonging to ten Brazilian breeds / genetic groups (locally adapted and specialized breeds) were analysed. Thirty-five different haplotypes belonging to 18 haplogroups were identified with 33 polymorphic sites. Haplotype diversity (varying from 0.20 to 0.96) and nucleotide diversity (varying from 0.0039 to 0.0239) was lower for locally adapted than for specialized breeds, with the same pattern observed for FST values. Haplogroups identified in Brazilian breeds are in agreement with previous findings in South American samples. The low variability observed mainly in locally adapted breeds, indicates that, to ensure conservation of these breeds, careful reproductive management is needed. Additional genetic characterization studies are required to support accurate decision-making.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits.

BACKGROUND: Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance. RESULTS: Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively. CONCLUSIONS: Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

Genome-wide analysis of the world's sheep breeds reveals high levels of historic mixture and strong recent selection.

Through their domestication and subsequent selection, sheep have been adapted to thrive in a diverse range of environments. To characterise the genetic consequence of both domestication and selection, we genotyped 49,034 SNP in 2,819 animals from a diverse collection of 74 sheep breeds. We find the majority of sheep populations contain high SNP diversity and have retained an effective population size much higher than most cattle or dog breeds, suggesting domestication occurred from a broad genetic base. Extensive haplotype sharing and generally low divergence time between breeds reveal frequent genetic exchange has occurred during the development of modern breeds. A scan of the genome for selection signals revealed 31 regions containing genes for coat pigmentation, skeletal morphology, body size, growth, and reproduction. We demonstrate the strongest selection signal has occurred in response to breeding for the absence of horns. The high density map of genetic variability provides an in-depth view of the genetic history for this important livestock species.

Identification of selection signatures in livestock species.

The identification of regions that have undergone selection is one of the principal goals of theoretical and applied evolutionary genetics. Such studies can also provide information about the evolutionary processes involved in shaping genomes, as well as physical and functional information about genes/genomic regions. Domestication followed by breed formation and selection schemes has allowed the formation of very diverse livestock breeds adapted to a wide variety of environments and with special characteristics. The advances in genomics in the last five years have enabled the development of several methods to detect selection signatures and have resulted in the publication of a considerable number of studies involving livestock species. The aims of this review are to describe the principal effects of natural/artificial selection on livestock genomes, to present the main methods used to detect selection signatures and to discuss some recent results in this area. This review should be useful also to research scientists working with wild animals/non-domesticated species and plant biologists working with breeding and evolutionary biology.

Validation of a microsatellite panel for parentage testing of locally adapted and commercial goats in Brazil.

Brazilian goats are generally kept in small herds and extensive rearing systems, mainly in the northeastern region of the country. Despite production improvement in recent years, the lack of pedigree control has affected genetic progress. This study aimed to validate a panel of 16 microsatellites for parentage testing in locally adapted and commercial goats breeds raised in Brazil, as well as to compare its efficiency with the panel recommended by the Brazilian Ministry of Agriculture, Livestock and Supplies (MAPA) in 2004. The number of alleles and expected heterozygosity (He) per marker ranged from four to 18, and from 0.051 to 0.831, respectively. Using all markers, 100% of parentage cases of the validation dataset were resolved with a strict confidence level of 95%. The 16 microsatellites panel showed adequate exclusion power (99.99%) and identity accuracy (99.99%). Suggestions for improvement of the marker panel endorsed by MAPA are provided.

Fine mapping and single nucleotide polymorphism effects estimation on pig chromosomes 1, 4, 7, 8, 17 and X.

Fine mapping of quantitative trait loci (QTL) from previous linkage studies was performed on pig chromosomes 1, 4, 7, 8, 17, and X which were known to harbor QTL. Traits were divided into: growth performance, carcass, internal organs, cut yields, and meat quality. Fifty families were used of a F2 population produced by crossing local Brazilian Piau boars with commercial sows. The linkage map consisted of 237 SNP and 37 microsatellite markers covering 866 centimorgans. QTL were identified by regression interval mapping using GridQTL. Individual marker effects were estimated by Bayesian LASSO regression using R. In total, 32 QTL affecting the evaluated traits were detected along the chromosomes studied. Seven of the QTL were known from previous studies using our F2 population, and 25 novel QTL resulted from the increased marker coverage. Six of the seven QTL that were significant at the 5% genome-wide level had SNPs within their confidence interval whose effects were among the 5% largest effects. The combined use of microsatellites along with SNP markers increased the saturation of the genome map and led to smaller confidence intervals of the QTL. The results showed that the tested models yield similar improvements in QTL mapping accuracy.

Multivariate analysis of heat tolerance characteristics in Santa Inês and crossbred lambs in the Federal District of Brazil.

Heat tolerance is an important trait for deciding on which genetic resource to use on farms in tropical regions. Three genetic groups of lambs were evaluated by physical and physiological parameters related to adaptation and multivariate analyses used to test if these traits were efficient in separating the genetic groups and which traits were more important in this separation. Crossbred lambs were more similar to each other than to Santa Inês, in general; but individually there was an overlap between Ile de France and Santa Inês lambs. Size traits had higher canonical correlation with physiological variables, while rectal temperature, heart rate, and hemoglobin were more highly correlated with physical traits. Error rates using only physiological traits were high, suggesting that these are inefficient in separating groups of lambs for heat tolerance, while separation with physical traits was more efficient.

Genetic parameters for type classification of Nelore cattle on central performance tests at pasture in Brazil.

The objective of this study was to characterize Nelore cattle on central performance tests in pasture, ranked by the visual classification method EPMURAS (structure, precocity, muscle, navel, breed, posture, and sexual characteristics), and to estimate genetic and phenotypic correlations between these parameters, including visual as well as production traits (initial and final weight on test, weight gain, and weight corrected for 550 days). The information used in the study was obtained on 21,032 Nelore bulls which were participants in the central performance test at pasture of the Brazilian Association for Zebu Breeders (ABCZ). Heritabilities obtained were from 0.19 to 0.50. Phenotypic correlations were positive from 0.70 to 0.97 between the weight traits, from 0.65 to 0.74 between visual characteristics, and from 0.29 to 0.47 between visual characteristics and weight traits. The genetic correlations were positive ranging from 0.80 to 0.98 between the characteristics of structure, precocity and musculature, from 0.13 to 0.64 between the growth characteristics, and from 0.41 to 0.97 between visual scores and weight gains. Heritability and genetic correlations indicate that the use of visual scores, along with the selection for growth characteristics, can bring positive results in selection of beef cattle for rearing on pasture.

Landscape Genetics for Brazilian Equines.

Optimization of DNA collection for National gene bank and conservation programs requires information on spatial and genetic distribution of animals countrywide. The relationship between genetic and geographic distances were examined in 8 Brazilian horse breeds (Baixadeiro, Crioulo, Campeiro, Lavradeiro, Marajoara, Mangalarga Marchador, Pantaneiro and Puruca) using Single Nucleotide Polymorphism markers and collection point locations. Mantel correlations, Genetic Landscape Shape Interpolation, Allelic Aggregation Index Analyses and Spatial autocorrelation tests indicated a nonrandom distribution of horses throughout the country. Minimum collection distances for the national Gene Bank should be 530km, with clear divisions seen in genetic structure of horse populations in both North/South and East/West directions. Comparing Pantaneiro and North/Northeastern breeds, physical distance is not necessarily the defining factor for genetic differentiation. This should be considered when sampling these local breeds. These data can help optimise GenBank collection routines and conservation strategies for these breeds.