Evaluation of Detection Efficiency for Trio Full Sibling Testing.

OBJECTIVES: To study the detection efficiency of trio full sibling with another known full sibling reference added under different number of autosomal STR typing systems. METHODS: Based on 43 detection systems consisting of 13 to 55 representative autosomal STR loci, 10 000 true families (full sibling group) and 10 000 false families (unrelated individual group) were randomly simulated. The full sibling index (FSI) was calculated based on the method of family reconstruction. The cumulative sibling relationship index (CFSI) of 0.000 1 and 10 000 were used as the evaluation thresholds, and the detection efficiency parameters were calculated and compared with the identification of the duo full sibling testing. RESULTS: With the increasing number of STR loci, the error rate and inability of judgement rate gradually decreased; the sensitivity, specificity, correct rate of judgment and other parameters gradually increased, and the system efficiency gradually improved. Under the same detection system, trio full sibling testing showed higher sensitivity, specificity, system efficiency and lower inability of judgement rate compared with duo full sibling testing. When the system efficiency was higher than 0.85 and inability of judgement rate was less than 0.01%, at least 20 STRs should be detected for trio full sibling testing, which was less than 29 STRs required by duo full sibling testing. CONCLUSIONS: The detection efficiency of trio full sibling testing is superior to that of duo full sibling testing with the same detection system, which is an effective identification scheme for laboratories with inadequate detection systems or for materials with limited conditions.

Calculation of Likelihood Ratios for Incest Cases Using IBD Patterns.

OBJECTIVES: To calculate the likelihood ratios of incest cases using identity by descent (IBD) patterns. METHODS: The unique IBD pattern was formed by denoting the alleles from the members in a pedigree with a same digital. The probability of each IBD pattern was obtained by multiplying the prior probability by the frequency of non-IBD alleles. The pedigree likelihoods of incest cases under different hypotheses were obtained by summing all IBD pattern probabilities, and the likelihood ratio(LR) was calculated by comparing the likelihoods of different pedigrees. RESULTS: The IBD patterns and the formulae of calculating LR for father-daughter incest and brother-sister incest were obtained. CONCLUSIONS: The calculations of LR for incest cases were illustrated based on IBD patterns.

Application of Microhaplotypes in Sibling Kinship Testing.

OBJECTIVES: To investigate the efficacy of different numbers of microhaplotype (MH) loci and the introduction of different reference samples on the identification of full sibling, half sibling and differentiation between full sibling and half sibling kinships, and to explore the effect of changing mutation rate on sibling testing. METHODS: First, a family map involving three generations was established, and four full sibling identification models, five half sibling identification models and five models distinguishing full and half siblings were constructed for different reference samples introduced. Based on the results of the previous study, two sets of nonbinary SNP-MH containing 34 and 54 loci were selected. Based on the above MH loci, 100 000 pairs of full sibling vs. unrelated individuals, 100 000 pairs of half sibling vs. unrelated individuals and 100 000 pairs of full sibling vs. half sibling were simulated based on the corresponding sibling kinship testing models, and the efficacy of each sibling kinship testing model was analyzed by the likelihood ratio algorithm under different thresholds. The mutant rate of 54 MH loci was changed to analyze the effect of mutation rate on sibling identification. RESULTS: In the same relationship testing model, the systematic efficacy of sibling testing was positively correlated with the number of MH loci detected. With the same number of MH loci, the efficacy of full sibling testing was better than that of uncle or grandfather when the reference sample introduced was a full sibling of A, but there was no significant difference in the identification efficacy of the four reference samples introduced for full sibling and half sibling differentiation testing. In addition, the mutation rate had a slight effect on the efficacy of sibling kinship testing. CONCLUSIONS: Increasing the number of MH loci and introducing reference samples of known relatives can increase the efficacy of full sibling testing, half sibling testing, and differentiation between full and half sibling kinships. The level of mutation rate in sibling testing by likelihood ratio method has a slight but insignificant effect on the efficacy.

Application of Microhaplotypes in Sibling Kinship Testing / 法医学杂志

