PCSK7 gene variation bridges atherogenic dyslipidemia with hepatic inflammation in NAFLD patients.

Dyslipidemia and altered iron metabolism are typical features of nonalcoholic fatty liver disease (NAFLD). Proprotein convertase subtilisin/kexin type 7 (PCSK7) gene variation has been associated with circulating lipids and liver damage during iron overload. The aim of this study was to examine the impact of the PCSK7 rs236918 variant on NAFLD-related traits in 1,801 individuals from the Liver Biopsy Cohort (LBC), 500,000 from the UK Biobank Cohort (UKBBC), and 4,580 from the Dallas Heart Study (DHS). The minor PCSK7 rs236918 C allele was associated with higher triglycerides, aminotransferases, and hepatic inflammation in the LBC (P < 0.05) and with hypercholesterolemia and liver disease in the UKBBC. In the DHS, PCSK7 missense variants were associated with circulating lipids. PCSK7 was expressed in hepatocytes and its hepatic expression correlated with that of lipogenic genes (P < 0.05). The rs236918 C allele was associated with upregulation of a new "intra-PCSK7" long noncoding RNA predicted to interact with the protein, higher hepatic and circulating PCSK7 protein (P < 0.01), which correlated with triglycerides (P = 0.04). In HepG2 cells, PCSK7 deletion reduced lipogenesis, fat accumulation, inflammation, transforming growth factor ß pathway activation, and fibrogenesis. In conclusion, PCSK7 gene variation is associated with dyslipidemia and more severe liver disease in high risk individuals, likely by modulating PCSK7 expression/activity.

Natural selection on TMPRSS6 associated with the blunted erythropoiesis and improved blood viscosity in Tibetan pigs.

Tibetan pigs, indigenous to Tibetan plateau, are well adapted to hypoxia. So far, there have been not any definitively described genes and functional sites responsible for hypoxia adaptation for the Tibetan pig. The whole genome-wide association studies in human suggested that genetic variations in TMPRSS6 was associated with hemoglobin concentration (HGB) and red cell counts (RBC). Here we conducted resequencing of the nearly entire genomic region (40.1 kb) of the candidate gene TMPRSS6 in 40 domestic pigs and 40 wild boars along continuous altitudes and identified 708 SNPs, in addition to an indel (CGTG/----) in the intron 10. We conduct the CGTG indel in 838 domestic pigs, both the CGTG deletion frequency and the pairwise r2 linkage disequilibrium showed an increase with elevated altitudes, suggesting that TMPRSS6 has been under Darwinian positive selection. As the conserved core sequence of hypoxia-response elements (HREs), the deletion of CGTG in Tibetan pigs decreased the expression levels of TMPRSS6 mRNA and protein in the liver revealed by real-time quantitative PCR and western blot, respectively. We compared domestic pigs and Tibetan pigs living continuous altitudes, found that the blood-related traits with the increase of altitude, however, the HGB did not increase with the elevation in Tibetan pigs. Genotype association analysis results dissected a genetic effect on reducing HGB by 13.25 g/L in Gongbo'gyamda Tibetan pigs, decreasing mean corpuscular volume (MCV) by 4.79 fl in Diqing Tibetan pigs. In conclusion, the CGTG deletion of TMPRSS6 resulted in lower HGB and smaller MCV, which could reflect a blunting erythropoiesis and improving blood viscosity as well as erythrocyte deformability. It remains to be determined whether a blunting of erythropoiesis for TMPRSS6 or others genetic effects are the physiological adaptations among Tibetan pigs.

QTL and QTL networks for cold tolerance at the reproductive stage detected using selective introgression in rice.

Low temperature stress is one of the major abiotic stresses limiting the productivity of Geng (japonica) rice grown the temperate regions as well as in tropical high lands worldwide. To develop rice varieties with improved cold tolerance (CT) at the reproductive stage, 84 BC2 CT introgression lines (ILs) were developed from five populations through backcross breeding. These CT ILs plus 310 random ILs from the same BC populations were used for dissecting genetic networks underlying CT in rice by detecting QTLs and functional genetic units (FGUs) contributing to CT. Seventeen major QTLs for CT were identified using five selective introgression populations and the method of segregation distortion. Of them, three QTLs were confirmed using the random populations and seven others locate in the regions with previously reported CT QTLs/genes. Using multi-locus probability tests and linkage disequilibrium (LD) analyses, 46 functional genetic units (FGUs) (37 single loci and 9 association groups or AGs) distributed in 37 bins (~20%) across the rice genome for CT were detected. Together, each of the CT loci (bins) was detected in 1.7 populations, including 18 loci detected in two or more populations. Putative genetic networks (multi-locus structures) underlying CT were constructed based on strong non-random associations between or among donor alleles at the unlinked CT loci/FGUs identified in the CT ILs, suggesting the presence of strong epistasis among the detected CT loci. Our results demonstrated the power and usefulness of using selective introgression for simultaneous improvement and genetic dissection of complex traits such as CT in rice.

