Use of microsatellite data and pedigree information in the genetic management of two long-term salmon conservation programs

In this paper, we describe the utility of microsatellite data and genetic pedigree information to guide the genetic management of two long-term conservation programs for endangered populations of salmon: Snake River Sockeye Salmon, Oncorhynchus nerka, and inner Bay of Fundy Atlantic Salmon, Salmo salar. Both programs are captive broodstock (live gene banking) programs for endangered populations of salmon. In order for these programs to be successful for recovery efforts, genetic change, including accumulation of inbreeding, loss of genetic variation, and adaptation to captivity, must be minimized. We provide an overview of each program, describe broodstock selection and pairing for spawning, and discuss how pedigree data are being used to evaluate different management practices. While there are inherent species and programmatic differences, both of these programs use widely accepted genetic conservation strategies (minimize mean kinship, reduce variance in family size, minimize inbreeding in the next generation, maintain large census and effective population size) to potentially mitigate some unintended side-effects associated with the rearing of small populations in captivity. These case studies highlight the benefits and practical limitations of applying these strategies in the genetic management of salmon, and may be used to inform other conservation programs.     

QTL linkage analysis of connected populations using ancestral marker and pedigree information

The common assumption in quantitative trait locus (QTL) linkage mapping studies that parents of multiple connected populations are unrelated is unrealistic for many plant breeding programs. We remove this assumption and propose a Bayesian approach that clusters the alleles of the parents of the current mapping populations from locus-specific identity by descent (IBD) matrices that capture ancestral marker and pedigree information. Moreover, we demonstrate how the parental IBD data can be incorporated into a QTL linkage analysis framework by using two approaches: a Threshold IBD model (TIBD) and a Latent Ancestral Allele Model (LAAM). The TIBD and LAAM models are empirically tested via numerical simulation based on the structure of a commercial maize breeding program. The simulations included a pilot dataset with closely linked QTL on a single linkage group and 100 replicated datasets with five linkage groups harboring four unlinked QTL. The simulation results show that including parental IBD data (similarly for TIBD and LAAM) significantly improves the power and particularly accuracy of QTL mapping, e.g., position, effect size and individuals’ genotype probability without significantly increasing computational demand.     

The use of geographical information systems in biodiversity exploration and conservation

We describe a method for applying geographical information systems (GIS) to exploring biodiversity in the wild relatives of crop species and illustrate its application to the wild common bean (Phaseolus vulgaris L.). We use the latitude, longitude and altitude of the location of origin of each accession in a germplasm collection of wild P. vulgaris, along with long-term monthly mean values of rainfall, temperature and diurnal temperature range for about 10000 stations throughout Latin America to produce maps indicating areas with bean-favouring climates. In a test case, these identify a new suitable area in Colombia where wild P. vulgaris has been reported in the literature, and two more areas which are strong candidates on other grounds. Dividing the bean-favouring climates into clusters identifies areas that have similar climates but are geographically remote, where we can expect to find wild beans with similar ecological adaptation. We discuss the implications of these results for conserving and improving the common bean, and the application of these methods to other species.     

Evaluating pedigree data. I. The estimation of pedigree error in the presence of marker mistyping

Pedigrees used in the analysis of genetic or medical data are usually ascertained from sources subject to a variety of errors including misidentification of individuals, faults in historical documents or record linkage, nonpaternity, and unidentified adoption. Genetic markers can be used to verify putative family and pedigree data through the search for inconsistencies, or genetic exclusions, between putative parents and offspring. The probability of observing an exclusion given the occurrence of an error depends upon the gene frequencies at the loci under study and the forms of error. In addition, inconsistencies can arise from laboratory errors in marker determination. Together, these problems make the proper statistical analysis of such data desirable. Here we give a model that specifies the combined effects of various kinds of pedigree error along with genetic marker error. This model allows the maximum-likelihood estimation of the rates of various forms of pedigree error and laboratory error from genetic marker data collected on putative families. The method is illustrated by applying it to data obtained from a South Pacific island population, Tokelau. From the observed distribution of genetic marker inconsistencies between the parents and offspring of putative families, derived from the extensive genealogy of this population, we are able to estimate that the error of a paternal link is 4%, the error of a maternal link is zero, and the overall system typing error is 1%.     

Efficiency of triple test cross for detecting epistasis with marker information