OBJECTIVES@#To investigate the efficacy of different numbers of microhaplotype (MH) loci and the introduction of different reference samples on the identification of full sibling, half sibling and differentiation between full sibling and half sibling kinships, and to explore the effect of changing mutation rate on sibling testing.@*METHODS@#First, a family map involving three generations was established, and four full sibling identification models, five half sibling identification models and five models distinguishing full and half siblings were constructed for different reference samples introduced. Based on the results of the previous study, two sets of nonbinary SNP-MH containing 34 and 54 loci were selected. Based on the above MH loci, 100 000 pairs of full sibling vs. unrelated individuals, 100 000 pairs of half sibling vs. unrelated individuals and 100 000 pairs of full sibling vs. half sibling were simulated based on the corresponding sibling kinship testing models, and the efficacy of each sibling kinship testing model was analyzed by the likelihood ratio algorithm under different thresholds. The mutant rate of 54 MH loci was changed to analyze the effect of mutation rate on sibling identification.@*RESULTS@#In the same relationship testing model, the systematic efficacy of sibling testing was positively correlated with the number of MH loci detected. With the same number of MH loci, the efficacy of full sibling testing was better than that of uncle or grandfather when the reference sample introduced was a full sibling of A, but there was no significant difference in the identification efficacy of the four reference samples introduced for full sibling and half sibling differentiation testing. In addition, the mutation rate had a slight effect on the efficacy of sibling kinship testing.@*CONCLUSIONS@#Increasing the number of MH loci and introducing reference samples of known relatives can increase the efficacy of full sibling testing, half sibling testing, and differentiation between full and half sibling kinships. The level of mutation rate in sibling testing by likelihood ratio method has a slight but insignificant effect on the efficacy.

Calculation of Likelihood Ratios for Incest Cases Using IBD Patterns / 法医学杂志

OBJECTIVES@#To calculate the likelihood ratios of incest cases using identity by descent (IBD) patterns.@*METHODS@#The unique IBD pattern was formed by denoting the alleles from the members in a pedigree with a same digital. The probability of each IBD pattern was obtained by multiplying the prior probability by the frequency of non-IBD alleles. The pedigree likelihoods of incest cases under different hypotheses were obtained by summing all IBD pattern probabilities, and the likelihood ratio(LR) was calculated by comparing the likelihoods of different pedigrees.@*RESULTS@#The IBD patterns and the formulae of calculating LR for father-daughter incest and brother-sister incest were obtained.@*CONCLUSIONS@#The calculations of LR for incest cases were illustrated based on IBD patterns.

Evaluation of Detection Efficiency for Trio Full Sibling Testing / 法医学杂志

OBJECTIVES@#To study the detection efficiency of trio full sibling with another known full sibling reference added under different number of autosomal STR typing systems.@*METHODS@#Based on 43 detection systems consisting of 13 to 55 representative autosomal STR loci, 10 000 true families (full sibling group) and 10 000 false families (unrelated individual group) were randomly simulated. The full sibling index (FSI) was calculated based on the method of family reconstruction. The cumulative sibling relationship index (CFSI) of 0.000 1 and 10 000 were used as the evaluation thresholds, and the detection efficiency parameters were calculated and compared with the identification of the duo full sibling testing.@*RESULTS@#With the increasing number of STR loci, the error rate and inability of judgement rate gradually decreased; the sensitivity, specificity, correct rate of judgment and other parameters gradually increased, and the system efficiency gradually improved. Under the same detection system, trio full sibling testing showed higher sensitivity, specificity, system efficiency and lower inability of judgement rate compared with duo full sibling testing. When the system efficiency was higher than 0.85 and inability of judgement rate was less than 0.01%, at least 20 STRs should be detected for trio full sibling testing, which was less than 29 STRs required by duo full sibling testing.@*CONCLUSIONS@#The detection efficiency of trio full sibling testing is superior to that of duo full sibling testing with the same detection system, which is an effective identification scheme for laboratories with inadequate detection systems or for materials with limited conditions.

Forensic efficiencies of individual identification, kinship testing and ancestral inference in three Yunnan groups based on a self-developed multiple DIP panel.

