Genomic Analysis of Inbreeding and Coancestry in Nordic Jersey and Holstein Dairy Cattle Populations.

In recent years, Genomic Selection (GS) has accelerated genetic gain in dairy cattle breeds worldwide. Despite the evident genetic progress, several dairy populations have also encountered challenges such as heightened inbreeding rates and reduced effective population sizes. The challenge has been to find a balance between achieving substantial genetic gain while managing genetic diversity within the population, thereby mitigating the negative effects of inbreeding depression. This study aims to elucidate the impact of GS on pedigree and genomic rates of inbreeding (ΔF) and coancestry (ΔC) in Nordic Jersey (NJ) and Holstein (NH) cattle populations. Furthermore, key genetic metrics including the generation interval (L), effective population size (Ne), and future effective population size (FNe) were assessed between 2 time periods, before and after GS, and across distinct animal cohorts in both breeds: females, bulls, and approved semen-producing bulls (AI-sires). Analysis of ΔF and ΔC revealed distinct trends across the studied periods and animal groups. Notably, there was a consistent increase in yearly ΔF for most animal groups in both breeds. An exception was observed in NH AI-sires, which demonstrated a slight decrease in yearly ΔF. Moreover, NJ displayed minimal changes in yearly ΔC between the periods, whereas NH exhibited elevated ΔC values across all animal groups. Particularly striking was the substantial increase in yearly ΔC within the NH female population, surging from 0.02% to 0.39% between the periods. Implementation of GS resulted in a reduction of the generation interval across all animal cohorts in both NJ and NH breeds. However, the extent of reduction was more pronounced in males compared with females. This reduction in generation interval influenced generational changes in ΔF and ΔC. Bulls and AI-sires of both breeds exhibited reduced generational ΔF between periods, in contrast to females that demonstrated an opposing pattern. Between the periods, NJ maintained a relatively stable Ne, 29.4 before and 30.3 after GS, while NH experienced a notable decline from 54.3 to 42.8. Female groups in both breeds displayed a negative Ne trend, while males demonstrated either neutral or positive Ne developments. Regarding FNe, NJ exhibited positive FNe development with an increase from 40.7 to 57.2. The opposite was observed in NH, where FNe decreased from 198.8 to 42.7. In summary, it was evident that the genomic methods could detect differences between the populations and changes in ΔF and ΔC more efficiently than pedigree methods. GS implementation yielded positive outcomes within the NJ population regarding the rate of coancestry but the opposite was observed with NH. Moreover, analysis of ΔC data hints at the potential to decrease future ΔF through informed mating strategies. Conversely, NH faces more pressing concerns, even though ΔF remains comparatively modest in contrast to what has been observed in other Holstein populations. These findings underscore the necessity of genomic control of inbreeding and coancestry with strategic changes in the Nordic breeding schemes for dairy to ensure long-term sustainability in the forthcoming years.

Fine-scale genetic structure in the orchid Gymnadenia conopsea is not associated with local density of flowering plants.

PREMISE: Density-dependent pollinator visitation can lead to density-dependent mating patterns and within-population genetic structure. In Gymnadenia conopsea, individuals in low-density patches receive more self pollen than individuals in high-density patches, suggesting higher relatedness at low density. Ongoing fragmentation is also expected to cause more local matings, potentially leading to biparental inbreeding depression. METHODS: To evaluate whether relatedness decreases with local density, we analyzed 1315 SNP loci in 113 individuals within two large populations. We quantified within-population genetic structure in one of the populations, recorded potential habitat barriers, and visualized gene flow using estimated effective migration surfaces (EEMS). We further estimated the magnitude of biparental inbreeding depression that would result from matings restricted to within 5 m. RESULTS: There was no significant relationship between local density and relatedness in any population. We detected significant fine-scale genetic structure consistent with isolation by distance, with positive kinship coefficients at distances below 10 m. Kinship coefficients were low, and predicted biparental inbreeding depression resulting from matings within the closest 5 m was a modest 1-3%. The EEMS suggested that rocks and bushes may act as barriers to gene flow within a population. CONCLUSIONS: The results suggest that increased self-pollen deposition in sparse patches does not necessarily cause higher selfing rates or that inbreeding depression results in low establishment success of inbred individuals. The modest relatedness suggests that biparental inbreeding depression is unlikely to be an immediate problem following fragmentation of large populations. The results further indicate that habitat structure may contribute to governing fine-scale genetic structure in G. conopsea.

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.

