Association between Polymorphisms in the TSHR Gene and Graves' Orbitopathy

Background

Graves'' orbitopathy (GO) as well as Graves'' disease (GD) hyperthyroidism originate from an autoimmune reaction against the common auto-antigen, thyroid-stimulating hormone receptor (TSHR). GO phenotype is associated with environmental risk factors, mainly nicotinism, as well as genetic risk factors which initiate an immunologic reaction. In some patients GO is observed before diagnosis of GD hyperthyroidism, while it can also be observed far after diagnosis. The intensity of GO symptoms varies greatly in these patients. Thus, the pathogenesis of GD and GO may correlate with different genetic backgrounds, which has been confirmed by studies of correlations between GO and polymorphisms in cytokines involved in orbit inflammation. The aim of our analysis was to assess genetic predisposition to GO in young patients (age of diagnosis ≤30 years of age), for whom environmental effects had less time to influence outcomes than in adults.

Methods

768 GD patients were included in the study. 359 of them had clinically evident orbitopathy (NOSPECS ≥2). Patients were stratified by age at diagnosis. Association analyses were performed for genes with a known influence on development of GD - TSHR, HLA-DRB1, cytotoxic T-lymphocyte antigen 4 (CTLA4) and lymphoid protein tyrosine phosphatase (PTPN22).

Results

The rs179247 TSHR polymorphism was associated with GO in young patients only. In young GO-free patients, allele A was statistically more frequent and homozygous carriers had a considerable lower risk of disease incidence than patients with AG or GG genotypes. Those differences were not found in either elderly patients or the group analyzed as a whole.

Conclusions

Allele A of the rs179247 polymorphism in the TSHR gene is associated with lower risk of GO in young GD patients.     

Differences in Gene-Gene Interactions in Graves’ Disease Patients Stratified by Age of Onset

Background

Graves’ disease (GD) is a complex disease in which genetic predisposition is modified by environmental factors. Each gene exerts limited effects on the development of autoimmune disease (OR = 1.2–1.5). An epidemiological study revealed that nearly 70% of the risk of developing inherited autoimmunological thyroid diseases (AITD) is the result of gene interactions. In the present study, we analyzed the effects of the interactions of multiple loci on the genetic predisposition to GD. The aim of our analyses was to identify pairs of genes that exhibit a multiplicative interaction effect.

Material and Methods

A total of 709 patients with GD were included in the study. The patients were stratified into more homogeneous groups depending on the age at time of GD onset: younger patients less than 30 years of age and older patients greater than 30 years of age. Association analyses were performed for genes that influence the development of GD: HLADRB1, PTPN22, CTLA4 and TSHR. The interactions among polymorphisms were analyzed using the multiple logistic regression and multifactor dimensionality reduction (MDR) methods.

Results

GD patients stratified by the age of onset differed in the allele frequencies of the HLADRB1*03 and 1858T polymorphisms of the PTPN22 gene (OR = 1.7, p = 0.003; OR = 1.49, p = 0.01, respectively). We evaluated the genetic interactions of four SNPs in a pairwise fashion with regard to disease risk. The coexistence of HLADRB1 with CTLA4 or HLADRB1 with PTPN22 exhibited interactions on more than additive levels (OR = 3.64, p = 0.002; OR = 4.20, p < 0.001, respectively). These results suggest that interactions between these pairs of genes contribute to the development of GD. MDR analysis confirmed these interactions.

Conclusion

In contrast to a single gene effect, we observed that interactions between the HLADRB1/PTPN22 and HLADRB1/CTLA4 genes more closely predicted the risk of GD onset in young patients.     

Fetal microchimeric cells in blood of women with an autoimmune thyroid disease

Context

Hashimoto''s thyroiditis (HT) and Graves'' disease (GD), two autoimmune thyroid diseases (AITD), occur more frequently in women than in men and show an increased incidence in the years following parturition. Persisting fetal cells could play a role in the development of these diseases.

Objective

Aim of this study was to detect and characterize fetal cells in blood of postpartum women with and without an AITD.

Participants

Eleven patients with an AITD and ten healthy volunteers, all given birth to a son maximum 5 years before analysis, and three women who never had been pregnant, were included. None of them had any other disease of the thyroid which could interfere with the results obtained.

Methods

Fluorescence in situ hybridization (FISH) and repeated FISH were used to count the number of male fetal cells. Furthermore, the fetal cells were further characterized.