Deletion/insertion polymorphism (DIP), as a short insertion/deletion sequence polymorphic genetic marker, has attracted the attention of forensic genetic scientist due to its lack of stutter, short amplicon and abundant ancestral information. In this study, based on a self-developed 43 autosomal deletion/insertion polymorphism (A-DIP) loci panel which could meet the forensic application purposes of individual identification, kinship testing and ancestral inference to some extent, we evaluated the forensic efficiencies of the above three forensic objectives in Chinese Yi, Hani and Miao groups of Yunnan province. The cumulative match probability (CPM) and combined probability of exclusion (CPE) of these three groups were 1.11433E-18, 8.24299E-19, 4.21721E-18; 0.999610217, 0.999629285 and 0.999582084, respectively. Average 96.65% full sibling pairs could be identified from unrelated individual pairs (as likelihood ratios > 1) using this DIP panel, whereas the average false positive rate was 3.69% in three target Yunnan groups. With the biogeographical ancestor prediction models constructed by extreme gradient boosting (XGBoost) and support vector machine (SVM) algorithms, 0.8239 (95% CI 0.7984, 0.8474) of the unrelated individuals could be correctly divided according to the continental origins based on the 43 A-DIPs which were large frequency distribution differentiations among different continental populations. The present results of principal component analysis (PCA), multidimensional scaling (MDS), neighbor joining (NJ) and maximum likelihood (ML) phylogenetic trees and STRUCTURE analyses indicated that these three Yunnan groups had relatively close genetic distances with East Asian populations.

Pairwise kinship testing with a combination of STR and SNP loci.

Individual commercially available kits exhibit limited discrimination power in full-sibling and second-degree kinship analysis, and therefore they are commonly combined with other kits to obtain more loci and a higher efficacy. However, few studies have systematically evaluated the discrimination power of combined loci. In this study, we combined the ForenSeq™ DNA Signature kit (containing 27 short tandem repeats [STRs] + 91 single nucleotide polymorphisms [SNPs]) with the AGCU NC 21 + 1 PCR amplification kit (containing 21 STRs) to obtain a non-overlapping set of 40 STR and 91 SNP markers. The discrimination power was evaluated for 74 full-sibling pairs, 114 uncle/aunt-nephew/niece pairs and 93 grandparent-grandson/granddaughter pairs. The results show that the efficacy of the 40 STR + 91 SNP combination is higher than the efficacy of either 27 STRs + 91 SNPs or 40 STRs alone. Both the sensitivity and specificity of the 40 STR + 91 SNP marker set achieved 100 % in full-sibling testing, with strong power to distinguish second-degree relatives from unrelated pairs. The 40 STR + 91 SNP set could also distinguish most full-sibling relatives from second-degree relatives but was insufficient to distinguish relatives who belong to the same autosomal kinship class. Our results suggest that ignoring linkage can lead to incorrect likelihood ratios for both related and unrelated pairs, while mutation had a relatively lower effect on the likelihood ratios. Moreover, linkage and mutation had a higher impact on full-sibling testing than on second-degree kinship testing. The discrimination power of the 40 STR and 91 SNP marker set could be strengthened by adding an additional relative.

A new framework to test parent-child and full sibling relationships with population substructure.

Familial search is a statistical approach that is used to infer genetic relationships between a forensic sample and individuals in a DNA database. Several authors have proposed likelihood ratio-based statistics for testing parent-child and full sibling relationships when population substructure exists. This paper proposes three new statistics and investigates performance of each statistic based on Type I error and power. Three statistics, defined by (1) the local allele frequency, (2) the Balding-Nichols approach and (3) the ratio between the maximum of the genotype probabilities over racial subgroup, were found to be good for testing these two types of familial relationships. Power analyses within racial groups are also included, with the power highest for African-American samples and lowest for Caucasian samples. Finally, simulation studies were done on both original and extended CODIS core loci. There were clear differences in power, with the power substantially higher for extended CODIS core loci.

Assessing thermal adaptation using family-based association and FST outlier tests in a threatened trout species.

