The VNTR in complex disorders: The forgotten polymorphisms? A functional way forward?

Mechanistic Characterization of De Novo Generation of Variable Number Tandem Repeats in Circular Plasmids during Site-Directed Mutagenesis and Optimization for Coding Gene Application

Site-directed mutagenesis for creating point mutations, sometimes, gives rise to plasmids carrying variable number tandem repeats (VNTRs) locally, which are arbitrarily regarded as polymerase chain reaction (PCR) related artifacts. Here, the alternative end-joining mechanism is reported rather than PCR artifacts accounts largely for that VNTRs formation and expansion. During generating a point mutation on GPLD1 gene, an unexpected formation of VNTRs employing the 31 bp mutagenesis primers is observed as the repeat unit in the pcDNA3.1-GPLD1 plasmid. The 31 bp VNTRs are formed in 24.75% of the resulting clones with copy number varied from 2 to 13. All repeat units are aligned with the same orientation as GPLD1 gene. 43.54% of the repeat junctions harbor nucleotide mutations while the rest don't. Their demonstrated short primers spanning the 3' part of the mutagenesis primers are essential for initial creation of the 2-copy tandem repeats (TRs) in circular plasmids. The dimerization of mutagenesis primers by the alternative end-joining in a correct orientation is required for further expansion of the 2-copy TRs. Lastly, a half-double priming strategy is established, verified the findings and offered a simple method for VNTRs creation on coding genes in circular plasmids without junction mutations.

Description of a New Simple and Cost-Effective Molecular Testing That Could Simplify MUC1 Variant Detection

Introduction

Patients with autosomal dominant tubulointerstitial kidney disease (ADTKD) usually present with nonspecific progressive chronic kidney disease (CKD) with mild to negative proteinuria and a family history. ADTKD-MUC1 leads to the formation of a frameshift protein that accumulates in the cytoplasm, leading to tubulointerstitial damage. ADTKD-MUC1 prevalence remains unclear because MUC1 variants are not routinely detected by standard next-generation sequencing (NGS) techniques.

Methods

We developed a bioinformatic counting script that can detect specific genetic sequences and count the number of occurrences. We used DNA samples from 27 patients for validation, 11 of them were patients from the Lille University Hospital in France and 16 were from the Wake Forest Hospital, NC. All patients from Lille were tested with an NGS gene panel with our script and all patients from Wake Forest Hospital were tested with the snapshot reference technique. Between January 2018 and February 2023, we collected data on all patients diagnosed with MUC1 variants with this script.

Results

A total of 27 samples were tested anonymously by the BROAD Institute reference technique for confirmation and we were able to get a 100% concordance for MUC1 diagnosis. Clinico-biologic characteristics in our cohort were similar to those previously described in ADTKD-MUC1.

Conclusion

We describe a new simple and cost-effective method for molecular testing of ADTKD-MUC1. Genetic analyses in our cohort suggest that MUC1 might be the first cause of ADTKD. Increasing the availability of MUC1 diagnosis tools will contribute to a better understanding of the disease and to the development of specific treatments.

Exploring the Relationship between CLPTM1L-MS2 Variants and Susceptibility to Bladder Cancer

CLPTM1L (Cleft Lip and Palate Transmembrane Protein 1-Like) has previously been implicated in tumorigenesis and drug resistance in cancer. However, the genetic link between CLPTM1L and bladder cancer remains uncertain. In this study, we investigated the genetic association of variable number of tandem repeats (VNTR; minisatellites, MS) regions within CLPTM1L with bladder cancer. We identified four CLPTM1L-MS regions (MS1~MS4) located in intron regions. To evaluate the VNTR polymorphic alleles, we analyzed 441 cancer-free controls and 181 bladder cancer patients. Our analysis revealed a higher frequency of specific repeat sizes within the MS2 region in bladder cancer cases compared to controls. Notably, 25 and 27 repeats were exclusively present in the bladder cancer group. Moreover, rare alleles within the medium-length repeat range (25-29 repeats) were associated with an elevated bladder cancer risk (odds ratio [OR] = 5.78, 95% confidence interval [CI]: 1.49-22.47, p = 0.004). We confirmed that all MS regions followed Mendelian inheritance, and demonstrated that MS2 alleles increased CLPTM1L promoter activity in the UM-UC3 bladder cancer cells through a luciferase assay. Our findings propose the utility of CLPTM1L-MS regions as DNA typing markers, particularly highlighting the potential of middle-length rare alleles within CLPTM1L-MS2 as predictive markers for bladder cancer risk.

Investigating the relationship between the VNTR variant of the interleukin-1 receptor antagonist gene and coronary in-stent restenosis

OBJECTIVE

This study aimed to examine the association between the interleukin-1 receptor antagonist gene (IL-1RN) and coronary in-stent restenosis (ISR) through the analysis of the VNTR variant based on the previously reported results.

MATERIALS AND METHODS

The samples were classified into two clearly defined groups: the case group, which comprised 45 patients diagnosed with in-stent restenosis (ISR+), and the control group, which included 60 patients without ISR (ISR-). Polymerase chain reaction (PCR) was performed to examine the 86-bp VNTR variant of the IL-1RN gene.

RESULTS

In the analysis of six identified groups consisting of variant alleles of 86 base pairs of VNTR of the IL-1RN gene statistically significant difference was observed for the presence of IL1RN*2 allele between cases and controls (p = 0.04, OR; 0.045).

CONCLUSION

Individuals with allele 2 of the IL-1Ra gene may be more predisposed to ISR. This could be due to an imbalance between IL-1Ra and IL-1β which is crucial in preventing the initiation or advancement of inflammatory diseases in specific organs. The observed phenomenon can be characterized by increased production of IL-1β and potential reduction of IL-1Ra as a result of functional VNTR variation in IL-RN gene.

Variable number tandem repeats (VNTRs) as modifiers of breast cancer risk in carriers of BRCA1 185delAG