Results

In patients with HT, 7 to 11 fetal cells per 1.000.000 maternal cells were detected, compared to 14 to 29 fetal cells in patients with GD (p = 0,0061). In patients with HT, mainly fetal CD8⁺ T cells were found, while in patients with GD, fetal B and CD4⁺ T cells were detected. In healthy volunteers with son, 0 to 5 fetal cells were observed, which was significantly less than the number observed in patients (p<0,05). In women who never had been pregnant, no male cells were detected.

Conclusion

This study shows a clear association between fetal microchimeric cells and autoimmune thyroid diseases.     

The Effect of Single Nucleotide Polymorphisms from Genome Wide Association Studies in Multiple Sclerosis on Gene Expression

Background

Multiple sclerosis (MS) is a complex neurological disorder. Its aetiology involves both environmental and genetic factors. Recent genome-wide association studies have identified a number of single nucleotide polymorphisms (SNPs) associated with susceptibility to (MS). We investigated whether these genetic variations were associated with alteration in gene expression.

Methods/Principal Findings

We used a database of mRNA expression and genetic variation derived from immortalised peripheral lymphocytes to investigate polymorphisms associated with MS for correlation with gene expression. Several SNPs were found to be associated with changes in expression: in particular two with HLA-DQA1, HLA-DQA2, HLA-DQB1, HLA-DRB1, HLA-DRB4 and HLA-DRB5, one with ZFP57, one with CD58, two with IL7 and FAM164A, and one with FAM119B, TSFM and KUB3. We found minimal cross-over with a recent whole genome expression study in MS patients.

Discussion

We have shown that many susceptibility loci in MS are associated with changes in gene expression using an unbiased expression database. Several of these findings suggest novel gene candidates underlying the effects of MS-associated genetic variation.     

Thyroid stimulating hormone receptor (TSHR) intron 1 variants are major risk factors for Graves' disease in three European Caucasian cohorts

Background

The thyroid stimulating hormone receptor (TSHR) gene is an established susceptibility locus for Graves'' disease (GD), with recent studies refining association to two single nucleotide polymorphisms (SNPs), rs179247 and rs12101255, within TSHR intron 1.

Methodology and Principal Findings

We aimed to validate association of rs179247 and rs12101255 in Polish and UK Caucasian GD case-control subjects, determine the mode of inheritance and to see if association correlates with specific GD clinical manifestations. We investigated three case-control populations; 558 GD patients and 520 controls from Warsaw, Poland, 196 GD patients and 198 controls from Gliwice, Poland and 2504 GD patients from the UK National collection and 2784 controls from the 1958 British Birth cohort. Both rs179247 (P = 1.2×10⁻²–6.2×10⁻¹⁵, OR = 1.38–1.45) and rs12101255 (P = 1.0×10⁻⁴–3.68×10⁻²¹, OR = 1.47–1.87) exhibited strong association with GD in all three cohorts. Logistic regression suggested association of rs179247 is secondary to rs12101255 in all cohorts. Inheritance modeling suggested a co-dominant mode of inheritance in all cohorts. Genotype-phenotype correlations provided no clear evidence of association with any specific clinical characteristics.

Conclusions

We have validated association of TSHR intron 1 SNPs with GD in three independent European cohorts and have demonstrated that the aetiological variant within the TSHR is likely to be in strong linkage disequilibrium with rs12101255. Fine mapping is now required to determine the exact location of the aetiological DNA variants within the TSHR.     

Dopaminergic Neuronal Imaging in Genetic Parkinson's Disease: Insights into Pathogenesis

Objectives

To compare the dopaminergic neuronal imaging features of different subtypes of genetic Parkinson''s Disease.

Methods

A retrospective study of genetic Parkinson''s diseases cases in which DaTSCAN (123I-FP-CIT) had been performed. Specific non-displaceable binding was calculated for bilateral caudate and putamen for each case. The right:left asymmetry index and striatal asymmetry index was calculated.

Results

Scans were available from 37 cases of monogenetic Parkinson''s disease (7 glucocerebrosidase (GBA) mutations, 8 alpha-synuclein, 3 LRRK2, 7 PINK1, 12 Parkin). The asymmetry of radioligand uptake for Parkinson''s disease with GBA or LRRK2 mutations was greater than that for Parkinson''s disease with alpha synuclein, PINK1 or Parkin mutations.

