Search

Harvest index is a measure of success in partitioning assimilated photosynthate. An improvement of harvest index means an increase in the economic portion of the plant. Our objective was to identify genetic markers associated with harvest index traits using 203 O. sativa accessions. The phenotyping for 14 traits was conducted in both temperate (Arkansas) and subtropical (Texas) climates and the genotyping used 154 SSRs and an indel marker. Heading, plant height and weight, and panicle length had negative correlations, while seed set and grain weight/panicle had positive correlations with harvest index across both locations. Subsequent genetic diversity and population structure analyses identified five groups in this collection, which corresponded to their geographic origins. Model comparisons revealed that different dimensions of principal components analysis (PCA) affected harvest index traits for mapping accuracy, and kinship did not help. In total, 36 markers in Arkansas and 28 markers in Texas were identified to be significantly associated with harvest index traits. Seven and two markers were consistently associated with two or more harvest index correlated traits in Arkansas and Texas, respectively. Additionally, four markers were constitutively identified at both locations, while 32 and 24 markers were identified specifically in Arkansas and Texas, respectively. Allelic analysis of four constitutive markers demonstrated that allele 253 bp of RM431 had significantly greater effect on decreasing plant height, and 390 bp of RM24011 had the greatest effect on decreasing panicle length across both locations. Many of these identified markers are located either nearby or flanking the regions where the QTLs for harvest index have been reported. Thus, the results from this association mapping study complement and enrich the information from linkage-based QTL studies and will be the basis for improving harvest index directly and indirectly in rice.

The advent of next-generation sequencing (NGS) has revolutionized genomic and transcriptomic approaches to biology. These new sequencing tools are also valuable for the discovery, validation and assessment of genetic markers in populations. Here we review and discuss best practices for several NGS methods for genome-wide genetic marker development and genotyping that use restriction enzyme digestion of target genomes to reduce the complexity of the target. These new methods -- which include reduced-representation sequencing using reduced-representation libraries (RRLs) or complexity reduction of polymorphic sequences (CRoPS), restriction-site-associated DNA sequencing (RAD-seq) and low coverage genotyping -- are applicable to both model organisms with high-quality reference genome sequences and, excitingly, to non-model species with no existing genomic data.

Genome organization in the family Pinaceae is complex and largely unknown. The main purpose of the present study was to develop and physically map species-diagnostic and species-specific molecular markers in pine and spruce. Five RAPD (random amplified polymorphic DNA) and one ISSR (inter-simple sequence repeat) species-diagnostic or species-specific markers for Picea mariana, Picea rubens, Pinus strobus, or Pinus monticola were identified, cloned, and sequenced. In situ hybridization of these sequences to spruce and pine chromosomes showed the sequences to be present in high copy number and evenly distributed throughout the genome. The analysis of centromeric and telomeric regions revealed the absence of significant clustering of species-diagnostic and species-specific sequences in all the chromosomes of the four species studied. Both RAPD and ISSR markers showed similar patterns.

The resolution of the chicken consensus linkage map has been dramatically improved in this study by genotyping 12,945 single nucleotide polymorphisms (SNPs) on three existing mapping populations in chicken: the Wageningen (WU), East Lansing (EL), and Uppsala (UPP) mapping populations. As many as 8599 SNPs could be included, bringing the total number of markers in the current consensus linkage map to 9268. The total length of the sex average map is 3228 cM, considerably smaller than previous estimates using the WU and EL populations, reflecting the higher quality of the new map. The current map consists of 34 linkage groups and covers at least 29 of the 38 autosomes. Sex-specific analysis and comparisons of the maps based on the three individual populations showed prominent heterogeneity in recombination rates between populations, but no significant heterogeneity between sexes. The recombination rates in the F(1) Red Jungle fowl/White Leghorn males and females were significantly lower compared with those in the WU broiler population, consistent with a higher recombination rate in purebred domestic animals under strong artificial selection. The recombination rate varied considerably among chromosomes as well as along individual chromosomes. An analysis of the sequence composition at recombination hot and cold spots revealed a strong positive correlation between GC-rich sequences and high recombination rates. The GC-rich cohesin binding sites in particular stood out from other GC-rich sequences with a 3.4-fold higher density at recombination hot spots versus cold spots, suggesting a functional relationship between recombination frequency and cohesin binding.

In genome-wide association studies (GWAS), single-marker analysis is usually employed to identify the most significant single nucleotide polymorphisms (SNPs). The trend test has been proposed for analysis of case-control association. Three trend tests, optimal for the recessive, additive and dominant models respectively, are available. When the underlying genetic model is unknown, the maximum of the three trend test results (MAX) has been shown to be robust against genetic model misspecification. Since the asymptotic distribution of MAX depends on the allele frequency of the SNP, using the P-value of MAX for ranking may be different from using the MAX statistic. Calculating the P-value of MAX for 300,000 (300 K) or more SNPs is computationally intensive and the software and program to obtain the P-value of MAX are not widely available. On the other hand, the MAX statistic is very easy to calculate without complex computer programs. Thus, we study whether or not one could use the MAX statistic instead of its P-value to rank SNPs in GWAS. The approaches using the MAX and its P-value to rank SNPs are referred to as MAX-rank and P-rank. By applying MAX-rank and P-rank to simulated and four real datasets from GWAS, we found the ranks of SNPs with true association are very similar using both approaches. Thus, we recommend to use MAX-rank for genome-wide scans. After the top-ranked SNPs are identified, their P-values based on MAX can be calculated and compared with the significance level.