Despite substantial efforts in identifying both rare and common variants affecting disease risk, in the majority of diseases, a large proportion of unexplained genetic risk remains. We propose that variable number tandem repeats (VNTRs) may explain a proportion of the missing genetic risk. Herein, in a pilot study with a retrospective cohort design, we tested whether VNTRs are causal modifiers of breast cancer risk in 347 female carriers of the BRCA1 185delAG pathogenic variant, an important group given their high risk of developing breast cancer. We performed targeted-capture to sequence VNTRs, called genotypes with adVNTR, tested the association of VNTRs and breast cancer risk using Cox regression models, and estimated the effect size using a retrospective likelihood approach. Of 303 VNTRs that passed quality control checks, 4 VNTRs were significantly associated with risk to develop breast cancer at false discovery rate [FDR] < 0.05 and an additional 4 VNTRs had FDR < 0.25. After determining the specific risk alleles, there was a significantly earlier age at diagnosis of breast cancer in carriers of the risk alleles compared to those without the risk alleles for seven of eight VNTRs. One example is a VNTR in exon 2 of LINC01973 with a per-allele hazard ratio of 1.58 (1.07-2.33) and 5.28 (2.79-9.99) for the homozygous risk-allele genotype. Results from this first systematic study of VNTRs demonstrate that VNTRs may explain a proportion of the unexplained genetic risk for breast cancer.

Detecting tandem repeat variants in coding regions using code-adVNTR

The human genome contains more than one million tandem repeats (TRs), DNA sequences containing multiple approximate copies of a motif repeated contiguously. TRs account for significant genetic variation, with 50 + diseases attributed to changes in motif number. A few diseases have been to be caused by small indels in variable number tandem repeats (VNTRs) including poly-cystic kidney disease type 1 (MCKD1) and monogenic type 1 diabetes. However, small indels in VNTRs are largely unexplored mainly due to the long and complex structure of VNTRs with multiple motifs. We developed a method, code-adVNTR, that utilizes multi-motif hidden Markov models to detect both, motif count variation and small indels, within VNTRs. In simulated data, code-adVNTR outperformed GATK-HaplotypeCaller in calling small indels within large VNTRs. We used code-adVNTR to characterize coding VNTRs in the 1000 genomes data identifying many population-specific variants, and to reliably call MUC1 mutations for MCKD1.

•

Detection of coding variants in tandem repeats is confounded by ambiguous mapping

•

Our method, code-adVNTR, detects variants in coding VNTRs using multi-motif HMMs

•

code-adVNTR outperforms GATK-HaplotypeCaller on indel detection in tandem repeats

•

A known frameshift variant within a VNTR in MUC1 gene was accurately detected

A hidden layer of structural variation in transposable elements reveals potential genetic modifiers in human disease-risk loci

Genome-wide association studies (GWAS) have been highly informative in discovering disease-associated loci but are not designed to capture all structural variations in the human genome. Using long-read sequencing data, we discovered widespread structural variation within SINE-VNTR-Alu (SVA) elements, a class of great ape-specific transposable elements with gene-regulatory roles, which represents a major source of structural variability in the human population. We highlight the presence of structurally variable SVAs (SV-SVAs) in neurological disease-associated loci, and we further associate SV-SVAs to disease-associated SNPs and differential gene expression using luciferase assays and expression quantitative trait loci data. Finally, we genetically deleted SV-SVAs in the BIN1 and CD2AP Alzheimer's disease-associated risk loci and in the BCKDK Parkinson's disease-associated risk locus and assessed multiple aspects of their gene-regulatory influence in a human neuronal context. Together, this study reveals a novel layer of genetic variation in transposable elements that may contribute to identification of the structural variants that are the actual drivers of disease associations of GWAS loci.

Chronic pain and delinquency partially explain the effect of the DRD4 gene polymorphism on adult substance use

Background: The dopamine receptor D4 [DRD4] has been reported to be associated with substance use. Yet, the roles that health conditions and behaviors may play in such association are understudied.Objective: This longitudinal study investigated the potential mediation effects of chronic pain and delinquency in adolescence on the association between the DRD4 2-repeat allele and substance use in adulthood. Sex, witnessing violence, and experiencing violence were also examined as potential moderators for the mediation pathways.Methods: We used the restricted and candidate gene data from the National Longitudinal Study of Adolescent to Adult Health (Waves I-IV) to conduct secondary analysis (N = 8,671; 47% male). A two-step approach was adopted to examine the mediation effects regarding four substance use outcomes in adulthood: number of lifetime alcohol use disorder symptoms, lifetime regular smoker status, past-month smoking, and lifetime "pain killer" misuse. The moderation effects were investigated using stratification and permutation.Results: The DRD4 2-repeat allele was associated with all adulthood substance use outcomes through adolescent chronic pain and delinquency (AORs/IRR range 1.08-3.78; all ps<0.01). The association between delinquency and smoking was higher among females. The association between delinquency and substance use was lower among the participants who witnessed violence in adolescence.Conclusions: This study identified modifiable mediators underlying the association between the DRD4 2-repeat allele and substance use behaviors, concluding that chronic pain and delinquency partially explain the effect of the DRD4 gene polymorphism on adult substance use.

Short Tandem Repeat Variation in the CNR1 Gene Associated With Analgesic Requirements of Opioids in Postoperative Pain Management

Short tandem repeats (STRs) and variable number of tandem repeats (VNTRs) that have been identified at approximately 0.7 and 0.5 million loci in the human genome, respectively, are highly multi-allelic variations rather than single-nucleotide polymorphisms. The number of repeats of more than a few thousand STRs was associated with the expression of nearby genes, indicating that STRs are influential genetic variations in human traits. Analgesics act on the central nervous system via their intrinsic receptors to produce analgesic effects. In the present study, we focused on STRs and VNTRs in the CNR1, GRIN2A, PENK, and PDYN genes and analyzed two peripheral pain sensation-related traits and seven analgesia-related traits in postoperative pain management. A total of 192 volunteers who underwent the peripheral pain sensation tests and 139 and 252 patients who underwent open abdominal and orthognathic cosmetic surgeries, respectively, were included in the study. None of the four STRs or VNTRs were associated with peripheral pain sensation. Short tandem repeats in the CNR1, GRIN2A, and PENK genes were associated with the frequency of fentanyl use, fentanyl dose, and visual analog scale pain scores 3 h after orthognathic cosmetic surgery (Spearman’s rank correlation coefficient ρ = 0.199, p = 0.002, ρ = 0.174, p = 0.006, and ρ = 0.135, p = 0.033, respectively), analgesic dose, including epidural analgesics after open abdominal surgery (ρ = −0.200, p = 0.018), and visual analog scale pain scores 24 h after orthognathic cosmetic surgery (ρ = 0.143, p = 0.023), respectively. The associations between STRs in the CNR1 gene and the frequency of fentanyl use and fentanyl dose after orthognathic cosmetic surgery were confirmed by Holm’s multiple-testing correction. These findings indicate that STRs in the CNR1 gene influence analgesia in the orofacial region.

Revisiting tandem repeats in psychiatric disorders from perspectives of genetics, physiology, and brain evolution