Conclusions

The asymmetry of radioligand uptake in Parkinsons disease associated with GBA or LRRK2 mutations suggests that interactions with additional genetic or environmental factors may be associated with dopaminergic neuronal loss.     

Genetic diversity of selected genes that are potentially economically important in feral sheep of New Zealand

Background

Feral sheep are considered to be a source of genetic variation that has been lost from their domestic counterparts through selection.

Methods

This study investigates variation in the genes KRTAP1-1, KRT33, ADRB3 and DQA2 in Merino-like feral sheep populations from New Zealand and its offshore islands. These genes have previously been shown to influence wool, lamb survival and animal health.

Results

All the genes were polymorphic, but no new allele was identified in the feral populations. In some of these populations, allele frequencies differed from those observed in commercial Merino sheep and other breeds found in New Zealand. Heterozygosity levels were comparable to those observed in other studies on feral sheep. Our results suggest that some of the feral populations may have been either inbred or outbred over the duration of their apparent isolation.

Conclusion

The variation described here allows us to draw some conclusions about the likely genetic origin of the populations and selective pressures that may have acted upon them, but they do not appear to be a source of new genetic material, at least for these four genes.     

CTLA4 Variants and Haplotype Contribute Genetic Susceptibility to Myasthenia Gravis in Northern Chinese Population

Background

Cytotoxic T lymphocyte-associated antigen-4 (CTLA4), a critical negative regulator of the T-cell response, has been considered a candidate for many autoimmune diseases. Evidence from Caucasians supported a genetic predisposition of CTLA4 to myasthenia gravis (MG), but the contribution in East Asians has not been established.

Objectives

To investigate the role of CTLA4 variants in the susceptibility to MG and the contribution to subtypes of MG.

Methods

Six autoimmune disease-related risk alleles of CTLA4 (rs1863800, rs733618, rs4553808, rs5742909, rs231775, and rs3087243) were investigated for MG in northern Chinese. 168 patients with MG (mean age 37.1±20.5 years, 64 men and 104 women) and 233 healthy controls (mean age 53.3±8.7 years, 96 men and 137 women) were screened, and the contribution of CTLA4 to the general risk of MG and each subgroup was explored.

Results

rs1863800*C, rs733618*C, and rs231775*G were significantly associated with the whole cohort of patients with MG after permutation correction for multiple-testing adjustment (P = 0.027, 0.001, and 0.032, respectively). A risk haplotype (CCACG) [odds ratio (OR) = 1.535, range = 1.150–2.059, P = 0.004)] was also identified. The stratified subtype analysis indicated that the positive contribution was possibly derived from early onset MG (EOMG), seropositive MG (SPMG), female patients, and MG without thymoma. No association was observed in juvenile MG/LOMG, and MG coupled with thymoma.

Conclusion

A predisposing effect of rs1863800*C, rs733618*C, and rs231775*G of CTLA4 gene to general risk of MG in Chinese was demonstrated for the first time, which was likely derived from EOMG, SPMG, MG without thymoma and the female patients.     

Rational design of a plasmid origin that replicates efficiently in both gram-positive and gram-negative bacteria

Background

Most plasmids replicate only within a particular genus or family.

Methodology/Principal Findings

Here we describe an engineered high copy number expression vector, pBAV1K-T5, that produces varying quantities of active reporter proteins in Escherichia coli, Acinetobacter baylyi ADP1, Agrobacterium tumefaciens, (all Gram-negative), Streptococcus pneumoniae, Leifsonia shinshuensis, Peanibacillus sp. S18-36 and Bacillus subtilis (Gram-positive).

Conclusions/Significance

Our results demonstrate the efficiency of pBAV1K-T5 replication in different bacterial species, thereby facilitating the study of proteins that don''t fold well in E. coli and pathogens not amenable to existing genetic tools.     

Non-HLA genes PTPN22, CDK6 and PADI4 are associated with specific autoantibodies in HLA-defined subgroups of rheumatoid arthritis

Introduction

Genetic susceptibility to complex diseases has been intensively studied during the last decade, yet only signals with small effect have been found leaving open the possibility that subgroups within complex traits show stronger association signals. In rheumatoid arthritis (RA), autoantibody production serves as a helpful discriminator in genetic studies and today anti-citrullinated cyclic peptide (anti-CCP) antibody positivity is employed for diagnosis of disease. The HLA-DRB1 locus is known as the most important genetic contributor for the risk of RA, but is not sufficient to drive autoimmunity and additional genetic and environmental factors are involved. Hence, we addressed the association of previously discovered RA loci with disease-specific autoantibody responses in RA patients stratified by HLA-DRB1*04.

