Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND: Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk. METHOD: We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation. RESULTS: Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05. CONCLUSION: This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk

Homozygous deletion mapping in myeloma samples identifies genes and an expression signature relevant to pathogenesis and outcome.

PURPOSE: Myeloma is a clonal malignancy of plasma cells. Poor-prognosis risk is currently identified by clinical and cytogenetic features. However, these indicators do not capture all prognostic information. Gene expression analysis can be used to identify poor-prognosis patients and this can be improved by combination with information about DNA-level changes. EXPERIMENTAL DESIGN: Using single nucleotide polymorphism-based gene mapping in combination with global gene expression analysis, we have identified homozygous deletions in genes and networks that are relevant to myeloma pathogenesis and outcome. RESULTS: We identified 170 genes with homozygous deletions and corresponding loss of expression. Deletion within the "cell death" network was overrepresented and cases with these deletions had impaired overall survival. From further analysis of these events, we have generated an expression-based signature associated with shorter survival in 258 patients and confirmed this signature in data from two independent groups totaling 800 patients. We defined a gene expression signature of 97 cell death genes that reflects prognosis and confirmed this in two independent data sets. CONCLUSIONS: We developed a simple 6-gene expression signature from the 97-gene signature that can be used to identify poor-prognosis myeloma in the clinical environment. This signature could form the basis of future trials aimed at improving the outcome of poor-prognosis myeloma.

Gene mapping and expression analysis of 16q loss of heterozygosity identifies WWOX and CYLD as being important in determining clinical outcome in multiple myeloma.

We performed fluorescent in situ hybridization (FISH) for 16q23 abnormalities in 861 patients with newly diagnosed multiple myeloma and identified deletion of 16q [del(16q)] in 19.5%. In 467 cases in which demographic and survival data were available, del(16q) was associated with a worse overall survival (OS). It was an independent prognostic marker and conferred additional adverse survival impact in cases with the known poor-risk cytogenetic factors t(4;14) and del(17p). Gene expression profiling and gene mapping using 500K single-nucleotide polymorphism (SNP) mapping arrays revealed loss of heterozygosity (LOH) involving 3 regions: the whole of 16q, a region centered on 16q12 (the location of CYLD), and a region centered on 16q23 (the location of the WW domain-containing oxidoreductase gene WWOX). CYLD is a negative regulator of the NF-kappaB pathway, and cases with low expression of CYLD were used to define a "low-CYLD signature." Cases with 16q LOH or t(14;16) had significantly reduced WWOX expression. WWOX, the site of the translocation breakpoint in t(14;16) cases, is a known tumor suppressor gene involved in apoptosis, and we were able to generate a "low-WWOX signature" defined by WWOX expression. These 2 genes and their corresponding pathways provide an important insight into the potential mechanisms by which 16q LOH confers poor prognosis.

Integration of global SNP-based mapping and expression arrays reveals key regions, mechanisms, and genes important in the pathogenesis of multiple myeloma.

Multiple myeloma is characterized by genomic alterations frequently involving gains and losses of chromosomes. Single nucleotide polymorphism (SNP)-based mapping arrays allow the identification of copy number changes at the sub-megabase level and the identification of loss of heterozygosity (LOH) due to monosomy and uniparental disomy (UPD). We have found that SNP-based mapping array data and fluorescence in situ hybridization (FISH) copy number data correlated well, making the technique robust as a tool to investigate myeloma genomics. The most frequently identified alterations are located at 1p, 1q, 6q, 8p, 13, and 16q. LOH is found in these large regions and also in smaller regions throughout the genome with a median size of 1 Mb. We have identified that UPD is prevalent in myeloma and occurs through a number of mechanisms including mitotic nondisjunction and mitotic recombination. For the first time in myeloma, integration of mapping and expression data has allowed us to reduce the complexity of standard gene expression data and identify candidate genes important in both the transition from normal to monoclonal gammopathy of unknown significance (MGUS) to myeloma and in different subgroups within myeloma. We have documented these genes, providing a focus for further studies to identify and characterize those that are key in the pathogenesis of myeloma.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