Process of Introduction of Australian Braford Cattle to South America: Configuration of Population Structure and Genetic Diversity Evolution.

This study analyzes the evolution of the population structure and genetic diversity of Braford cattle in South America from 1949 to 2019 to suggest effective strategies for breeding in the future. The percentage of bulls historically increased. The average generational interval decreased to 11.78 years for the current population. Average inbreeding (F) and coancestry (C) are low and show a historically increasing trend (0.001% to 0.002%, respectively). The degree of nonrandom mating (α) increased from -0.0001 to 0.0001 denoting a change in the trend to mate similar individuals. The average relatedness coefficient (ΔR) increased in the current period from 0.002% to 0.004%. A single ancestor explained 4.55% to 7.22% of the population's gene pool. While the effective population size based on the individual inbreeding rate (NeFi) was 462.963, when based on the individual coancestry rate (NeCi), it was 420.168. Genetic diversity loss is small and mainly ascribed to bottlenecks (0.12%) and to unequal contributions of the founders (0.02%). Even if adequate levels of diversity can be found, practices that consider the overuse of individual bulls (conditioned by nature or not), could lead to a long-term reduction in diversity. The present results permit tailoring genetic management strategies that are perfectly adapted to the needs that the population demands internationally.

Evolution of inbreeding: a gaze into five Italian beef cattle breeds history.

In the last decades, intensive selection programs have led to sustained increases of inbreeding in dairy cattle, a feature that might have adverse consequences on the viability and phenotypic performance of their offspring. This study aimed to determine the evolution of inbreeding of five Italian beef cattle breeds (Marchigiana, Chianina, Romagnola, Maremmana, and Podolica) during a period of almost 20 years (2002-2019). The estimates of Ho, He, Fhat2 , and Fped averaged across years (2002-2019) in the studied breeds fluctuated between 0.340-0.401, 0.348-0.392, -0.121-0.072, and 0.000-0.068, respectively. Moreover, annual rates of increase of the estimated inbreeding coefficients have been very low (Fhat2 = 0.01-0.02%; Fped = 0.003-0.004%). The use of a high number of bulls combined with strategies implemented by the Association of Italian Beef Cattle Breeders ANABIC to minimize inbreeding might explain these results. Despite the fact that diversity and inbreeding have remained quite stable during the last two decades, we have detected a sustained decrease of the population effective size of these five breeds. Such results should be interpreted with caution due to the inherent difficulty of estimating Ne from SNPs data in a reliable manner.

Genomic inbreeding and population structure of northern pike (Esox lucius) in Xinjiang, China.

Northern pike (Esox lucius) was widely distributed in the high latitudes of the northern hemisphere. In China, northern pike was originally distributed only in the upper reaches of the Irtysh River in Xinjiang and has appeared in many water bodies outside the Irtysh River Basin in Northern Xinjiang. A total of four populations were collected from north to south in Xinjiang, including Irtysh River (RIR), Ulungu Lake (LUL), a small lake nearby Ulungu River (LJD), and Bosten Lake (LBO). We estimated population genomic parameters, performed gene flow analysis, and estimated the effective population size of each population. The proportion of individuals with high inbreeding coefficient (F ≥ 0.0625) accounted for 36.4% (44/121) of all sequenced individuals, approximately 4.5% (1/22) in LUL, 25.9% (7/27) in LBO, 42.9% (18/42) in RIR, and 60% (18/30) in LJD. RIR had the highest mean of genomic relatedness (coancestry coefficient = 0.025 ± 0.040, IBD = 0.036 ± 0.078). Gene flow results showed that the population spreading was from RIR into two branches, one was LBO, and the other continued to split into LUL and LJD, and migration signal from LBO to LUL was detected. Our results suggested that the extinction risk of northern pike was very low in Xinjiang of China, and the controlled capture fishery of northern pike could be developed reasonably.

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.