A genome-wide association study of early-maturation traits in upland cotton based on the CottonSNP80K array.

Genome-wide association studies (GWASs) efficiently identify genetic loci controlling traits at a relatively high resolution. In this study, variations in major early-maturation traits, including seedling period (SP), bud period (BP), flower and boll period (FBP), and growth period (GP), of 169 upland cotton accessions were investigated, and a GWAS of early maturation was performed based on a CottonSNP80K array. A total of 49,650 high-quality single-nucleotide polymorphisms (SNPs) were screened, and 29 significant SNPs located on chromosomes A6, A7, A8, D1, D2, and D9, were repeatedly identified as associated with early-maturation traits, in at least two environments or two algorithms. Of these 29 significant SNPs, 1, 12, 11, and 5 were related to SP, BP, FBP, and GP, respectively. Six peak SNPs, TM47967, TM13732, TM20937, TM28428, TM50283, and TM72552, exhibited phenotypic contributions of approximately 10%, which could allow them to be used for marker-assisted selection. One of these, TM72552, as well as four other SNPs, TM72554, TM72555, TM72558, and TM72559, corresponded to the quantitative trait loci previously reported. In total, 274 candidate genes were identified from the genome sequences of upland cotton and were categorized based on their functional annotations. Finally, our studies identified Gh_D01G0340 and Gh_D01G0341 as potential candidate genes for improving cotton early maturity.

Minor diplotypes of FMO3 might protect children and adolescents from obesity and insulin resistance.

In order to assess whether flavin-containing monooxygenase-3 (FMO3) might be involved in early cardiovascular risk, we assessed adiposity and traditional metabolic variables in children/adolescents grouped according to their genotypes in two FMO3 exonic polymorphisms, rs2266782 (E158K) and rs2266780 (E308G), which are in linkage disequilibrium and have been associated with decreased FMO3 activity. Among 776 children/adolescents (10.8 ± 2.2 years) recruited from the general population (452) and from our obesity outpatient clinic (324), the 68 carrying either the 158K-308G/158K-308E or the 158K-308G/158K-308G diplotype had lower mean z-BMI and prevalence of obesity compared to their 708 peers carrying any of the other diplotypes (0.39 vs 0.80, p = 0.01; OR = 0.39[0.17-0.87], p = 0.018, respectively), and to the sub-sample of 303 children carrying the major diplotype (158E-308E/158E-308E) (0.39 vs 0.87, p = 0.008; OR = 0.35[0.16-0.81], p = 0.014, respectively). They also had lower z-BMI-adjusted lnHOMA-IR compared to all the other children (0.75 vs 0.97, p = 0.001) and those carrying the major diplotype, (0.75 vs 0.98, p = 0.03), as well as lower z-BMI-adjusted iln-triglycerides compared to all the other children (3.98 vs 4.17, p = 0.037). These associations provide the first evidence that FMO3 may be involved in early body weight, insulin sensitivity, and lipid regulation in humans.

Effects of the Ordering of Natural Selection and Population Regulation Mechanisms on Wright-Fisher Models.

We explore the effect of different mechanisms of natural selection on the evolution of populations for one- and two-locus systems. We compare the effect of viability and fecundity selection in the context of the Wright-Fisher model with selection under the assumption of multiplicative fitness. We show that these two modes of natural selection correspond to different orderings of the processes of population regulation and natural selection in the Wright-Fisher model. We find that under the Wright-Fisher model these two different orderings can affect the distribution of trajectories of haplotype frequencies evolving with genetic recombination. However, the difference in the distribution of trajectories is only appreciable when the population is in significant linkage disequilibrium. We find that as linkage disequilibrium decays the trajectories for the two different models rapidly become indistinguishable. We discuss the significance of these findings in terms of biological examples of viability and fecundity selection, and speculate that the effect may be significant when factors such as gene migration maintain a degree of linkage disequilibrium.

