Genetic Mapping and QTL Analysis of Growth-Related Traits in Pinctada fucata Using Restriction-Site Associated DNA Sequencing

The pearl oyster, Pinctada fucata (P. fucata), is one of the marine bivalves that is predominantly cultured for pearl production. To obtain more genetic information for breeding purposes, we constructed a high-density linkage map of P. fucata and identified quantitative trait loci (QTL) for growth-related traits. One F1 family, which included the two parents, 48 largest progeny and 50 smallest progeny, was sampled to construct a linkage map using restriction site-associated DNA sequencing (RAD-Seq). With low coverage data, 1956.53 million clean reads and 86,342 candidate RAD loci were generated. A total of 1373 segregating SNPs were used to construct a sex-average linkage map. This spanned 1091.81 centimorgans (cM), with 14 linkage groups and an average marker interval of 1.41 cM. The genetic linkage map coverage, Coa, was 97.24%. Thirty-nine QTL-peak loci, for seven growth-related traits, were identified using the single-marker analysis, nonparametric mapping Kruskal-Wallis (KW) test. Parameters included three for shell height, six for shell length, five for shell width, four for hinge length, 11 for total weight, eight for soft tissue weight and two for shell weight. The QTL peak loci for shell height, shell length and shell weight were all located in linkage group 6. The genotype frequencies of most QTL peak loci showed significant differences between the large subpopulation and the small subpopulation (P<0.05). These results highlight the effectiveness of RAD-Seq as a tool for generation of QTL-targeted and genome-wide marker data in the non-model animal, P. fucata, and its possible utility in marker-assisted selection (MAS).     

Ethnicity and human genetic linkage maps

Human genetic linkage maps are based on rates of recombination across the genome. These rates in humans vary by the sex of the parent from whom alleles are inherited, by chromosomal position, and by genomic features, such as GC content and repeat density. We have examined--for the first time, to our knowledge--racial/ethnic differences in genetic maps of humans. We constructed genetic maps based on 353 microsatellite markers in four racial/ethnic groups: whites, African Americans, Mexican Americans, and East Asians (Chinese and Japanese). These maps were generated using 9,291 subjects from 2,900 nuclear families who participated in the National Heart, Lung, and Blood Institute-funded Family Blood Pressure Program, the largest sample used for map construction to date. Although the maps for the different groups are generally similar, we did find regional and genomewide differences across ethnic groups, including a longer genomewide map for African Americans than for other populations. Some of this variation was explained by genotyping artifacts--namely, null alleles (i.e., alleles with null phenotypes) at a number of loci--and by ethnic differences in null-allele frequencies. In particular, null alleles appear to be the likely explanation for the excess map length in African Americans. We also found that nonrandom missing data biases map results. However, we found regions on chromosome 8p and telomeric segments with significant ethnic differences and a suggestive interval on chromosome 12q that were not due to genotype artifacts. The difference on chromosome 8p is likely due to a polymorphic inversion in the region. The results of our investigation have implications for inferences of possible genetic influences on human recombination as well as for future linkage studies, especially those involving populations of nonwhite ethnicity.     

Evaluating the results of genomewide linkage scans of complex traits by locus counting

The evaluation of results from primary genomewide linkage scans of complex human traits remains an area of importance and considerable debate. Apart from the usual assessment of statistical significance by use of asymptotic and empirical calculations, an additional means of evaluation--based on counting the number of distinct regions showing evidence of linkage--is possible. We have explored the characteristics of such a locus-counting method over a range of experimental conditions typically encountered during genomewide scans for complex trait loci. Under the null hypothesis, factors that have an impact on the informativeness of the data--such as map density, availability of parental data, and completeness of genotyping--are seen to markedly influence the number of regions of excess allele sharing and the empirically derived genomewide significance of the associated LOD score thresholds. In some circumstances, the expected number of regions is less than one-quarter of that predicted under the assumption of a dense map and complete extraction of inheritance information. We have applied this method to a previously analyzed data set--the Warren 2 genome scan for type 2-diabetes susceptibility--and demonstrate that more regions showing evidence for linkage were observed in the primary genome scan than would be expected by chance, across the whole range of LOD scores, even though no single linkage result achieved empirical genomewide statistical significance. Locus counting may be useful in assessing the results from genome scans for complex traits in general, especially because relatively few scans generate evidence for linkage reaching genomewide significance by dense-map criteria. By taking account of the effects of reduced data informativeness on the expected number of regions showing evidence for linkage, a more meaningful, and less conservative, evaluation of the results from such linkage studies is possible.     

