Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation.

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges.

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

Digital genotyping of macrosatellites and multicopy genes reveals novel biological functions associated with copy number variation of large tandem repeats.

Tandem repeats are common in eukaryotic genomes, but due to difficulties in assaying them remain poorly studied. Here, we demonstrate the utility of Nanostring technology as a targeted approach to perform accurate measurement of tandem repeats even at extremely high copy number, and apply this technology to genotype 165 HapMap samples from three different populations and five species of non-human primates. We observed extreme variability in copy number of tandemly repeated genes, with many loci showing 5-10 fold variation in copy number among humans. Many of these loci show hallmarks of genome assembly errors, and the true copy number of many large tandem repeats is significantly under-represented even in the high quality 'finished' human reference assembly. Importantly, we demonstrate that most large tandem repeat variations are not tagged by nearby SNPs, and are therefore essentially invisible to SNP-based GWAS approaches. Using association analysis we identify many cis correlations of large tandem repeat variants with nearby gene expression and DNA methylation levels, indicating that variations of tandem repeat length are associated with functional effects on the local genomic environment. This includes an example where expansion of a macrosatellite repeat is associated with increased DNA methylation and suppression of nearby gene expression, suggesting a mechanism termed "repeat induced gene silencing", which has previously been observed only in transgenic organisms. We also observed multiple signatures consistent with altered selective pressures at tandemly repeated loci, suggesting important biological functions. Our studies show that tandemly repeated loci represent a highly variable fraction of the genome that have been systematically ignored by most previous studies, copy number variation of which can exert functionally significant effects. We suggest that future studies of tandem repeat loci will lead to many novel insights into their role in modulating both genomic and phenotypic diversity.

Electron density shape analysis of a family of through-space and through-bond interactions.

A family of styrene derivatives has been used to study the effects of through-space and through-bond interactions on the local and global shapes of electron densities of complete molecules and a set of substituents on their central rings. Shape analysis methods which have been used extensively in the past for the study of molecular property-molecular shape correlations have shown that in these molecules a complementary role is played by the through-space and through-bond interactions. For each specific example, the dominance of either one of the two interactions can be identified and interpreted in terms of local shapes and the typical reactivities of the various substituents. Three levels of quantum chemical computational methods have been applied for these structures, including the B3LYP/cc-pVTZ level of density functional methodology, and the essential conclusions are the same for all three levels. The general approach is suggested as a tool for the identification of specific interaction types which are able to modify molecular electron densities. By separately influencing the through-space and through-bond components using polar groups and groups capable of conjugation, some fine-tuning of the overall effects becomes possible. The method described may contribute to an improved understanding and control of molecular properties involving complex interactions with a possible role in the emerging field of molecular design.

Noncoding variants alter GATA2 expression in rhombomere 4 motor neurons and cause dominant hereditary congenital facial paresis.

Hereditary congenital facial paresis type 1 (HCFP1) is an autosomal dominant disorder of absent or limited facial movement that maps to chromosome 3q21-q22 and is hypothesized to result from facial branchial motor neuron (FBMN) maldevelopment. In the present study, we report that HCFP1 results from heterozygous duplications within a neuron-specific GATA2 regulatory region that includes two enhancers and one silencer, and from noncoding single-nucleotide variants (SNVs) within the silencer. Some SNVs impair binding of NR2F1 to the silencer in vitro and in vivo and attenuate in vivo enhancer reporter expression in FBMNs. Gata2 and its effector Gata3 are essential for inner-ear efferent neuron (IEE) but not FBMN development. A humanized HCFP1 mouse model extends Gata2 expression, favors the formation of IEEs over FBMNs and is rescued by conditional loss of Gata3. These findings highlight the importance of temporal gene regulation in development and of noncoding variation in rare mendelian disease.

Formation of novel CENP-A domains on tandem repetitive DNA and across chromosome breakpoints on human chromosome 8q21 neocentromeres.