AFLP-based genetic mapping of the " bud-flowering" trait in heather (Calluna vulgaris)

Background: Calluna vulgaris is one of the most important landscaping plants produced in Germany. Its enormous economic success is due to the prolonged flower attractiveness of mutants in flower morphology, the so-called bud-bloomers. In this study, we present the first genetic linkage map of C. vulgaris in which we mapped a locus of the economically highly desired trait " flower type" .Results: The map was constructed in JoinMap 4.1. using 535 AFLP markers from a single mapping population. A large fraction (40%) of markers showed distorted segregation. To test the effect of segregation distortion on linkage estimation, these markers were sorted regarding their segregation ratio and added in groups to the data set. The plausibility of group formation was evaluated by comparison of the " two-way pseudo-testcross" and the " integrated" mapping approach. Furthermore, regression mapping was compared to the multipoint-likelihood algorithm. The majority of maps constructed by different combinations of these methods consisted of eight linkage groups corresponding to the chromosome number of C. vulgaris.Conclusions: All maps confirmed the independent inheritance of the most important horticultural traits " flower type" , " flower colour" , and " leaf colour". An AFLP marker for the most important breeding target " flower type" was identified. The presented genetic map of C. vulgaris can now serve as a basis for further molecular marker selection and map-based cloning of the candidate gene encoding the unique flower architecture of C. vulgaris bud-bloomers. © 2013 Behrend et al.; licensee BioMed Central Ltd.

Towards in silico IBV vaccine design: defining the role of polymorphism in viral attenuation

The gammacoronavirus, Infectious Bronchitis Virus (IBV), is a respiratory pathogen of chickens. IBV is a constant threat to poultry production as established vaccines are often ineffective against emerging strains. This requires constant and rapid vaccine production by a process of viral attenuation by egg passage, but the essential forces leading to attenuation in the virus have not yet been characterised. Knowledge of these factors will lead to the development of more effective, rationally attenuated, live vaccines and reduction of the mortality and morbidity caused by this pathogen. M41 CK strain was egg passaged four times many years ago at Houghton Poultry Research Station and stored as M41-CK EP4 (stock virus at The Pirbright Institute since 1992). It was the first egg passage to have its genome pyrosequenced and was therefore used as the baseline reference. The overall aim of this project was to analyse deep sequence data obtained from four IBV isolates (called A, A1, C and D) each originating from the common M41-CK EP4 (ep4) and independently passaged multiple times in embryonated chicken eggs (figure 1.1). Highly polymorphic encoding regions of the IBV genome were then identified which are likely involved in the attenuation process through the formation of independent SNPs and/or SNP clusters. This was then used to direct targeted investigation of SNPs during the attenuation process of the four IBV passages. A previously generated deep sequence dataset was used as a preliminary map of attenuation for one virulent strain of IBV. This investigation showed the nucleocapsid and spike as two highly polymorphic encoding regions within the IBV genome with the highest proportion of SNPs compared to encoding region size. This analysis then led to more focussed studies of the nucleocapsid and spike encoding region with the ultimate aim of mapping key attenuating regions and nucleotide positions. The 454 pyrosequencing data and further investigation of nucleocapsid and spike encoding regions have identified the SNPs present at the same nucleotide positions within analysed A, A1, C and D isolates. These SNPs probably play a crucial role in viral attenuation and universal vaccine production but it is not clear if independent SNPs are also involved in loss of virulence. The majority of SNPs accumulated at different nucleotide positions without further continuation in Sanger sequenced egg passages presenting S2 subunit (spike) and nucleocapsid as polymorphic encoding regions which in nature remain highly conserved.

Genetic Predictors Of Long-term Response To Growth Hormone (gh) Therapy In Children With Gh Deficiency And Turner Syndrome: The Influence Of A Socs2 Polymorphism.

