Parallel Evolution of Polydactyly Traits in Chinese and European Chickens

Polydactyly is one of the most common hereditary congenital limb malformations in chickens and other vertebrates. The zone of polarizing activity regulatory sequence (ZRS) is critical for the development of polydactyly. The causative mutation of polydactyly in the Silkie chicken has been mapped to the ZRS; however, the causative mutations of other chicken breeds are yet to be established. To understand whether the same mutation decides the polydactyly phenotype in other chicken breeds, we detected the single-nucleotide polymorphism in 26 different chicken breeds, specifically, 24 Chinese indigenous breeds and 2 European breeds. The mutation was found to have fully penetrated chickens with polydactyly in China, indicating that it is causative for polydactyly in Chinese indigenous chickens. In comparison, the mutation showed no association with polydactyly in Houdan chickens, which originate from France, Europe. Based on the different morphology of polydactyly in Chinese and European breeds, we assumed that the trait might be attributable to different genetic foundations. Therefore, we subsequently performed genome-wide association analysis (GWAS) to locate the region associated with polydactyly. As a result, a ~0.39 Mb genomic region on GGA2p was identified. The region contains six candidate genes, with the causative mutation found in Chinese indigenous breeds also being located in this region. Our results demonstrate that polydactyly in chickens from China and Europe is caused by two independent mutation events that are closely located in the chicken genome.     

Comparative genomic analysis of the false killer whale (Pseudorca crassidens) LMBR1 locus

The sequencing and comparative genomic analysis of LMBR1 loci in mammals or other species, including human, would be very important in understanding evolutionary genetic changes underlying the evolution of limb development. In this regard, comparative genomic annotation of the false killer whale LMBR1 locus could shed new light on the evolution of limb development. We sequenced two false killer whale BAC clones, corresponding to 156 kb and 144 kb, respectively, harboring the tightly linked RNF32, LMBR1, and NOM1 genes. Our annotation of the false killer whale LMBR1 gene showed that it consists of 17 exons (1473 bp), in contrast to 18 exons (1596 bp) in human, and it displays 93.1% and 95.6% nucleotide and amino acid sequence similarity, respectively, compared with the human gene. In particular, we discovered that exon 10, deleted in the false killer whale LMBR1 gene, is present only in primates, and this fact strongly implies that exon 10 might be crucial in determining primate-specific limb development. ZRS and TFBS sequences have been well conserved across 11 species, suggesting that these regions could be involved in an important function of limb development and limb patterning. The neighboring gene RNF32 showed several lineage-conserved exons, such as exons 2 through 9 conserved in eutherian mammals, exons 3 through 9 conserved in mammals, and exons 5 through 9 conserved in vertebrates. The other neighboring gene, NOM1, had undergone a substitution (ATG→GTA) at the start codon, giving rise to a 36 bp shorter N-terminal sequence compared with the human sequence. Our comparative analysis of the false killer whale LMBR1 genomic locus provides important clues regarding the genetic regions that may play crucial roles in limb development and patterning.     

A single C to T transition in intron 5 of LMBR1 gene is associated with triphalangeal thumb-polysyndactyly syndrome in a Chinese family

Triphalangeal thumb-polysyndactyly syndrome (TPT-PS) is a type of human hand-foot malformation. In this study, we collected data from a Chinese family with TPT-PS and mapped the disease region to chromosome 7q36. By using a fine mapping study and a haplotype analysis, we narrowed the affected region to 1.7cM between markers D7S2465 and D7S2423, which contains four candidate genes: HLXB9, LMBR1, NOM1, and RNF32. By sequence analysis, we found no sequence alterations, which are specific to the patients in the transcribed regions and in the intron-exon boundaries among the four genes. After closely examining intron 5 of the LMBR1 gene, we discovered a single C to T transition in the affected TPT-PS individuals of the Chinese subject family. The position of this C to T transition is located close to other sequence alterations involved in several preaxial polydactyly (PPD) families, supporting the notion that intron 5 of LMBR1 contains a cis-acting regulator of limb-specific Sonic Hedgehog (SHH). We postulate that the disruption of this cis-regulator via a single C to T transition results in the dysregulation of SHH, which leads to the TPT-PS found in this case.     

The expression of preaxial polydactyly is influenced by modifying genetic elements and is not maintained by chromosomal inversion in an avian biomedical model