Discovering genetic markers associated with phenotypic or ecological characteristics can improve our understanding of adaptation and guide conservation of key evolutionary traits. The Lahontan cutthroat trout (Oncorhynchus clarkii henshawi) of the northern Great Basin Desert, USA, demonstrated exceptional tolerance to high temperatures in the desert lakes where it resided historically. This trait is central to a conservation hatchery effort to protect the genetic legacy of the nearly extinct lake ecotype. We genotyped full-sibling families from this conservation broodstock and samples from the only two remaining, thermally distinct, native lake populations at 4,644 new single nucleotide polymorphisms (SNPs). Family-based genome-wide association testing of the broodstock identified nine and 26 SNPs associated with thermal tolerance (p < 0.05 and p < 0.1), measured in a previous thermal challenge experiment. Genes near the associated SNPs had complex functions related to immunity, growth, metabolism and ion homeostasis. Principal component analysis using the thermotolerance-related SNPs showed unexpected divergence between the conservation broodstock and the native lake populations at these loci. FST outlier tests on the native lake populations identified 18 loci shared between two or more of the tests, with two SNPs identified by all three tests (p < 0.01); none overlapped with loci identified by association testing in the broodstock. A recent history of isolation and the complex genetic and demographic backgrounds of Lahontan cutthroat trout probably limited our ability to find shared thermal tolerance loci. Our study extends the still relatively rare application of genomic tools testing for markers associated with important phenotypic or environmental characteristics in species of conservation concern.

Equation Derivation of the Probability Distribution of IBS Score among Full Sibling Pairs.

ABSTRACT: Objective To derive the general equation of the probability distribution of identity by state （IBS） score among biological full sibling pairs by calculating STR allele frequency. Methods Based on the Mendelian genetics law and the hypothesis that parents of biological full siblings （FS） were unrelated individuals, the IBS score and corresponding probability of different genotype combinations in the offspring when unrelated individuals of different genotype combinations give birth to two offsprings were derived. Results Given fi （i=1, 2, …, m） as the frequency of the ith allele of a STR locus, the probability of sharing 2 alleles （p2FS）, 1 allele （p1FS） or 0 allele （p0FS） with biological full sibling pairs on the locus can be respectively expressed as follows： (see the text). The sum of p2FS, p1FS and p0FS must be 1. As for the multiple genotyping system that contained n STR loci, IBS scores between biological full sibling pairs conform to binomial distribution： IBS~B（2n, π1）. The population rate π1, can be given by the formula： (see the text). Conclusion The alternative hypothesis in biological full sibling testing is that two appraised individuals are biological full siblings. The probability of the corresponding alternative hypothesis of any STR locus combination or IBS score can be directly calculated by the equations presented in this study, and the calculation results are the basis for explanations of the evidence.

Formula Derivation for the Probability Distribution of IBS Score in Unrelated Individual Pairs / 法医学杂志

@#Objective To derive the general equation of the probability distribution of identity by state （IBS） score among biological full sibling pairs by calculating STR allele frequency. Methods Based on the Mendelian genetics law and the hypothesis that parents of biological full siblings（FS） were unrelated individuals, the IBS score and corresponding probability of different genotype combinations in the offspring when unrelated individuals of different genotype combinations give birth to two offsprings were derived. Results Given fi （i=1, 2, …, m） as the frequency of the ith allele of a STR locus, the probability of sharing 2 alleles（p2FS）, 1 allele（p1FS） or 0 allele （p0FS） with biological full sibling pairs on the locus can be respectively expressed as follows： p2FS= 14 ×[1 + 2∑im= 1 fi2 + 2（∑im= 1 fi2）2 -∑im= 1 fi4] , p1FS= 12 ×[1 +∑im= 1 fi2 - 2（∑im= 1 fi2）2 - 2∑im= 1 fi3 + 2∑im= 1 fi4] and p0FS= 14 ×[1 - 4∑im= 1 fi2 + 2（∑im= 1 fi2）2 + 4∑im= 1 fi3 -3∑im= 1 fi4] . The sum of p2FS, p1FS and p0FS must be 1. As for the multiple genotyping system that contained n STR loci, IBS scores between biological full sibling pairs conform to binomial distribution：IBS~B（2n, π1）. The population rate π1, can be given by the formula：π1= 1n∑ln= 1 p2FSl + 21n∑ln= 1 p1FSl . Conclusion The alternative hypothesis in biological full sibling testing is that two appraised individuals are biological full siblings. The probability of the corresponding alternative hypothesis of any STR locus combination or IBS score can be directly calculated by the equations presented in this study, and the calculation results are the basis for explanations of the evidence.