Methods

We investigated 2178 patients from three RA cohorts from Sweden and Spain for 41 genetic variants and four autoantibodies, including the generic anti-CCP as well as specific responses towards citrullinated peptides from vimentin, alpha-enolase and type II collagen.

Results

Our data demonstrated different genetic associations of autoantibody-positive disease subgroups in relation to the presence of DRB1*04. Two specific subgroups of autoantibody-positive RA were identified. The SNP in PTPN22 was associated with presence of anti-citrullinated enolase peptide antibodies in carriers of HLA-DRB1*04 (Cochran-Mantel-Haenszel test P = 0.0001, P_corrected <0.05), whereas SNPs in CDK6 and PADI4 were associated with anti-CCP status in DRB1*04 negative patients (Cochran-Mantel-Haenszel test P = 0.0004, P_corrected <0.05 for both markers). Additionally we see allelic correlation with autoantibody titers for PTPN22 SNP rs2476601 and anti-citrullinated enolase peptide antibodies in carriers of HLA-DRB1*04 (Mann Whitney test P = 0.02) and between CDK6 SNP rs42041 and anti-CCP in non-carriers of HLA-DRB1*04 (Mann Whitney test P = 0.02).

Conclusion

These data point to alternative pathways for disease development in clinically similar RA subgroups and suggest an approach for study of genetic complexity of disease with strong contribution of HLA.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-014-0414-3) contains supplementary material, which is available to authorized users.     

Comprehensive genotyping in two homogeneous Graves' disease samples reveals major and novel HLA association alleles

Background

Graves'' disease (GD) is the leading cause of hyperthyroidism and thyroid eye disease inherited as a complex trait. Although geoepidemiology studies showed relatively higher prevalence of GD in Asians than in Caucasians, previous genetic studies were contradictory concerning whether and/or which human leukocyte antigen (HLA) alleles are associated with GD in Asians.

Methodology/Principal Findings

We conducted a case-control association study (499 unrelated GD cases and 504 controls) and a replication in an independent family sample (419 GD individuals and their 282 relatives in 165 families). To minimize genetic and phenotypic heterogeneity, we included only ethnic Chinese Han population in Taiwan and excluded subjects with hypothyroidism. We performed direct and comprehensive genotyping of six classical HLA loci (HLA-A, -B, -C, -DPB1, -DQB1 and -DRB1) to 4-digit resolution. Combining the data of two sample populations, we found that B*46:01 (odds ratio under dominant model [OR]  = 1.33, Bonferroni corrected combined P [P_Bc]  = 1.17×10⁻²), DPB1*05:01 (OR  = 2.34, P_Bc = 2.58×10⁻¹⁰), DQB1*03:02 (OR  = 0.62, P_Bc  = 1.97×10⁻²), DRB1*15:01 (OR  = 1.68, P_Bc = 1.22×10⁻²) and DRB1*16:02 (OR  = 2.63, P_Bc  = 1.46×10⁻⁵) were associated with GD. HLA-DPB1*05:01 is the major gene of GD in our population and singly accounts for 48.4% of population-attributable risk.

Conclusions/Significance

These GD-associated alleles we identified in ethnic Chinese Hans, and those identified in other Asian studies, are totally distinct from the known associated alleles in Caucasians. Identification of population-specific association alleles is the critical first step for individualized medicine. Furthermore, comparison between different susceptibility/protective alleles across populations could facilitate generation of novel hypothesis about GD pathophysiology and indicate a new direction for future investigation.     

Genetic Variants in Meiotic Program Initiation Pathway Genes Are Associated with Spermatogenic Impairment in a Han Chinese Population

Background

The meiotic program initiation pathway genes (CYP26B1, NANOS1 and STRA8) have been proposed to play key roles in spermatogenesis.

Objective

To elucidate the exact role of the genetic variants of the meiosis initiation genes in spermatogenesis, we genotyped the potential functional genetic variants of CYP26B1, NANOS1 and STRA8 genes, and evaluated their effects on spermatogenesis in our study population.

Design, Setting, and Participants