Genome-wide association analysis for heat tolerance at flowering detected a large set of genes involved in adaptation to thermal and other stresses.

Fertilization sensitivity to heat in rice is a major issue within climate change scenarios in the tropics. A panel of 167 indica landraces and improved varieties was phenotyped for spikelet sterility (SPKST) under 38°C during anthesis and for several secondary traits potentially affecting panicle micro-climate and thus the fertilization process. The panel was genotyped with an average density of one marker per 29 kb using genotyping by sequencing. Genome-wide association analyses (GWAS) were conducted using three methods based on single marker regression, haplotype regression and simultaneous fitting of all markers, respectively. Fourteen loci significantly associated with SPKST under at least two GWAS methods were detected. A large number of associations was also detected for the secondary traits. Analysis of co-localization of SPKST associated loci with QTLs detected in progenies of bi-parental crosses reported in the literature allowed to narrow -down the position of eight of those QTLs, including the most documented one, qHTSF4.1. Gene families underlying loci associated with SPKST corresponded to functions ranging from sensing abiotic stresses and regulating plant response, such as wall-associated kinases and heat shock proteins, to cell division and gametophyte development. Analysis of diversity at the vicinity of loci associated with SPKST within the rice three thousand genomes, revealed widespread distribution of the favourable alleles across O. sativa genetic groups. However, few accessions assembled the favourable alleles at all loci. Effective donors included the heat tolerant variety N22 and some Indian and Taiwanese varieties. These results provide a basis for breeding for heat tolerance during anthesis and for functional validation of major loci governing this trait.

Reducing GWAS Complexity.

Genome-wide association studies (GWAS) have revealed numerous genomic 'hits' associated with complex phenotypes. In most cases these hits, along with surrogate genetic variation as measure by numerous single nucleotide polymorphisms (SNPs) that are in linkage disequilibrium, are not in coding genes making assignment of functionality or causality intractable. Here we propose that fine-mapping along with the matching of risk SNPs at chromatin biofeatures lessen this complexity by reducing the number of candidate functional/causal SNPs. For example, we show here that only on average 2 SNPs per prostate cancer risk locus are likely candidates for functionality/causality; we further propose that this manageable number should be taken forward in mechanistic studies. The candidate SNPs can be looked up for each prostate cancer risk region in 2 recent publications in 2015 (1,2) from our groups.

Associative transcriptomics study dissects the genetic architecture of seed glucosinolate content in Brassica napus.

Breeding new varieties with low seed glucosinolate (GS) concentrations has long been a prime target in Brassica napus. In this study, a novel association mapping methodology termed 'associative transcriptomics' (AT) was applied to a panel of 101 B. napus lines to define genetic regions and also candidate genes controlling total seed GS contents. Over 100,000 informative single-nucleotide polymorphisms (SNPs) and gene expression markers (GEMs) were developed for AT analysis, which led to the identification of 10 SNP and 7 GEM association peaks. Within these peaks, 26 genes were inferred to be involved in GS biosynthesis. A weighted gene co-expression network analysis provided additional 40 candidate genes. The transcript abundance in leaves of two candidate genes, BnaA.GTR2a located on chromosome A2 and BnaC.HAG3b on C9, was correlated with seed GS content, explaining 18.8 and 16.8% of phenotypic variation, respectively. Resequencing of genomic regions revealed six new SNPs in BnaA.GTR2a and four insertions or deletions in BnaC.HAG3b. These deletion polymorphisms were then successfully converted into polymerase chain reaction-based diagnostic markers that can, due to high linkage disequilibrium observed in these regions of the genome, be used for marker-assisted breeding for low seed GS lines.

Admixture mapping of prostate cancer in African Americans participating in the North Carolina-Louisiana Prostate Cancer Project (PCaP).

BACKGROUND: Few genetic risk factors have been uncovered that contribute specifically to the racial disparity in prostate cancer (CaP) observed in African Americans (AA). With the advent of ancestry informative marker (AIM) single nucleotide polymorphism (SNP) panels and powerful genetic strategies such as mapping by admixture linkage disequilibrium (MALD) it is possible to discover genes that underlie ethnic variation in disease risk. METHODS: One thousand one hundred thirty AA CaP cases enrolled in the North Carolina-Louisiana Prostate Cancer Project (PCaP) were genotyped using a 1,509 AIM SNP panel. MALD was performed using ADMIXMAP to test for linkage between CaP risk and ancestry estimates at each AIM SNP. RESULTS: The largest increase of African ancestry was observed at marker rs12543473 (P = 0.0011), located on chromosome 8q24.21, and the greatest excess of European ancestry was observed at marker rs10768140 (P = 0.0004) at chromosome 11p13. CONCLUSIONS: The study confirmed the 8q24 risk loci and identified a novel genomic region on 11p13 that is associated with CaP risk. These findings should be replicated in larger AA populations and combined with fine mapping data to further refine the novel 11p13 CaP risk loci.