Equation Derivation of the Probability Distribution of IBS Score among Full Sibling Pairs / 法医学杂志

Objective To derive the general equation of the probability distribution of identity by state （IBS） score among biological full sibling pairs by calculating STR allele frequency. Methods Based on the Mendelian genetics law and the hypothesis that parents of biological full siblings （FS） were unrelated individuals, the IBS score and corresponding probability of different genotype combinations in the offspring when unrelated individuals of different genotype combinations give birth to two offsprings were derived. Results Given fi （i=1, 2, …, m） as the frequency of the ith allele of a STR locus, the probability of sharing 2 alleles （p2FS）, 1 allele （p1FS） or 0 allele （p0FS） with biological full sibling pairs on the locus can be respectively expressed as follows： (see the text). The sum of p2FS, p1FS and p0FS must be 1. As for the multiple genotyping system that contained n STR loci, IBS scores between biological full sibling pairs conform to binomial distribution： IBS~B（2n, π1）. The population rate π1, can be given by the formula： (see the text). Conclusion The alternative hypothesis in biological full sibling testing is that two appraised individuals are biological full siblings. The probability of the corresponding alternative hypothesis of any STR locus combination or IBS score can be directly calculated by the equations presented in this study, and the calculation results are the basis for explanations of the evidence.

Formula Derivation for the Probability Distribution of IBS Score in Unrelated Individual Pairs.

OBJECTIVES: To derive the probability equation given by STR allele frequencies of identity by state （IBS） score shared by unrelated individual pairs. METHODS: By comparing the STR genotypes of two unrelated individuals, three mutually exclusive combinations could be obtained： （1） sharing 2 identical alleles, a2=1, otherwise a2=0; （2） sharing 1 identical allele, a1=1, otherwise a1=0; （3） sharing 0 identical allele, a0=1, otherwise a0=0. And the IBS score of the one STR locus in this unrelated individual pair could be given by the formula： ibs=2a2+a1. The probability of a2=1 （p2）, a1=1 （p1） and a0=1 （p0） were derived and expressed in powers of the allele frequencies. Subsequently, for a genotyping system including n independent STR loci, the characteristics of binomial distribution of IBS score shared by a pair of unrelated individuals could be given by p2l and p1l （l=1, 2, …, n）. RESULTS: All the general equations of p2, p1 and p0 were derived from the basic conceptions of a2, a1 and a0, respectively. Given fi （i=1, 2, …, m） as the ith allele frequency of a STR locus, the general equations of p2, p1 and p0 could be respectively expressed in powers of fi： [Formula: see text],[Formula: see text] and [Formula: see text]. The sum of p2, p1 and p0 must be equal to 1. Then, the binomial distribution of IBS score shared by unrelated individual pairs genotyped with n independently STR loci could be written by： IBS~B（2n, π）, and the general probability, π, could be given by the formula： [Formula: see text]. CONCLUSIONS: In the biological full sibling identification, the probability of null hypothesis corresponding to any specific IBS score can be directly calculated by the general equations presented in this study, which is the basement of the evidence explanation.

Study of spontaneous mutations in the transmission of poplar chloroplast genomes from mother to offspring.

