A general and efficient method for estimating continuous IBD functions for use in genome scans for QTL

<p>Abstract</p> <p>Background</p> <p>Identity by descent (IBD) matrix estimation is a central component in mapping of Quantitative Trait Loci (QTL) using variance component models. A large number of algorithms have been developed for estimation of IBD between individuals in populations at discrete locations in the genome for use in genome scans to detect QTL affecting various traits of interest in experimental animal, human and agricultural pedigrees. Here, we propose a new approach to estimate IBD as continuous functions rather than as discrete values.</p> <p>Results</p> <p>Estimation of IBD functions improved the computational efficiency and memory usage in genome scanning for QTL. We have explored two approaches to obtain continuous marker-bracket IBD-functions. By re-implementing an existing and fast deterministic IBD-estimation method, we show that this approach results in IBD functions that produces the exact same IBD as the original algorithm, but with a greater than 2-fold improvement of the computational efficiency and a considerably lower memory requirement for storing the resulting genome-wide IBD. By developing a general IBD function approximation algorithm, we show that it is possible to estimate marker-bracket IBD functions from IBD matrices estimated at marker locations by any existing IBD estimation algorithm. The general algorithm provides approximations that lead to QTL variance component estimates that even in worst-case scenarios are very similar to the true values. The approach of storing IBD as polynomial IBD-function was also shown to reduce the amount of memory required in genome scans for QTL.</p> <p>Conclusion</p> <p>In addition to direct improvements in computational and memory efficiency, estimation of IBD-functions is a fundamental step needed to develop and implement new efficient optimization algorithms for high precision localization of QTL. Here, we discuss and test two approaches for estimating IBD functions based on existing IBD estimation algorithms. Our approaches provide immediately useful techniques for use in single QTL analyses in the variance component QTL mapping framework. They will, however, be particularly useful in genome scans for multiple interacting QTL, where the improvements in both computational and memory efficiency are the key for successful development of efficient optimization algorithms to allow widespread use of this methodology.</p

Modelling dominance in a flexible intercross analysis

<p>Abstract</p> <p>Background</p> <p>The aim of this paper is to develop a flexible model for analysis of quantitative trait loci (QTL) in outbred line crosses, which includes both additive and dominance effects. Our flexible intercross analysis (FIA) model accounts for QTL that are not fixed within founder lines and is based on the variance component framework. Genome scans with FIA are performed using a score statistic, which does not require variance component estimation.</p> <p>Results</p> <p>Simulations of a pedigree with 800 <it>F</it>₂individuals showed that the power of FIA including both additive and dominance effects was almost 50% for a QTL with equal allele frequencies in both lines with complete dominance and a moderate effect, whereas the power of a traditional regression model was equal to the chosen significance value of 5%. The power of FIA without dominance effects included in the model was close to those obtained for FIA with dominance for all simulated cases except for QTL with overdominant effects. A genome-wide linkage analysis of experimental data from an <it>F</it>₂intercross between Red Jungle Fowl and White Leghorn was performed with both additive and dominance effects included in FIA. The score values for chicken body weight at 200 days of age were similar to those obtained in FIA analysis without dominance.</p> <p>Conclusion</p> <p>We have extended FIA to include QTL dominance effects. The power of FIA was superior, or similar, to standard regression methods for QTL effects with dominance. The difference in power for FIA with or without dominance is expected to be small as long as the QTL effects are not overdominant. We suggest that FIA with only additive effects should be the standard model to be used, especially since it is more computationally efficient.</p

Mutations in or near the Transmembrane Domain Alter PMEL Amyloid Formation from Functional to Pathogenic

PMEL is a pigment cell-specific protein that forms physiological amyloid fibrils upon which melanins ultimately deposit in the lumen of the pigment organelle, the melanosome. Whereas hypomorphic PMEL mutations in several species result in a mild pigment dilution that is inherited in a recessive manner, PMEL alleles found in the Dominant white (DW) chicken and Silver horse (HoSi)—which bear mutations that alter the PMEL transmembrane domain (TMD) and that are thus outside the amyloid core—are associated with a striking loss of pigmentation that is inherited in a dominant fashion. Here we show that the DW and HoSi mutations alter PMEL TMD oligomerization and/or association with membranes, with consequent formation of aberrantly packed fibrils. The aberrant fibrils are associated with a loss of pigmentation in cultured melanocytes, suggesting that they inhibit melanin production and/or melanosome integrity. A secondary mutation in the Smoky chicken, which reverts the dominant DW phenotype, prevents the accumulation of PMEL in fibrillogenic compartments and thus averts DW–associated pigment loss; a secondary mutation found in the Dun chicken likely dampens a HoSi–like dominant mutation in a similar manner. We propose that the DW and HoSi mutations alter the normally benign amyloid to a pathogenic form that antagonizes melanosome function, and that the secondary mutations found in the Smoky and Dun chickens revert or dampen pathogenicity by functioning as null alleles, thus preventing the formation of aberrant fibrils. We speculate that PMEL mutations can model the conversion between physiological and pathological amyloid

