Prediction of multilocus identity-by-descent

Previous studies have enabled exact prediction of probabilities of identity-by-descent (IBD) in randommating populations for a few loci (up to four or so), with extension to more using approximate regression methods. Here we present a precise predictor of multiple-locus IBD using simple formulas based on exact results for two loci. In particular, the probability of non-IBD X ABC at each of ordered loci A, B, and C can be well approximated by XABC = XABXBC/XB and generalizes to X123. . .k = X12X23. . .Xk-1,k/ Xk-2, where X is the probability of non-IBD at each locus. Predictions from this chain rule are very precise with population bottlenecks and migration, but are rather poorer in the presence of mutation. From these coefficients, the probabilities of multilocus IBD and non-IBD can also be computed for genomic regions as functions of population size, time, and map distances. An approximate but simple recurrence formula is also developed, which generally is less accurate than the chain rule but is more robust with mutation. Used together with the chain rule it leads to explicit equations for non-IBD in a region. The results can be applied to detection of quantitative trait loci (QTL) by computing the probability of IBD at candidate loci in terms of identity-by-state at neighboring markers.

Identity-by-Descent Mapping to Detect Rare Variants Conferring Susceptibility to Multiple Sclerosis

Genome-wide association studies (GWAS) have identified around 60 common variants associated with multiple sclerosis (MS), but these loci only explain a fraction of the heritability of MS. Some missing heritability may be caused by rare variants that have been suggested to play an important role in the aetiology of complex diseases such as MS. However current genetic and statistical methods for detecting rare variants are expensive and time consuming. 'Population-based linkage analysis' (PBLA) or so called identity-by-descent (IBD) mapping is a novel way to detect rare variants in extant GWAS datasets. We employed BEAGLE fastIBD to search for rare MS variants utilising IBD mapping in a large GWAS dataset of 3,543 cases and 5,898 controls. We identified a genome-wide significant linkage signal on chromosome 19 (LOD = 4.65; p = 1.9×10-6). Network analysis of cases and controls sharing haplotypes on chromosome 19 further strengthened the association as there are more large networks of cases sharing haplotypes than controls. This linkage region includes a cluster of zinc finger genes of unknown function. Analysis of genome wide transcriptome data suggests that genes in this zinc finger cluster may be involved in very early developmental regulation of the CNS. Our study also indicates that BEAGLE fastIBD allowed identification of rare variants in large unrelated population with moderate computational intensity. Even with the development of whole-genome sequencing, IBD mapping still may be a promising way to narrow down the region of interest for sequencing priority. © 2013 Lin et al.

Identity by descent and candidate gene mapping of Richieri-Costa and Pereira syndrome

The Richieri-Costa-Pereira syndrome is a rare autosomal recessive disorder characterized by short stature, Robin sequence, cleft mandible, pre/postaxial anomalies and clubfoot. of 15 families reported with this disorder 14 are from Brazil suggesting a founder effect. We studied 15 families using identity-by-descent as a hypothesis to attempt gene localization We have examined through linkage analysis 497 polymorphicmarkers and also performed direct sequencing of exons for 10 candidate genes selected on the basis of their expression in the developing mandible and limb. No evidence for allele sharing at any locus tested or mutations in candidate genes was found. Additional higher resolution mapping, new families and other candidate genes might improve future chances of gene identification. (C) 2003 Wiley-Liss, Inc.

Comparison of identity by descent and identity by state for detecting genetic regions under selection in a sorghum pedigree breeding program

Seventy sorghum inbred lines which formed part of the Queensland Department of Primary Industries (QDPI) sorghum breeding program were screened with 104 previously mapped RFLP markers. The lines were related by pedigree and consisted of ancestral source lines, intermediate lines and recent releases from the program. We compared the effect of defining marker alleles using either identity by state (IBS) or identity by descent (IBD) on our capacity to trace markers through the pedigree and detect evidence of selection for particular alleles. Allelic identities defined using IBD were much more sensitive for detecting non-Mendelian segregation in this pedigree. Only one marker allele showed significant evidence of selection when IBS was used compared with ten regions with particular allelic identities when IBD was used. Regions under selection were compared with the location of QTLs for agronomic traits known to be under selection in the breeding program. Only two of the ten regions were associated with known QTLs that matched with knowledge of the agronomic characteristics of the ancestral lines. Some of the other regions were hypothesised to be associated with genes for particular traits based on the properties of the ancestral source lines.

Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within- family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.

