Genome-wide association study of intraocular pressure identifies the GLCCI1/ICA1 region as a glaucoma susceptibility locus

To discover quantitative trait loci for intraocular pressure, amajor risk factor for glaucoma and the only modifiable one,weperformed agenome-wide association studyonadiscoverycohort of2175individualsfromSydney, Australia. We found a novel association between intraocular pressure and a common variant at 7p21 near to GLCCI1 and ICA1. The findings in this region were confirmed through two UK replication cohorts totalling 4866 individuals (rs59072263, Pcombined = 1.10 × 10-8). A copy of the G allele at this SNP is associated with an increase in mean IOP of 0.45 mmHg (95%CI = 0.30-0.61 mmHg). These results lend support to the implication of vesicle trafficking and glucocorticoid inducibility pathways in the determination of intraocular pressure and in the pathogenesis of primary open-angle glaucoma.

Detecting quantitative trait loci affecting beef tenderness on bovine chromosome 7 near calpastatin and lysyl oxidase

From a study of 3 large half-sib families of cattle, we describe linkage between DNA polymorphisms on bovine chromosome 7 and meat tenderness. Quantitative trait loci (QTL) for Longissimus lumborum peak force (LLPF) and Semitendonosis adhesion (STADH) were located to this map of DNA markers, which includes the calpastatin ( CAST) and lysyl oxidase (LOX) genes. The LLPF QTL has a maximum lodscore of 4.9 and allele substitution of approximately 0.80 of a phenotypic standard deviation, and the peak is located over the CAST gene. The STADH QTL has a maximum lodscore of 3.5 and an allele substitution of approximately 0.37 of a phenotypic standard deviation, and the peak is located over the LOX gene. This suggests 2 separate likelihood peaks on the chromosome. Further analyses of meat tenderness measures in the Longissimus lumborum, LLPF and LL compression (LLC), in which outlier individuals or kill groups are removed, demonstrate large shifts in the location of LLPF QTL, as well as confirming that there are indeed 2 QTL on bovine chromosome 7. We found that both QTL are reflected in both LLPF and LLC measurements, suggesting that both these components of tenderness, myofibrillar and connective tissue, are detected by both measurements in this muscle.

Simple sequence repeat markers associated with three quantitative trait loci for black point resistance can be used to enrich selection populations in bread wheat

Black point in wheat has the potential to cost the Australian industry $A30.4 million a year. It is difficult and expensive to screen for resistance, so the aim of this study was to validate 3 previously identified quantitative trait loci (QTLs) for black point resistance on chromosomes 2B, 4A, and 3D of the wheat variety Sunco. Black point resistance data and simple sequence repeat (SSR) markers, linked to the resistance QTLs and suited to high-throughput assay, were analysed in the doubled haploid population, Batavia (susceptible) × Pelsart (resistant). Sunco and Pelsart both have Cook in their pedigree and both have the Triticum timopheevii translocation on 2B. SSR markers identified for the 3 genetic regions were gwm319 (2B, T. timopheevii translocation), wmc048 (4AS), and gwm341 (3DS). Gwm319 and wmc048 were associated with black point resistance in the validation population. Gwm341 may have an epistatic influence on the trait because when resistance alleles were present at both gwm319 and wmc048, the Batavia-derived allele at gwm341 was associated with a higher proportion of resistant lines. Data are presented showing the level of enrichment achieved for black point resistance, using 1, 2, or 3 of these molecular markers, and the number of associated discarded resistant lines. The level of population enrichment was found to be 1.83-fold with 6 of 17 resistant lines discarded when gwm319 and wmc048 were both used for selection. Interactions among the 3 QTLs appear complex and other genetic and epigenetic factors influence susceptibility to black point. Polymorphism was assessed for these markers within potential breeding material. This indicated that alternative markers to wmc048 may be required for some parental combinations. Based on these results, marker-assisted selection for the major black point resistance QTLs can increase the rate of genetic gain by improving the selection efficiency and may facilitate stacking of black point resistances from different sources.