Polydactyly (Po) is a common mutation found in many vertebrates. The UCD-Po.003 congenic chicken line was previously characterized for Po inheritance (autosomal dominant) and the mutation was mapped to chromosome 2p. Here, we describe for the first time the range and variability of the phenotype in this congenic line. Further, we studied the hypothesis that a chromosomal inversion was responsible for the maintenance of a large (6.3 Mb) candidate gene region. Fluorescence in situ hybridization employing BACs encompassing a 10.7-Mb region of GGA2p showed that the Po chromosome was normal, i.e. exhibits the wild-type BAC order. Continued fine-mapping along with a change in breeding strategy reduced the size of the causative region to 1.43 Mb. Recent research indicates that the cause of preaxial Po resides within a 794-bp highly conserved zone of polarizing activity regulatory sequence (ZRS) element located in intron 5 of the LMBR1 gene; however, the ZRS polymorphism of interest is found in some but not all breeds of polydactylous chicken. Therefore, we sequenced the ZRS in 101 heterozygous and 30 unaffected (wild-type) individuals to establish the relevance of this region to the Po condition in the UCD-Po.003 congenic line. A single point mutation (C/A at coordinate GGA2p: 8,414,121) within the ZRS segregated with carrier status. The polydactylous UCD-Silkie line also maintains this SNP in addition to a single base deletion. An inheritance analysis of the phenotypic variation in UCD-Po.003 suggests recessive epistasis as the mode of inheritance for the additional modifying genetic elements, residing outside the ZRS, to impact the preaxial polydactyl phenotype. These results contribute to our understanding of the cause of Po in an important vertebrate model.     

Haplotype combination of SREBP-1c gene sequence variants is associated with growth traits in cattle

The aim of this study was to examine the association of the SREBP-1c polymorphism with growth traits in cattle breeds. Five sequence variants (SVs) were identified within the bovine sterol regulatory element-binding protein-1c gene (SREBP-1c), using DNA sequencing, PCR, PCR–RFLP, and forced PCR–RFLP methods. These polymorphisms include three missense mutations (SV1, SV4, and SV5) in exons 7, 9, and 12, a silent mutation (SV3) in exon 9, and a large deletion (SV2) in intron 7. Overall, we report the validation of polymorphisms within the bovine SREBP-1c gene, and the haplotype variability and extent of linkage disequilibrium (LD) in 1061 individuals representing the five main cattle breeds from China. We also investigated haplotype frequencies and LD coefficients for five SVs in all study populations. LD and haplotype structure of SREBP-1c were different between breeds. The result of haplotype analysis of five SVs showed that 27 different haplotypes were identified by all breeds. Two haplotypes (Hap1 and Hap2) shared by all five populations accounted for 42.75%, 35.68%, 36.44%, 25.43%, and 96.26% of all haplotypes observed in the cattle breeds Nanyang, Qinchuan, Jiaxian, Jinnan, and Chinese Holstein, respectively. The statistical analyses indicated that one single SV and 38 combined haplotypes were significantly associated with growth traits in the Nanyang cattle population (P < 0.05 or P < 0.01). The results of this study suggest that the SREBP-1c gene possibly is a strong candidate gene that affects growth traits in the Chinese beef cattle breeding program.     

Linkage disequilibrium mapping in domestic dog breeds narrows the progressive rod-cone degeneration interval and identifies ancestral disease-transmitting chromosome

Canine progressive rod-cone degeneration (prcd) is a retinal disease previously mapped to a broad, gene-rich centromeric region of canine chromosome 9. As allelic disorders are present in multiple breeds, we used linkage disequilibrium (LD) to narrow the approximately 6.4-Mb interval candidate region. Multiple dog breeds, each representing genetically isolated populations, were typed for SNPs and other polymorphisms identified from BACs. The candidate region was initially localized to a 1.5-Mb zero recombination interval between growth factor receptor-bound protein 2 (GRB2) and SEC14-like 1 (SEC14L). A fine-scale haplotype of the region was developed, which reduced the LD interval to 106 kb and identified a conserved haplotype of 98 polymorphisms present in all prcd-affected chromosomes from 14 different dog breeds. The findings strongly suggest that a common ancestor transmitted the prcd disease allele to many of the modern dog breeds and demonstrate the power of the LD approach in the canine model.     

Prospects for whole genome linkage disequilibrium mapping in domestic dog breeds

