Genotype matrix mapping: searching for quantitative trait loci interactions in genetic variation in complex traits.

In order to reveal quantitative trait loci (QTL) interactions and the relationship between various interactions in complex traits, we have developed a new QTL mapping approach, named genotype matrix mapping (GMM), which searches for QTL interactions in genetic variation. The central approach in GMM is the following. (1) Each tested marker is given a virtual matrix, named a genotype matrix (GM), containing intersecting lines and rows equal to the total allele number for that marker in the population analyzed. (2) QTL interactions are then estimated and compared through virtual networks among the GMs. To evaluate the contribution of marker combinations to a quantitative phenotype, the GMM method divides the samples into two non-overlapping subclasses, S(0) and S(1); the former contains the samples that have a specific genotype pattern to be evaluated, and the latter contains samples that do not. Based on this division, the F-measure is calculated as an index of significance. With the GMM method, we extracted significant marker combinations consisting of one to three interacting markers. The results indicated there were multiple QTL interactions affecting the phenotype (flowering date). GMM will be a valuable approach to identify QTL interactions in genetic variation of a complex trait within a variety of organisms.     

Microsatellite loci for genetic mapping in the turkey (Meleagris gallopavo)

New microsatellite loci for the turkey (Meleagris gallopavo) were developed from two small insert DNA libraries. Polymorphism at these new loci was examined in domestic birds and two resource populations designed for genetic linkage mapping. The majority of loci (152 of 168) was polymorphic in domestic turkeys and informative in two mapping resource populations and thus will be useful for genetic linkage mapping.     

Use of Mutant-Assisted Gene Identification and Characterization (MAGIC) to identify novel genetic loci that modify the maize hypersensitive response

The partially dominant, autoactive maize disease resistance gene Rp1-D21 causes hypersensitive response (HR) lesions to form spontaneously on leaves and stems in the absence of pathogen recognition. The maize nested association mapping (NAM) population consists of 25 200-line subpopulations each derived from a cross between the maize line B73 and one of 25 diverse inbred lines. By crossing a line carrying the Rp1-D21 gene with lines from three of these subpopulations and assessing the F(1) progeny, we were able to map several novel loci that modify the maize HR, using both single-population quantitative trait locus (QTL) and joint analysis of all three populations. Joint analysis detected QTL in greater number and with greater confidence and precision than did single population analysis. In particular, QTL were detected in bins 1.02, 4.04, 9.03, and 10.03. We have previously termed this technique, in which a mutant phenotype is used as a "reporter" for a trait of interest, Mutant-Assisted Gene Identification and Characterization (MAGIC).     

High genetic diversity among fossorial lizard populations (Anniella pulchra) in a rapidly developing landscape (Central California)

California is a biodiversity hotspot facing unbridled human population growth, especially in Central California. One of the poorly known, sensitive species in this area is the California legless lizard (Anniella pulchra), a fossorial worm-like reptile. We report mt and nuDNA sequences from 69 museum-vouchered samples of Anniella (A. pulchra and its sister species A. geronimensis) from 48 localities. Our genetic survey reveals substantially more genetic diversity within A. pulchra than previously reported. Our two independently evolving markers (mt and nuDNA) reveal five major lineages of A. pulchra. Two of the five major lineages of A. pulchra correspond to a north-south split found in other widespread California reptiles. These northern and southern clades also correspond to a previous study showing variation in chromosomal number. Unlike most other Californian reptiles, A. pulchra has major genetic lineages that are endemic to Central California including two that are endemic to the San Joaquin Valley and Carrizo Plain. Although A. pulchra is threatened throughout its range, the distinct San Joaquin lineages are seriously imperiled by urban sprawl. Some of the localities for the newly recognized genetic lineages have already been destroyed by development.     

Allozyme variation and genetic relationships among species in the Carex willdenowii complex (Cyperaceae)