CACTA-like transposable element in ZmCCT attenuated photoperiod sensitivity and accelerated the postdomestication spread of maize.

The postdomestication adaptation of maize to longer days required reduced photoperiod sensitivity to optimize flowering time. We performed a genome-wide association study and confirmed that ZmCCT, encoding a CCT domain-containing protein, is associated with the photoperiod response. In early-flowering maize we detected a CACTA-like transposable element (TE) within the ZmCCT promoter that dramatically reduced flowering time. TE insertion likely occurred after domestication and was selected as maize adapted to temperate zones. This process resulted in a strong selective sweep within the TE-related block of linkage disequilibrium. Functional validations indicated that the TE represses ZmCCT expression to reduce photoperiod sensitivity, thus accelerating maize spread to long-day environments.

Graphics processing unit-accelerated quantitative trait Loci detection.

Mapping quantitative trait loci (QTL) using genetic marker information is a time-consuming analysis that has interested the mapping community in recent decades. The increasing amount of genetic marker data allows one to consider ever more precise QTL analyses while increasing the demand for computation. Part of the difficulty of detecting QTLs resides in finding appropriate critical values or threshold values, above which a QTL effect is considered significant. Different approaches exist to determine these thresholds, using either empirical methods or algebraic approximations. In this article, we present a new implementation of existing software, QTLMap, which takes advantage of the data parallel nature of the problem by offsetting heavy computations to a graphics processing unit (GPU). Developments on the GPU were implemented using Cuda technology. This new implementation performs up to 75 times faster than the previous multicore implementation, while maintaining the same results and level of precision (Double Precision) and computing both QTL values and thresholds. This speedup allows one to perform more complex analyses, such as linkage disequilibrium linkage analyses (LDLA) and multiQTL analyses, in a reasonable time frame.

Common exonic missense variants in the C2 domain of the human KIBRA protein modify lipid binding and cognitive performance.

The human KIBRA gene has been linked to human cognition through a lead intronic single-nucleotide polymorphism (SNP; rs17070145) that is associated with episodic memory performance and the risk to develop Alzheimer's disease. However, it remains unknown how this relates to the function of the KIBRA protein. Here, we identified two common missense SNPs (rs3822660G/T [M734I], rs3822659T/G [S735A]) in exon 15 of the human KIBRA gene to affect cognitive performance, and to be in almost complete linkage disequilibrium with rs17070145. The identified SNPs encode variants of the KIBRA C2 domain with distinct Ca(2+) dependent binding preferences for monophosphorylated phosphatidylinositols likely due to differences in the dynamics and folding of the lipid-binding pocket. Our results further implicate the KIBRA protein in higher brain function and provide direction to the cellular pathways involved.

Genome-wide associations for feed utilisation complex in primiparous Holstein-Friesian dairy cows from experimental research herds in four European countries.

Genome-wide association studies for difficult-to-measure traits are generally limited by the sample size with accurate phenotypic data. The objective of this study was to utilise data on primiparous Holstein­Friesian cows from experimental farms in Ireland, the United Kingdom, the Netherlands and Sweden to identify genomic regions associated with the feed utilisation complex: fat and protein corrected milk yield (FPCM), dry matter intake (DMI), body condition score (BCS) and live-weight (LW). Phenotypic data and 37 590 single nucleotide polymorphisms (SNPs) were available on up to 1629 animals. Genetic parameters of the traits were estimated using a linear animal model with pedigree information, and univariate genome-wide association analyses were undertaken using Bayesian stochastic search variable selection performed using Gibbs sampling. The variation in the phenotypes explained by the SNPs on each chromosome was related to the size of the chromosome and was relatively consistent for each trait with the possible exceptions of BTA4 for BCS, BTA7, BTA13, BTA14, BTA18 for LW and BTA27 for DMI. For LW, BCS, DMI and FPCM, 266, 178, 206 and 254 SNPs had a Bayes factor .3, respectively. Olfactory genes and genes involved in the sensory smell process were overrepresented in a 500 kbp window around the significant SNPs. Potential candidate genes were involved with functions linked to insulin, epidermal growth factor and tryptophan.