Linkage disequilibrium (LD) mapping is commonly used as a fine mapping tool in human genome mapping and has been used with some success for initial disease gene isolation in certain isolated inbred human populations. An understanding of the population history of domestic dog breeds suggests that LD mapping could be routinely utilized in this species for initial genome-wide scans. Such an approach offers significant advantages over traditional linkage analysis. Here, we demonstrate, using canine copper toxicosis in the Bedlington terrier as the model, that LD mapping could be reasonably expected to be a useful strategy in low-resolution, genome-wide scans in pure-bred dogs. Significant LD was demonstrated over distances up to 33.3 cM. It is very unlikely, for a number of reasons discussed, that this result could be extrapolated to the rest of the genome. It is, however, consistent with the expectation given the population structure of canine breeds and, in this breed at least, with the hypothesis that it may be possible to utilize LD in a genome-wide scan. In this study, LD mapping confirmed the location of the copper toxicosis in Bedlington terrier gene (CT-BT) and was able to do so in a population that was refractory to traditional linkage analysis.     

Single-nucleotide polymorphisms and haplotype LD analysis of the 29-kb IGF2 region on chromosome 11p15.5 in the Korean population

OBJECTIVE: We investigated sequence variations of the 29-kb insulin-like growth factor 2 (IGF2) region in human chromosome region 11p15.5 in the Korean population. This region consists of IGF2, insulin-like growth factor 2 antisense (IGF2AS), and the insulin gene, all important candidate genes for various diseases, including cancer, obesity, diabetes, and coronary disease. While single nucleotide polymorphisms (SNPs) have been identified for this region and used in association studies, ethnic differences in genetic variation at this site have not been addressed. To date, SNPs for the entire 29-kb region in the Korean population have not been reported. METHODS: We surveyed a population of 108 Koreans for SNPs in the 29-kb IGF2 region. RESULTS: We identified 62 SNPs, consisting of 6 SNPs in the promoter region, 17 in the untranslated region, 19 in introns, and 20 in the intergenic region. We also analyzed linkage disequilibrium (LD) patterns and haplotypes using 36 high-frequency (> 5%)SNPs and found a well-defined LD block spanning about 13 kb that includes 8 kb of the IGF2AS gene, with two hot-spot regions flanking the LD block. CONCLUSION: These SNPs may be useful as genetic markers in disease association studies in the Korean population.     

Intriguing balancing selection on the intron 5 region of LMBR1 in human population

Background

The intron 5 of gene LMBR1 is the cis-acting regulatory module for the sonic hedgehog (SHH) gene. Mutation in this non-coding region is associated with preaxial polydactyly, and may play crucial roles in the evolution of limb and skeletal system.

Methodology/Principal Findings

We sequenced a region of the LMBR1 gene intron 5 in East Asian human population, and found a significant deviation of Tajima''s D statistics from neutrality taking human population growth into account. Data from HapMap also demonstrated extended linkage disequilibrium in the region in East Asian and European population, and significantly low degree of genetic differentiation among human populations.

Conclusion/Significance

We proposed that the intron 5 of LMBR1 was presumably subject to balancing selection during the evolution of modern human.     

Haplotypes within genes of β-chemokines in 17q11 are associated with multiple sclerosis: a second phase study

We previously defined haplotypes of single nucleotide polymorphisms (SNP) with possible relevance to multiple sclerosis (MS) in 2 CC chemokine ligand (CCL) clusters in chromosome 17q11. The 17q11 region was also identified as a susceptibility locus by a meta-analysis of linkage studies. To confirm and refine the previous finding in a second, high resolution SNP scan in a new set of families. We genotyped 232 SNPs in 1369 individuals in 361 MS families. Transmission of marker alleles and haplotypes from unaffected parents to affected offspring was tested by using the pedigree disequilibrium test, the TRANSMIT 2.5 program, and the family and haplotype based association tests. Distribution of linkage disequilibrium (LD) was assessed by ldmax. In consensus with observations in the first scan, the present study identified haplotypes within CCL3 and CCL15 in the telomeric CCL cluster. There was also an overlap in the findings in the centromeric CCL cluster. Strong and extensive LD was detected both within the centromeric and telomeric CCL gene clusters. The present study replicates our previous findings and further suggests the existence of MS associated haplotypes within genes of CCL3 and CCL15. Haplotypes of interest are also present within the centromeric gene cluster (including CCL2, CCL7, CCL11, CCL8, and CCL13), but extensive LD prevents further refinement of these haplotypes by using the methods applied. Sequencing of the identified chromosomal segments and their flanking regions will be necessary to define specific variants with direct relevance to MS pathogenesis.     

Linkage disequilibrium across six prion gene regions spanning 20 kbp in U.S. sheep