There is great interindividual variability in the response to GH therapy. Ascertaining genetic factors can improve the accuracy of growth response predictions. Suppressor of cytokine signaling (SOCS)-2 is an intracellular negative regulator of GH receptor (GHR) signaling. The objective of the study was to assess the influence of a SOCS2 polymorphism (rs3782415) and its interactive effect with GHR exon 3 and -202 A/C IGFBP3 (rs2854744) polymorphisms on adult height of patients treated with recombinant human GH (rhGH). Genotypes were correlated with adult height data of 65 Turner syndrome (TS) and 47 GH deficiency (GHD) patients treated with rhGH, by multiple linear regressions. Generalized multifactor dimensionality reduction was used to evaluate gene-gene interactions. Baseline clinical data were indistinguishable among patients with different genotypes. Adult height SD scores of patients with at least one SOCS2 single-nucleotide polymorphism rs3782415-C were 0.7 higher than those homozygous for the T allele (P < .001). SOCS2 (P = .003), GHR-exon 3 (P= .016) and -202 A/C IGFBP3 (P = .013) polymorphisms, together with clinical factors accounted for 58% of the variability in adult height and 82% of the total height SD score gain. Patients harboring any two negative genotypes in these three different loci (homozygosity for SOCS2 T allele; the GHR exon 3 full-length allele and/or the -202C-IGFBP3 allele) were more likely to achieve an adult height at the lower quartile (odds ratio of 13.3; 95% confidence interval of 3.2-54.2, P = .0001). The SOCS2 polymorphism (rs3782415) has an influence on the adult height of children with TS and GHD after long-term rhGH therapy. Polymorphisms located in GHR, IGFBP3, and SOCS2 loci have an influence on the growth outcomes of TS and GHD patients treated with rhGH. The use of these genetic markers could identify among rhGH-treated patients those who are genetically predisposed to have less favorable outcomes.

Rh, Kell, Duffy, Kidd And Diego Blood Group System Polymorphism In Brazilian Japanese Descendants.

Polymorphisms of Rh, Kell, Duffy, Kidd and Diego blood group systems were studied in 209 unrelated Brazilian Japanese descendants from South of Brazil. The methods used were multiplex-PCR, AS-PCR and RFLP-PCR. The differences in frequencies among the populations were evaluated using chi-square test. The frequencies for Rh, Kell, Kidd and Diego system were similar to those of the Japanese. RHCE(*)CC, RHCE(*)EE genotypes and FY(*)01 allele were lower and FY(*)01N.01 was higher than Japanese. These differences in the frequencies between Brazilian Japanese descendants and Japanese could indicate a gene flow in Brazilian population and reinforce the importance of this knowledge to achieve safe red blood cells.

The Ccr5Δ32 Polymorphism In Brazilian Patients With Sickle Cell Disease.

Previous studies on the role of inflammation in the pathophysiology of sickle cell disease (SCD) suggested that the CCR5Δ32 allele, which is responsible for the production of truncated C-C chemokine receptor type 5 (CCR5), could confer a selective advantage on patients with SCD because it leads to a less efficient Th1 response. We determined the frequency of the CCR5Δ32 polymorphism in 795 Afro-Brazilian SCD patients followed up at the Pernambuco Hematology and Hemotherapy Center, in Northeastern Brazil, divided into a pediatric group (3 months-17 years, n = 483) and an adult group (18-70 years, n = 312). The adult patients were also compared to a healthy control group (blood donors, 18-61 years, n = 247). The CCR5/CCR5Δ32 polymorphism was determined by allele-specific PCR. No homozygous patient for the CCR5Δ32 allele was detected. The frequency of heterozygotes in the study population (patients and controls) was 5.8%, in the total SCD patients 5.1%, in the children 5.4%, in the adults with SCD 4.8%, and in the adult controls 8.1%. These differences did not reach statistical significance. Our findings failed to demonstrate an important role of the CCR5Δ32 allele in the population sample studied here.