Isolation and characterization of cross-amplification microsatellite panels for species of Procapra (Bovidae; Antilopinae).

The three Procapra species, Tibetan gazelle (P. picticaudata), Mongolian gazelle (P. gutturosa) and Przewalski's gazelle (P. przewalskii) are endemic to Asia. Several intraspecific genetic issues have been studied with species-specific microsatellite loci in these Asian gazelles. However, cross-species microsatellite panels are absent, which inhibits comparative conservation and evolutionary studies of the Procapra. In this study, we isolated 20 cross-species microsatellite loci for Procapra from both related species and the genomic library of P. przewalskii. Fifty-three samples of the three gazelles were used to characterize the markers. Allele numbers ranged from three to 20, with a mean of 7.93 per locus. Observed heterozygosity (H(O)) averaged 0.680 and expected heterozygosity (H(E)) 0.767. The mean polymorphic information content (PIC) was 0.757 for P. picticaudata, 0.803 for P. gutturosa and 0.590 for P. przewalskii. Nine loci were significantly deviated from Hardy-Weinberg (H-W) equilibrium in the three species. Significant linkage disequilibrium was detected in four pairs of loci in P. przewalskii, five pairs in P. gutturosa and 51 pairs in P. picticaudata. Considering the abundance of published loci and their high success rates of cross-amplification, testing and utilization of loci from related species is efficient for wild species of Bovidae. The cross-species microsatellite loci we developed will facilitate further interspecies genetic studies in Procapra.

La genética de la depresión: ¿qué información aportan las nuevas aproximaciones metodológicas? / The genetics of depression: What information can new methodologic approaches provide?

El trastorno depresivo mayor forma parte del grupo de enfermedades denominadas genéticamente complejas en cuya base se encuentran involucrados una serie de genes de efecto menor o susceptibilidad cuya expresión podría estar modulada por un gran número de factores ambientales. Desde los primeros estudios clásicos de ligamiento hasta las nuevas metodologías de estudio de genoma completo se ha puesto de manifiesto la dificultad para comprender las bases genéticas de la depresión mayor. En muchos estudios se han identificado una o varias regiones génicas cuya variabilidad confiere un riesgo pequeño para desarrollar un trastorno depresivo; es decir, dichas variantes explicarían un porcentaje muy pequeño del componente genético total de la enfermedad en la población y, por tanto, poseerían un valor predictivo bajo. Aunque los resultados obtenidos hasta el momento no son concluyentes, las nuevas aproximaciones basadas en estudios de interacción gen-ambiente así como los análisis de vías biológicas (basados en los estudios GWAS) abren nuevas y prometedoras perspectivas en la investigación de las bases genéticas y biológicas de la depresión mayor (AU)

Major depressive disorder is a genetically complex disease involving several minor, or susceptibility, genes whose expression may be modulated by many environmental factors. From the classical early linkage studies to the complete genome-wide association (GWA) study methodologies, it is evident that it is difficult to understand the genetic bases of this mental disorder. Many studies have identified one or more gene regions whose variability confers a small risk for developing depressive disorder, explaining a small percentage of the total genetic component of disease with low predictive value. Although the results to date are inconclusive, new approaches based on gene-environment interactions and biological pathway analysis (based on GWA studies) open new and promising perspectives in the investigation of the genetic and biological basis of major depression (AU)

Invited review: Genomic analysis of data from physiological studies.

Physiology deals with the functions of living organisms and their systems, and its scientific endeavors can be viewed as having temporally occurred in 3 phases. The first phase of physiology studies focused on determining the functions of particular organs and tissues and their functional differences according to physiological status. The second phase of studies focused on characterizing differences in these functions according to the environment, or productivity. The third phase of studies focuses on determining the physiological causes of differences in productivity. Distinguishing cause from effect in physiological systems of inter-related processes is problematic, such that science has struggled to identify the root physiological mechanisms and their role in the network of genes leading to differences in productivity. Genomics is the study of the entire genome and provides powerful new tools that will accelerate third-phase discoveries of causal physiological processes. That research exploits information on DNA polymorphisms known as markers, complete DNA sequence, RNA sequence, and RNA expression in particular tissues at specific life stages. Physiologists can determine the genetic cause of mutant animals, identify genetic differences between cases and controls, and identify genes responsible for differences in performance between average and above-average animals. In some species, physiologists can leverage genomic data being used to predict genetic merit in elite seedstock populations, as a starting point to identify genes that will then motivate detailed physiological studies in the organs or tissues and stages of life in which those genes are expressed. Such work will increase our knowledge of biology and may lead to novel approaches to manipulate animal performance.