Single nucleotide polymorphisms (SNPs) and haplotype alleles within the prion gene (PRNP) coding sequence of domestic sheep (Ovis aries) are associated with genetic predisposition to scrapie, a transmissible spongiform encephalopathy disease of sheep. This report describes regions of linkage disequilibrium (LD) throughout the PRNP gene region in U.S. sheep and provides a genetic framework for identifying additional PRNP determinants associated with scrapie resistance. Four sequence tagged sites (i.e., STS or amplicons) totaling 3869 bp and spanning 20 kbp of genomic PRNP sequence were sequenced in a diverse panel of 90 sires representing ten popular U.S. breeds of sheep. Analysis of these sequences identified 36 previously unreported polymorphisms. In combination with two previously characterized STS, 62 polymorphisms were analyzed in a 20-kbp PRNP region in this panel of U.S. sheep. Two regions of strong LD and ten common haplotypes were identified. The haplotype encoding amino acid residues A, R, and Q at codons 136, 154, and 171, respectively, was observed on nine larger haplotypes spanning PRNP from the promoter region to the 3′ untranslated region. The haplotype encoding VRQ was observed on two larger haplotypes, whereas ARR, ARH, and AHQ were each present on a single haplotype. The existence of multiple haplotypes encoding ARQ raises the question of whether sheep bearing these different haplotypes are equally susceptible to scrapie. The haplotype structure within the 20-kbp region of PRNP identified in this study is important for higher-resolution analysis of genetics contributions to scrapie susceptibility. Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the accession number DQ077504.     

Assessment of the functionality of genome‐wide canine SNP arrays and implications for canine disease association studies

Domestic dogs share a wide range of important disease conditions with humans, including cancers, diabetes and epilepsy. Many of these conditions have similar or identical underlying pathologies to their human counterparts and thus dogs represent physiologically relevant natural models of human disorders. Comparative genomic approaches whereby disease genes can be identified in dog diseases and then mapped onto the human genome are now recognized as a valid method and are increasing in popularity. The majority of dog breeds have been created over the past few hundred years and, as a consequence, the dog genome is characterized by extensive linkage disequilibrium (LD), extending usually from hundreds of kilobases to several megabases within a breed, rather than tens of kilobases observed in the human genome. Genome‐wide canine SNP arrays have been developed, and increasing success of using these arrays to map disease loci in dogs is emerging. No equivalent of the human HapMap currently exists for different canine breeds, and the LD structure for such breeds is far less understood than for humans. This study is a dedicated large‐scale assessment of the functionalities (LD and SNP tagging performance) of canine genome‐wide SNP arrays in multiple domestic dog breeds. We have used genotype data from 18 breeds as well as wolves and coyotes genotyped by the Illumina 22K canine SNP array and Affymetrix 50K canine SNP array. As expected, high tagging performance was observed with most of the breeds using both Illumina and Affymetrix arrays when multi‐marker tagging was applied. In contrast, however, large differences in population structure, LD coverage and pairwise tagging performance were found between breeds, suggesting that study designs should be carefully assessed for individual breeds before undertaking genome‐wide association studies (GWAS).     

Nucleotide variability and linkage disequilibrium patterns in the porcine MUC4 gene

ABSTRACT: BACKGROUND: MUC4 is a type of membrane anchored glycoprotein and serves as the major constituent of mucus that covers epithelial surfaces of many tissues such as trachea, colon and cervix. MUC4 plays important roles in the lubrication and protection of the surface epithelium, cell proliferation and differentiation, immune response, cell adhesion and cancer development. To gain insights into the evolution of the porcine MUC4 gene, we surveyed the nucleotide variability and linkage disequilibrium (LD) within this gene in Chinese indigenous breeds and Western commercial breeds. RESULTS: A total of 53 SNPs covering the MUC4 gene were genotyped on 5 wild boars and 307 domestic pigs representing 11 Chinese breeds and 3 Western breeds. The nucleotide variability, haplotype phylogeny and LD extent of MUC4 were analyzed in these breeds. Both Chinese and Western breeds had considerable nucleotide diversity at the MUC4 locus. Western pig breeds like Duroc and Large White have comparable nucleotide diversity as many of Chinese breeds, thus artificial selection for lean pork production have not reduced the genetic variability of MUC4 in Western commercial breeds. Haplotype phylogeny analyses indicated that MUC4 had evolved divergently in Chinese and Western pigs. The dendrogram of genetic differentiation between breeds generally reflected demographic history and geographical distribution of these breeds. LD patterns were unexpectedly similar between Chinese and Western breeds, in which LD usually extended less than 20 kb. This is different from the presumed high LD extent (more than 100 kb) in Western commercial breeds. The significant positive Tajima'D, and Fu and Li's D statistics in a few Chinese and Western breeds implied that MUC4 might undergo balancing selection in domestic breeds. Nevertheless, we cautioned that the significant statistics could be upward biased by SNP ascertainment process. CONCLUSIONS: Chinese and Western breeds have similar nucleotide diversity but evolve divergently in the MUC4 region. Western breeds exhibited unusual low LD extent at the MUC4 locus, reflecting the complexity of nucleotide variability of pig genome. The finding suggests that high density (e.g. 1SNP/10 kb) markers are required to capture the underlying causal variants at such regions.     