Endogenous human centromeres form on megabase-sized arrays of tandemly repeated alpha satellite DNA. Human neocentromeres form epigenetically at ectopic sites devoid of alpha satellite DNA and permit analysis of centromeric DNA and chromatin organization. In this study, we present molecular cytogenetic and CENP-A chromatin immunoprecipitation (ChIP) on CHIP analyses of two neocentromeres that have formed in chromosome band 8q21 each with a unique DNA and CENP-A chromatin configuration. The first neocentromere was found on a neodicentric chromosome 8 with an inactivated endogenous centromere, where the centromeric activity and CENP-A domain were repositioned to band 8q21 on a large tandemly repeated DNA. This is the first example of a neocentromere forming on repetitive DNA, as all other mapped neocentromeres have formed on single copy DNA. Quantitative fluorescent in situ hybridization (FISH) analysis showed a 60% reduction in the alpha satellite array size at the inactive centromere compared to the active centromere on the normal chromosome 8. This neodicentric chromosome may provide insight into centromere inactivation and the role of tandem DNA in centromere structure. The second neocentromere was found on a neocentric ring chromosome that contained the 8q21 tandemly repeated DNA, although the neocentromere was localized to a different genomic region. Interestingly, this neocentromere is composed of two distinct CENP-A domains in bands 8q21 and 8q24, which are brought into closer proximity on the ring chromosome. This neocentromere suggests that chromosomal rearrangement and DNA breakage may be involved in neocentromere formation. These novel examples provide insight into the formation and structure of human neocentromeres.

Large inverted repeats within Xp11.2 are present at the breakpoints of isodicentric X chromosomes in Turner syndrome.

Turner syndrome (TS) results from whole or partial monosomy X and is mediated by haploinsufficiency of genes that normally escape X-inactivation. Although a 45,X karyotype is observed in half of all TS cases, the most frequent variant TS karyotype includes the isodicentric X chromosome alone [46,X,idic(X)(p11)] or as a mosaic [46,X,idic(X)(p11)/45,X]. Given the mechanism of idic(X)(p11) rearrangement is poorly understood and breakpoint sequence information is unknown, this study sought to investigate the molecular mechanism of idic(X)(p11) formation by determining their precise breakpoint intervals. Karyotype analysis and fluorescence in situ hybridization mapping of eight idic(X)(p11) cell lines and three unbalanced Xp11.2 translocation lines identified the majority of breakpoints within a 5 Mb region, from approximately 53 to 58 Mb, in Xp11.1-p11.22, clustering into four regions. To further refine the breakpoints, a high-resolution oligonucleotide microarray (average of approximately 350 bp) was designed and array-based comparative genomic hybridization (aCGH) was performed on all 11 idic(X)(p11) and Xp11.2 translocation lines. aCGH analyses identified all breakpoint regions, including an idic(X)(p11) line with two potential breakpoints, one breakpoint shared between two idic(X)(p11) lines and two Xp translocations that shared breakpoints with idic(X)(p11) lines. Four of the breakpoint regions included large inverted repeats composed of repetitive gene clusters and segmental duplications, which corresponded to regions of copy-number variation. These data indicate that the rearrangement sites on Xp11.2 that lead to isodicentric chromosome formation and translocations are probably not random and suggest that the complex repetitive architecture of this region predisposes it to rearrangements, some of which are recurrent.

Tetrasomy 15q26: a distinct syndrome or Shprintzen-Goldberg syndrome phenocopy?

PURPOSE: The aim of this study was to characterize the clinical phenotype of patients with tetrasomy of the distal 15q chromosome in the form of a neocentric marker chromosome and to evaluate whether the phenotype represents a new clinical syndrome or is a phenocopy of Shprintzen-Goldberg syndrome. METHODS: We carried out comprehensive clinical evaluation of four patients who were identified with a supernumerary marker chromosome. The marker chromosome was characterized by G-banding, fluorescence in situ hybridization, single nucleotide polymorphism oligonucleotide microarray analysis, and immunofluorescence with antibodies to centromere protein C. RESULTS: The marker chromosomes were categorized as being neocentric with all showing tetrasomy for regions distal to 15q25 and the common region of overlap being 15q26→qter. CONCLUSION: Tetrasomy of 15q26 likely results in a distinct syndrome as the patients with tetrasomy 15q26 share a strikingly more consistent phenotype than do the patients with Shprintzen-Goldberg syndrome, who show remarkable clinical variation.