BACKGROUND: Chloroplasts have their own genomes, independent from nuclear genomes, that play vital roles in growth, which is a major targeted trait for genetic improvement in Populus. Angiosperm chloroplast genomes are maternally inherited, but the chloroplast' variation pattern of poplar at the single-base level during the transmission from mother to offspring remains unknown. RESULTS: Here, we constructed high-quality and almost complete chloroplast genomes for three poplar clones, 'NL895' and its parents, 'I69' and 'I45', from the short-read datasets using multi-pass sequencing (15-16 times per clone) and ultra-high coverage (at least 8500× per clone), with the four-step strategy of Simulation-Assembly-Merging-Correction. Each of the three resulting chloroplast assemblies contained contigs covering > 99% of Populus trichocarpa chloroplast DNA as a reference. A total of 401 variant loci were identified by a hybrid strategy of genome comparison-based and mapping-based single nucleotide polymorphism calling. The genotypes of 94 variant loci were different among the three poplar clones. However, only 1 of the 94 loci was a missense mutation, which was located in the exon region of rpoC1 encoding the ß' subunit of plastid-encoded RNA polymerase. The genotype of the loci in NL895 and its female parent (I69) was different from that of its male parent (I45). CONCLUSIONS: This research provides resources for further chloroplast genomic studies of a F1 full-sibling family derived from a cross between I69 and I45, and will improve the application of chloroplast genomic information in modern Populus breeding programs.

Formula Derivation for the Probability Distribution of IBS Score in Unrelated Individual Pairs / 法医学杂志

OBJECTIVES@#To derive the probability equation given by STR allele frequencies of identity by state （IBS） score shared by unrelated individual pairs.@*METHODS@#By comparing the STR genotypes of two unrelated individuals, three mutually exclusive combinations could be obtained： （1） sharing 2 identical alleles, a₂=1, otherwise a₂=0; （2） sharing 1 identical allele, a₁=1, otherwise a₁=0; （3） sharing 0 identical allele, a₀=1, otherwise a₀=0. And the IBS score of the one STR locus in this unrelated individual pair could be given by the formula： ibs=2a₂+a₁. The probability of a₂=1 （p₂）, a₁=1 （p₁） and a₀=1 （p₀） were derived and expressed in powers of the allele frequencies. Subsequently, for a genotyping system including n independent STR loci, the characteristics of binomial distribution of IBS score shared by a pair of unrelated individuals could be given by p₂l and p₁l （l=1, 2, …, n）.@*RESULTS@#All the general equations of p₂, p₁ and p₀ were derived from the basic conceptions of a₂, a₁ and a₀, respectively. Given fi （i=1, 2, …, m） as the ith allele frequency of a STR locus, the general equations of p₂, p₁ and p₀ could be respectively expressed in powers of fi： [Formula: see text],[Formula: see text] and [Formula: see text]. The sum of p₂, p₁ and p₀ must be equal to 1. Then, the binomial distribution of IBS score shared by unrelated individual pairs genotyped with n independently STR loci could be written by： IBS~B（2n, π）, and the general probability, π, could be given by the formula： [Formula: see text].@*CONCLUSIONS@#In the biological full sibling identification, the probability of null hypothesis corresponding to any specific IBS score can be directly calculated by the general equations presented in this study, which is the basement of the evidence explanation.

Development of F1 hybrid population and the high-density linkage map for European aspen (Populus tremula L.) using RADseq technology.

BACKGROUND: Restriction-site associated DNA sequencing (RADseq) technology was recently employed to identify a large number of single nucleotide polymorphisms (SNP) for linkage mapping of a North American and Eastern Asian Populus species. However, there is also the need for high-density genetic linkage maps for the European aspen (P. tremula) as a tool for further mapping of quantitative trait loci (QTLs) and marker-assisted selection of the Populus species native to Europe. RESULTS: We established a hybrid F1 population from the cross of two aspen parental genotypes diverged in their phenological and morphological traits. We performed RADseq of 122 F1 progenies and two parents yielding 15,732 high-quality SNPs that were successfully identified using the reference genome of P. trichocarpa. 2055 SNPs were employed for the construction of maternal and paternal linkage maps. The maternal linkage map was assembled with 1000 SNPs, containing 19 linkage groups and spanning 3054.9 cM of the genome, with an average distance of 3.05 cM between adjacent markers. The paternal map consisted of 1055 SNPs and the same number of linkage groups with a total length of 3090.56 cM and average interval distance of 2.93 cM. The linkage maps were employed for QTL mapping of one-year-old seedlings height variation. The most significant QTL (LOD = 5.73) was localized to LG5 (96.94 cM) of the male linkage map, explaining 18% of the phenotypic variation. CONCLUSIONS: The set of 15,732 SNPs polymorphic in aspen and high-density genetic linkage maps constructed for the P. tremula intra-specific cross will provide a valuable source for QTL mapping and identification of candidate genes facilitating marker-assisted selection in European aspen.