How to make a zone of polarizing activity: insights into limb development via the abnormality preaxial polydactyly

Early in vertebrate limb development, a program initiates that polarizes the limb along the antero-posterior axis. The mesenchyme at the posterior margin is ultimately responsible for the asymmetry due to a region called the zone of polarizing activity (ZPA). The ZPA produces and secretes the molecule SHH, which coordinates the patterning of the resulting digits. Preaxial polydactyly (PPD) is a commonly occurring limb abnormality; investigating the genetic basis of this defect has provided insights into our understanding of digit patterning. PPD disrupts limb asymmetry by producing an ectopic ZPA at the opposite margin of the limb bud. Mutations in the long-range, limb-specific regulatory element of the Shh gene are responsible for the defect. Genetic analysis of this limb abnormality provides an important approach in understanding the mechanisms that control digit patterning.     

Large-scale SNP genotyping with canine buccal swab DNA

The dog is an attractive model for genetic studies of complex disease. With drafts of the canine genome complete, a large number of single-nucleotide polymorphisms (SNPs) that are potentially useful for gene-mapping studies and empirical estimations of canine diversity and linkage disequilibrium (LD) are now available. Unfortunately, most canine SNPs remain uncharacterized, and the amount and quality of DNA available from population-based samples are limited. We assessed how these real-world challenges influence automated SNP genotyping methods such as Illumina's GoldenGate assay. We examined 384 SNPs on canine chromosome 9 and successfully genotyped a minimum of 217 and a maximum of 275 SNPs using buccal swab samples for 181 dogs (86 beagles, 76 border collies, and 15 Australian shepherds). Call rates per SNP and sample averaged 97%, with reproducibility within and between analyses averaging 98%. The majority of these SNPs were polymorphic across all 3 breeds. We observed extensive LD, albeit less than reported for surveys using fewer dogs, consistent between breeds. Analyses of population substructure indicated that beagles are distinct from border collies and Australian shepherds. These results demonstrate the suitability of amplified canine buccal samples for high-throughput multiplex genotyping and confirm extensive LD in the dog.     

Combined Haplotypes of CaSR Gene Sequence Variants and Their Associations with Growth Traits in Cattle

The calcium-sensing receptor (CaSR) is a Class C G-protein coupled receptor that regulates food intake and assimilation. However, studies on the relationship between CaSR gene and growth traits in cattle are deficient. The aim of this study was to examine the association of the CaSR polymorphism with growth traits in cattle breeds. Four novel single nucleotide polymorphisms (SNPs) and one previously reported SNP (NC_007299.5: g.67630865T>C, 67638409G>C, 67660395G>C, 67661546C>G, and 67661892A>C) were identified in the bovine CaSR gene using DNA sequencing and PCR-SSCP methods in 520 individuals from three representative breeds. The three SNP P4_2, P7_1, and P7_4 in LX, QC, and JX cattle populations belonged to intermediate genetic diversity (0.25?相似文献   

Identification of genetic variation and haplotype structure of the canine ABCA4 gene for retinal disease association studies

Over 200 mutations in the retina specific member of the ATP-binding cassette transporter superfamily (ABCA4) have been associated with a diverse group of human retinal diseases. The disease mechanisms, and genotype–phenotype associations, nonetheless, remain elusive in many cases. As orthologous genes are commonly mutated in canine models of human blinding disorders, canine ABCA4 appears to be an ideal candidate gene to identify and study sequence changes in dogs affected by various forms of inherited retinal degeneration. However, the size of the gene and lack of haplotype assignment significantly limit targeted association and/or linkage approaches. This study assessed the naturally observed sequence diversity of ABCA4 in the dog, identifying 80% of novel variations. While none of the observed polymorphisms have been associated with blinding disorders to date, breed and potentially disease specific haplotypes have been identified. Moreover, a tag SNP map of 17 (15) markers has been established that accurately predicts common ABCA4 haplotypes (frequency > 5%) explaining >85% (>80%) of the observed genetic diversity and will considerably advance future studies. Our sequence analysis of the complete canine ABCA4 coding region will clearly provide a baseline and tools for future association studies and comparative genomics to further delineate the role of ABCA4 in canine blinding disorders.