A quantitative-trait genome-wide association study of alcoholism risk in the community: findings and implications

BACKGROUND Given moderately strong genetic contributions to variation in alcoholism and heaviness of drinking (50% to 60% heritability) with high correlation of genetic influences, we have conducted a quantitative trait genome-wide association study (GWAS) for phenotypes related to alcohol use and dependence. METHODS Diagnostic interview and blood/buccal samples were obtained from sibships ascertained through the Australian Twin Registry. Genome-wide single nucleotide polymorphism (SNP) genotyping was performed with 8754 individuals (2062 alcohol-dependent cases) selected for informativeness for alcohol use disorder and associated quantitative traits. Family-based association tests were performed for alcohol dependence, dependence factor score, and heaviness of drinking factor score, with confirmatory case-population control comparisons using an unassessed population control series of 3393 Australians with genome-wide SNP data. RESULTS No findings reached genome-wide significance (p = 8.4 x 10(-8) for this study), with lowest p value for primary phenotypes of 1.2 x 10(-7). Convergent findings for quantitative consumption and diagnostic and quantitative dependence measures suggest possible roles for a transmembrane protein gene (TMEM108) and for ANKS1A. The major finding, however, was small effect sizes estimated for individual SNPs, suggesting that hundreds of genetic variants make modest contributions (1/4% of variance or less) to alcohol dependence risk. CONCLUSIONS We conclude that: - 1) meta-analyses of consumption data may contribute usefully to gene discovery; - 2) translation of human alcoholism GWAS results to drug discovery or clinically useful prediction of risk will be challenging, and; - 3) through accumulation across studies, GWAS data may become valuable for improved genetic risk differentiation in research in biological psychiatry (e.g., prospective high-risk or resilience studies).

Quantitative Trait Loci for CD4:CD8 Lymphocyte Ratio Are Associated with Risk of Type 1 Diabetes and HIV-1 Immune Control

Abnormal expansion or depletion of particular lymphocyte subsets is associated with clinical manifestations such as HIV progression to AIDS and autoimmune disease. We sought to identify genetic predictors of lymphocyte levels and reasoned that these may play a role in immune-related diseases. We tested 2.3 million variants for association with five lymphocyte subsets, measured in 2538 individuals from the general population, including CD4+ T cells, CD8+ T cells, CD56+ natural killer (NK) cells, and the derived measure CD4:CD8 ratio. We identified two regions of strong association. The first was located in the major histocompatibility complex (MHC), with multiple SNPs strongly associated with CD4:CD8 ratio (rs2524054, p = 2.1 × 10−28). The second region was centered within a cluster of genes from the Schlafen family and was associated with NK cell levels (rs1838149, p = 6.1 × 10−14). The MHC association with CD4:CD8 replicated convincingly (p = 1.4 × 10−9) in an independent panel of 988 individuals. Conditional analyses indicate that there are two major independent quantitative trait loci (QTL) in the MHC region that regulate CD4:CD8 ratio: one is located in the class I cluster and influences CD8 levels, whereas the second is located in the class II cluster and regulates CD4 levels. Jointly, both QTL explained 8% of the variance in CD4:CD8 ratio. The class I variants are also strongly associated with durable host control of HIV, and class II variants are associated with type-1 diabetes, suggesting that genetic variation at the MHC may predispose one to immune-related diseases partly through disregulation of T cell homeostasis.

Quantitative trait loci detection and benefits from marker-assisted selection in dairy cattle