A taxonomic study by Naczi, Reznicek, and Ford (American Journal of Botany, 85, 434-447, 1998) has determined that three species (Carex willdenowii, C. basiantha, and C. superata) can be recognized within the C. willdenowii complex. To determine the amount of genetic divergence within and between these species, allozyme analyses were conducted on 14 populations distributed from Pennsylvania to eastern Texas. Seventeen loci were surveyed, 13 of which were polymorphic, with all populations being polymorphic at one or more loci. Interspecific genetic identities ranged from 0.560 (C. willdenowii and C. basiantha) to 0.807 (C. basiantha and C. superata). Alleles for the isozymes Aat-1, Dia-1, Idh-2, Mdh-2, Per-1, Pgm-1, and Pgm-2 served to distinguish C. willdenowii from C. basiantha and C. superata. Carex basiantha and C. superata were recognized by alleles of Mdh-2, Pgm-1, and Tpi-2. The genetic identities of populations within species were high and exceeded 0.957. A caespitose growth habit and perigynia in close proximity to the staminate flowers suggest adaptations for selfing and therefore low levels of heterozygosity. Paradoxically, the values for expected heterozygosities (Hexp) were always lower than those obtained by direct count (Hobs): F values were highly negative, indicating heterozygous excess. Disassortative mating and selection are discussed as possible mechanisms for maintaining heterozygous excess within populations.     

On the use of multifactor dimensionality reduction (MDR) and classification and regression tree (CART) to identify haplotype-haplotype interactions in genetic studies

Haplotype-based approaches may have greater power than single-locus analyses when the SNPs are in strong linkage disequilibrium with the risk locus. To overcome potential complexities owing to large numbers of haplotypes in genetic studies, we evaluated two data mining approaches, multifactor dimensionality reduction (MDR) and classification and regression tree (CART), with the concept of haplotypes considering their haplotype uncertainty to detect haplotype-haplotype (HH) interactions. In evaluation of performance for detecting HH interactions, MDR had higher power than CART, but MDR gave a slightly higher type I error. Additionally, we performed an HH interaction analysis with a publicly available dataset of Parkinson's disease and confirmed previous findings that the RET proto-oncogene is associated with the disease. In this study, we showed that using HH interaction analysis is possible to assist researchers in gaining more insight into identifying genetic risk factors for complex diseases.     

Diversity at two genetic loci associated with the major histocompatibility complex in the golden snub-nosed monkey (Rhinopithecus roxellana)

Habitat loss via human activity has fragmented populations of the golden snub-nosed monkey (Rhinopithecus roxellana), and thus affected patterns of gene flow. We investigated in-depth a single troop in the Qinling Mountains, central China, two major histocompatibility complex (MHC) class II loci, DQA1 and DQB1, and compared the resulting data with data from troops from the wider Qinling Mountains region. We found that a novel DQB1 allele was only present in the study troop and relatively few divergent alleles at the DQA1 and DQB1 loci were present compared with the wider population. The inbreeding coefficient (F_is) at the MHC region was lower than previous measurements, which may have reflected different sampling strategies. However, R. roxellana has relatively high diversity in MHC genes, even though it has probably experienced serious past population bottlenecks and reduced gene flow between populations. We also found that some alleles present in the wider population had been lost in the study troop, and suggest that conservation management strategies be implemented to increase gene flow between troops in order to increase genetic variation.     

GM1-gangliosidosis (genetic beta-galactosidase deficiency): identification of four mutations in different clinical phenotypes among Japanese patients

GM1-gangliosidosis is a genetic neurological disorder caused by mutations in the lysosomal acid beta-galactosidase gene. While its phenotypic expression is complex, it is usually classified as being of infantile, juvenile, or adult form, on the basis of age at onset, the rate of symptomatic progression, and severity of central nervous system involvement. We have analyzed the acid beta-galactosidase gene in 12 Japanese patients from nine families. The aim was to identify mutations in individual patients and then to examine possible correlation between the mutations and the clinical phenotypes. Northern blotting studies with a full-length human beta-galactosidase cDNA showed that the mRNA ranged from undetectable to substantially decreased in the infantile patients but was normal in quantity and size in all juvenile and adult patients. Four distinct missense mutations have been identified, each limited to the respective clinical forms within our small-size samples. In the infantile patient with decreased but detectable mRNA, a point mutation was found resulting in Arg49----Cys. In the infantile patient with nearly undetectable mRNA, mutation Arg457----Ter was identified. The mutation Arg201----Cys was found in all four of the juvenile patients, while all six adult patients were homozygous for the point mutation Ile51----Thr. The mutations found in the juvenile and adult patients alter restriction sites in the normal gene and thus are amendable to quick screening. The prediction that these mutations are responsible for the clinical disease was confirmed by no expression of the catalytic activity of the mutant proteins in the COS-I cell expression system.(ABSTRACT TRUNCATED AT 250 WORDS)     

