Variable number tandem repeats - Their emerging role in sickness and health

A comprehensive tandem repeat catalog of the human genome

With the increasing availability of long-read sequencing data, high-quality human genome assemblies, and software for fully characterizing tandem repeats, genome-wide genotyping of tandem repeat loci on a population scale becomes more feasible. Such efforts not only expand our knowledge of the tandem repeat landscape in the human genome but also enhance our ability to differentiate pathogenic tandem repeat mutations from benign polymorphisms. To this end, we analyzed 272 genomes assembled using datasets from three public initiatives that employed different long-read sequencing technologies. Here, we report a catalog of over 18 million tandem repeat loci, many of which were previously unannotated. Some of these loci are highly polymorphic, and many of them reside within coding sequences.

Association between PDCD6-VNTR polymorphism and urinary cancer susceptibility

BACKGROUND

Programmed cell death 6 (PDCD6) is known to be involved in apoptosis and tumorigenesis. Given the reported association with urinary cancer susceptibility through SNP analysis, we further analyzed the entire genomic structure of PDCD6.

METHODS

Three VNTR regions (MS1-MS3) were identified through the analysis of the genomic structure of PDCD6. To investigate the association between these VNTR regions and urinary cancer susceptibility, genomic DNA was extracted from 413 cancer-free male controls, 267 bladder cancer patients, and 331 prostate cancer patients. Polymerase chain reaction (PCR) was performed to analyze the PDCD6-MS regions. Statistical analysis was performed to determine the association between specific genotypes and cancer risk. In addition, the effect of specific VNTRs on PDCD6 expression was also confirmed using a reporter vector.

RESULTS

Among the three VNTR regions, MS1 and MS2 exhibited monomorphism, while the MS3 region represented polymorphism, with its transmission to subsequent generations through meiosis substantiating its utility as a DNA typing marker. In a case-control study, the presence of rare alleles within PDCD6-MS3 exhibited significant associations with both bladder cancer (OR = 2.37, 95% CI: 1.33-4.95, P = 0.019) and prostate cancer (OR = 2.11, 95% CI: 1.03-4.36, P = 0.038). Furthermore, through luciferase assays, we validated the impact of the MS3 region on modulating PDCD6 expression.

CONCLUSIONS

This study suggests that the PDCD6-MS3 region could serve as a prognostic marker for urinary cancers, specifically bladder cancer and prostate cancer. Moreover, the subdued influence exerted by PDCD6-MS3 on the expression of PDCD6 offers another insight concerning the progression of urinary cancer.

Deciphering the role of a SINE-VNTR-Alu retrotransposon polymorphism as a biomarker of Parkinson's disease progression

SINE-VNTR-Alu (SVA) retrotransposons are transposable elements which represent a source of genetic variation. We previously demonstrated that the presence/absence of a human-specific SVA, termed SVA_67, correlated with the progression of Parkinson's disease (PD). In the present study, we demonstrate that SVA_67 acts as expression quantitative trait loci, thereby exhibiting a strong regulatory effect across the genome using whole genome and transcriptomic data from the Parkinson's progression markers initiative cohort. We further show that SVA_67 is polymorphic for its variable number tandem repeat domain which correlates with both regulatory properties in a luciferase reporter gene assay in vitro and differential expression of multiple genes in vivo. Additionally, this variation's utility as a biomarker is reflected in a correlation with a number of PD progression markers. These experiments highlight the plethora of transcriptomic and phenotypic changes associated with SVA_67 polymorphism which should be considered when investigating the missing heritability of neurodegenerative diseases.

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.

Modification of Huntington's disease by short tandem repeats