Artificial selection has created myriad breeds of domestic animals, each characterized by unique phenotypes pertaining to behavior, morphology, physiology, and disease. Most domestic animal populations share features with isolated founder populations, making them well suited for positional cloning. Genome sequences are now available for most domestic species, and with them a panoply of tools including high-density single-nucleotide polymorphism panels. As a result, domestic animal populations are becoming invaluable resources for studying the molecular architecture of complex traits and of adaptation. Here we review recent progress and issues in the positional identification of genes underlying complex traits in domestic animals. As many phenotypes studied in animals are quantitative, we focus on mapping, fine mapping, and cloning of quantitative trait loci.

A directed search for QTL affecting carcass traits was carried out in the region of growth differentiation factor 8 (GDF8, also known as myostatin) on ovine chromosome 2 in seven Texel-sired half-sib families totaling 927 progeny. Weights were recorded at birth, weaning, ultrasound scanning, and slaughter. Ultrasonic measures of LM cross-sectional dimensions and s.c. fat above the LM were made, with the same measurements made on the LM after slaughter. Following slaughter, linear measurements of carcass length and width were made on all carcasses, and legs and loins from 540 lambs were dissected. Genotyping was carried out using eight microsatellite markers from FCB128 to RM356 on OAR 2 and analyzed using Haley-Knott regression. There was no evidence for QTL for growth rates or linear carcass traits. There was some evidence for QTL affecting LM dimensions segregating in some sire families, although it was not consistent between ultrasound and carcass measures of the same traits. There was strong and consistent evidence for a QTL affecting muscle and fat traits in the leg that mapped between markers BM81124 and BULGE20 for the four sires that were heterozygous in this region, but not for the three sires that were homozygous. The size of the effect varied across the four sires, ranging from 0.5 to 0.9 of an adjusted SD for weight-adjusted leg muscle traits, and ranging from 0.6 to 1.2 of an adjusted SD for weight-adjusted leg fat traits. The clearest effect shown was for multivariate analysis combining all leg muscle and fat traits analyzed across sires, where the -log(10) probability was 14. Animals carrying the favorable haplotype had 3.3% more muscle and 9.9% less fat in the leg relative to animals carrying other haplotypes. There was evidence for a second peak in the region of marker TEXAN2 for one sire group. It seems that a QTL affecting muscle and fat traits exists within the New Zealand Texel population, and it maps to the region of GDF8 on OAR2.

To understand the genetic basis of pathways involved in the control of breathing, a large scale, high-throughput study using chromosomal substitution strains of rats is underway. Eight new consomic rat stains (SS-2(BN), SS-4(BN), SS-6(BN), SS-7(BN), SS-8(BN), SS-11(BN), SS-12(BN), SS-14(BN), SS-Y(BN)), containing one homozygous BN/NHsdMcwi (BN) chromosome on a background of SS/JrHsdMcwi (SS), were created by PhysGen (http://pga.mcw.edu) Program for Genomic Applications. Male and female rats were studied using standard plethysmography under control conditions and during acute hypoxia (inspired oxygen fraction = 0.12) and hypercapnia (inspired CO(2) fraction = 0.07). The rats were also studied during treadmill exercise. Both male and female BN rats had a significantly lower ventilatory response during 7% CO(2) compared with SS rats of the same gender. SS-6(BN) female rats had a significantly reduced ventilatory response, similar to BN rats due primarily to a reduced tidal volume. Male SS-6(BN) rats had a significantly reduced tidal volume response to hypercapnia but a slightly increased frequency response during hypercapnia. Gene(s) on the Y chromosome may play a role in this increased frequency response in the male rats because the SS-Y(BN) hypercapnic ventilatory response involves a significantly increased frequency response. Several chromosomal substitutions slightly altered the ventilatory responses to hypoxia and exercise. However, genes on chromosomes 6 and Y of those studied are of primary importance in aspects of ventilatory control currently studied.

Pooling DNA from subjects within a group and comparing the pooled DNA across groups for a dense map of DNA markers offers a solution to the conundrum that linkage is systematic but not powerful whereas allelic association is powerful but not systematic. We used DNA pooling to screen 66 markers on chromosome 22 in original and replication samples of children of high general cognitive ability (g) and controls of average g. Although none of these markers survived our three-stage screening design (original pooling, replication pooling, individual genotyping), the results of DNA pooling were largely confirmed by individual genotyping. We can therefore exclude associations of major effect size on chromosome 22 for g, a key variable for cognitive neuroscience research on learning and memory.

Myostatin belongs to the transforming growth factor-beta superfamily, and is expressed specifically in developing and mature skeletal muscle. Myostatin appears to act as a negative regulator of muscle development, since mice with targeted disruption of this gene display a large increase in muscle mass. In this study, the porcine myostatin gene was mapped to chromosome 15q2.3 by fluorescence in situ hybridization. Myostatin was also positioned within the chromosome 15 linkage group using both a polymorphism located in the second intron and an associated microsatellite. The development of highly polymorphic markers associated with myostatin will support population studies to identify alleles of this gene that affects muscle mass and/or fat deposition in swine.

Radiation hybrid mapping has become an established tool for building physical maps. It represents a powerful way of constructing YAC contigs and high-resolution maps for positional cloning experiments. Ideally, radiation hybrids should not only provide support for the true order of the markers, but also accurate estimates of the physical distances between them. Statistical analysis of radiation hybrids has proved difficult because of the number of parameters (representing the fragment retention probabilities) that must be estimated, and simplifying assumptions are needed to analyze large numbers of markers simultaneously. The ramifications of these assumptions for the calculation of physical distances are investigated. A simple two-locus model is presented to demonstrate that variation in marker retention can lead to distortions in the estimates of distance. Multilocus simulations show that, when marker retention is constant across the chromosome, good estimates of physical distance can be derived using simple models of retention. However, further simulations exploring variable retention schemes demonstrate that significant errors in the estimates of map distances can occur. Ways of minimizing these distortions are discussed.