Thirty-three microsatellite loci for noninvasive genetic studies of the giant panda (Ailuropoda melanoleuca)

Limited microsatellite markers useable in noninvasive genetic methods have hampered the studies of dispersal patterns and mating systems of giant pandas. Therefore, we describe in this paper the characterization of 15 novel microsatellite loci from genomic DNA-enriched libraries and 18 redesigned microsatellite loci from published papers on the giant panda. The number of alleles per locus in 60 individuals ranged from 2 to 13, the average observed heterozygosity per locus from 0.168 to 0.800, and the average expected heterozygosity per locus from 0.152 to 0.882. All loci followed Hardy-Weinberg expectations. Four pairs of significant linkage association were found among all these loci. Moreover, the 33 microsatellite loci showed high amplification successes rate in noninvasive samples, which indicated that these loci will be of use in studying dispersal patterns and mating systems of giant pandas using noninvasive genetic methods.     

Spatial and genetic structure of a Melittis sarmatica Klok. (Lamiaceae) metapopulation located in the Volkovysk Upland (Belarus)

The spatial and genetic structure of the metapopulation of Melittis sarmatica Klok. from the Volkovysk Upland (Grodno region, Belarus) has been studied. The use of the RAPD method for genetic analysis of the similarity of populations revealed the exchange of genetic information between metapopulation components. We also revealed the most stable connections within the metapopulation, the key and problem elements of its structure, and some genetically separate groups. A negative linear dependence between the level of genetic similarity of populations and the distance between them was observed.     

Genomic and genetic relationships among species of Leymus (Poaceae: Triticeae) inferred from 18S-26S ribosomal genes

The 18S-26S ribosomal genes in three closely related species of Leymus (Poaceae: Triticeae) were examined using fluorescence in situ hybridization (FISH) and restriction fragment length polymorphism (RFLP). Both approaches revealed a close relationship between L. arenarius (8x = 56, northern European) and L. racemosus (4x = 28, central Eurasian), whereas L. mollis (4x = 28, northern American/Pacific) was distinct. Each species had three homologous pairs of major rDNA loci: a1, a2, and a3 for L. arenarius; m1, m2, and m3 for L. mollis; and r1, r2, and r3 for L. racemosus. Leymus arenarius had in addition three minor loci, a4, a5, and a6. The major loci of L. arenarius and L. racemosus were identical, indicating that the former species could have originated from the latter, via interspecific hybridization and/or polyploidy. The rDNA-RFLPs further indicated relationships of these species to other species of Leymus (L. karellini, 8x = 56 and L. angustus, 12x = 84) and Psathyrostachys (P. fragilis, P. huashanica, P. juncea, and P. lanuginosa, which are all diploids). A phenogram constructed from 20 BamHI, EcoRI, and DraI rDNA fragments revealed closer relationship between the two genera, Leymus and Psathyrostachys, than that among species within a genus.     

Phylogenetic relationships and genetic diversity among members of theFestuca-Lolium complex (Poaceae) based on ITS sequence data