Expansions of glutamine-coding CAG trinucleotide repeats cause a number of neurodegenerative diseases, including Huntington's disease and several of spinocerebellar ataxias. In general, age-at-onset of the polyglutamine diseases is inversely correlated with the size of the respective inherited expanded CAG repeat. Expanded CAG repeats are also somatically unstable in certain tissues, and age-at-onset of Huntington's disease corrected for individual HTT CAG repeat length (i.e. residual age-at-onset), is modified by repeat instability-related DNA maintenance/repair genes as demonstrated by recent genome-wide association studies. Modification of one polyglutamine disease (e.g. Huntington's disease) by the repeat length of another (e.g. ATXN3, CAG expansions in which cause spinocerebellar ataxia 3) has also been hypothesized. Consequently, we determined whether age-at-onset in Huntington's disease is modified by the CAG repeats of other polyglutamine disease genes. We found that the CAG measured repeat sizes of other polyglutamine disease genes that were polymorphic in Huntington's disease participants but did not influence Huntington's disease age-at-onset. Additional analysis focusing specifically on ATXN3 in a larger sample set (n = 1388) confirmed the lack of association between Huntington's disease residual age-at-onset and ATXN3 CAG repeat length. Additionally, neither our Huntington's disease onset modifier genome-wide association studies single nucleotide polymorphism data nor imputed short tandem repeat data supported the involvement of other polyglutamine disease genes in modifying Huntington's disease. By contrast, our genome-wide association studies based on imputed short tandem repeats revealed significant modification signals for other genomic regions. Together, our short tandem repeat genome-wide association studies show that modification of Huntington's disease is associated with short tandem repeats that do not involve other polyglutamine disease-causing genes, refining the landscape of Huntington's disease modification and highlighting the importance of rigorous data analysis, especially in genetic studies testing candidate modifiers.

Genetic modifiers of repeat expansion disorders

Repeat expansion disorders (REDs) are monogenic diseases caused by a sequence of repetitive DNA expanding above a pathogenic threshold. A common feature of the REDs is a strong genotype-phenotype correlation in which a major determinant of age at onset (AAO) and disease progression is the length of the inherited repeat tract. Over a disease-gene carrier's life, the length of the repeat can expand in somatic cells, through the process of somatic expansion which is hypothesised to drive disease progression. Despite being monogenic, individual REDs are phenotypically variable, and exploring what genetic modifying factors drive this phenotypic variability has illuminated key pathogenic mechanisms that are common to this group of diseases. Disease phenotypes are affected by the cognate gene in which the expansion is found, the location of the repeat sequence in coding or non-coding regions and by the presence of repeat sequence interruptions. Human genetic data, mouse models and in vitro models have implicated the disease-modifying effect of DNA repair pathways via the mechanisms of somatic mutation of the repeat tract. As such, developing an understanding of these pathways in the context of expanded repeats could lead to future disease-modifying therapies for REDs.

Characterisation of retrotransposon insertion polymorphisms in whole genome sequencing data from individuals with amyotrophic lateral sclerosis

The genetics of an individual is a crucial factor in understanding the risk of developing the neurodegenerative disease amyotrophic lateral sclerosis (ALS). There is still a large proportion of the heritability of ALS, particularly in sporadic cases, to be understood. Among others, active transposable elements drive inter-individual variability, and in humans long interspersed element 1 (LINE1, L1), Alu and SINE-VNTR-Alu (SVA) retrotransposons are a source of polymorphic insertions in the population. We undertook a pilot study to characterise the landscape of non-reference retrotransposon insertion polymorphisms (non-ref RIPs) in 15 control and 15 ALS individuals' whole genomes from Project MinE, an international project to identify potential genetic causes of ALS. The combination of two bioinformatics tools (mobile element locator tool (MELT) and TEBreak) identified on average 1250 Alu, 232 L1 and 77 SVA non-ref RIPs per genome across the 30 analysed. Further PCR validation of individual polymorphic retrotransposon insertions showed a similar level of accuracy for MELT and TEBreak. Our preliminary study did not identify a specific RIP or a significant difference in the total number of non-ref RIPs in ALS compared to control genomes. The use of multiple bioinformatic tools improved the accuracy of non-ref RIP detection and our study highlights the potential importance of studying these elements further in ALS.

Mismatch repair is a double-edged sword in the battle against microsatellite instability

Roughly 3% of the human genome consists of microsatellites or tracts of short tandem repeats (STRs). These STRs are often unstable, undergoing high-frequency expansions (increases) or contractions (decreases) in the number of repeat units. Some microsatellite instability (MSI) is seen at multiple STRs within a single cell and is associated with certain types of cancer. A second form of MSI is characterised by expansion of a single gene-specific STR and such expansions are responsible for a group of 40+ human genetic disorders known as the repeat expansion diseases (REDs). While the mismatch repair (MMR) pathway prevents genome-wide MSI, emerging evidence suggests that some MMR factors are directly involved in generating expansions in the REDs. Thus, MMR suppresses some forms of expansion while some MMR factors promote expansion in other contexts. This review will cover what is known about the paradoxical effect of MMR on microsatellite expansion in mammalian cells.

A polymorphic transcriptional regulatory domain in the amyotrophic lateral sclerosis risk gene CFAP410 correlates with differential isoform expression

We describe the characterisation of a variable number tandem repeat (VNTR) domain within intron 1 of the amyotrophic lateral sclerosis (ALS) risk gene CFAP410 (Cilia and flagella associated protein 410) (previously known as C21orf2), providing insight into how this domain could support differential gene expression and thus be a modulator of ALS progression or risk. We demonstrated the VNTR was functional in a reporter gene assay in the HEK293 cell line, exhibiting both the properties of an activator domain and a transcriptional start site, and that the differential expression was directed by distinct repeat number in the VNTR. These properties embedded in the VNTR demonstrated the potential for this VNTR to modulate CFAP410 expression. We extrapolated these findings in silico by utilisation of tagging SNPs for the two most common VNTR alleles to establish a correlation with endogenous gene expression. Consistent with in vitro data, CFAP410 isoform expression was found to be variable in the brain. Furthermore, although the number of matched controls was low, there was evidence for one specific isoform being correlated with lower expression in those with ALS. To address if the genotype of the VNTR was associated with ALS risk, we characterised the variation of the CFAP410 VNTR in ALS cases and matched controls by PCR analysis of the VNTR length, defining eight alleles of the VNTR. No significant difference was observed between cases and controls, we noted, however, the cohort was unlikely to contain sufficient power to enable any firm conclusion to be drawn from this analysis. This data demonstrated that the VNTR domain has the potential to modulate CFAP410 expression as a regulatory element that could play a role in its tissue-specific and stimulus-inducible regulation that could impact the mechanism by which CFAP410 is involved in ALS.

Non-reference genome transposable elements (TEs) have a significant impact on the progression of the Parkinson's disease

The pathophysiology of Parkinson's disease (PD) is a complex process of the interaction between genetic and environmental factors. Studies on the genetic component of PD have predominantly focused on single nucleotide polymorphisms (SNPs) using a cross-sectional case-control design in large genome-wide association studies. This approach while giving insight into a significant portion of the genetics of PD does not fully account for all the genetic components resulting in missing heritability. In this study, we approached this problem by focusing on the non-reference genome transposable elements (TEs) and their impact on the progression of PD using a longitudinal study design within the Parkinson's progression markers initiative (PPMI) cohort. We analyzed 2886 Alu repeats, 360 LINE1 and 128 SINE-VNTR-Alus (SVAs) that were called from the whole-genome sequence data which are not within the reference genome. The presence or absence of these non-reference TE variants is known as a retrotransposon insertion polymorphism, and measuring this polymorphism describes the impact of TEs on the traits. The variations for the presence or absence of the non-reference TE elements were modeled to align with the changes in the 114 outcome measures during the five-year follow-up period of the PPMI cohort. Linear mixed-effects models were used, and many TEs were found to have a highly significant effect on the longitudinal changes in the clinically important PD outcomes such as UPDRS subscale II, UPDRS total scores, and modified Schwab and England ADL scale. In addition, the progression of several imaging and functional measures, including the Caudate/Putamen ratio and levodopa equivalent daily dose (LEDD) were also significantly affected by the TEs. In conclusion, this study identified the overwhelming effect of the non-reference TEs on the progression of PD and is a good example of the impact the variations in the "junk DNA" have on complex diseases.

TRIB3 promoter 33 bp VNTR is associated with the risk of cerebrovascular disease in type 2 diabetic patients

Previous studies have demonstrated that TRIB3 is closely related to insulin resistance, metabolic disorders and vascular diseases. Recently, it was reported that a 33 bp variable number of tandem repeats (VNTR) located in the TRIB3 promoter could considerably alter its transcriptional activity. Nonetheless, whether the shift of TRIB3 transcriptional activity has the effect of inducing diabetic vascular complications is still unclear. Therefore, in our study, we aimed to explore the relationship between the TRIB3 33bp VNTR and diabetic vascular complications. The TRIB3 33bp VNTR polymorphisms were determined by PCR and Sanger sequencing, a total of 798 eligible Chinese patients with type 2 diabetes (T2DM) were included in our study and then evaluated with clinical data. After adjusting for age, gender, BMI, smoking history, drinking history and duration of diabetes, we found that the high number of 33 bp tandem repeats (repeats>8) was significantly associated with an increase in the risk of cerebrovascular diseases compared with the low number of 33 bp tandem repeats (repeats≤6) in patients with T2DM(OR 2.66, 95% CI 1.29–5.47, p = 0.008). The intermediate number of 33bp tandem repeats (6 < repeat≤8) was markedly associated with a decreased risk of diabetic retinopathy compared with the low number of tandem repeats (OR 0.65, 95% CI 0.46–0.91, p = 0.012). Adjusting for gender, age and BMI, there was a significant difference in DBP levels among patients with the number of different 33 bp tandem repeats (Low vs. Intermediate vs. High, 81.6 ± 12.8 vs. 79.8 ± 12.4 vs. 78.7 ± 12.6 mmHg; p = 0.045). Subgroup analysis found that TRIB3 VNTR was significantly correlated with the difference in systolic blood pressure (SBP) in T2DM patients taking ACEI/ARB drugs (Low vs. Intermediate vs. High, 146.27 ± 18.23 vs. 140.01 ± 19.91 vs. 140.77 ± 18.64 mmHg; p = 0.018). Our results indicated that TRIB3 promoter 33bp VNTR is related to vascular diseases in T2DM patients, and may serve as a new biomarker for individualized prevention and therapy of T2DM.

SARS-CoV-2 and human retroelements: a case for molecular mimicry?

Background

The factors driving the late phase of COVID-19 are still poorly understood. However, autoimmunity is an evolving theme in COVID-19’s pathogenesis. Additionally, deregulation of human retroelements (RE) is found in many viral infections, and has also been reported in COVID-19.

Results

Unexpectedly, coronaviruses (CoV) – including SARS-CoV-2 – harbour many RE-identical sequences (up to 35 base pairs), and some of these sequences are part of SARS-CoV-2 epitopes associated to COVID-19 severity. Furthermore, RE are expressed in healthy controls and human cells and become deregulated after SARS-CoV-2 infection, showing mainly changes in long interspersed nuclear element (LINE1) expression, but also in endogenous retroviruses.

Conclusion

CoV and human RE share coding sequences, which are targeted by antibodies in COVID-19 and thus could induce an autoimmune loop by molecular mimicry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01040-2.

Expression Quantitative Trait Loci (eQTLs) Associated with Retrotransposons Demonstrate their Modulatory Effect on the Transcriptome

Transposable elements (TEs) are repetitive elements that belong to a variety of functional classes and have an important role in shaping genome evolution. Around 50% of the human genome contains TEs, and they have been termed the "dark matter" of the genome because relatively little is known about their function. While TEs have been shown to participate in aberrant gene regulation and the pathogenesis of diseases, only a few studies have explored the systemic effect of TEs on gene expression. In the present study, we analysed whole genome sequences and blood whole transcriptome data from 570 individuals within the Parkinson's Progressive Markers Initiative (PPMI) cohort to identify expression quantitative trait loci (eQTL) regulating genome-wide gene expression associated with TEs. We identified 2132 reference TEs that were polymorphic for their presence or absence in our study cohort. The presence or absence of the TE element could change the expression of the gene or gene clusters from zero to tens of thousands of copies of RNA. The main finding is that many TEs possess very strong regulatory effects, and they have the potential to modulate large genetic networks with hundreds of target genes over the genome. We illustrate the plethora of regulatory mechanisms using examples of their action at the HLA gene cluster and data showing different TEs' convergence to modulate WFS1 gene expression. In conclusion, the presence or absence of polymorphisms of TEs has an eminent genome-wide regulatory function with large effect size at the level of the whole transcriptome. The role of TEs in explaining, in part, the missing heritability for complex traits is convincing and should be considered.