Genomic phylogeography of the White-crowned Manakin Pseudopipra pipra (Aves: Pipridae) illuminates a continental-scale radiation out of the Andes.

The complex landscape history of the Neotropics has generated opportunities for population isolation and diversification that place this region among the most species-rich in the world. Detailed phylogeographic studies are required to uncover the biogeographic histories of Neotropical taxa, to identify evolutionary correlates of diversity, and to reveal patterns of genetic connectivity, disjunction, and potential differentiation among lineages from different areas of endemism. The White-crowned Manakin (Pseudopipra pipra) is a small suboscine passerine bird that is broadly distributed through the subtropical rainforests of Central America, the lower montane cloud forests of the Andes from Colombia to central Peru, the lowlands of Amazonia and the Guianas, and the Atlantic forest of southeast Brazil. Pseudopipra is currently recognized as a single, polytypic biological species. We studied the effect of the Neotropical landscape on genetic and phenotypic differentiation within this species using genomic data derived from double digest restriction site associated DNA sequencing (ddRAD), and mitochondrial DNA. Most of the genetic breakpoints we identify among populations coincide with physical barriers to gene flow previously associated with avian areas of endemism. The phylogenetic relationships among these populations imply a novel pattern of Andean origination for this group, with subsequent diversification into the Amazonian lowlands. Our analysis of genomic admixture and gene flow reveals a complex history of introgression between some western Amazonian populations. These reticulate processes confound our application of standard concatenated and coalescent phylogenetic methods and raise the question of whether a lineage in the western Napo area of endemism should be considered a hybrid species. Lastly, analysis of variation in vocal and plumage phenotypes in the context of our phylogeny supports the hypothesis that Pseudopipra is a species-complex composed of at least 8, and perhaps up to 17 distinct species which have arisen in the last ~2.5 Ma.

The ghosts of propagation past: haplotype information clarifies the relative influence of stocking history and phylogeographic processes on contemporary population structure of walleye (Sander vitreus).

Stocking of fish is an important tool for maintaining fisheries but can also significantly alter population genetic structure and erode the portfolio of within-species diversity that is important for promoting resilience and adaptability. Walleye (Sander vitreus) are a highly valued sportfish in the midwestern United States, a region characterized by postglacial recolonization from multiple lineages and an extensive history of stocking. We leveraged genomic data and recently developed analytical approaches to explore the population structure of walleye from two midwestern states, Minnesota and Wisconsin. We genotyped 954 walleye from 23 populations at ~20,000 loci using genotyping by sequencing and tested for patterns of population structure with single-SNP and microhaplotype data. Populations from Minnesota and Wisconsin were highly differentiated from each other, with additional substructure found in each state. Population structure did not consistently adhere to drainage boundaries, as cases of high intra-drainage and low inter-drainage differentiation were observed. Low genetic structure was observed between populations from the upper Wisconsin and upper Chippewa river watersheds, which are found as few as 50 km apart and were likely homogenized through historical stocking. Nevertheless, we were able to differentiate these populations using microhaplotype-based co-ancestry analysis, providing increased resolution over previous microsatellite studies and our other single SNP-based analyses. Although our results illustrate that walleye population structure has been influenced by past stocking practices, native ancestry still exists in most populations and walleye populations may be able to purge non-native alleles and haplotypes in the absence of stocking. Our study is one of the first to use genomic tools to investigate the influence of stocking on population structure in a nonsalmonid fish and outlines a workflow leveraging recently developed analytical methods to improve resolution of complex population structure that will be highly applicable in many species and systems.

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.