Genome-wide association studies (GWASs) have revealed substantial genetic components comprised of single nucleotide polymorphisms (SNPs) in the heritable risk of psychiatric disorders. However, genetic risk factors not covered by GWAS also play pivotal roles in these illnesses. Tandem repeats, which are likely functional but frequently overlooked by GWAS, may account for an important proportion in the "missing heritability" of psychiatric disorders. Despite difficulties in characterizing and quantifying tandem repeats in the genome, studies have been carried out in an attempt to describe impact of tandem repeats on gene regulation and human phenotypes. In this review, we have introduced recent research progress regarding the genomic distribution and regulatory mechanisms of tandem repeats. We have also summarized the current knowledge of the genetic architecture and biological underpinnings of psychiatric disorders brought by studies of tandem repeats. These findings suggest that tandem repeats, in candidate psychiatric risk genes or in different levels of linkage disequilibrium (LD) with psychiatric GWAS SNPs and haplotypes, may modulate biological phenotypes related to psychiatric disorders (e.g., cognitive function and brain physiology) through regulating alternative splicing, promoter activity, enhancer activity and so on. In addition, many tandem repeats undergo tight natural selection in the human lineage, and likely exert crucial roles in human brain evolution. Taken together, the putative roles of tandem repeats in the pathogenesis of psychiatric disorders is strongly implicated, and using examples from previous literatures, we wish to call for further attention to tandem repeats in the post-GWAS era of psychiatric disorders.

A multicenter case-control study of the effect of e-nos VNTR polymorphism on upper gastrointestinal hemorrhage in NSAID users