Set up of cutoff thresholds for kinship determination using SNP loci.

The usefulness of single nucleotide polymorphism (SNP) loci for kinship testing has been demonstrated in many case works, and suggested as a promising marker for relationship identification. For interpreting results based on the calculation of the likelihood ratio (LR) in kinship testing, it is important to prepare cutoffs for respective relatives which are dependent on genetic relatedness. For this, analysis using true pedigree data is significant and reliable as it reflects the actual frequencies of markers in the population. In this study, the kinship index was explored through 1209 parent-child pairs, 1373 full sibling pairs, and 247 uncle-nephew pairs using 136 SNP loci. The cutoffs for LR were set up using different numbers of SNP loci with accuracy, sensitivity, and specificity. It is expected that this study can support the application of SNP loci-based kinship testing for various relationships.

Bayesian pedigree inference with small numbers of single nucleotide polymorphisms via a factor-graph representation.

We develop a computational framework for addressing pedigree inference problems using small numbers (80-400) of single nucleotide polymorphisms (SNPs). Our approach relaxes the assumptions, which are commonly made, that sampling is complete with respect to the pedigree and that there is no genotyping error. It relies on representing the inferred pedigree as a factor graph and invoking the Sum-Product algorithm to compute and store quantities that allow the joint probability of the data to be rapidly computed under a large class of rearrangements of the pedigree structure. This allows efficient MCMC sampling over the space of pedigrees, and, hence, Bayesian inference of pedigree structure. In this paper we restrict ourselves to inference of pedigrees without loops using SNPs assumed to be unlinked. We present the methodology in general for multigenerational inference, and we illustrate the method by applying it to the inference of full sibling groups in a large sample (n=1157) of Chinook salmon typed at 95 SNPs. The results show that our method provides a better point estimate and estimate of uncertainty than the currently best-available maximum-likelihood sibling reconstruction method. Extensions of this work to more complex scenarios are briefly discussed.

Differential immune gene expression profiles in susceptible and resistant full-sibling families of Atlantic salmon (Salmo salar) challenged with infectious pancreatic necrosis virus (IPNV).

This study aims to identify at the expression level the immune-related genes associated with IPN-susceptible and resistant phenotypes in Atlantic salmon full-sibling families. We have analyzed thirty full-sibling families infected by immersion with IPNV and then classified as resistant or susceptible using a multivariate survival analysis based on a gamma-Cox frailty model and the Kaplan-Meier mortality curves. In four families within each group head kidneys were pooled for real-time PCR and one-color salmon-specific oligonucleotide microarray (21K) analysis at day 1 and 5 post-infection. Transcripts involved in innate response (IL-6, IFN-α), antigen presentation (HSP-70, HSP-90, MHC-I), TH1 response (IL-12, IFN-γ, CRFB6), immunosuppression (IL-10, TGF-ß1) and leukocyte activation and migration (CCL-19, CD18) showed a differential expression pattern between both phenotypes, except in IL-6. In susceptible families, except for IFN-γ, the expressions dropped to basal values at day 5 post-infection. In resistant families, unlike susceptible families, levels remained high or increased (except for IL-6) at day 5. Transcriptomic analysis showed that both families have a clear differential expression pattern, resulting in a marked down-regulation in immune related genes involved in innate response, complement system, antigen recognition and activation of immune response in IPN-resistant. Down-regulation of genes, mainly related to tissue differentiation and protein degradation metabolism, was also observed in resistant families. We have identified an immune-related gene patterns associated with susceptibility and resistance to IPNV infection of Atlantic salmon. This suggests that a limited immune response is associated with resistant fish phenotype to IPNV challenge while a highly inflammatory but short response is associated with susceptibility.