There have been few DNA sequencebased studies of phylogenetic relationships within theFestuca-Lolium complex. Here we infer the phylogeny of 31Festuca-Lolium taxa with a dataset of 116 ITS sequences. The results are consistent with previous studies that resolved two majorFestuca clades: one clade of fine fescues and another clade that containsLolium and associatedFestuca species. This study is unique in suggesting a third, basalFestuca clade, but support for the basal position of this group is low. Extensive sampling permitted investigation of the effects of lineage sorting and reticulate events on the evolution of the complex. Roughly half of the taxa show evidence of lineage sorting or reticulation, and the monophyly ofLolium has likely been obscured by reticulate events. Overall, polyploid species harbor higher levels of ITS sequence diversity than diploids; ITS sequence variants may provide clues to the identity of allopolyploid parents.     

Polymorphism Interaction Analysis (PIA): a method for investigating complex gene-gene interactions

Background  

The risk of common diseases is likely determined by the complex interplay between environmental and genetic factors, including single nucleotide polymorphisms (SNPs). Traditional methods of data analysis are poorly suited for detecting complex interactions due to sparseness of data in high dimensions, which often occurs when data are available for a large number of SNPs for a relatively small number of samples. Validation of associations observed using multiple methods should be implemented to minimize likelihood of false-positive associations. Moreover, high-throughput genotyping methods allow investigators to genotype thousands of SNPs at one time. Investigating associations for each individual SNP or interactions between SNPs using traditional approaches is inefficient and prone to false positives.     

Signatures of adaptation and genetic structure among the mainland populations of Pinus radiata (D. Don) inferred from SNP loci

Insights into the relative contributions of locus specific and genome-wide effects on population genetic diversity can be gained through separation of their resulting genetic signals. Here we explore patterns of adaptive and neutral genetic diversity in the disjunct natural populations of Pinus radiata (D. Don) from mainland California. A first-generation common garden of 447 individuals revealed significant differentiation of wood phenotypes among populations (P _ST), possibly reflecting local adaptation in response to environment. We subsequently screened all trees for genetic diversity at 149 candidate gene single nucleotide polymorphism (SNP) loci for signatures of adaptation. Ten loci were identified as being possible targets of diversifying selection following F _ST outlier tests. Multivariate canonical correlation performed on a data set of 444 individuals identified significant covariance between environment, adaptive phenotypes and outlier SNP diversity, lending support to the case for local adaptation suggested from F _ST and P _ST tests. Covariation among discrete sets of outlier SNPs and adaptive phenotypes (inferred from multivariate loadings) with environment are supported by existing studies of candidate gene function and genotype–phenotype association. Canonical analyses failed to detect significant correlations between environment and 139 non-outlier SNP loci, which were applied to estimate neutral patterns of genetic differentiation among populations (F _ST 4.3 %). Using this data set, significant hierarchical structure was detected, indicating three populations on the mainland. The hierarchical relationships based on neutral SNP markers (and SSR) were in contrast with those inferred from putatively adaptive loci, potentially highlighting the independent action of selection and demography in shaping genetic structure in this species.     

Multiple loci and epistases control genetic variation for seed dormancy in weedy rice (Oryza sativa)

Weedy rice has much stronger seed dormancy than cultivated rice. A wild-like weedy strain SS18-2 was selected to investigate the genetic architecture underlying seed dormancy, a critical adaptive trait in plants. A framework genetic map covering the rice genome was constructed on the basis of 156 BC(1) [EM93-1 (nondormant breeding line)//EM93-1/SS18-2] individuals. The mapping population was replicated using a split-tiller technique to control and better estimate the environmental variation. Dormancy was determined by germination of seeds after 1, 11, and 21 days of after-ripening (DAR). Six dormancy QTL, designated as qSD(S)-4, -6, -7-1, -7-2, -8, and -12, were identified. The locus qSD(S)-7-1 was tightly linked to the red pericarp color gene Rc. A QTL x DAR interaction was detected for qSD(S)-12, the locus with the largest main effect at 1, 11, and 21 DAR (R(2) = 0.14, 0.24, and 0.20, respectively). Two, three, and four orders of epistases were detected with four, six, and six QTL, respectively. The higher-order epistases strongly suggest the presence of genetically complex networks in the regulation of variation for seed dormancy in natural populations and make it critical to select for a favorable combination of alleles at multiple loci in positional cloning of a target dormancy gene.