Bleeding in non-steroidal anti-inflammatory drug (NSAID) users limited their prescription. This first multicenter full case–control study (325 cases and 744 controls), explored the association of e-NOS intron 4 variable number tandem repeat (VNTR) polymorphism with upper gastrointestinal hemorrhage (UGIH) in NSAID exposed and unexposed populations and assessed any interaction between this polymorphism and NSAIDs. NSAID users carrying e-NOS intron 4 wild type genotype or VNTR polymorphism have higher odds of UGIH than those unexposed to NSAIDs [Odds Ratio (OR): 6.62 (95% Confidence Interval (CI): 4.24, 10.36) and OR: 5.41 (95% CI 2.62, 11.51), respectively], with no effect modification from VNTR polymorphism-NSAIDs interaction [Relative Excess Risk due to Interaction (RERI): −1.35 (95% CI −5.73, 3.03); Synergism Index (S): 0.77 (95% CI 0.31, 1.94)]. Similar findings were obtained for aspirin exposure. Non-aspirin NSAID users who carry e-NOS intron 4 VNTR polymorphism have lower odds of UGIH [OR: 4.02 (95% CI 1.85, 8.75) than those users with wild type genotype [OR: 6.52 (95% CI 4.09, 10.38)]; though the interaction estimates are not statistically significant [RERI: −2.68 (95% CI −6.67, 1.31); S: 0.53 (95% CI 0.18, 1.55)]. This exploratory study suggests that the odds of UGIH in NSAID or aspirin users does not modify according to patient´s e-NOS intron 4 genotype.

The position of single-base deletions in the VNTR sequence of the carboxyl ester lipase (CEL) gene determines proteotoxicity

Variable number of tandem repeat (VNTR) sequences in the genome can have functional consequences that contribute to human disease. This is the case for the CEL gene, which is specifically expressed in pancreatic acinar cells and encodes the digestive enzyme carboxyl ester lipase. Rare single-base deletions (DELs) within the first (DEL1) or fourth (DEL4) VNTR segment of CEL cause maturity-onset diabetes of the young, type 8 (MODY8), an inherited disorder characterized by exocrine pancreatic dysfunction and diabetes. Studies on the DEL1 variant have suggested that MODY8 is initiated by CEL protein misfolding and aggregation. However, it is unclear how the position of single-base deletions within the CEL VNTR affects pathogenic properties of the protein. Here, we investigated four naturally occurring CEL variants, arising from single-base deletions in different VNTR segments (DEL1, DEL4, DEL9, and DEL13). When the four variants were expressed in human embryonic kidney 293 cells, only DEL1 and DEL4 led to significantly reduced secretion, increased intracellular aggregation, and increased endoplasmic reticulum stress compared with normal CEL protein. The level of O-glycosylation was affected in all DEL variants. Moreover, all variants had enzymatic activity comparable with that of normal CEL. We conclude that the longest aberrant protein tails, resulting from single-base deletions in the proximal VNTR segments, have highest pathogenic potential, explaining why DEL1 and DEL4 but not DEL9 and DEL13 have been observed in patients with MODY8. These findings further support the view that CEL mutations cause pancreatic disease through protein misfolding and proteotoxicity, leading to endoplasmic reticulum stress and activation of the unfolded protein response.

Variable number tandem repeats mediate the expression of proximal genes

Variable number tandem repeats (VNTRs) account for significant genetic variation in many organisms. In humans, VNTRs have been implicated in both Mendelian and complex disorders, but are largely ignored by genomic pipelines due to the complexity of genotyping and the computational expense. We describe adVNTR-NN, a method that uses shallow neural networks to genotype a VNTR in 18 seconds on 55X whole genome data, while maintaining high accuracy. We use adVNTR-NN to genotype 10,264 VNTRs in 652 GTEx individuals. Associating VNTR length with gene expression in 46 tissues, we identify 163 "eVNTRs". Of the 22 eVNTRs in blood where independent data is available, 21 (95%) are replicated in terms of significance and direction of association. 49% of the eVNTR loci show a strong and likely causal impact on the expression of genes and 80% have maximum effect size at least 0.3. The impacted genes are involved in diseases including Alzheimer's, obesity and familial cancers, highlighting the importance of VNTRs for understanding the genetic basis of complex diseases.

DRD4 genotyping may differentiate symptoms of attention-deficit/hyperactivity disorder and sluggish cognitive tempo

OBJECTIVE

Studies to reduce the heterogeneity of attention-deficit/hyperactivity disorder (ADHD) have increased interest in the concept of sluggish cognitive tempo (SCT). The aim of this study was to investigate if the prevalence of two variable-number tandem repeats (VNTRs) located within the 3'-untranslated region of the DAT1 gene and in exon 3 of the dopamine D4 receptor (DRD4) gene differ among four groups (31 subjects with SCT but no ADHD, 146 individuals with ADHD but no SCT, 67 subjects with SCT + ADHD, and 92 healthy controls).

METHODS

We compared the sociodemographic profiles, neurocognitive domains, and prevalence of two VNTRs in SCT and ADHD subjects versus typically developing (TD) controls.

RESULTS

The SCT without ADHD group had a higher proportion of females and lower parental educational attainment. Subjects in this group performed worse on neuropsychological tests, except for psychomotor speed and commission errors, compared to controls. However, the ADHD without SCT group performed significantly worse on all neuropsychological domains than controls. We found that 4R homozygosity for the DRD4 gene was most prevalent in the ADHD without SCT group. The SCT without ADHD group had the highest 7R allele frequency, differing significantly from the ADHD without SCT group.

CONCLUSION

The 7R allele of DRD4 gene was found to be significantly more prevalent in SCT cases than in ADHD cases. No substantial neuropsychological differences were found between SCT and ADHD subjects.

Extreme enrichment of VNTR-associated polymorphicity in human subtelomeres: genes with most VNTRs are predominantly expressed in the brain

The human genome harbors numerous structural variants (SVs) which, due to their repetitive nature, are currently underexplored in short-read whole-genome sequencing approaches. Using single-molecule, real-time (SMRT) long-read sequencing technology in combination with FALCON-Unzip, we generated a de novo assembly of the diploid genome of a 115-year-old Dutch cognitively healthy woman. We combined this assembly with two previously published haploid assemblies (CHM1 and CHM13) and the GRCh38 reference genome to create a compendium of SVs that occur across five independent human haplotypes using the graph-based multi-genome aligner REVEAL. Across these five haplotypes, we detected 31,680 euchromatic SVs (>50 bp). Of these, ~62% were comprised of repetitive sequences with 'variable number tandem repeats' (VNTRs), ~10% were mobile elements (Alu, L1, and SVA), while the remaining variants were inversions and indels. We observed that VNTRs with GC-content >60% and repeat patterns longer than 15 bp were 21-fold enriched in the subtelomeric regions (within 5 Mb of the ends of chromosome arms). VNTR lengths can expand to exceed a critical length which is associated with impaired gene transcription. The genes that contained most VNTRs, of which PTPRN2 and DLGAP2 are the most prominent examples, were found to be predominantly expressed in the brain and associated with a wide variety of neurological disorders. Repeat-induced variation represents a sizeable fraction of the genetic variation in human genomes and should be included in investigations of genetic factors associated with phenotypic traits, specifically those associated with neurological disorders. We make available the long and short-read sequence data of the supercentenarian genome, and a compendium of SVs as identified across 5 human haplotypes.

Alzheimer Disease Pathology-Associated Polymorphism in a Complex Variable Number of Tandem Repeat Region Within the MUC6 Gene, Near the AP2A2 Gene

We found evidence of late-onset Alzheimer disease (LOAD)-associated genetic polymorphism within an exon of Mucin 6 (MUC6) and immediately downstream from another gene: Adaptor Related Protein Complex 2 Subunit Alpha 2 (AP2A2). PCR analyses on genomic DNA samples confirmed that the size of the MUC6 variable number tandem repeat (VNTR) region was highly polymorphic. In a cohort of autopsied subjects with quantitative digital pathology data (n = 119), the size of the polymorphic region was associated with the severity of pTau pathology in neocortex. In a separate replication cohort of autopsied subjects (n = 173), more pTau pathology was again observed in subjects with longer VNTR regions (p = 0.031). Unlike MUC6, AP2A2 is highly expressed in human brain. AP2A2 expression was lower in a subset analysis of brain samples from persons with longer versus shorter VNTR regions (p = 0.014 normalizing with AP2B1 expression). Double-label immunofluorescence studies showed that AP2A2 protein often colocalized with neurofibrillary tangles in LOAD but was not colocalized with pTau proteinopathy in progressive supranuclear palsy, or with TDP-43 proteinopathy. In summary, polymorphism in a repeat-rich region near AP2A2 was associated with neocortical pTau proteinopathy (because of the unique repeats, prior genome-wide association studies were probably unable to detect this association), and AP2A2 was often colocalized with neurofibrillary tangles in LOAD.

Determining the association between polymorphisms of the DAT1 and DRD4 genes with attention deficit hyperactivity disorder in children from Java Island

Attention deficit hyperactivity disorder (ADHD) is one of the most common neurobehavioural in the children. Genetic factor is known one of the factors which contributed in ADHD development. VNTR polymorphism in 3'UTR exon 15 of DAT1 gene and exon 3 of DRD4 gene are reported to be associated in ADHD. In this study we examine the association of ADHD with VNTR polymorphism of DAT1 and DRD4 gene in Indonesian children. Sixty-five ADHD children and 70 normal children (6-13 years of age), were included in the study, we matched by age and gender. ADHD was diagnosed by DSM-IV. We performed a casecontrol study to found the association between ADHD and VNTR polymorphism of DAT1 and DRD4 genes. The 10-repeat allele of DAT1 and 2-repeat allele of DRD4 were higher in Indonesian children. Although the frequency of these allele was higher, but it was similar both in ADHD and control groups. Neither DAT1 nor DRD4 gene showed showed significant difference in genotype distribution and frequency allele between both groups (p > 0.05). No association between ADHD and VNTR polymorphism of DAT1 and DRD4 genes found in Indonesian children. This data suggest that DAT1 and DRD4 do not contribute to etiology of ADHD in Indonesian children. Further studies are needed to clarify association between VNTR polymorphism of DAT1 and DRD4 genetic with ADHD of Indonesian children in larger sample size and family based study.

The MUC6/AP2A2 Locus and Its Relevance to Alzheimer's Disease: A Review

We recently reported evidence of Alzheimer's disease (AD)-linked genetic variation within the mucin 6 (MUC6) gene on chromosome 11p, nearby the adaptor-related protein complex 2 subunit alpha 2 (AP2A2) gene. This locus has interesting features related to human genomics and clinical research. MUC6 gene variants have been reported to potentially influence viral-including herpesvirus-immunity and the gut microbiome. Within the MUC6 gene is a unique variable number of tandem repeat (VNTR) region. We discovered an association between MUC6 VNTR repeat expansion and AD pathologic severity, particularly tau proteinopathy. Here, we review the relevant literature. The AD-linked VNTR polymorphism may also influence AP2A2 gene expression. AP2A2 encodes a polypeptide component of the adaptor protein complex, AP-2, which is involved in clathrin-coated vesicle function and was previously implicated in AD pathogenesis. To provide background information, we describe some key knowledge gaps in AD genetics research. The "missing/hidden heritability problem" of AD is highlighted. Extensive portions of the human genome, including the MUC6 VNTR, have not been thoroughly evaluated due to limitations of existing high-throughput sequencing technology. We present and discuss additional data, along with cautionary considerations, relevant to the hypothesis that MUC6 repeat expansion influences AD pathogenesis.

A human minisatellite hosts an alternative transcription start site for NPRL3 driving its expression in a repeat number-dependent manner

Minisatellites, also called variable number of tandem repeats (VNTRs), are a class of repetitive elements that may affect gene expression at multiple levels and have been correlated to disease. Their identification and role as expression quantitative trait loci (eQTL) have been limited by their absence in comparative genomic hybridization and single nucleotide polymorphisms arrays. By taking advantage of cap analysis of gene expression (CAGE), we describe a new example of a minisatellite hosting a transcription start site (TSS) which expression is dependent on the repeat number. It is located in the third intron of the gene nitrogen permease regulator like protein 3 (NPRL3). NPRL3 is a component of the GAP activity toward rags 1 protein complex that inhibits mammalian target of rapamycin complex 1 (mTORC1) activity and it is found mutated in familial focal cortical dysplasia and familial focal epilepsy. CAGE tags represent an alternative TSS identifying TAGNPRL3 messenger RNAs (mRNAs). TAGNPRL3 is expressed in red blood cells both at mRNA and protein levels, it interacts with its protein partner NPRL2 and its overexpression inhibits cell proliferation. This study provides an example of a minisatellite that is both a TSS and an eQTL as well as identifies a new VNTR that may modify mTORC1 activity.

NanoSatellite: accurate characterization of expanded tandem repeat length and sequence through whole genome long-read sequencing on PromethION

Technological limitations have hindered the large-scale genetic investigation of tandem repeats in disease. We show that long-read sequencing with a single Oxford Nanopore Technologies PromethION flow cell per individual achieves 30× human genome coverage and enables accurate assessment of tandem repeats including the 10,000-bp Alzheimer's disease-associated ABCA7 VNTR. The Guppy "flip-flop" base caller and tandem-genotypes tandem repeat caller are efficient for large-scale tandem repeat assessment, but base calling and alignment challenges persist. We present NanoSatellite, which analyzes tandem repeats directly on electric current data and improves calling of GC-rich tandem repeats, expanded alleles, and motif interruptions.

The heterogeneous pathogenesis of type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM) results from the destruction of pancreatic β-cells that is mediated by the immune system. Multiple genetic and environmental factors found in variable combinations in individual patients are involved in the development of T1DM. Genetic risk is defined by the presence of particular allele combinations, which in the major susceptibility locus (the HLA region) affect T cell recognition and tolerance to foreign and autologous molecules. Multiple other loci also regulate and affect features of specific immune responses and modify the vulnerability of β-cells to inflammatory mediators. Compared with the genetic factors, environmental factors that affect the development of T1DM are less well characterized but contact with particular microorganisms is emerging as an important factor. Certain infections might affect immune regulation, and the role of commensal microorganisms, such as the gut microbiota, are important in the education of the developing immune system. Some evidence also suggests that nutritional factors are important. Multiple islet-specific autoantibodies are found in the circulation from a few weeks to up to 20 years before the onset of clinical disease and this prediabetic phase provides a potential opportunity to manipulate the islet-specific immune response to prevent or postpone β-cell loss. The latest developments in understanding the heterogeneity of T1DM and characterization of major disease subtypes might help in the development of preventive treatments.