Duplication of the ZIC2 gene is not associated with holoprosencephaly.

Cytogenetic testing using genomic microarrays presents a clinical challenge when data regarding the phenotypic consequences of the genomic alteration are not available. We describe a chromosome 13q32.3 duplication discovered by microarray testing in a fetus with a prenatally detected apparently balanced de novo translocation 46,XY,t(2;13)(q37;q32). Microarray analysis on the fetal DNA showed duplications of 384 and 564 kb at the breakpoint regions on chromosomes 2q37.3 and 13q32.3, respectively. There were no disease-associated genes in the duplicated region on chromosome 2q37. The duplicated region on chromosome 13q contains the ZIC2 gene. Haploinsufficiency of ZIC2 is known to cause holoprosencephaly and other brain malformations. Studies in the mouse models have suggested that over expression of ZIC2 may also lead to brain malformations. Fetal MRI of the brain was normal and the family elected to continue the pregnancy. An apparently normal baby was born at term. At 3 months of age a physical exam showed no abnormalities and no developmental delay. This report shows that duplication of ZIC2 is not necessarily associated with brain malformations. We also describe the phenotype from four additional patients with duplications of the region of chromosome 13 containing ZIC2 and three previously described patients with supernumerary marker chromosomes derived from distal chromosome 13. None of the eight patients had holoprosencephaly or brain malformations, indicating that duplication of ZIC2 is not associated with brain anomalies. This information will be useful for counseling in other occurrences of this duplication identified by microarray.

Atoms and bonds in molecules from radial densities.

Radial densities are explored as an alternative method for partitioning the molecular density into atomic regions and bonding regions. The radial densities for atoms in molecules are similar to those of an isolated atom. The method may also provide an alternative to Bragg-Slater radii.

Analysis of transposon interruptions suggests selection for L1 elements on the X chromosome.

It has been hypothesised that the massive accumulation of L1 transposable elements on the X chromosome is due to their function in X inactivation, and that the accumulation of Alu elements near genes is adaptive. We tested the possible selective advantage of these two transposable element (TE) families with a novel method, interruption analysis. In mammalian genomes, a large number of TEs interrupt other TEs due to the high overall abundance and age of repeats, and these interruptions can be used to test whether TEs are selectively neutral. Interruptions of TEs, which are beneficial for the host, are expected to be deleterious and underrepresented compared with neutral ones. We found that L1 elements in the regions of the X chromosome that contain the majority of the inactivated genes are significantly less frequently interrupted than on the autosomes, while L1s near genes that escape inactivation are interrupted with higher frequency, supporting the hypothesis that L1s on the X chromosome play a role in its inactivation. In addition, we show that TEs are less frequently interrupted in introns than in intergenic regions, probably due to selection against the expansion of introns, but the insertion pattern of Alus is comparable to other repeats.

Biased distributions and decay of long interspersed nuclear elements in the chicken genome.