Preference Mapping Of Dulce De Leche Commercialized In Brazilian Markets.

Dulce de leche samples available in the Brazilian market were submitted to sensory profiling by quantitative descriptive analysis and acceptance test, as well sensory evaluation using the just-about-right scale and purchase intent. External preference mapping and the ideal sensory characteristics of dulce de leche were determined. The results were also evaluated by principal component analysis, hierarchical cluster analysis, partial least squares regression, artificial neural networks, and logistic regression. Overall, significant product acceptance was related to intermediate scores of the sensory attributes in the descriptive test, and this trend was observed even after consumer segmentation. The results obtained by sensometric techniques showed that optimizing an ideal dulce de leche from the sensory standpoint is a multidimensional process, with necessary adjustments on the appearance, aroma, taste, and texture attributes of the product for better consumer acceptance and purchase. The optimum dulce de leche was characterized by high scores for the attributes sweet taste, caramel taste, brightness, color, and caramel aroma in accordance with the preference mapping findings. In industrial terms, this means changing the parameters used in the thermal treatment and quantitative changes in the ingredients used in formulations.

Characterization Of Cumulative Joint Damage Patterns In Patients With Rheumatoid Arthritis: A Clinical, Serological, And Gene Polymorphism Perspective.

To characterize cumulative joint damage (CJD) patterns in rheumatoid arthritis (RA) and determine their associations with demographic/clinical features and HLA-DRB1 gene polymorphism. Hand and foot radiographs were obtained from 404 patients with RA. CJD patterns were determined by 3 derivations from Sharp/van der Heijde scores, obtained by the mathematical division of scores for hands/feet (Sharp-h/f score), fingers/wrists (Sharp-f/w score), and erosion/space narrowing (Sharp-e/sn score), respectively. DNA and serum were obtained for determination of HLA-DRB1 polymorphism, rheumatoid factor (RF), and anticitrullinated protein antibodies (ACPA). Patients with wrist-dominant CJD pattern were more likely to have severe RA than those with finger-dominant pattern (68.4% vs 46.0%; p = 0.036) as were those with foot-dominant vs hand-dominant CJD pattern (76.5% vs 56.4%; p = 0.044). HLA-DRB1 shared epitope (SE) alleles were associated with erosion-dominant CJD pattern (p = 0.021). Patients with erosion-dominant CJD pattern had higher levels of RF and ACPA than those with space-narrowing-dominant CJD pattern (median RF 71.35 U/ml vs 22.05 U/ml, respectively; p = 0.003; median ACPA 187.9 U/ml vs 143.2 U/ml, respectively; p < 0.001). The majority of triple-positive patients (SE+, RF+, ACPA+) had erosion-dominant CJD pattern (62.3%) while the majority of triple-negative patients (SE-, FR-, ACPA-) had space narrowing-dominant CJD pattern (75%; p = 0.017). ACPA was associated with HLA-DRB1 SE alleles (p < 0.05). Patients with foot-dominant CJD pattern were taller than those with hand-dominant CJD pattern (p = 0.002); those with erosion-dominant CJD pattern had higher weight and body mass index than those with space narrowing-dominant CJD pattern (p = 0.014, p = 0.001). CJD patterns were associated with disease severity, HLA-DRB1 SE status, presence and titer of ACPA and RF, and morphometric features.

Phylogenetic And Phylogeographic Mapping Of The Avian Coronavirus Spike Protein-encoding Gene In Wild And Synanthropic Birds.