No Support for Historical Candidate Gene or Candidate Gene-by-Interaction Hypotheses for Major Depression Across Multiple Large Samples

OBJECTIVE

Interest in candidate gene and candidate gene-by-environment interaction hypotheses regarding major depressive disorder remains strong despite controversy surrounding the validity of previous findings. In response to this controversy, the present investigation empirically identified 18 candidate genes for depression that have been studied 10 or more times and examined evidence for their relevance to depression phenotypes.

METHODS

Utilizing data from large population-based and case-control samples (Ns ranging from 62,138 to 443,264 across subsamples), the authors conducted a series of preregistered analyses examining candidate gene polymorphism main effects, polymorphism-by-environment interactions, and gene-level effects across a number of operational definitions of depression (e.g., lifetime diagnosis, current severity, episode recurrence) and environmental moderators (e.g., sexual or physical abuse during childhood, socioeconomic adversity).

RESULTS

No clear evidence was found for any candidate gene polymorphism associations with depression phenotypes or any polymorphism-by-environment moderator effects. As a set, depression candidate genes were no more associated with depression phenotypes than noncandidate genes. The authors demonstrate that phenotypic measurement error is unlikely to account for these null findings.

CONCLUSIONS

The study results do not support previous depression candidate gene findings, in which large genetic effects are frequently reported in samples orders of magnitude smaller than those examined here. Instead, the results suggest that early hypotheses about depression candidate genes were incorrect and that the large number of associations reported in the depression candidate gene literature are likely to be false positives.

Imputation of behavioral candidate gene repeat variants in 486,551 publicly-available UK Biobank individuals

Some of the most widely studied variants in psychiatric genetics include variable number tandem repeat variants (VNTRs) in SLC6A3, DRD4, SLC6A4, and MAOA. While initial findings suggested large effects, their importance with respect to psychiatric phenotypes is the subject of much debate with broadly conflicting results. Despite broad interest, these loci remain absent from the largest available samples, such as the UK Biobank, limiting researchers' ability to test these contentious hypotheses rigorously in large samples. Here, using two independent reference datasets, we report out-of-sample imputation accuracy estimates of >0.96 for all four VNTR variants and one modifying SNP, depending on the reference and target dataset. We describe the imputation procedures of these candidate variants in 486,551 UK Biobank individuals, and have made the imputed variant data available to UK Biobank researchers. This resource, provided to the scientific community, will allow the most rigorous tests to-date of the roles of these variants in behavioral and psychiatric phenotypes.

Perspective: DNA Copy Number Variations in Cardiovascular Diseases