The genomes of birds are much smaller than mammalian genomes, and transposable elements (TEs) make up only 10% of the chicken genome, compared with the 45% of the human genome. To study the mechanisms that constrain the copy numbers of TEs, and as a consequence the genome size of birds, we analyzed the distributions of LINEs (CR1's) and SINEs (MIRs) on the chicken autosomes and Z chromosome. We show that (1) CR1 repeats are longest on the Z chromosome and their length is negatively correlated with the local GC content; (2) the decay of CR1 elements is highly biased, and the 5'-ends of the insertions are lost much faster than their 3'-ends; (3) the GC distribution of CR1 repeats shows a bimodal pattern with repeats enriched in both AT-rich and GC-rich regions of the genome, but the CR1 families show large differences in their GC distribution; and (4) the few MIRs in the chicken are most abundant in regions with intermediate GC content. Our results indicate that the primary mechanism that removes repeats from the chicken genome is ectopic exchange and that the low abundance of repeats in avian genomes is likely to be the consequence of their high recombination rates.

Analysis of the largest tandemly repeated DNA families in the human genome.

BACKGROUND: Tandemly Repeated DNA represents a large portion of the human genome, and accounts for a significant amount of copy number variation. Here we present a genome wide analysis of the largest tandem repeats found in the human genome sequence. RESULTS: Using Tandem Repeats Finder (TRF), tandem repeat arrays greater than 10 kb in total size were identified, and classified into simple sequence e.g. GAATG, classical satellites e.g. alpha satellite DNA, and locus specific VNTR arrays. Analysis of these large sequenced regions revealed that several "simple sequence" arrays actually showed complex domain and/or higher order repeat organization. Using additional methods, we further identified a total of 96 additional arrays with tandem repeat units greater than 2 kb (the detection limit of TRF), 53 of which contained genes or repeated exons. The overall size of an array of tandem 12 kb repeats which spanned a gap on chromosome 8 was found to be 600 kb to 1.7 Mbp in size, representing one of the largest non-centromeric arrays characterized. Several novel megasatellite tandem DNA families were observed that are characterized by repeating patterns of interspersed transposable elements that have expanded presumably by unequal crossing over. One of these families is found on 11 different chromosomes in >25 arrays, and represents one of the largest most widespread megasatellite DNA families. CONCLUSION: This study represents the most comprehensive genome wide analysis of large tandem repeats in the human genome, and will serve as an important resource towards understanding the organization and copy number variation of these complex DNA families.

Evolutionary history of mammalian transposons determined by genome-wide defragmentation.

The constant bombardment of mammalian genomes by transposable elements (TEs) has resulted in TEs comprising at least 45% of the human genome. Because of their great age and abundance, TEs are important in comparative phylogenomics. However, estimates of TE age were previously based on divergence from derived consensus sequences or phylogenetic analysis, which can be unreliable, especially for older more diverged elements. Therefore, a novel genome-wide analysis of TE organization and fragmentation was performed to estimate TE age independently of sequence composition and divergence or the assumption of a constant molecular clock. Analysis of TEs in the human genome revealed approximately 600,000 examples where TEs have transposed into and fragmented other TEs, covering >40% of all TEs or approximately 542 Mbp of genomic sequence. The relative age of these TEs over evolutionary time is implicit in their organization, because newer TEs have necessarily transposed into older TEs that were already present. A matrix of the number of times that each TE has transposed into every other TE was constructed, and a novel objective function was developed that derived the chronological order and relative ages of human TEs spanning >100 million years. This method has been used to infer the relative ages across all four major TE classes, including the oldest, most diverged elements. Analysis of DNA transposons over the history of the human genome has revealed the early activity of some MER2 transposons, and the relatively recent activity of MER1 transposons during primate lineages. The TEs from six additional mammalian genomes were defragmented and analyzed. Pairwise comparison of the independent chronological orders of TEs in these mammalian genomes revealed species phylogeny, the fact that transposons shared between genomes are older than species-specific transposons, and a subset of TEs that were potentially active during periods of speciation.

DNA modification and functional delivery into human cells using Escherichia coli DH10B.

The availability of almost the complete human genome as cloned BAC libraries represents a valuable resource for functional genomic analysis, which, however, has been somewhat limited by the ability to modify and transfer this DNA into mammalian cells intact. Here we report a novel comprehensive Escherichia coli-based vector system for the modification, propagation and delivery of large human genomic BAC clones into mammalian cells. The GET recombination inducible homologous recombination system was used in the BAC host strain E.coli DH10B to precisely insert an EGFPneo cassette into the vector portion of a approximately 200 kb human BAC clone, providing a relatively simple method to directly convert available BAC clones into suitable vectors for mammalian cells. GET recombination was also used for the targeted deletion of the asd gene from the E.coli chromosome, resulting in defective cell wall synthesis and diaminopimelic acid auxotrophy. Transfer of the Yersinia pseudotuberculosis invasin gene into E.coli DH10B asd(-) rendered it competent to invade HeLa cells and deliver DNA, as judged by transient expression of green fluorescent protein and stable neomycin-resistant colonies. The efficiency of DNA transfer and survival of HeLa cells has been optimized for incubation time and multiplicity of infection of invasive E.coli with HeLa cells. This combination of E.coli-based homologous recombination and invasion technologies using BAC host strain E.coli DH10B will greatly improve the utility of the available BAC libraries from the human and other genomes for gene expression and functional genomic studies.

CP248, a derivative of exisulind, causes growth inhibition, mitotic arrest, and abnormalities in microtubule polymerization in glioma cells.

Exisulind (sulindac sulfone) and two potent derivatives, CP248 and CP461, have been shown previously to cause growth inhibition and apoptosis in several types of human carcinoma cell lines. These and related compounds have not been previously studied with respect to glioma cell lines. In the present study, we found that these three compounds caused marked growth inhibition in four rat glioma and eight human glioma cell lines, with IC50 values of 150, 1, and 0.075 microm, respectively. When studied at these concentrations exisulind and CP461 had no significant effect on the cell cycle profile of glioma cells, but CP248 caused marked arrest in mitosis. Detailed studies of CP248 in the 9L rat gliosarcoma cell line indicated that treatment with 0.075 microM CP248 caused abnormalities in the spindle apparatus and activation of the spindle assembly check point. In interphase glioma cells, CP248 stabilized microtubules (MTs) at low concentrations (0.075 microM) and depolymerized MTs at higher concentrations (0.2-0.4 microM). In NIH 3T3 fibroblasts, 0.1 microM CP248 caused extensive MT depolymerization. CP248 also caused MT depolymerization when added to assembled MTs in vitro, which indicated that it can directly affect MTs, perhaps because it shares certain structural similarities with Colcemid. In glioma cells, the effects of CP248 on MTs were independent of the previously reported effects of this compound on activation of protein kinase G. Therefore, CP248 is a novel MT-active agent that may be useful in the treatment of glioblastoma, and possibly other types of cancer, because of its dual effects on protein kinase G and MTs.

Identification of a neocentromere in a rearranged Y chromosome with no detectable DYZ3 centromeric sequence.

An 18-year-old woman was evaluated because of primary amenorrhea and hypogonadism. Chromosome analysis from peripheral blood lymphocytes revealed a nonmosaic 46,X,+mar constitution. The marker was shown to be a rearranged Y chromosome consisting of an inverted duplication of the long arm: rea(Y)(qter-q11::q11-qter). Deletion mapping analysis with Y-specific STS showed that the marker lacked Yp and Y-centromeric (DYZ3) sequences, but it was positive for Yq sequences tested. Fluorescence in situ hybridization analysis with Y and X chromosome centromeric and pancentromeric probes showed no hybridization signals. The marker chromosome is present in 100% of the cells; therefore, it is mitotically stable despite the absence of DYZ3 centromeric sequence. Hybridization with CENP-A and CENP-C specific antibodies localized a neocentromere close to the breakpoint.

Unique case of mosaicism involving two morphologically similar marker chromosomes of different centric origin in a patient with developmental delay.

A five-year-old Caucasian male presented with developmental delay, minor dysmorphic features, and hyperactivity. Cytogenetic analysis showed the presence of a marker chromosome in the majority of cells analyzed. Fluorescence in situ hybridization (FISH) analyses using several alpha satellite probes, including D13Z1/D21Z1, did not reveal any signal on the marker chromosome. Subsequent multicolor FISH (M-FISH) indicated the marker to be derived from chromosome 13, and FISH with a chromosome 13 paint confirmed this finding. The absence of D13Z1/D21Z1 signal on the marker suggested that it was analphoid in nature. Comparative genomic hybridization (CGH) was utilized to further characterize the region of chromosome 13 from which the marker originated, and unexpectedly revealed a gain of chromosomal material at both the centromeric regions of chromosomes 3 and 13. In view of the CGH results, extensive FISH studies with D3Z1 and D13Z1/D21Z1 were performed and revealed the presence of four cell lines comprising one normal cell line (50.5%), a cell line with a chromosome 3 derived marker (19%), a cell line containing a marker derived from chromosome 13 (20%), and a cell line with both markers (10.5%). As the two markers appeared morphologically similar by GTG banding, all 47,XY metaphases in the initial analysis were thought to comprise only a single marker. This is the first report, to our knowledge, of the presence of a chromosome 3 and a chromosome 13 marker in mosaic condition in a congenital disorder. In light of our experience, we urge caution in interpreting karyotypes with marker chromosomes. Our case illustrates the limitations of fluorescent DNA probes and sampling errors.

Chromosome 13q neocentromeres: molecular cytogenetic characterization of three additional cases and clinical spectrum.

We report three new cases of chromosome 13 derived marker chromosomes, found in unrelated patients with dysmorphisms and/or developmental delay. Molecular cytogenetic analysis was performed using fluorescence in situ hybridization (FISH) with chromosome-specific painting probes, alpha satellite probes, and physically mapped probes from chromosome 13q, as well as comparative genomic hybridization (CGH). This analysis demonstrated that these markers consisted of inversion duplications of distal portions of chromosome 13q that have separated from the endogenous chromosome 13 centromere and contain no detectable alpha satellite DNA. The presence of a functional neocentromere on these marker chromosomes was confirmed by immunofluorescence with antibodies to centromere protein-C (CENP-C). The cytogenetic location of a neocentromere in band 13q32 was confirmed by simultaneous FISH with physically mapped YACs from 13q32 and immunofluorescence with anti-CENP-C. The addition of these three new cases brings the total number of described inv dup 13q neocentic chromosomes to 11, representing 21% (11/52) of the current overall total of 52 described cases of human neocentric chromosomes. This higher than expected frequency suggests that chromosome 13q may have an increased propensity for neocentromere formation. The clinical spectrum of all 11 cases is presented, representing a unique collection of polysomy for different portions of chromosome 13q without aneuploidies for additional chromosomal regions. The complexity and variability of the phenotypes seen in these patients does not support a simple reductionist view of phenotype/genotype correlation with polysomy for certain chromosomal regions.

Expansion of a 12-kb VNTR containing the REXO1L1 gene cluster underlies the microscopically visible euchromatic variant of 8q21.2.

Copy number variants visible with the light microscope have been described as euchromatic variants (EVs) and EVs with extra G-light material at 8q21.2 have been reported only once before. We report four further patients with EVs of 8q21.2 ascertained for clinical (3) or reproductive reasons (1). Enhanced signal strength from two overlapping bacterial artificial chromosomes (BACs) and microarray analysis mapped the EV to a 284-kb interval in the reference genome. This interval consists of a sequence gap flanked by segmental duplications that contain the 12-kb components of one of the largest Variable Number Tandem Repeat arrays in the human genome. Using digital NanoString technology with a custom probe for the RNA exonuclease 1 homologue (S. cerevisiae)-like 1 (REXO1L1) gene within each 12-kb repeat, significantly enhanced diploid copy numbers of 270 and 265 were found in an EV family and a median diploid copy number of 166 copies in 216 controls. These 8q21.2 EVs are not thought to have clinical consequences as the phenotypes of the probands were inconsistent, those referred for reproductive reasons were otherwise phenotypically normal and the REXO1L1 gene has no known disease association. This EV was found in 4/3078 (1 in 770) consecutive referrals for chromosome analysis and needs to be distinguished from pathogenic imbalances of medial 8q. The REXO1L1 gene product is a marker of hepatitis C virus (HCV) infection and a possible association between REXO1L1 copy number and susceptibility to HCV infection, progression or response to treatment has not yet been excluded.