Casework applications of probabilistic genotyping methods for DNA mixtures that allow relationships between contributors.

In both criminal cases and civil cases there is an increasing demand for the analysis of DNA mixtures involving relationships. The goal might be, for example, to identify the contributors to a DNA mixture where the unknown donors may be related, or to infer the relationship between individuals based on a DNA mixture. This paper applies a new approach to modelling and computation for DNA mixtures involving contributors with arbitrarily complex relationships to two real cases from the Spanish Forensic Police.

Improved dairy cattle mating plans at herd level using genomic information.

From 2012 to 2018, 223 180 Montbéliarde females were genotyped in France and the number of newly genotyped females increased at a rate of about 33% each year. With female genotyping information, farmers have access to the genomic estimated breeding values of the females in their herd and to their carrier status for genetic defects or major genes segregating in the breed. This information, combined with genomic coancestry, can be used when planning matings in order to maximize the expected on-farm profit of future female offspring. We compared different mating allocation approaches for their capacity to maximize the expected genetic gain while limiting expected progeny inbreeding and the probability to conceive an offspring homozygous for a lethal recessive allele. Three mate allocation strategies (random mating (RAND), sequential mating (gSEQ€) and linear programing mating (gLP€)) were compared on 160 actual Montbéliarde herds using male and female genomic information. Then, we assessed the benefit of using female genomic information by comparing matings planned using only female pedigree information with the equivalent strategy using genomic information. We measured the benefit of adding genomic expected inbreeding and risk of conception of an offspring homozygous for a lethal recessive allele to Net merit in mating plans. The influence of three constraints was tested: by relaxing the constraint on availability of a particular semen type (sexed or conventional) for bulls, by adding an upper limit of 8.5% coancestry between mate pairs or by using a more stringent maximum use of a bull in a herd (5% vs 10%). The use of genomic information instead of pedigree information improved the mate allocation method in terms of progeny expected genetic merit, genetic diversity and risk to conceive an offspring homozygous for a lethal recessive allele. Optimizing mate allocation using linear programming and constraining coancestry to a maximum of 8.5% per mate pair reduced the average coancestry with a small impact on expected Net Merit. In summary, for male and female selection pathways, using genomic information is more efficient than using pedigree information to maximize genetic gain while constraining the expected inbreeding of the progeny and the risk to conceive an offspring homozygous for a lethal recessive allele. This study also underlines the key role of semen type (sexed vs conventional) and the associated constraints on the mate allocation algorithm to maximize genetic gain while maintaining genetic diversity and limiting the risk to conceive an offspring homozygous for a lethal recessive allele.

A unified linear mixed model for familial relatedness and population structure in genetic association studies.

Familial relatedness (FR) and population structure (PS) are two major sources for genetic correlation. In the human population, both FR and PS can further break down into additive and dominant components to account for potential additive and dominant genetic effects. In this study, besides the classical additive genomic relationship matrix, a dominant genomic relationship matrix is introduced. A link between the additive/dominant genomic relationship matrices and the coancestry (or kinship)/double coancestry coefficients is also established. In addition, a way to separate the FR and PS correlations based on the estimates of coancestry and double coancestry coefficients from the genomic relationship matrices is proposed. A unified linear mixed model is also developed, which can account for both the additive and dominance effects of FR and PS correlations as well as their possible random interactions. Finally, this unified linear mixed model is applied to analyze two study cohorts from UK Biobank.

Inbreeding, effective population size, and coancestry in the Latxa dairy sheep breed.