Development of a linkage map and QTL scan for growth traits in North American bison

PCR protocols incorporating fluorescently labeled multiplexed primer combinations were developed to produce a linkage map for bison. Three hundred fifty eight microsatellite loci spanning all 29 autosomes were genotyped via 83 PCR multiplexes and nine individual amplifications. A total of 292 markers were integrated into an autosomal linkage map for bison. The sex averaged bison map (2,647 cM) was approximately 9% longer than the corresponding USDA MARC map, which covered 2,415 cM. Utilizing weaning, yearling and 17-month weights from two private bison herds, a QTL scan was conducted using the developed linkage map. LOD peaks suggestive of QTL were identified on chromosomes 2, 7, 15, and 24 for weaning weight, chromosomes 4, 14, and 15 for yearling weight and chromosomes 8, 14, and 25 for 17-month weight. Four of the identified chromosomes have conserved synteny with regions harboring growth QTL in cattle.     

Efficient strategies for genomic searching using the affected-pedigree-member method of linkage analysis.

The affected-pedigree-member (APM) method of linkage analysis is a nonparametric statistic that tests for nonrandom cosegregation of a disease and marker loci. The APM statistic is based on the observation that if a marker locus is near a disease-susceptibility locus, then affected individuals within a family should be more similar at the marker locus than is expected by chance. The APM statistic measures marker similarity in terms of identity by state (IBS) of marker alleles; that is, two alleles are IBS if they are the same, regardless of their ancestral origin. Since the APM statistic measures increased marker similarity, it makes no assumptions concerning how the disease is inherited; this can be an advantage when dealing with complex diseases for which the mode of inheritance is difficult to determine. We investigate here the power of the APM statistic to detect linkage in the context of a genomewide search. In such a search, the APM statistic is evaluated at a grid of markers. Then regions with high APM statistics are investigated more thoroughly by typing more markers in the region. Using simulated data, we investigate various search strategies and recommend an optimal search strategy that maximizes the power to detect linkage while minimizing the false-positive rate and number of markers. We determine an optimal series of three increasing cut-points and an independent criterion for significance.     

猪的雌雄连锁图谱的比较研究

用统计的方法,对以一个商品猪群为参考家系,采用163个微卫星标记和3个I-型分子标记（RYR1、PRKAG3、PIT1）构建的猪常染色体雌、雄连锁图谱的长度进行了比较。结果表明,常染色体的雌性连锁图谱的总长度为2625.9 cM,雄性连锁图谱的总长度为2259.7 cM,二者比率为1.16 ：1;除了1号和14号染色体以外,其余染色体的雌性连锁图谱的长度均比雄性连锁图谱长。1、3、5、6、7、8、10、11、12、13、14、16、17、18号染色体的雌雄连锁图谱的长度差异极显著（P<0.01）;9号染色体的雌雄连锁图谱的长度差异为显著（P<0.05）;2、4、12、15号染色体的雌雄图谱的长度差异不显著。 Abstract：The difference between the length of female- and male-linkage map, which was created with a reference pedigree based on a commercial porcine population and using 163 microsatellite markers as well as 3 type-I markers (RYR1, PRKAG3, PIT1), was statistic analyzed. The results showed that the total length of female linkage map of autosomes is 2625.9 cM and the total length of the male linkage map is 2259.7 cM; the ratio between the total length of the female- and male-linkage maps is 1.16 ：1; except for the chromosomes 1 and 14, the female linkage maps of the other chromosomes are longer than the male linkage maps. The difference between the length of female- and male-linkage maps of chromosomes 1, 3, 5, 6, 7, 8, 10, 11, 12, 13, 14, 16, 17 and 18 is very significant (P<0.01) and the difference of chromosome 9 is significant (P<0.05); but there is no significance on chromosomes 2, 4, 12 and 15.     