The triple test cross (TTC) is an experimental design for detecting epistasis and estimating the components of genetic variance for quantitative traits. In this paper, we extend the analysis to include molecular information. The statistical power of the mating design was assessed under a model assuming that a finite number of loci affect the trait in question. Formulae are developed for the analysis with or without marker information relating to the recombination fraction between loci, the genetical properties of quantitative trait controlled by the quantitative trait loci (QTL), the linkage phases of the parents and population size. Application of these formulae showed that the recombination fraction between genes and the magnitude and the types of epistasis have important interactions in their effects on power. The results demonstrate that the TTC may have increased power to detect epistasis when marker information is present. However, the simulation experiments show that the standard deviation of the estimated expected mean square was higher with one marker than that with two, whereas the corresponding value without marker information was the lowest. In addition, we demonstrate that the relative position of QTL and markers and the number of markers can both affect the power of epistatic detection.     

Efficiency and robustness of pedigree segregation analysis.

Different pedigree structures and likelihoods are examined to determine their efficiency for parameter estimation under one-locus models. For the cases simulated, family size has little effect; estimates based on unconditional likelihoods are generally more efficient than those based on conditional likelihoods. The proposed method of pedigree analysis under a one-locus model is found to be robust in the analysis of nuclear families: skewness of the data and polygenic inheritance will not lead to the spurious detection of major loci unless they occur simultaneously, and together with a moderate amount of environmental correlation among sibs.     

Relationships among durum wheat accessions. II. A comparison of molecular and pedigree information.

The study of direct ancestry relationships provides information with which to determine essential derivation. SSR profiles were used to determine the pattern of relatedness among 134 durum wheat accessions, representing the most important modern durum wheat gene pools. Simple sequence repeat (SSR)- and amplified fragment length polymorphism (AFLP)-based genetic similarities among cultivars with accurate pedigrees were compared with pedigree-based coefficients of parentage. Sizeable departures of molecular similarities from the expected ones were observed, indicating the unreliability of inferring the pattern of genetic relatedness from the coefficient of parentage. Case studies consisting of parent-progeny cultivar trios and pairs, identified on the basis of their registered pedigree, were studied to evaluate the probability of ancestry of each progeny cultivar, compared with all the remaining accessions. Rare alleles and haplotype sharing were also explored. When the results did not agree with the registered parentages, SSR markers provided information with which to identify the most probable parents (or the corresponding "breeding lineages") in the collection.     

Haplotype inference in crossbred populations without pedigree information

Background

Current methods for haplotype inference without pedigree information assume random mating populations. In animal and plant breeding, however, mating is often not random. A particular form of nonrandom mating occurs when parental individuals of opposite sex originate from distinct populations. In animal breeding this is called crossbreeding and hybridization in plant breeding. In these situations, association between marker and putative gene alleles might differ between the founding populations and origin of alleles should be accounted for in studies which estimate breeding values with marker data. The sequence of alleles from one parent constitutes one haplotype of an individual. Haplotypes thus reveal allele origin in data of crossbred individuals.

Results

We introduce a new method for haplotype inference without pedigree that allows nonrandom mating and that can use genotype data of the parental populations and of a crossbred population. The aim of the method is to estimate line origin of alleles. The method has a Bayesian set up with a Dirichlet Process as prior for the haplotypes in the two parental populations. The basic idea is that only a subset of the complete set of possible haplotypes is present in the population.

Conclusion

Line origin of approximately 95% of the alleles at heterozygous sites was assessed correctly in both simulated and real data. Comparing accuracy of haplotype frequencies inferred with the new algorithm to the accuracy of haplotype frequencies inferred with PHASE, an existing algorithm for haplotype inference, showed that the DP algorithm outperformed PHASE in situations of crossbreeding and that PHASE performed better in situations of random mating.     

A unified approach to adjusting association tests for population admixture with arbitrary pedigree structure and arbitrary missing marker information

A general approach to family-based examinations of association between marker alleles and traits is proposed. The approach is based on computing p values by comparing test statistics for association to their conditional distributions given the minimal sufficient statistic under the null hypothesis for the genetic model, sampling plan and population admixture. The approach can be applied with any test statistic, so any kind of phenotype and multi-allelic markers may be examined, and covariates may be included in analyses. By virtue of the conditioning, the approach results in correct type I error probabilities regardless of population admixture, the true genetic model and the sampling strategy. An algorithm for computing the conditional distributions is described, and the results of the algorithm for configurations of nuclear families are presented. The algorithm is applicable with all pedigree structures and all patterns of missing marker allele information.     

NABIC marker database: A molecular markers information network of agricultural crops