Knowing the chromosomal areas or actual genes affecting the traits under selection would add more information to be used in the selection decisions which would potentially lead to higher genetic response. The first objective of this study was to map quantitative trait loci (QTL) affecting economically important traits in the Finnish Ayrshire population. The second objective was to investigate the effects of using QTL information in marker-assisted selection (MAS) on the genetic response and the linkage disequilibrium between the different parts of the genome. Whole genome scans were carried out on a grand-daughter design with 12 half-sib families and a total of 493 sons. Twelve different traits were studied: milk yield, protein yield, protein content, fat yield, fat content, somatic cell score (SCS), mastitis treatments, other veterinary treatments, days open, fertility treatments, non-return rate, and calf mortality. The average spacing of the typed markers was 20 cM with 2 to 14 markers per chromosome. Associations between markers and traits were analyzed with multiple marker regression. Significance was determined by permutation and genome-wise P-values obtained by Bonferroni correction. The benefits from MAS were investigated by simulation: a conventional progeny testing scheme was compared to a scheme where QTL information was used within families to select among full-sibs in the male path. Two QTL on different chromosomes were modelled. The effects of different starting frequencies of the favourable alleles and different size of the QTL effects were evaluated. A large number of QTL, 48 in total, were detected at 5% or higher chromosome-wise significance. QTL for milk production were found on 8 chromosomes, for SCS on 6, for mastitis treatments on 1, for other veterinary treatments on 5, for days open on 7, for fertility treatments on 7, for calf mortality on 6, and for non-return rate on 2 chromosomes. In the simulation study the total genetic response was faster with MAS than with conventional selection and the advantage of MAS persisted over the studied generations. The rate of response and the difference between the selection schemes reflected clearly the changes in allele frequencies of the favourable QTL. The disequilibrium between the polygenes and QTL was always negative and it was larger with larger QTL size. The disequilibrium between the two QTL was larger with QTL of large effect and it was somewhat larger with MAS for scenarios with starting frequencies below 0.5 for QTL of moderate size and below 0.3 for large QTL. In conclusion, several QTL affecting economically important traits of dairy cattle were detected. Further studies are needed to verify these QTL, check their presence in the present breeding population, look for pleiotropy and fine map the most interesting QTL regions. The results of the simulation studies show that using MAS together with embryo transfer to pre-select young bulls within families is a useful approach to increase the genetic merit of the AI-bulls compared to conventional selection.

Identification of quantitative trait loci for growth-related traits in the Pacific abalone Haliotis discus hannai Ino

The locations and effects of quantitative trait loci (QTL) were estimated for nine characters for growth-related traits in the Pacific abalone (Haliotis discus hannai Ino) using a randomly amplified polymorphic DNA (RAPD), amplification fragment length polymorphism (AFLP) and SSR genetic linkage map. Twenty-eight putatively significant QTLs (LOD > 2.4) were detected for nine traits (shell length, shell width, total weight, shell weight, weight of soft part, muscle weight, gonad and digestive gland weight, mantle weight and gill weight). The percentage of phenotypic variation explained by a single QTL ranged from 8.0% to 35.9%. The significant correlations (P < 0.001) were found among all the growth-related traits, and Pearson's correlation coefficients were more than 0.81. For the female map, the QTL for growth were concentrated on groups 1 and 4 linkage maps. On the male map, the QTL that influenced growth-related traits gathered on the groups 1 and 9 linkage maps. Genetic linkage map construction and QTL analysis for growth-related traits are the basis for the marker-assisted selection and will eventually improve production and quality of the Pacific abalone.

Genetic and genomic analyses of musculoskeletal differences between BEH and BEL strains

Berlin high (BEH) and Berlin low (BEL) strains selected for divergent growth differ 3-fold in body weight. We aimed at examining muscle mass, which is a major contributor to body weight, by exploring anatomical characteristics of the soleus muscle, its fiber numbers and their cross sectional area (CSA), by analysing transcriptome of the gastrocnemius and by initiating quantitative trait locus (QTL) mapping. BEH muscles were 4-to-8 times larger compared to BEL strain. In sub-strain BEH+/+, mutant myostatin was replaced with a wild type allele, however, BEH+/+muscles still were 2-to-4 times larger compared to the BEL strain. BEH soleus contained 2-times more (P<0.0001) and 2-times larger in CSA (P<0.0001) fibers compared to BEL strain. In addition, soleus femoral attachment anomaly (SFAA) was observed in all BEL mice. One significant (chromosome 1) and four suggestive (chromosomes 3, 4, 6 and 9) muscle weight QTLs were mapped in 21-day old F2 intercross (n=296) between BEH and BEL strains. The frequency of SFAA incidence in the F2 and in the backcross to BEL strain (BCL) suggested the presence of more than one causative gene. Two suggestive SFAA QTLs were mapped in BCL, however, their peak markers were not associated with the phenotype in F2. RNA-Seq analysis revealed 2,148 differentially expressed (P<0.1) genes and 45,673 SNPs and >2,000 indels between BEH+/+ and BEL males. In conclusion, contrasting muscle traits, genomic and gene expression differences between BEH and BEL strains provide a promising model for the search of genes involved in muscle growth and musculoskeletal morphogenesis.