Genetic variation in odorant receptors contributes to variation in olfactory behavior in a natural population of Drosophila melanogaster.

Chemoreception is a principle modality by which organisms gain information from their environment, and extensive variation in odor-mediated behavior has been documented within and among species. To examine the mechanisms by which sensory systems mediate these responses, we ask to what extent variation in Drosophila melanogaster odorant receptor genes contributes to variation in odor-mediated behavior. Significant differences in behavioral responses to structurally similar odorants, methyl hexanoate and ethyl hexanoate, were found in a natural population. Polymorphisms in 3 genomic regions (Or22a/Or22b, Or35a, and Or47a) were identified and associated with variation in behavior to these esters. Overall similarity in association profiles for both odorants was observed, except for Or47a in which polymorphisms were associated solely with variation in responses to ethyl hexanoate. Our analyses were then extended to examine polymorphisms in 3 odorant receptors previously reported to contribute to variation in olfactory behavior for the chemically distinct odorants benzaldehyde and acetophenone. Two Or10a polymorphisms were associated with variation in response to ethyl hexanoate. Finally, differences in Or35a and Or47a expression were associated with variation in responses to ethyl hexanoate. These results demonstrate that the genetic variation at the peripheral sensory stage plays a role in mediating differences in odor-mediated behavior.

Genetic comparisons between seed bank and Stipa krylovii plant populations.

The soil seed bank represents the potential plant population since it is the source for population replacement. The genetic structure of a Stipa krylovii (Roshev.) plant population and its soil seed bank was investigated in the Xilinguole Steppe of Inner Mongolia using random amplified polymorphic DNA (RAPD) analyses. The population was sampled at two sites that were in close proximity to each other (0.5 km apart). Thirty plants and 18 seed bank samples were taken from each site to determine the genetic diversity between sites and between sources (plant or seed). The material was analyzed using 13 primers to produce 92 loci. Eighty-six were multi-loci, of which 23 loci (26.74%) of allele frequencies showed significant differences (P < or = 0.05). The genetic similarity between two seed bank sites was 0.9843 while the genetic similarity between two plant sites was 0.9619. Their similarities were all greater than that between the seed bank and plant populations. An analysis of their genetic structure showed that 87.86% of total variation was derived by two-loci. Genetic structures between plant and soil seed bank populations in S. krylovii were different due to the variance of mean gametic disequilibria and mean gene diversity. AMOVA results showed that the majority of variance (88.62%) occurred within sites, 12.75% was from between-groups. Further research is needed to investigate the selective function in maintaining the genetic diversity of Stipa krylovii plant populations.

Nucleotide diversity and linkage disequilibrium of nine genes with putative effects on flowering time in perennial ryegrass (Lolium perenne L.).

Optimization of flowering is an important breeding goal in forage and turf grasses, such as perennial ryegrass (Lolium perenne L.). Nine floral control genes including Lolium perenne CONSTANS (LpCO), SISTER OF FLOWERING LOCUS T (LpSFT), TERMINAL FLOWER1 (LpTFL1), VERNALIZATION1 (LpVRN1, identical to LpMADS1) and five additional MADS-box genes, were analyzed for nucleotide diversity and linkage disequilibrium (LD). For each gene, about 1 kb genomic fragments were isolated from 10 to 20 genotypes of perennial ryegrass of diverse origin. Four to twelve haplotypes per gene were observed. On average, one single nucleotide polymorphism (SNP) was present per 127 bp between two randomly sampled sequences for the nine genes (π = 0.00790). Two MADS-box genes, LpMADS1 and LpMADS10, involved in timing of flowering showed high nucleotide diversity and rapid LD decay, whereas MADS-box genes involved in floral organ identity were found to be highly conserved and showed extended LD. For LpMADS4, LpMADS5, LpCO, LpSFT and LpTFL1, LD extended over the entire region analyzed. The results are compared to previously published results on resistance genes within the same collection of genotypes and the prospects for association mapping of floral control in perennial ryegrass are discussed.