Traditionally, breeding programs have estimated and managed inbreeding based on pedigree information. The availability of genomic marker panels has made possible new alternatives to achieve more precise estimates, for example in case of missing pedigree. The objective of the present study was to assess and compare, different estimation methods (pedigree-based methodologies, single SNP-based approach (homozygosity) and runs of homozygosity-based method) to analyze the evolution of genetic diversity measured as inbreeding or as coancestry of 3 selected populations of Latxa dairy sheep (Latxa Cara Rubia and Latxa Cara Negra from Euskadi and Navarre). Genomic data came from 972 artificial insemination rams genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA) whose genealogy consisted of 4,484 animals. Inbreeding estimates based on molecular data were more similar between them than compared with those based on pedigree information. However, the SNP-based approach estimations of effective population size differed more, reflecting the sensitivity of effective population size to small changes in the evolution of inbreeding. The 2 Latxa Cara Negra populations showed increases of inbreeding rates with time and effective population sizes between 64 and 103 animals, depending on breed and methodology used. The Latxa Cara Rubia population did not show an increase in inbreeding rate, mainly due to semen importation from the related French population of Manech Tête Rousse. The effective size estimates based on coancestry increase show a higher variability and they are more sensitive to the source of information and the data structure considered. Realized effective population size based on individual increase in inbreeding were in agreement with the previous estimates. Coancestry evolution analysis based on DNA information showed an increase on coancestry during the last 10 yr in all breeds, as a consequence of the selection process. Moreover, the increase on coancestry between Latxa Cara Rubia and Manech Tête Rousse was more noticeable between than within each of those breeds.

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.

Paternity analysis, pollen flow, and spatial genetic structure of a natural population of Euterpe precatoria in the Brazilian Amazon.

Euterpe precatoria, known as açaí do Amazonas, is a regionally important palm of the Amazon rainforest for the fruit production through extractive agriculture. Little information is available with regard to genetic diversity, gene flow, and spatial genetic structure (SGS) of açaí populations, which are essential for the use, management, and conservation of genetic resources of the species. This research aimed to assess the genetic diversity, inbreeding level, SGS, and gene flow in four ontogenetic stages of a natural E. precatoria population in the Brazilian Amazon, based on 18 microsatellite loci. The study was carried out in a natural population dispersed in an area of about 10 ha. Leaf tissues of 248 plants were mapped and sampled and classified into four ontogenetic stages: reproductive (59), immature (70), young (60), and seedling (59). Genetic diversity indices were high for all ontogenetic stages. The fixation index (F) for all ontogenetic stages was not significantly different from zero, indicating the absence of inbreeding. A significant SGS was found for all ontogenetic stages (68-110 m), indicating seed dispersal over short distances. Paternity analysis detected pollen immigration of 39.1%, a selfing rate of 4.2%, and a mean pollen dispersal distance within the population of 531 m. The results indicate substantial allele input in the population via pollen immigration, contributing to the maintenance of the genetic diversity of the population. However, within a population, the renewal with new progenies selected from seed plants spaced at least 110 m apart is important to avoid collecting seeds from related plants.

Estimating Genetic Relatedness in Admixed Populations.

Estimating genetic relatedness, and inbreeding coefficients is important to the fields of quantitative genetics, conservation, genome-wide association studies (GWAS), and population genetics. Traditional estimators of genetic relatedness assume an underlying model of population structure. Each individual is assigned to a population, depending on a priori assumptions about geographical location of sampling, proximity, or genetic similarity. But often, this population assignment is unknown and assumptions about assignment can lead to erroneous estimates of genetic relatedness. I develop a generalized method of estimating relatedness in admixed populations, to account for (1) multi-allelic genomic data, (2) including all nine Identity By Descent (IBD) states, and implement a maximum likelihood based estimator of pairwise genetic relatedness in structured populations, part of the software, InRelate. Replicated estimations of genetic relatedness between admixed full sib (FS), half sib (HS), first cousin (FC), parent-offspring (PO) and unrelated (UR) dyads in simulated and empirical data from the HGDP-CEPH panel show considerably low bias and error while using InRelate, compared to several previously developed methods. I also propose a bootstrap scheme, and a series of Wald Tests to assign relatedness categories to pairs of individuals.