Inactivation of Pmel Alters Melanosome Shape But Has Only a Subtle Effect on Visible Pigmentation

PMEL is an amyloidogenic protein that appears to be exclusively expressed in pigment cells and forms intralumenal fibrils within early stage melanosomes upon which eumelanins deposit in later stages. PMEL is well conserved among vertebrates, and allelic variants in several species are associated with reduced levels of eumelanin in epidermal tissues. However, in most of these cases it is not clear whether the allelic variants reflect gain-of-function or loss-of-function, and no complete PMEL loss-of-function has been reported in a mammal. Here, we have created a mouse line in which the Pmel gene has been inactivated (Pmel−/−). These mice are fully viable, fertile, and display no obvious developmental defects. Melanosomes within Pmel−/− melanocytes are spherical in contrast to the oblong shape present in wild-type animals. This feature was documented in primary cultures of skin-derived melanocytes as well as in retinal pigment epithelium cells and in uveal melanocytes. Inactivation of Pmel has only a mild effect on the coat color phenotype in four different genetic backgrounds, with the clearest effect in mice also carrying the brown/Tyrp1 mutation. This phenotype, which is similar to that observed with the spontaneous silver mutation in mice, strongly suggests that other previously described alleles in vertebrates with more striking effects on pigmentation are dominant-negative mutations. Despite a mild effect on visible pigmentation, inactivation of Pmel led to a substantial reduction in eumelanin content in hair, which demonstrates that PMEL has a critical role for maintaining efficient epidermal pigmentation

Targeted analysis of four breeds narrows equine Multiple Congenital Ocular Anomalies locus to 208 kilobases

The syndrome Multiple Congenital Ocular Anomalies (MCOA) is the collective name ascribed to heritable congenital eye defects in horses. Individuals homozygous for the disease allele (MCOA phenotype) have a wide range of eye anomalies, while heterozygous horses (Cyst phenotype) predominantly have cysts that originate from the temporal ciliary body, iris, and/or peripheral retina. MCOA syndrome is highly prevalent in the Rocky Mountain Horse but the disease is not limited to this breed. Affected horses most often have a Silver coat color; however, a pleiotropic link between these phenotypes is yet to be proven. Locating and possibly isolating these traits would provide invaluable knowledge to scientists and breeders. This would favor maintenance of a desirable coat color while addressing the health concerns of the affected breeds, and would also provide insight into the genetic basis of the disease. Identical-by-descent mapping was used to narrow the previous 4.6-Mb region to a 264-kb interval for the MCOA locus. One haplotype common to four breeds showed complete association to the disease (Cyst phenotype, n = 246; MCOA phenotype, n = 83). Candidate genes from the interval, SMARCC2 and IKZF4, were screened for polymorphisms and genotyped, and segregation analysis allowed the MCOA syndrome region to be shortened to 208 kb. This interval also harbors PMEL17, the gene causative for Silver coat color. However, by shortening the MCOA locus by a factor of 20, 176 other genes have been unlinked from the disease and only 15 genes remain

Detecting parent of origin and dominant QTL in a two-generation commercial poultry pedigree using variance component methodology

<p>Abstract</p> <p>Introduction</p> <p>Variance component QTL methodology was used to analyse three candidate regions on chicken chromosomes 1, 4 and 5 for dominant and parent-of-origin QTL effects. Data were available for bodyweight and conformation score measured at 40 days from a two-generation commercial broiler dam line. One hundred dams were nested in 46 sires with phenotypes and genotypes on 2708 offspring. Linear models were constructed to simultaneously estimate fixed, polygenic and QTL effects. Different genetic models were compared using likelihood ratio test statistics derived from the comparison of full with reduced or null models. Empirical thresholds were derived by permutation analysis.</p> <p>Results</p> <p>Dominant QTL were found for bodyweight on chicken chromosome 4 and for bodyweight and conformation score on chicken chromosome 5. Suggestive evidence for a maternally expressed QTL for bodyweight and conformation score was found on chromosome 1 in a region corresponding to orthologous imprinted regions in the human and mouse.</p> <p>Conclusion</p> <p>Initial results suggest that variance component analysis can be applied within commercial populations for the direct detection of segregating dominant and parent of origin effects.</p