In this study, all subjects were volunteers from the affiliated hospitals of Nanjing Medical University between March 2004 and July 2009 (NJMU Infertile Study). Total 719 idiopathic infertile cases were recruited and divided into three groups according to WHO semen parameters: 201 azoospermia patients (no sperm in the ejaculate even after centrifugation), 155 oligozoospermia patients (sperm counts <20×10⁶/ml) and 363 infertility/normozoospermia subjects (sperm counts >20×10⁶/ml). The control group consisted of 383 subjects with normal semen parameters, all of which had fathered at least one child without assisted reproductive technologies.

Measurements

Eight single nucleotide polymorphisms (SNPs) in CYP26B1, NANOS1 and STRA8 genes were determined by TaqMan allelic discrimination assay in 719 idiopathic infertile men and 383 healthy controls.

Results and Limitations

The genetic variant rs10269148 of STRA8 gene showed higher risk of spermatogenic impairment in the groups of abnormospermia (including azoospermia subgroup and oligozoospermia subgroup) and azoospermia than the controls with odds ratios and 95% confidence intervals of 2.52 (1.29–4.94) and 2.92 (1.41–6.06), respectively (P = 0.006, 0.002 respective). Notably, larger sample size studies and in vivo or in vitro functional studies are needed to substantiate the biological roles of these variants.

Conclusions

Our results provided epidemiological evidence supporting the involvement of genetic polymorphisms of the meiotic program initiation genes in modifying the risk of azoospermia and oligozoospermia in a Han-Chinese population.     

Credibility Analysis of Putative Disease-Causing Genes Using Bioinformatics

Background

Genetic studies are challenging in many complex diseases, particularly those with limited diagnostic certainty, low prevalence or of old age. The result is that genes may be reported as disease-causing with varying levels of evidence, and in some cases, the data may be so limited as to be indistinguishable from chance findings. When there are large numbers of such genes, an objective method for ranking the evidence is useful. Using the neurodegenerative and complex disease amyotrophic lateral sclerosis (ALS) as a model, and the disease-specific database ALSoD, the objective is to develop a method using publicly available data to generate a credibility score for putative disease-causing genes.

Methods

Genes with at least one publication suggesting involvement in adult onset familial ALS were collated following an exhaustive literature search. SQL was used to generate a score by extracting information from the publications and combined with a pathogenicity analysis using bioinformatics tools. The resulting score allowed us to rank genes in order of credibility. To validate the method, we compared the objective ranking with a rank generated by ALS genetics experts. Spearman''s Rho was used to compare rankings generated by the different methods.

Results

The automated method ranked ALS genes in the following order: SOD1, TARDBP, FUS, ANG, SPG11, NEFH, OPTN, ALS2, SETX, FIG4, VAPB, DCTN1, TAF15, VCP, DAO. This compared very well to the ranking of ALS genetics experts, with Spearman''s Rho of 0.69 (P = 0.009).

Conclusion

We have presented an automated method for scoring the level of evidence for a gene being disease-causing. In developing the method we have used the model disease ALS, but it could equally be applied to any disease in which there is genotypic uncertainty.     

Mutation Prevalence of Cerebral Cavernous Malformation Genes in Spanish Patients

Objective

To study the molecular genetic and clinical features of cerebral cavernous malformations (CCM) in a cohort of Spanish patients.

Methods

We analyzed the CCM1, CCM2, and CCM3 genes by MLPA and direct sequencing of exons and intronic boundaries in 94 familial forms and 41 sporadic cases of CCM patients of Spanish extraction. When available, RNA studies were performed seeking for alternative or cryptic splicing.

Results

A total of 26 pathogenic mutations, 22 of which predict truncated proteins, were identified in 29 familial forms and in three sporadic cases. The repertoire includes six novel non-sense and frameshift mutations in CCM1 and CCM3. We also found four missense mutations, one of them located at the third NPXY motif of CCM1 and another one that leads to cryptic splicing of CCM1 exon 6. We found four genomic deletions with the loss of the whole CCM2 gene in one patient and a partial loss of CCM1and CCM2 genes in three other patients. Four families had mutations in CCM3. The results include a high frequency of intronic variants, although most of them localize out of consensus splicing sequences. The main symptoms associated to clinical debut consisted of cerebral haemorrhage, migraines and epileptic seizures. The rare co-occurrence of CCM with Noonan and Chiari syndromes and delayed menarche is reported.