A genetic linkage map for hexaploid,cultivated oat (Avena sativa L.) based on an intraspecific cross 'Ogle/MAM17-5'

Genetic research and breeding of oat ( Avena sativa L.) would be aided by development of a genetic linkage map for a breeding population. Such a map could be used for localization of qualitative and quantitative trait loci, marker-assisted selection and other genetic analysis in an adapted, agronomically useful background. The objectives of this research were to develop a genetic linkage map of hexaploid cultivated oat, to identify homoeologous relationships of linkage groups, and to compare homologous linkage groups between this map and the previously published hexaploid oat map from the cross 'Kanota/Ogle' (KO). A total of 510 markers, including 172 restriction fragment length polymorphisms (RFLP), 324 amplified fragment length polymorphisms (AFLP) and 14 simple sequence repeats (SSR), were assessed on a recombinant inbred population of 152 F(5:6) lines derived from the cross, 'Ogle/MAM17-5' (OM). Twenty eight linkage groups of 5 cM or longer were formed using 476 of the markers, while 34 markers remained either unlinked or in small fragments less than 5 cM. The 28 linkage groups contained from 3 to 33 markers, and varied in size from 5.2 to 123.0 cM, representing a total map length of 1,396.7 cM. Three putative homoeologous groups (OM7, OM8 and OM18; OM2 and OM23; OM13 and OM16) were identified. Comparison with the published KO map indicated that nine OM linkage groups could be determined to be homologous to linkage groups in the KO map. Further comparison of the homologous linkage groups revealed that residual differences in genomic rearrangements existed between the two hexaploid oat populations. Some linkage groups were significantly extended compared with the KO map. Since the OM mapping population is segregating for a number of agronomically important traits, this genetic map will provide a useful tool for identification of qualitative and quantitative loci for these traits.     

Generalized T2 test for genome association studies

Recent progress in the development of single-nucleotide polymorphism (SNP) maps within genes and across the genome provides a valuable tool for fine-mapping and has led to the suggestion of genomewide association studies to search for susceptibility loci for complex traits. Test statistics for genome association studies that consider a single marker at a time, ignoring the linkage disequilibrium between markers, are inefficient. In this study, we present a generalized T2 statistic for association studies of complex traits, which can utilize multiple SNP markers simultaneously and considers the effects of multiple disease-susceptibility loci. This generalized T2 statistic is a corollary to that originally developed for multivariate analysis and has a close relationship to discriminant analysis and common measure of genetic distance. We evaluate the power of the generalized T2 statistic and show that power to be greater than or equal to those of the traditional chi2 test of association and a similar haplotype-test statistic. Finally, examples are given to evaluate the performance of the proposed T2 statistic for association studies using simulated and real data.     

The X chromosome in quantitative trait locus mapping

The X chromosome requires special treatment in the mapping of quantitative trait loci (QTL). However, most QTL mapping methods, and most computer programs for QTL mapping, have focused exclusively on autosomal loci. We describe a method for appropriate treatment of the X chromosome for QTL mapping in experimental crosses. We address the important issue of formulating the null hypothesis of no linkage appropriately. If the X chromosome is treated like an autosome, a sex difference in the phenotype can lead to spurious linkage on the X chromosome. Further, the number of degrees of freedom for the linkage test may be different for the X chromosome than for autosomes, and so an X chromosome-specific significance threshold is required. To address this issue, we propose a general procedure to obtain chromosome-specific significance thresholds that controls the genomewide false positive rate at the desired level. We apply our methods to data on gut length in a large intercross of mice carrying the Sox10Dom mutation, a model of Hirschsprung disease. We identified QTL contributing to variation in gut length on chromosomes 5 and 18. We found suggestive evidence of linkage to the X chromosome, which would be viewed as strong evidence of linkage if the X chromosome was treated as an autosome. Our methods have been implemented in the package R/qtl.     