The Crest Phenotype in Chicken Is Associated with Ectopic Expression of HOXC8 in Cranial Skin

The Crest phenotype is characterised by a tuft of elongated feathers atop the head. A similar phenotype is also seen in several wild bird species. Crest shows an autosomal incompletely dominant mode of inheritance and is associated with cerebral hernia. Here we show, using linkage analysis and genome-wide association, that Crest is located on the E22C19W28 linkage group and that it shows complete association to the HOXC-cluster on this chromosome. Expression analysis of tissues from Crested and non-crested chickens, representing 26 different breeds, revealed that HOXC8, but not HOXC12 or HOXC13, showed ectopic expression in cranial skin during embryonic development. We propose that Crest is caused by a cis-acting regulatory mutation underlying the ectopic expression of HOXC8. However, the identification of the causative mutation(s) has to await until a method becomes available for assembling this chromosomal region. Crest is unfortunately located in a genomic region that has so far defied all attempts to establish a contiguous sequence

A Genome-Wide SNP Scan Reveals Novel Loci for Egg Production and Quality Traits in White Leghorn and Brown-Egg Dwarf Layers

Availability of the complete genome sequence as well as high-density SNP genotyping platforms allows genome-wide association studies (GWAS) in chickens. A high-density SNP array containing 57,636 markers was employed herein to identify associated variants underlying egg production and quality traits within two lines of chickens, i.e., White Leghorn and brown-egg dwarf layers. For each individual, age at first egg (AFE), first egg weight (FEW), and number of eggs (EN) from 21 to 56 weeks of age were recorded, and egg quality traits including egg weight (EW), eggshell weight (ESW), yolk weight (YW), eggshell thickness (EST), eggshell strength (ESS), albumen height(AH) and Haugh unit(HU) were measured at 40 and 60 weeks of age. A total of 385 White Leghorn females and 361 brown-egg dwarf dams were selected to be genotyped. The genome-wide scan revealed 8 SNPs showing genome-wise significant (P<1.51E-06, Bonferroni correction) association with egg production and quality traits under the Fisher's combined probability method. Some significant SNPs are located in known genes including GRB14 and GALNT1 that can impact development and function of ovary, but more are located in genes with unclear functions in layers, and need to be studied further. Many chromosome-wise significant SNPs were also detected in this study and some of them are located in previously reported QTL regions. Most of loci detected in this study are novel and the follow-up replication studies may be needed to further confirm the functional significance for these newly identified SNPs

Genome-Wide Association Study Identified a Narrow Chromosome 1 Region Associated with Chicken Growth Traits

Chicken growth traits are important economic traits in broilers. A large number of studies are available on finding genetic factors affecting chicken growth. However, most of these studies identified chromosome regions containing putative quantitative trait loci and finding causal mutations is still a challenge. In this genome-wide association study (GWAS), we identified a narrow 1.5 Mb region (173.5–175 Mb) of chicken (Gallus gallus) chromosome (GGA) 1 to be strongly associated with chicken growth using 47,678 SNPs and 489 F2 chickens. The growth traits included aggregate body weight (BW) at 0–90 d of age measured weekly, biweekly average daily gains (ADG) derived from weekly body weight, and breast muscle weight (BMW), leg muscle weight (LMW) and wing weight (WW) at 90 d of age. Five SNPs in the 1.5 Mb KPNA3-FOXO1A region at GGA1 had the highest significant effects for all growth traits in this study, including a SNP at 8.9 Kb upstream of FOXO1A for BW at 22–48 d and 70 d, a SNP at 1.9 Kb downstream of FOXO1A for WW, a SNP at 20.9 Kb downstream of ENSGALG00000022732 for ADG at 29–42 d, a SNP in INTS6 for BW at 90 d, and a SNP in KPNA3 for BMW and LMW. The 1.5 Mb KPNA3-FOXO1A region contained two microRNA genes that could bind to messenger ribonucleic acid (mRNA) of IGF1, FOXO1A and KPNA3. It was further indicated that the 1.5 Mb GGA1 region had the strongest effects on chicken growth during 22–42 d