Systems Genetics of Metabolism: The Use of the BXD Murine Reference Panel for Multiscalar Integration of Traits.

Metabolic homeostasis is achieved by complex molecular and cellular networks that differ significantly among individuals and are difficult to model with genetically engineered lines of mice optimized to study single gene function. Here, we systematically acquired metabolic phenotypes by using the EUMODIC EMPReSS protocols across a large panel of isogenic but diverse strains of mice (BXD type) to study the genetic control of metabolism. We generated and analyzed 140 classical phenotypes and deposited these in an open-access web service for systems genetics (www.genenetwork.org). Heritability, influence of sex, and genetic modifiers of traits were examined singly and jointly by using quantitative-trait locus (QTL) and expression QTL-mapping methods. Traits and networks were linked to loci encompassing both known variants and novel candidate genes, including alkaline phosphatase (ALPL), here linked to hypophosphatasia. The assembled and curated phenotypes provide key resources and exemplars that can be used to dissect complex metabolic traits and disorders.

A functional genomics approach reveals novel quantitative trait loci associated with platelet signalling pathways

Platelet response to activation varies widely between individuals but shows interindividual consistency and strong heritability. The genetic basis of this variation has not been properly explored. We therefore systematically measured the effect on function of sequence variation in 97 candidate genes in the collagen and adenosine-diphosphate (ADP) signaling pathways. Resequencing of the genes in 48 European DNA samples nearly doubled the number of known single nucleotide polymorphisms (SNPs) and informed the selection of 1327 SNPs for genotyping in 500 healthy Northern European subjects with known platelet responses to collagen-related peptide (CRP-XL) and ADP. This identified 17 novel associations with platelet function (P < .005) accounting for approximately 46% of the variation in response. Further investigations with platelets of known genotype explored the mechanisms behind some of the associations. SNPs in PEAR1 associated with increased platelet response to CRP-XL and increased PEAR1 protein expression after platelet degranulation. The minor allele of a 3' untranslated region (UTR) SNP (rs2769668) in VAV3 was associated with higher protein expression (P = .03) and increased P-selectin exposure after ADP activation (P = .004). Furthermore the minor allele of the intronic SNP rs17786144 in ITPR1 modified Ca2+ levels after activation with ADP (P < .004). These data provide novel insights into key hubs within platelet signaling networks.

Adaptation of tree growth to elevated CO2: quantitative trait loci for biomass in Populus

Information on the genetic variation of plant response to elevated CO2 (e[CO2]) is needed to understand plant adaptation and to pinpoint likely evolutionary response to future high atmospheric CO2 concentrations.• Here, quantitative trait loci (QTL) for above- and below-ground tree growth were determined in a pedigree – an F2 hybrid of poplar (Populus trichocarpa and Populus deltoides), following season-long exposure to either current day ambient CO2 (a[CO2]) or e[CO2] at 600 µl l−1, and genotype by environment interactions investigated.• In the F2 generation, both above- and below-ground growth showed a significant increase in e[CO2]. Three areas of the genome on linkage groups I, IX and XII were identified as important in determining above-ground growth response to e[CO2], while an additional three areas of the genome on linkage groups IV, XVI and XIX appeared important in determining root growth response to e[CO2].• These results quantify and identify genetic variation in response to e[CO2] and provide an insight into genomic response to the changing environment