Human genome contains many variations, often called mutations, which are difficult to detect and have remained a challenge for years. A substantial part of the genome encompasses repeats and when such repeats are in the coding region they may lead to change in the gene expression profile followed by pathological conditions. Structural variants are alterations which change one or more sequence feature in the chromosome such as change in the copy number, rearrangements, and translocations of a sequence and can be balanced or unbalanced. Copy number variants (CNVs) may increase or decrease the copies of a given region and have a pivotal role in the onset of many diseases including cardiovascular disorders. Cardiovascular disorders have a magnitude of well-established risk factors and etiology, but their correlation with CNVs is still being studied. In this article, we have discussed history of CNVs and a summary on the diseases associated with CNVs. To detect such variations, we shed light on the number of techniques introduced so far and their limitations. The lack of studies on cardiovascular diseases to determine the frequency of such variants needs clinical studies with larger cohorts. This review is a compilation of articles suggesting the importance of CNVs in multitude of cardiovascular anomalies. Finally, future perspectives for better understanding of CNVs and cardiovascular disorders have also been discussed.

Targeted genotyping of variable number tandem repeats with adVNTR

Whole-genome sequencing is increasingly used to identify Mendelian variants in clinical pipelines. These pipelines focus on single-nucleotide variants (SNVs) and also structural variants, while ignoring more complex repeat sequence variants. Here, we consider the problem of genotyping Variable Number Tandem Repeats (VNTRs), composed of inexact tandem duplications of short (6–100 bp) repeating units. VNTRs span 3% of the human genome, are frequently present in coding regions, and have been implicated in multiple Mendelian disorders. Although existing tools recognize VNTR carrying sequence, genotyping VNTRs (determining repeat unit count and sequence variation) from whole-genome sequencing reads remains challenging. We describe a method, adVNTR, that uses hidden Markov models to model each VNTR, count repeat units, and detect sequence variation. adVNTR models can be developed for short-read (Illumina) and single-molecule (Pacific Biosciences [PacBio]) whole-genome and whole-exome sequencing, and show good results on multiple simulated and real data sets.

A Census of Tandemly Repeated Polymorphic Loci in Genic Regions Through the Comparative Integration of Human Genome Assemblies

Polymorphic Tandem Repeat (PTR) is a common form of polymorphism in the human genome. A PTR consists in a variation found in an individual (or in a population) of the number of repeating units of a Tandem Repeat (TR) locus of the genome with respect to the reference genome. Several phenotypic traits and diseases have been discovered to be strongly associated with or caused by specific PTR loci. PTR are further distinguished in two main classes: Short Tandem Repeats (STR) when the repeating unit has size up to 6 base pairs, and Variable Number Tandem Repeats (VNTR) for repeating units of size above 6 base pairs. As larger and larger populations are screened via high throughput sequencing projects, it becomes technically feasible and desirable to explore the association between PTR and a panoply of such traits and conditions. In order to facilitate these studies, we have devised a method for compiling catalogs of PTR from assembled genomes, and we have produced a catalog of PTR for genic regions (exons, introns, UTR and adjacent regions) of the human genome (GRCh38). We applied four different TR discovery software tools to uncover in the first phase 55,223,485 TR (after duplicate removal) in GRCh38, of which 373,173 were determined to be PTR in the second phase by comparison with five assembled human genomes. Of these, 263,266 are not included by state-of-the-art PTR catalogs. The new methodology is mainly based on a hierarchical and systematic application of alignment-based sequence comparisons to identify and measure the polymorphism of TR. While previous catalogs focus on the class of STR of small total size, we remove any size restrictions, aiming at the more general class of PTR, and we also target fuzzy TR by using specific detection tools. Similarly to other previous catalogs of human polymorphic loci, we focus our catalog toward applications in the discovery of disease-associated loci. Validation by cross-referencing with existing catalogs on common clinically-relevant loci shows good concordance. Overall, this proposed census of human PTR in genic regions is a shared resource (web accessible), complementary to existing catalogs, facilitating future genome-wide studies involving PTR.

The association of insertions/deletions (INDELs) and variable number tandem repeats (VNTRs) with obesity and its related traits and complications

BACKGROUND

Despite the fact that insertions/deletions (INDELs) are the second most common type of genetic variations and variable number tandem repeats (VNTRs) represent a large portion of the human genome, they have received far less attention than single nucleotide polymorphisms (SNPs) and larger forms of structural variation like copy number variations (CNVs), especially in genome-wide association studies (GWAS) of complex diseases like polygenic obesity. This is exemplified by the vast amount of review papers on the role of SNPs and CNVs in obesity, its related traits (like anthropometric measurements, biochemical variables, and eating behavior), and its related complications (like hypertension, hypertriglyceridemia, hypercholesterolemia, and insulin resistance-collectively known as metabolic syndrome). Hence, this paper reviews the types of INDELs and VNTRs that have been studied for association with obesity and its related traits and complications. These INDELs and VNTRs could be found in the obesity loci or genes from the earliest GWAS and candidate gene association studies, like FTO, genes in the leptin-proopiomelanocortin pathway, and UCP2/3. Given the important role of the brain serotonergic and dopaminergic reward system in obesity susceptibility, the association of INDELs and VNTRs in these neurotransmitters' metabolism and transport genes with obesity is also reviewed. Next, the role of INS VNTR in obesity and its related traits is questionable, since recent large-scale studies failed to replicate the earlier positive associations. As obesity results in chronic low-grade inflammation of the adipose tissue, the proinflammatory cytokine gene IL1RA and anti-inflammatory cytokine gene IL4 have VNTRs that are implicated in obesity. A systemic proinflammatory state in combination with activation of the renin-angiotensin system and decreased nitric oxide bioavailability as found in obesity leads to endothelial dysfunction. This explains why VNTR and INDEL in eNOS and ACE, respectively, could be predisposing factors of obesity. Finally, two novel genes, DOCK5 and PER3, which are involved in the regulation of the Akt/MAPK pathway and circadian rhythm, respectively, have VNTRs and INDEL that might be associated with obesity.

SHORT CONCLUSION

In conclusion, INDELs and VNTRs could have important functional consequences in the pathophysiology of obesity, and research on them should be continued to facilitate obesity prediction, prevention, and treatment.

Dynamic interaction between fetal adversity and a genetic score reflecting dopamine function on developmental outcomes at 36 months

Background

Fetal adversity, evidenced by poor fetal growth for instance, is associated with increased risk for several diseases later in life. Classical cut-offs to characterize small (SGA) and large for gestational age (LGA) newborns are used to define long term vulnerability. We aimed at exploring the possible dynamism of different birth weight cut-offs in defining vulnerability in developmental outcomes (through the Bayley Scales of Infant and Toddler Development), using the example of a gene vs. fetal adversity interaction considering gene choices based on functional relevance to the studied outcome.

Methods