The evolution and population dynamics of avian coronaviruses (AvCoVs) remain underexplored. In the present study, in-depth phylogenetic and Bayesian phylogeographic studies were conducted to investigate the evolutionary dynamics of AvCoVs detected in wild and synanthropic birds. A total of 500 samples, including tracheal and cloacal swabs collected from 312 wild birds belonging to 42 species, were analysed using molecular assays. A total of 65 samples (13%) from 22 bird species were positive for AvCoV. Molecular evolution analyses revealed that the sequences from samples collected in Brazil did not cluster with any of the AvCoV S1 gene sequences deposited in the GenBank database. Bayesian framework analysis estimated an AvCoV strain from Sweden (1999) as the most recent common ancestor of the AvCoVs detected in this study. Furthermore, the analysis inferred an increase in the AvCoV dynamic demographic population in different wild and synanthropic bird species, suggesting that birds may be potential new hosts responsible for spreading this virus.

Microsatellite Markers For Urochloa Humidicola (poaceae) And Their Transferability To Other Urochloa Species.

Urochloa humidicola is a warm-season grass commonly used as forage in the tropics and is recognized for its tolerance to seasonal flooding. This grass is an important forage species for the Cerrado and Amazon regions of Brazil. U. humidicola is a polyploid species with variable ploidy (6X-9X) and facultative apomixis with high phenotypic plasticity. However, this apomixis and ploidy, as well as the limited knowledge of the genetic basis of the germplasm collection, have constrained genetic breeding activities, yet microsatellite markers may enable a better understanding of the species' genetic composition. This study aimed to develop and characterize new polymorphic microsatellite molecular markers in U. humidicola and to evaluate their transferability to other Urochloa species. A set of microsatellite markers for U. humidicola was identified from two new enriched genomic DNA libraries: the first library was constructed from a single sexual genotype and the second from a pool of eight apomictic genotypes selected on the basis of previous results. Of the 114 loci developed, 72 primer pairs presented a good amplification product, and 64 were polymorphic among the 34 genotypes tested. The number of bands per simple sequence repeat (SSR) locus ranged from 1 to 29, with a mean of 9.6 bands per locus. The mean polymorphism information content (PIC) of all loci was 0.77, and the mean discrimination power (DP) was 0.87. STRUCTURE analysis revealed differences among U. humidicola accessions, hybrids, and other Urochloa accessions. The transferability of these microsatellites was evaluated in four species of the genus, U. brizantha, U. decumbens, U. ruziziensis, and U. dictyoneura, and the percentage of transferability ranged from 58.33% to 69.44% depending on the species. This work reports new polymorphic microsatellite markers for U. humidicola that can be used for breeding programs of this and other Urochloa species, including genetic linkage mapping, quantitative trait loci identification, and marker-assisted selection.

Polymorphism In Lep And Lepr May Modify Leptin Levels And Represent Risk Factors For Thyroid Cancer.

Purpose. To understand the role of polymorphisms in the LEP (rs7799039 and rs2167270) and LEPR (rs1137101 and rs1137100) genes in DTC susceptibility and their effect on leptin levels. Methods. We studied 153 patients with DTC and 234 controls through TaqMan SNP Genotyping and ELISA, comparing these data to the clinicopathological data of patients with DTC. Results. Patients with AA genotype of rs7799039 had higher levels of serum leptin (9.22 ± 0.98 ng/mL) than those with AG genotype (10.07 ± 0.60 ng/mL; P = 0.005). Individuals with AG genotype of rs2167270 also produced higher serum leptin levels (10.05 ± 0.59 ng/mL) than the subjects with GG genotype (9.52 ± 0.79 ng/mL; P < 0.05). A multivariate logistic regression adjusted for gender, age, and BMI showed that the AG genotype of rs7799039 was an independent risk for DTC (OR, 11.689; P = 0.0183; 95% CI, 1.516-90.119). Similarly, AG and GG genotypes of rs1137101 increased the susceptibility to DTC (OR, 3.747; P = 0.027; 95% CI, 1.161-12.092 and OR, 5.437; P = 0.013; 95% CI, 1.426-20.729). Conclusions. We demonstrated that rs7799039 and rs2167270 polymorphisms modify the serum leptin concentrations in patients with DTC. Furthermore, polymorphisms rs7799039 and rs1137101 increase the risk of DTC development, although they do not correlate with tumor aggressiveness.