Genetical and comparative genomics of Brassica under altered Ca supply identifies Arabidopsis Ca-transporter orthologs

Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca2+ transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca2+ transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.

Identification of quantitative trait loci affecting resistance to gastrointestinal parasites in a double backcross population of Red Maasai and Dorper sheep

A genome-wide scan for quantitative trait loci (QTL) affecting gastrointestinal nematode resistance in sheep was completed using a double backcross population derived from Red Maasai and Dorper ewes bred to F1 rams. This design provided an opportunity to map potentially unique genetic variation associated with a parasite-tolerant breed like Red Maasai, a breed developed to survive East African grazing conditions. Parasite indicator phenotypes (blood packed cell volume PCV and faecal egg count FEC) were collected on a weekly basis from 1064 lambs during a single 3-month post-weaning grazing challenge on infected pastures. The averages of last measurements for FEC (AVFEC) and PCV (AVPCV), along with decline in PCV from challenge start to end (PCVD), were used to select lambs (N = 371) for genotyping that represented the tails (10% threshold) of the phenotypic distributions. Marker genotypes for 172 microsatellite loci covering 25 of 26 autosomes (1560.7 cm) were scored and corrected by Genoprob prior to qxpak analysis that included BoxCox transformed AVFEC and arcsine transformed PCV statistics. Significant QTL for AVFEC and AVPCV were detected on four chromosomes, and this included a novel AVFEC QTL on chromosome 6 that would have remained undetected without BoxCox transformation methods. The most significant P-values for AVFEC, AVPCV and PCVD overlapped the same marker interval on chromosome 22, suggesting the potential for a single causative mutation, which remains unknown. In all cases, the favourable QTL allele was always contributed from Red Maasai, providing support for the idea that future marker-assisted selection for genetic improvement of production in East Africa will rely on markers in linkage disequilibrium with these QTL.

A Bayesian approach for constructing genetic maps when markers are miscoded

The advent of molecular markers has created opportunities for a better understanding of quantitative inheritance and for developing novel strategies for genetic improvement of agricultural species, using information on quantitative trait loci (QTL). A QTL analysis relies on accurate genetic marker maps. At present, most statistical methods used for map construction ignore the fact that molecular data may be read with error. Often, however, there is ambiguity about some marker genotypes. A Bayesian MCMC approach for inferences about a genetic marker map when random miscoding of genotypes occurs is presented, and simulated and real data sets are analyzed. The results suggest that unless there is strong reason to believe that genotypes are ascertained without error, the proposed approach provides more reliable inference on the genetic map.

Factors that cause genotype by environment interaction and use of a multiple-trait herd-cluster model for milk yield of Holstein cattle from Brazil and Colombia

Descriptive herd variables (DVHE) were used to explain genotype by environment interactions (G x E) for milk yield (MY) in Brazilian and Colombian production environments and to develop a herd-cluster model to estimate covariance components and genetic parameters for each herd environment group. Data consisted of 180,522 lactation records of 94,558 Holstein cows from 937 Brazilian and 400 Colombian herds. Herds in both countries were jointly grouped in thirds according to 8 DVHE: production level, phenotypic variability, age at first calving, calving interval, percentage of imported semen, lactation length, and herd size. For each DVHE, REML bivariate animal model analyses were used to estimate genetic correlations for MY between upper and lower thirds of the data. Based on estimates of genetic correlations, weights were assigned to each DVHE to group herds in a cluster analysis using the FASTCLUS procedure in SAS. Three clusters were defined, and genetic and residual variance components were heterogeneous among herd clusters. Estimates of heritability in clusters 1 and 3 were 0.28 and 0.29, respectively, but the estimate was larger (0.39) in Cluster 2. The genetic correlations of MY from different clusters ranged from 0.89 to 0.97. The herd-cluster model based on DVHE properly takes into account G x E by grouping similar environments accordingly and seems to be an alternative to simply considering country borders to distinguish between environments.