Genomewide linkage analyses of bipolar disorder: a new sample of 250 pedigrees from the National Institute of Mental Health Genetics Initiative

We conducted genomewide linkage analyses on 1,152 individuals from 250 families segregating for bipolar disorder and related affective illnesses. These pedigrees were ascertained at 10 sites in the United States, through a proband with bipolar I affective disorder and a sibling with bipolar I or schizoaffective disorder, bipolar type. Uniform methods of ascertainment and assessment were used at all sites. A 9-cM screen was performed by use of 391 markers, with an average heterozygosity of 0.76. Multipoint, nonparametric linkage analyses were conducted in affected relative pairs. Additionally, simulation analyses were performed to determine genomewide significance levels for this study. Three hierarchical models of affection were analyzed. Significant evidence for linkage (genomewide P<.05) was found on chromosome 17q, with a peak maximum LOD score of 3.63, at the marker D17S928, and on chromosome 6q, with a peak maximum LOD score of 3.61, near the marker D6S1021. These loci met both standard and simulation-based criteria for genomewide significance. Suggestive evidence of linkage was observed in three other regions (genomewide P<.10), on chromosomes 2p, 3q, and 8q. This study, which is based on the largest linkage sample for bipolar disorder analyzed to date, indicates that several genes contribute to bipolar disorder.     

The PiGMaP consortium linkage map of the pig (Sus scrofa)

A linkage map of the porcine genome has been developed by segregation analysis of 239 genetic markers. Eighty-one of these markers correspond to known genes. Linkage groups have been assigned to all 18 autosomes plus the X Chromosome (Chr). As 69 of the markers on the linkage map have also been mapped physically (by others), there is significant integration of linkage and physical map data. Six informative markers failed to show linkage to these maps. As in other species, the genetic map of the heterogametic sex (male) was significantly shorter (16.5 Morgans) than the genetic map of the homogametic sex (female) (21.5 Morgans). The sex-averaged genetic map of the pig was estimated to be 18 Morgans in length. Mapping information for 61 Type I loci (genes) enhances the contribution of the pig gene map to comparative gene mapping. Because the linkage map incorporates both highly polymorphic Type II loci, predominantly microsatellites, and Type I loci, it will be useful both for large experiments to map quantitative trait loci and for the subsequent isolation of trait genes following a comparative and candidate gene approach.     

A Genomewide Screen for Generalized Vitiligo: Confirmation of AIS1 on Chromosome 1p31 and Evidence for Additional Susceptibility Loci

Generalized vitiligo is a common autoimmune disorder characterized by the development of white patches of skin and overlying hair due to loss of pigment-forming melanocytes from the involved areas. Family clustering of cases is not uncommon, in a pattern suggestive of multifactorial, polygenic inheritance, and there is strong association between vitiligo and other autoimmune diseases. To map genetic loci that confer susceptibility to generalized vitiligo and perhaps other autoimmune diseases, we performed a genomewide linkage scan in 71 white multiplex families with vitiligo from North America and the United Kingdom. Linkage was assessed by multipoint nonparametric linkage analyses. One linkage signal, AIS1, located at 1p31, met genomewide criteria for highly significant linkage (nonparametric LOD 5.56; P=.000000282), establishing its importance as a major vitiligo susceptibility locus. An additional seven signals, on chromosomes 1, 7, 8, 11, 19, and 22, met genomewide criteria for "suggestive linkage," and will thus be of particular importance for follow-up studies.     

Chromosome-level homeology in paleopolyploid soybean (Glycine max) revealed through integration of genetic and chromosome maps

Soybean has 20 chromosome pairs that are derived from at least two rounds of genomewide duplication or polyploidy events although, cytogenetically, soybean behaves like a diploid and has disomic inheritance for most loci. Genetically anchored genomic clones were used as probes for fluorescence in situ hybridization (FISH) to determine the level of postpolyploid chromosomal rearrangements and to integrate the genetic and physical maps to (1) assign linkage groups to specific chromosomes, (2) assess chromosomal structure, and (3) determine the distribution of recombination along the length of a chromosome. FISH mapping of seven putatively gene-rich BACs from linkage group L (chromosome 19) revealed that most of the genetic map correlates to the highly euchromatic long arm and that there is extensive homeology with another chromosome pair, although colinearity of some loci does appear to be disrupted. Moreover, mapping of BACs containing high-copy sequences revealed sequestration of high-copy repeats to the pericentromeric regions of this chromosome. Taken together, these data present a model of chromosome structure in a highly duplicated but diploidized eukaryote, soybean.     