36-month-old children from an established prospective birth cohort (Maternal Adversity, Vulnerability, and Neurodevelopment) were classified according to birth weight ratio (BWR) (SGA ≤0.85, LGA >1.15, exploring a wide range of other cut-offs) and genotyped for polymorphisms associated with dopamine signaling (TaqIA-A1 allele, DRD2-141C Ins/Ins, DRD4 7-repeat, DAT1-10- repeat, Met/Met-COMT), composing a score based on the described function, in which hypofunctional variants received lower scores.

Results

There were 251 children (123 girls and 128 boys). Using the classic cut-offs (0.85 and 1.15), there were no statistically significant interactions between the neonatal groups and the dopamine genetic score. However, when changing the cut-offs, it is possible to see ranges of BWR that could be associated with vulnerability to poorer development according to the variation in the dopamine function.

Conclusion

The classic birth weight cut-offs to define SGA and LGA newborns should be seen with caution, as depending on the outcome in question, the protocols for long-term follow up could be either too inclusive—therefore most costly, or unable to screen true vulnerabilities—and therefore ineffective to establish early interventions and primary prevention.

Genetic sequence-based prediction of long-range chromatin interactions suggests a potential role of short tandem repeat sequences in genome organization

Background

Knowing the three-dimensional (3D) structure of the chromatin is important for obtaining a complete picture of the regulatory landscape. Changes in the 3D structure have been implicated in diseases. While there exist approaches that attempt to predict the long-range chromatin interactions, they focus only on interactions between specific genomic regions — the promoters and enhancers, neglecting other possibilities, for instance, the so-called structural interactions involving intervening chromatin.

Results

We present a method that can be trained on 5C data using the genetic sequence of the candidate loci to predict potential genome-wide interaction partners of a particular locus of interest. We have built locus-specific support vector machine (SVM)-based predictors using the oligomer distance histograms (ODH) representation. The method shows good performance with a mean test AUC (area under the receiver operating characteristic (ROC) curve) of 0.7 or higher for various regions across cell lines GM12878, K562 and HeLa-S3. In cases where any locus did not have sufficient candidate interaction partners for model training, we employed multitask learning to share knowledge between models of different loci. In this scenario, across the three cell lines, the method attained an average performance increase of 0.09 in the AUC. Performance evaluation of the models trained on 5C data regarding prediction on an independent high-resolution Hi-C dataset (which is a rather hard problem) shows 0.56 AUC, on average. Additionally, we have developed new, intuitive visualization methods that enable interpretation of sequence signals that contributed towards prediction of locus-specific interaction partners. The analysis of these sequence signals suggests a potential general role of short tandem repeat sequences in genome organization.

Conclusions

We demonstrated how our approach can 1) provide insights into sequence features of locus-specific interaction partners, and 2) also identify their cell-line specificity. That our models deem short tandem repeat sequences as discriminative for prediction of potential interaction partners, suggests that they could play a larger role in genome organization. Thus, our approach can (a) be beneficial to broadly understand, at the sequence-level, chromatin interactions and higher-order structures like (meta-) topologically associating domains (TADs); (b) study regions omitted from existing prediction approaches using various information sources (e.g., epigenetic information); and (c) improve methods that predict the 3D structure of the chromatin.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1624-x) contains supplementary material, which is available to authorized users.

New PAH gene promoter KLF1 and 3'-region C/EBPalpha motifs influence transcription in vitro

Phenylketonuria (PKU) is a metabolic disease caused by mutations in the phenylalanine hydroxylase (PAH) gene. Although the PAH genotype remains the main determinant of PKU phenotype severity, genotype-phenotype inconsistencies have been reported. In this study, we focused on unanalysed sequences in non-coding PAH gene regions to assess their possible influence on the PKU phenotype. We transiently transfected HepG2 cells with various chloramphenicol acetyl transferase (CAT) reporter constructs which included PAH gene non-coding regions. Selected non-coding regions were indicated by in silico prediction to contain transcription factor binding sites. Furthermore, electrophoretic mobility shift assay (EMSA) and supershift assays were performed to identify which transcriptional factors were engaged in the interaction. We found novel KLF1 motif in the PAH promoter, which decreases CAT activity by 50 % in comparison to basal transcription in vitro. The cytosine at the c.-170 promoter position creates an additional binding site for the protein complex involving KLF1 transcription factor. Moreover, we assessed for the first time the role of a multivariant variable number tandem repeat (VNTR) region located in the 3'-region of the PAH gene. We found that the VNTR3, VNTR7 and VNTR8 constructs had approximately 60 % of CAT activity. The regulation is mediated by the C/EBPalpha transcription factor, present in protein complex binding to VNTR3. Our study highlighted two novel promoter KLF1 and 3'-region C/EBPalpha motifs in the PAH gene which decrease transcription in vitro and, thus, could be considered as PAH expression modifiers. New transcription motifs in non-coding regions will contribute to better understanding of the PKU phenotype complexity and may become important for the optimisation of PKU treatment.

Monogenic Diabetes: What It Teaches Us on the Common Forms of Type 1 and Type 2 Diabetes

To date, more than 30 genes have been linked to monogenic diabetes. Candidate gene and genome-wide association studies have identified > 50 susceptibility loci for common type 1 diabetes (T1D) and approximately 100 susceptibility loci for type 2 diabetes (T2D). About 1-5% of all cases of diabetes result from single-gene mutations and are called monogenic diabetes. Here, we review the pathophysiological basis of the role of monogenic diabetes genes that have also been found to be associated with common T1D and/or T2D. Variants of approximately one-third of monogenic diabetes genes are associated with T2D, but not T1D. Two of the T2D-associated monogenic diabetes genes-potassium inward-rectifying channel, subfamily J, member 11 (KCNJ11), which controls glucose-stimulated insulin secretion in the β-cell; and peroxisome proliferator-activated receptor γ (PPARG), which impacts multiple tissue targets in relation to inflammation and insulin sensitivity-have been developed as major antidiabetic drug targets. Another monogenic diabetes gene, the preproinsulin gene (INS), is unique in that INS mutations can cause hyperinsulinemia, hyperproinsulinemia, neonatal diabetes mellitus, one type of maturity-onset diabetes of the young (MODY10), and autoantibody-negative T1D. Dominant heterozygous INS mutations are the second most common cause of permanent neonatal diabetes. Moreover, INS gene variants are strongly associated with common T1D (type 1a), but inconsistently with T2D. Variants of the monogenic diabetes gene Gli-similar 3 (GLIS3) are associated with both T1D and T2D. GLIS3 is a key transcription factor in insulin production and β-cell differentiation during embryonic development, which perturbation forms the basis of monogenic diabetes as well as its association with T1D. GLIS3 is also required for compensatory β-cell proliferation in adults; impairment of this function predisposes to T2D. Thus, monogenic forms of diabetes are invaluable "human models" that have contributed to our understanding of the pathophysiological basis of common T1D and T2D.