Using identity-by-de scent information in affected sib pairs to increase the efficiency of genetic association studies

T he aim of this study was to determine whether identity-by-descent (IBD) information for affected sib pairs (ASPs) can be used to select a sample of cases for a genetic case-control study which will provide more power for detecting association with loci in a known linkage region. By modeling the expected frequency of the disease allele in ASPs showing IBD sharing of 0, 1, or 2 alleles, and considering additive, recessive, and dominant disease models, we show that cases selected from IBD 2 families are best for this purpose, followed by those selected from IBD 1 families; least useful are cases selected from IBD 0 families.

Effect of IGF1, GH, and PIT1 markers on the genetic parameters of growth and reproduction traits in Canchim cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Myoepithelial molecular markers in human breast carcinoma PMC42-LA cells are induced by extracellular matrix and stromal cells

The microenvironment plays a key role in the cellular differentiation of the two main cell lineages of the human breast, luminal epithelial, and myoepithelial. It is not clear, however, how the components of the microenvironment control the development of these cell lineages. To investigate how lineage development is regulated by 3-D culture and microenvironment components, we used the PMC42-LA human breast carcinoma cell line, which possesses stem cell characteristics. When cultured on a two-dimensional glass substrate, PMC42-LA cells formed a monolayer and expressed predominantly luminal epithelial markers, including cytokeratins 8, 18, and 19; E-cadherin; and sialomucin. The key myoepithelial-specific proteins α-smooth muscle actin and cytokeratin 14 were not expressed. When cultured within Engelbreth-Holm- Swarm sarcoma-derived basement membrane matrix (EHS matrix), PMC42-LA cells formed organoids in which the expression of luminal markers was reduced and the expression of other myoepithelial-specific markers (cytokeratin 17 and P-cadherin) was promoted. The presence of primary human mammary gland fibroblasts within the EHS matrix induced expression of the key myoepithelial-specific markers, α-smooth muscle actin and cytokeratin 14. Immortalized human skin fibroblasts were less effective in inducing expression of these key myoepithelial-specific markers. Confocal dual-labeling showed that individual cells expressed luminal or myoepithelial proteins, but not both. Conditioned medium from the mammary fibroblasts was equally effective in inducing myoepithelial marker expression. The results indicate that the myoepithelial lineage is promoted by the extracellular matrix, in conjunction with products secreted by breast-specific fibroblasts. Our results demonstrate a key role for the breast microenvironment in the regulation of breast lineage development.

Marker selection by Akaike information criterion and Bayesian information criterion

We carried out a discriminant analysis with identity by descent (IBD) at each marker as inputs, and the sib pair type (affected-affected versus affected-unaffected) as the output. Using simple logistic regression for this discriminant analysis, we illustrate the importance of comparing models with different number of parameters. Such model comparisons are best carried out using either the Akaike information criterion (AIC) or the Bayesian information criterion (BIC). When AIC (or BIC) stepwise variable selection was applied to the German Asthma data set, a group of markers were selected which provide the best fit to the data (assuming an additive effect). Interestingly, these 25-26 markers were not identical to those with the highest (in magnitude) single-locus lod scores.

Investigation of mode characteristics for microdisk resonators by S-matrix and three-dimensional finite-difference time-domain technique