Construction of a genetic linkage map in the wild plant Mimulus using RAPD and isozyme markers.

As a first step to mapping quantitative trait loci for mating system differences, a genetic linkage map was generated from an interspecific backcross between Mimulus guttatus and Mimulus platycalyx. The linkage map consists of 99 RAPD and two isozyme markers. Eighty-one of these markers were mapped to 15 linkage groups, spanning 1437 contiguous centiMorgans, and covering 58% of the estimated genome. The genome length of Mimulus is estimated at 2474 +/- 35 cM; bootstrapping indicates that only ca. 40 markers are needed to give an accurate estimate of genome length. Further statistical analyses indicate that many RAPD markers cannot be ordered with certainty and that uncertain linkage groups tend to map nonlinearly even under commonly used mapping functions. Strategies for speeding up the mapping process for a wild species and possible applications of a partial linkage map in evolutionary studies are discussed. Key words : linkage map, mating system, Mimulus, RAPD.     

Applying novel genome-wide linkage strategies to search for loci influencing type 2 diabetes and adult height in American Samoa

Type 2 diabetes mellitus (T2DM) is a common complex phenotype that by the year 2010 is predicted to affect 221 million people globally. In the present study we performed a genome-wide linkage scan using the allele-sharing statistic Sall implemented in Allegro and a novel two-dimensional genome-wide strategy implemented in Merloc that searches for pairwise interaction between genetic markers located on different chromosomes linked to T2DM. In addition, we used a robust score statistic from the newly developed QTL-ALL software to search for linkage to variation in adult height. The strategies were applied to a study sample consisting of 238 sib-pairs affected with T2DM from American Samoa. We did not detect any genome-wide significant susceptibility loci for T2DM. However, our two-dimensional linkage investigation detected several loci pairs of interest, including 11q22 and 21q21, 9q21 and 11q22, 1p22-p21 and 4p15, and 4p15 and 15q11-q14, with a two-loci maximum LOD score (MLS) greater than 2.00. Most detected individual loci have previously been identified as susceptibility loci for diabetes-related traits. Our two-dimensional linkage results may facilitate the selection of potential candidate genes and molecular pathways for further diabetes studies because these results, besides providing candidate loci, also demonstrate that polygenic effects may play an important role in T2DM. Linkage was detected (p value of 0.005) for variation in adult height on chromosome 9q31, which was reported previously in other populations. Our finding suggests that the 9q31 region may be a strong quantitative trait locus for adult height, which is likely to be of importance across populations.     

Mapping unexplored genomes: a genetic linkage map of the woody sonchus alliance (Asteraceae: Sonchinae) in the Macaronesian Islands

As an initial step to mapping quantitative trait loci for species differences and adaptive radiation of insular endemics in Macaronesia, a genetic linkage map was constructed from an intergeneric backcross between Lactucosonchus webbii and Sonchus radicatus, core members of the tree lettuces in the Macaronesian Islands. A total of 152 amplified fragment length polymorphism markers were mapped into 10 major and 3 minor linkage groups for a total map length of 644 cM with an average distance of 4.53 cM for the 10 major groups. The genetic linkage map length is considerably less than the estimated, and this may reflect incomplete genomic coverage in the current study or reduced recombination, which is a common feature of maps for hybrids of divergent taxa. Segregation distortion occurred in 34% of the mapped markers, and they were located primarily in 4 linkage groups. Segregation distortion in the current BC(1) intergeneric population is slightly lower than average (40%) for BC(1) interspecific populations. This level of segregation distortion implies that unlike what we normally assume no to few reproductive barriers, oceanic island plant taxa do exhibit some degree of postmating reproductive isolation.     