Urban rat races: spatial population genomics of brown rats (Rattus norvegicus) compared across multiple cities.

Urbanization often substantially influences animal movement and gene flow. However, few studies to date have examined gene flow of the same species across multiple cities. In this study, we examine brown rats (Rattus norvegicus) to test hypotheses about the repeatability of neutral evolution across four cities: Salvador, Brazil; New Orleans, USA; Vancouver, Canada; and New York City, USA. At least 150 rats were sampled from each city and genotyped for a minimum of 15 000 genome-wide single nucleotide polymorphisms. Levels of genome-wide diversity were similar across cities, but varied across neighbourhoods within cities. All four populations exhibited high spatial autocorrelation at the shortest distance classes (less than 500 m) owing to limited dispersal. Coancestry and evolutionary clustering analyses identified genetic discontinuities within each city that coincided with a resource desert in New York City, major waterways in New Orleans, and roads in Salvador and Vancouver. Such replicated studies are crucial to assessing the generality of predictions from urban evolution, and have practical applications for pest management and public health. Future studies should include a range of global cities in different biomes, incorporate multiple species, and examine the impact of specific characteristics of the built environment and human socioeconomics on gene flow.

Pedigree analysis and inbreeding effects over morphological traits in Campolina horse population.

Genetic improvement, without control of inbreeding, can go to loss of genetic variability, reducing the potential for genetic gains in the domestic populations. The aim of this study was to analyze the population structure and the inbreeding depression in Campolina horses. Phenotype information from 43 465 individuals was analyzed, data provided by the Campolina Breeders Association. A pedigree file containing 107 951 horses was used to connected the phenotyped individuals. The inbreeding coefficient was performed by use of the diagonal of the relationship matrix and the genealogical parameters were computed using proper softwares. The effective population size was estimated based on the rate of inbreeding and census information, and the stratification of the population was verified by the average relationship coefficient between animals born in different regions of Brazil. The effects of inbreeding on morphological traits were made by the use of inbreeding coefficient as a covariate in the model of random regression. The inbreeding coefficient increased from 1990 on, impacting effective population size and, consequently, shrinking genetic variability. The paternal inbreeding was greater than maternal, which may be attributed to the preference for inbred animals in reproduction. The average genetic relationship coefficient of animals born in different states was lower than individuals born within the same state. The increase in the inbreeding coefficient was negatively associated with all studied traits, showing the importance to avoid genetic losses in the long term. Although results do not indicate a severe narrowing of the population until the present date, the average relationship coefficient shows signs of increase, which could cause a drastic reduction in genetic variability if inbred mating is not successfully controlled in the Campolina horse population.

Evolving gene banks: improving diverse populations of crop and exotic germplasm with optimal contribution selection.

We simulated pre-breeding in evolving gene banks - populations of exotic and crop types undergoing optimal contribution selection for long-term genetic gain and management of population genetic diversity. The founder population was based on crosses between elite crop varieties and exotic lines of field pea (Pisum sativum) from the primary genepool, and was subjected to 30 cycles of recurrent selection for an economic index composed of four traits with low heritability: black spot resistance, flowering time and stem strength (measured on single plants), and grain yield (measured on whole plots). We compared a small population with low selection pressure, a large population with high selection pressure, and a large population with moderate selection pressure. Single seed descent was compared with S0-derived recurrent selection. Optimal contribution selection achieved higher index and lower population coancestry than truncation selection, which reached a plateau in index improvement after 40 years in the large population with high selection pressure. With optimal contribution selection, index doubled in 38 years in the small population with low selection pressure and 27-28 years in the large population with moderate selection pressure. Single seed descent increased the rate of improvement in index per cycle but also increased cycle time.