The mode characteristics of a three-dimensional (3D) microdisk with a vertical refractive index distribution of n(2)/3.4/n(2) are investigated by the S-matrix method and 3D finite-difference time-domain (FDTD) technique. For the microdisk with a thickness of 0.2 mu m. and a radius of 1 mu m, the mode wavelengths and quality factors for the HE7,1 mode obtained by 3D FDTD simulation and the S-matrix method are in good agreement as n(2) increases from 1.0 to 2.6. But the Q factor obtained by the 3D FDTD rapidly decreases from 1.12 X 10(4) to 379 as n2 increases from 2.65 to 2.8 owing to the vertical radiation losses, which cannot be predicted by the proposed S-matrix method. The comparisons also show that quality factors obtained from the analytical solution of two-dimensional microdisks under the effective index approximation are five to seven times smaller than those of the 3D FDTD as n(2) = 1 and R = 1 mu m. (c) 2006 Optical Society of America.

Kingship by Descent or Kingship by Election? The Contested Title of James VI and I

Throughout the reign of Elizabeth I, a steady stream of tracts appeared in English print to vindicate the succession of the most prominent contenders, Mary and James Stuart of Scotland. This article offers a comprehensive account of the polemical battle between the supporters and opponents of the Stuarts, and further identifies various theories of English kingship, most notably the theory of corporate kingship, developed by the Stuart polemicists to defend the Scottish succession. James's accession to the English throne in March 1603 marked the protracted end of the debate over the succession. The article concludes by suggesting that, while powerfully renouncing the opposition to his succession, over the course of his attempt to unify his two kingdoms, James and his supporters ultimately departed from the polemic of corporate kingship, for a more assertive language of kingship by natural and divine law.

Novel Statistical Methods in Quantitative Genetics : Modeling Genetic Variance for Quantitative Trait Loci Mapping and Genomic Evaluation

This thesis develops and evaluates statistical methods for different types of genetic analyses, including quantitative trait loci (QTL) analysis, genome-wide association study (GWAS), and genomic evaluation. The main contribution of the thesis is to provide novel insights in modeling genetic variance, especially via random effects models. In variance component QTL analysis, a full likelihood model accounting for uncertainty in the identity-by-descent (IBD) matrix was developed. It was found to be able to correctly adjust the bias in genetic variance component estimation and gain power in QTL mapping in terms of precision.  Double hierarchical generalized linear models, and a non-iterative simplified version, were implemented and applied to fit data of an entire genome. These whole genome models were shown to have good performance in both QTL mapping and genomic prediction. A re-analysis of a publicly available GWAS data set identified significant loci in Arabidopsis that control phenotypic variance instead of mean, which validated the idea of variance-controlling genes.  The works in the thesis are accompanied by R packages available online, including a general statistical tool for fitting random effects models (hglm), an efficient generalized ridge regression for high-dimensional data (bigRR), a double-layer mixed model for genomic data analysis (iQTL), a stochastic IBD matrix calculator (MCIBD), a computational interface for QTL mapping (qtl.outbred), and a GWAS analysis tool for mapping variance-controlling loci (vGWAS).

Acheiropodia is caused by a genomic deletion in C7orf2, the human orthologue of the Lmbr1 gene

Acheiropodia is an autosomal recessive developmental disorder presenting with bilateral congenital amputations of the upper and lower extremities and aplasia of the hands and feet. This severely handicapping condition appears to affect only the extremities, with no other systemic manifestations reported. Recently, a locus for acheiropodia was mapped on chromosome 7q36. Herein we report the narrowing of the critical region for the acheiropodia gene and the subsequent identification of a common mutation in C7orf2-the human orthologue of the mouse Lmbr1 gene-that is responsible for the disease. Analysis of five families with acheiropodia, by means of 15 polymorphic markers, narrowed the critical region to 1.3 cM, on the basis of identity by descent, and to <0.5 Mb, on the basis of physical mapping. Analysis of C7orf2, the human orthologue of the mouse Lmbr1 gene, identified a deletion in all five families, thus identifying a common acheiropodia mutation. The deletion was identified at both the genomic-DNA and mRNA level. It leads to the production of a C7orf2 transcript lacking exon 4 and introduces a premature stop codon downstream of exon 3. Given the nature of the acheiropodia phenotype, it appears likely that the Lmbr1 gene plays an important role in limb development.