Toward a saturated linkage map in tomato based on isozymes and random cDNA sequences

A linkage map in tomato has been developed based on isozyme and random cDNA clones derived from mRNA. Interspecific backcross and F₂ populations of Lycopersicon esculentum and L. pennellii were employed in the linkage analysis. Allelic differences in cDNA markers were based on restriction fragment length polymorphisms detected through Southern analysis. A total of 57 unique cDNA clones have been analyzed. The majority of cDNA markers correspond to single loci and are dispersed throughtout the genome. Of those clones that hybridize to two or more loci, most show genetic independence (ie., they are unlinked). The combination of isozyme, cDNA and previously mapped DNA markers total 112 loci. It is estimated that approximately 92% of the genome can be monitored during segregation with these markers. Molecular maps, such as the one being constructed in tomato, may allow genetic and breeding experiments that previously were not possible.     

A genetic linkage map of the baboon (Papio hamadryas) genome based on human microsatellite polymorphisms

A first-generation genetic linkage map of the baboon (Papio hamadryas) genome was developed for use in biomedical and evolutionary genetics. Pedigreed baboons (n = 694) were selected from the breeding colony maintained by the Southwest Foundation for Biomedical Research. To facilitate comparison with the human genome, the baboon linkage map consists primarily of human microsatellite loci amplified using published human PCR primers. Genotypes for 325 human microsatellites and 6 novel baboon microsatellites were used in linkage analyses performed with the MultiMap expert system. The resulting sex-averaged meiotic recombination map covers all 20 baboon autosomes, with average spacing among loci of 7.2 cM. Direct comparison among homologous (orthologous) loci reveals that, for 7 human autosomes, locus order is conserved between humans and baboons. For the other 15 autosomes, one or more rearrangements distinguish the two genomes. The total centimorgan distances among homologous markers are 28.0% longer in the human genome than in the baboon, suggesting that rates of recombination may be higher in humans. This baboon linkage map is the first reported for any nonhuman primate species and creates opportunities for mapping quantitative trait loci in baboons, as well as for comparative evolutionary analyses of genome structure.     

Twenty-eight loci form a continuous linkage map of markers for human chromosome 1

We have used a combination of 30 serological, protein electromorphic, and DNA markers defining 28 loci to construct a linkage map of chromosome 1. These markers form a continuous linkage group of 320 cM in males and 608 cM in females; female genetic distances were on average twofold higher than those of males across the map. Among the DNA markers are 10 highly polymorphic markers reflecting loci that contain a variable number of tandem repeats, well distributed over the length of the chromosome, that will be highly efficient anchor points for application of this map to studies of human genetic disease.     

Linkage maps for 20 enzyme loci in Aedes triseriatus

Twenty enzyme loci were mapped on the three linkage groups of Aedes triseriatus using intraspecific and interspecific matings. Large numbers of single-pair forced matings were made among field-collected A. triseriatus. Parents with appropriate isozyme linkage genotypes were identified and the progeny analyzed using standard electrophoretic procedures. Interspecific data were obtained by performing single-pair forced matings between A. triseriatus and either A. hendersoni or A. brelandi and then backcrossing to one of the parental species. Interspecific recombination values were adjusted to compensate for reduced chiasmata (and crossovers) in progeny of interspecific crosses. Four loci--Aat2, Me, Idh 1, and Mpi-- were associated with sex on linkage group (LG)I. The LG I map was about 24% longer than the predicted length of 62 map units. Eleven loci--Gpi, Hk4, Odh, Est2, Pgm, Sod1, Gpd, Had, Aco2, Idh2, and Est5--were assigned to LG II and spanned approximately 60 map units. Five loci--Mdh2, Pgd, Aat1, Gapd, and Fum--were assigned to LG III, but exact positions and distances of loci were not definitely established. The linkage relationships of enzyme loci of A. triseriatus were compared to maps of five other Aedes species in four subgenera. Map differences indicated several major inversions and translocations that separated the subgenera. In addition, several linkage groups appeared to have been conserved during Aedes subgeneric divergence.