A novel polymorphic repeat in the upstream regulatory region of the estrogen-induced gene EIG121 is not associated with the risk of developing breast or endometrial cancer

BACKGROUND

The estrogen-induced gene 121 (EIG121) has been associated with breast and endometrial cancers, but its mechanism of action remains unknown. In a genome-wide search for tandem repeats, we found that EIG121 contains a short tandem repeat (STR) in its upstream regulatory region which has the potential to alter gene expression. The presence of this STR has not previously been analysed in relation to breast or endometrial cancer risk.

RESULTS

In this study, the lengths of this STR were determined by PCR, fragment analysis and sequencing using DNA from 223 breast cancer patients, 204 endometrial cancer patients and 220 healthy controls to determine if they were associated with the risk of developing breast or endometrial cancer. We found this repeat to be highly variable with the number of copies of the AG motif ranging from 27 to 72 and having a bimodal distribution. No statistically significant association was identified between the length of this STR and the risk of developing breast or endometrial cancer or age at diagnosis.

CONCLUSIONS

The STR in the upstream regulatory region of EIG121 is highly polymorphic, but is not associated with the risk of developing breast or endometrial cancer in the cohorts analysed here. While this polymorphic STR in the regulatory region of EIG121 appears to have no impact on the risk of developing breast or endometrial cancer, its association with disease recurrence or overall survival remains to be determined.

Differences of Variable Number Tandem Repeats in XRCC5 Promoter Are Associated with Increased or Decreased Risk of Breast Cancer in BRCA Gene Mutation Carriers

Ku80 is a subunit of the Ku heterodimer that binds to DNA double-strand break ends as part of the non-homologous end joining (NHEJ) pathway. Ku80 is also involved in homologous recombination (HR) via its interaction with BRCA1. Ku80 is encoded by the XRCC5 gene that contains a variable number tandem repeat (VNTR) insertion in its promoter region. Different VNTR genotypes can alter XRCC5 expression and affect Ku80 production, thereby affecting NHEJ and HR pathways. VNTR polymorphism is associated with multiple types of sporadic cancer. In this study, we investigated its potential association with familial breast cancer at the germline level. Using PCR, PAGE, Sanger sequencing, and statistical analyses, we compared VNTR genotypes in the XRCC5 promoter between healthy individuals and three types of familial breast cancer cases: mutated BRCA1 (BRCA1⁺), mutated BRCA2 (BRCA2⁺), and wild-type BRCA1/BRCA2 (BRCAx). We observed significant differences of VNTR genotypes between control and BRCA1⁺ group (P < 0.0001) and BRCA2⁺ group (P = 0.0042) but not BRCAx group (P = 0.2185), and the differences were significant between control and cancer-affected BRCA1⁺ cases (P < 0.0001) and BRCA2⁺ cases (P = 0.0092) but not cancer-affected BRCAx cases (P = 0.4251). Further analysis indicated that 2R/2R (OR = 1.94, 95%CI = 1.26–2.95, P = 0.0096) and 2R/1R (OR = 1.58, 95%CI = 1.11–2.26, P = 0.0388) were associated with increased risk but 1R/1R (OR = 0.55, 95%CI = 0.35–0.84, P = 0.0196) and 1R/0R (OR = 0, 95%CI = 0–0.29, P = 0.0012) were associated with decreased risk in cancer-affected BRCA1⁺ group; 2R/1R (OR = 1.94, 95%CI = 1.14–3.32, P = 0.0242) was associated with increased risk in cancer-affected BRCA2⁺ group. No correlation was observed for the altered risk between cancer-affected or -unaffected carriers and between different age of cancer diagnosis in cancer-affected carriers. The frequently observed VNTR association with in BRCA1⁺ and BRCA2⁺ breast cancer group indicates that VNTR polymorphism in the XRCC5 promoter is associated with altered risk of breast cancer in BRCA1⁺ and BRCA2⁺ carriers.

Complete sequencing and characterization of equine aggrecan

OBJECTIVES

To fully sequence and characterize equine aggrecan and confirm conservation of major aggrecanase, calpain and matrix metalloproteinase (MMP) cleavage sites.

METHODS

Reverse transcription-polymerase chain reaction and rapid amplification of cDNA ends were used to generate clones that encompassed the complete equine aggrecan sequence. Clones were sequenced and compared with the equine genome database to determine intron-exon boundaries.

RESULTS

The aggrecan gene spans over 61 kb on chromosome 1 and is encoded by 17 exons. Two major variants of aggrecan were cloned; one containing 8187 bp (2728 amino acids) and a second sequence of 8061 nucleotides (2686 amino acids). The variation was due to a CS1 domain polymorphism. Both sequences are substantially larger than predicted by the genomic database; 11 CS1 repeat elements are absent in the database sequence. The equine amino acid sequence was compared with human, bovine and murine sequences. Globular domains 1, 2 and 3 are highly conserved (overall identity over 80%). Equine CS1 is considerably larger than in other species and, therefore, is the least conserved domain (an overall amino acid identity of 22%). Previously defined aggrecanase, calpain and MMP cleavage sites were identified. Western blotting of chondrocyte culture samples showed complex post-secretion processing.

CLINICAL SIGNIFICANCE

The complete equine aggrecan sequence will support more in-depth research on aggrecan processing and degradation in equine articular cartilage and other musculoskeletal tissues.

Diversification of 13S globulins, allergenic seed storage proteins, of common buckwheat

The α polypeptide of the 13S globulin subunit of common buckwheat is the counterpart of the major allergenic β polypeptide. Trypsin digestibility varies between variants of the α polypeptide with and without a tandem repeat insert. To evaluate the intra-species diversity of 13S globulin, the comprehensive screening of a genomic DNA library was performed, resulting in the isolation of 14 and 3 genes for Met-poor and Met-rich subunits, respectively. Although most tandem repeat units were 45 bp in length, the two-repeat gene Glb2B and all one-repeat genes contained an additional 3 bp. Secondary structure predictions and polyacrylamide gel electrophoresis demonstrated that the sense strand of Glb2B-CCG, the additional 3 bp-deletion clone of Glb2B, formed a more rigid secondary structure than that of the wild-type. Thus, the large intra-species variation of 13S globulin revealed in this study and its diversification might be attributable to the unique nature of the tandem repeat sequences.