In 2013, National Agricultural Biotechnology Information Center (NABIC) reconstructs a molecular marker database for useful genetic resources. The web-based marker database consists of three major functional categories: map viewer, RSN marker and gene annotation. It provides 7250 marker locations, 3301 RSN marker property, 3280 molecular marker annotation information in agricultural plants. The individual molecular marker provides information such as marker name, expressed sequence tag number, gene definition and general marker information. This updated marker-based database provides useful information through a user-friendly web interface that assisted in tracing any new structures of the chromosomes and gene positional functions using specific molecular markers.

Availability

The database is available for free at http://nabic.rda.go.kr/gere/rice/molecularMarkers/     

Identification and potential use of a molecular marker for rust resistance in common bean

Summary The Up ₂ gene of common bean (Phaseolus Vulgaris L.) is an important source of dominant genetic resistance to the bean rust pathogen [Uromyces appendiculatus (Pers. ex Pers.) Unger var appendiculatus [syn U. Phaseoli (Reben) Wint.]. Up ₂ in combination with other rust resistance genes may be used to obtain potentially stable genetic resistance. It is difficult, however, to combine rust resistance genes effective against a single race due to epistatic interactions that frequently occur between them. A strategy that employed bulked DNA samples formed separately from the DNA of three BC₆F₂ individuals with Up ₂ and three without Up ₂ as contrasting near-isogenic lines (NILs) was used to identify random amplified polymorphic DNA fragments (RAPDs) tightly linked to the Up ₂ locus. Only 1 of 931 fragments amplified by 167 10-mer primers of arbitrary sequence in the polymerase chain reaction (PCR) was polymorphic. The RAPD marker (OA14₁₁₀₀) amplified by the 5-TCTGTGCTGG-3 primer was repeatable and its presence and absence easy to score. No recombination was observed between OA14₁₁₀₀ and the dominant Up ₂ allele within a segregating BC₆F₂ population of 84 individuals. This result suggests that OA14₁₁₀₀ and Up ₂ are tightly linked. Andean and Mesoamerican bean germ plasm, with and without the Up ₂ allele, were assayed for the presence of OA14₁₁₀₀. Apparently, the marker is of Andean origin because all Andean lines, with or without the Up ₂ allele, contained the marker, and the marker was absent in all Mesoamerican germ plasm except the lines to which Up-2 had been purposely transferred. These results suggest that OA14₁₁₀₀ will be most useful for pyramiding Up ₂ with other rust resistance genes into germ plasm of Mesoamerican origin where the marker does not traditionally exist. The use of bulked DNA samples may have concentrated resources toward the identification of RAPDs that were tightly linked to the target locus. Marker-based selection may provide an alternative to the time-consuming testcrosses required to pyramid bean rust resistance genes that exhibit epistasis.Research supported by the Michigan Agricultural Research Station and the USDA-ARS. Mention of a trademark or proprietary product does not constitute a guarantee or warranty of the product by the USDA and does not imply its approval to the exclusion of other products that may also be suitable     

Determining informativity of marker typing for genetic counseling in a pedigree

Summary For the situation of a Mendelian disease linked to a genetic marker, a new method is described that allows evaluating for genetic counseling the information potentially available from the linked marker before the marker data are actually obtained, that is, prior to drawing blood for marker typing. For a consultand in a family pedigree, the method determines the risk distribution (small families) or an approximation to it (larger families) and calculates the probability that the risk will deviate beyond certain limits from its a priori value, which exists without marker data, for example, that the risk will be smaller than 0.10 or larger than 0.90. The method was applied here to a pedigree of 15 individuals for which analytical calculations would be difficult to carry out.     

Estimating relatedness between individuals in general populations with a focus on their use in conservation programs

Relatedness estimators are widely used in genetic studies, but effects of population structure on performance of estimators, criteria to evaluate estimators, and benefits of using such estimators in conservation programs have to date received little attention. In this article we present new estimators, based on the relationship between coancestry and molecular similarity between individuals, and compare them with existing estimators using Monte Carlo simulation of populations, either panmictic or structured. Estimators were evaluated using statistical criteria and a diversity criterion that minimized relatedness. Results show that ranking of estimators depends on the population structure. An existing estimator based on two-gene and four-gene coefficients of identity performs best in panmictic populations, whereas a new estimator based on coancestry performs best in structured populations. The number of marker alleles and loci did not affect ranking of estimators. Statistical criteria were insufficient to evaluate estimators for their use in conservation programs. The regression coefficient of pedigree relatedness on estimated relatedness (beta2) was substantially lower than unity for all estimators, causing overestimation of the diversity conserved. A simple correction to achieve beta2 = 1 improves both existing and new estimators. Using relatedness estimates with correction considerably increased diversity in structured populations, but did not do so or even decreased diversity in panmictic populations.     

Genetic parameters for growth and faecal worm egg count following Haemonchus contortus experimental infestations using pedigree and molecular information

Background

Haemonchosis is a parasitic disease that causes severe economic losses in sheep industry. In recent years, the increasing resistance of the parasite to anthelmintics has raised the need for alternative control strategies. Genetic selection is a promising alternative but its efficacy depends on the availability of genetic variation and on the occurrence of favourable genetic correlations between the traits included in the breeding goal. The objective of this study was twofold. First, to estimate both the heritability of and the genetic correlations between growth traits and parasite resistance traits, using bivariate linear mixed animal models, from the phenotypes and genotypes of 1004 backcross lambs (considered as a single population), which underwent two subsequent experimental infestations protocols with Haemonchus contortus. Second, to compare the precision of the estimates when using two different relationship matrices: including pedigree information only or including also SNP (single nucleotide polymorphism) information.

Results

Heritabilities were low for average daily gain before infestation (0.10 to 0.15) and average daily gain during the first infestation (0.11 to 0.16), moderate for faecal egg counts during the first infestation (0.21 to 0.38) and faecal egg counts during the second infestation (0.48 to 0.55). Genetic correlations between both growth traits and faecal egg count during the naïve infestation were equal to zero but the genetic correlation between faecal egg count during the second infestation and growth was positive in a Haemonchus contortus free environment and negative in a contaminated environment. The standard errors of the estimates obtained by including SNP information were smaller than those obtained by including pedigree information only.

Conclusions

The genetic parameters estimates suggest that growth performance can be selected for independently of selection on resistance to naïve infestation. Selection for increased growth in a non-contaminated environment could lead to more susceptible animals with long-term exposure to the infestation but it could be possible to select for increased growth in a contaminated environment while also increasing resistance to the long-term exposure to the parasite. The use of molecular information increases the precision of the estimates.     

Breast cancer risk assessment: use of complete pedigree information and the effect of misspecified ages at diagnosis of affected relatives

Received: 10 July 1997 / Accepted: 16 November 1997     

Genetic adaptation to captivity in species conservation programs

As wild environments are often inhospitable, many species have to be captive-bred to save them from extinction. In captivity, species adapt genetically to the captive environment and these genetic adaptations are overwhelmingly deleterious when populations are returned to wild environments. I review empirical evidence on (i) the genetic basis of adaptive changes in captivity, (ii) factors affecting the extent of genetic adaptation to captivity, and (iii) means for minimizing its deleterious impacts. Genetic adaptation to captivity is primarily due to rare alleles that in the wild were deleterious and partially recessive. The extent of adaptation to captivity depends upon selection intensity, genetic diversity, effective population size and number of generation in captivity, as predicted by quantitative genetic theory. Minimizing generations in captivity provides a highly effective means for minimizing genetic adaptation to captivity, but is not a practical option for most animal species. Population fragmentation and crossing replicate captive populations provide practical means for minimizing the deleterious effects of genetic adaptation to captivity upon populations reintroduced into the wild. Surprisingly, equalization of family sizes reduces the rate of genetic adaptation, but not the deleterious impacts upon reintroduced populations. Genetic adaptation to captivity is expected to have major effects on reintroduction success for species that have spent many generations in captivity. This issue deserves a much higher priority than it is currently receiving.     

The role of pedigree information in combined linkage disequilibrium and linkage mapping of quantitative trait loci in a general complex pedigree

Combined linkage disequilibrium and linkage (LDL) mapping can exploit historical as well as recent and observed recombinations in a recorded pedigree. We investigated the role of pedigree information in LDL mapping and the performance of LDL mapping in general complex pedigrees. We compared using complete and incomplete genotypic data, spanning 5 or 10 generations of known pedigree, and we used bi- or multiallelic markers that were positioned at 1- or 5-cM intervals. Analyses carried out with or without pedigree information were compared. Results were compared with linkage mapping in some of the data sets. Linkage mapping or LDL mapping with sparse marker spacing ( approximately 5 cM) gave a poorer mapping resolution without considering pedigree information compared to that with considering pedigree information. The difference was bigger in a pedigree of more generations. However, LDL mapping with closely linked markers ( approximately 1 cM) gave a much higher mapping resolution regardless of using pedigree information. This study shows that when marker spacing is dense and there is considerable linkage disequilibrium generated from historical recombinations between flanking markers and QTL, the loss of power due to ignoring pedigree information is negligible and mapping resolution is very high.