Conclusions

Analysis of CCM genes by sequencing and MLPA has detected mutations in almost 35% of a Spanish cohort (36% of familial cases and 10% of sporadic patients). The results include 13 new mutations of CCM genes and the main clinical symptoms that deserves consideration in molecular diagnosis and genetic counselling of cerebral cavernous malformations.     

Shared epitope alleles remain a risk factor for anti-citrullinated proteins antibody (ACPA)--positive rheumatoid arthritis in three Asian ethnic groups

Background

To investigate the associations between HLA-DRB1 shared epitope (SE) alleles and rheumatoid arthritis in subsets of rheumatoid arthritis defined by autoantibodies in three Asian populations from Malaysia.

Methods

1,079 rheumatoid arthritis patients and 1,470 healthy controls were included in the study. Levels of antibodies to citrullinated proteins (ACPA) and rheumatoid factors were assessed and the PCR-SSO method was used for HLA-DRB1 genotyping.

Results

The proportion of ACPA positivity among Malay, Chinese and Indian rheumatoid arthritis patients were 62.9%, 65.2% and 68.6%, respectively. An increased frequency of SE alleles was observed in ACPA-positive rheumatoid arthritis among the three Asian ethnic groups. HLA-DRB1*10 was highly associated with rheumatoid arthritis susceptibility in these Asian populations. HLA-DRB1*0405 was significantly associated with susceptibility to rheumatoid arthritis in Malays and Chinese, but not in Indians. HLA-DRB1*01 did not show any independent effect as a risk factor for rheumatoid arthritis in this study and HLA-DRB1*1202 was protective in Malays and Chinese. There was no association between SE alleles and ACPA- negative rheumatoid arthritis in any of the three Asian ethnic groups.

Conclusion

The HLA-DRB1 SE alleles increase the risk of ACPA-positive rheumatoid arthritis in all three Asian populations from Malaysia.     

BRCA1 and BRCA2 germline mutations in Malaysian women with early-onset breast cancer without a family history

Background

In Asia, breast cancer is characterised by an early age of onset: In Malaysia, approximately 50% of cases occur in women under the age of 50 years. A proportion of these cases may be attributable, at least in part, to genetic components, but to date, the contribution of genetic components to breast cancer in many of Malaysia''s ethnic groups has not been well-characterised.

Methodology

Given that hereditary breast carcinoma is primarily due to germline mutations in one of two breast cancer susceptibility genes, BRCA1 and BRCA2, we have characterised the spectrum of BRCA mutations in a cohort of 37 individuals with early-onset disease (≤40 years) and no reported family history. Mutational analysis of BRCA1 and BRCA2 was conducted by full sequencing of all exons and intron-exon junctions.

Conclusions

Here, we report a total of 14 BRCA1 and 17 BRCA2 sequence alterations, of which eight are novel (3 BRCA1 and 5 BRCA2). One deleterious BRCA1 mutation and 2 deleterious BRCA2 mutations, all of which are novel mutations, were identified in 3 of 37 individuals. This represents a prevalence of 2.7% and 5.4% respectively, which is consistent with other studies in other Asian ethnic groups (4–9%).     

Anthocyanin biosynthetic genes in Brassica rapa

Background

Anthocyanins are a group of flavonoid compounds. As a group of important secondary metabolites, they perform several key biological functions in plants. Anthocyanins also play beneficial health roles as potentially protective factors against cancer and heart disease. To elucidate the anthocyanin biosynthetic pathway in Brassica rapa, we conducted comparative genomic analyses between Arabidopsis thaliana and B. rapa on a genome-wide level.

Results

In total, we identified 73 genes in B. rapa as orthologs of 41 anthocyanin biosynthetic genes in A. thaliana. In B. rapa, the anthocyanin biosynthetic genes (ABGs) have expanded and most genes exist in more than one copy. The anthocyanin biosynthetic structural genes have expanded through whole genome and tandem duplication in B. rapa. More structural genes located upstream of the anthocyanin biosynthetic pathway have been retained than downstream. More negative regulatory genes are retained in the anthocyanin biosynthesis regulatory system of B. rapa.

Conclusions

These results will promote an understanding of the genetic mechanism of anthocyanin biosynthesis, as well as help the improvement of the nutritional quality of B. rapa through the breeding of high anthocyanin content varieties.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-426) contains supplementary material, which is available to authorized users.