A rigorous approach to facilitate and guarantee the correctness of the genetic testing management in human genome information systems

Abstract Background Recent medical and biological technology advances have stimulated the development of new testing systems that have been providing huge, varied amounts of molecular and clinical data. Growing data volumes pose significant challenges for information processing systems in research centers. Additionally, the routines of genomics laboratory are typically characterized by high parallelism in testing and constant procedure changes. Results This paper describes a formal approach to address this challenge through the implementation of a genetic testing management system applied to human genome laboratory. We introduced the Human Genome Research Center Information System (CEGH) in Brazil, a system that is able to support constant changes in human genome testing and can provide patients updated results based on the most recent and validated genetic knowledge. Our approach uses a common repository for process planning to ensure reusability, specification, instantiation, monitoring, and execution of processes, which are defined using a relational database and rigorous control flow specifications based on process algebra (ACP). The main difference between our approach and related works is that we were able to join two important aspects: 1) process scalability achieved through relational database implementation, and 2) correctness of processes using process algebra. Furthermore, the software allows end users to define genetic testing without requiring any knowledge about business process notation or process algebra. Conclusions This paper presents the CEGH information system that is a Laboratory Information Management System (LIMS) based on a formal framework to support genetic testing management for Mendelian disorder studies. We have proved the feasibility and showed usability benefits of a rigorous approach that is able to specify, validate, and perform genetic testing using easy end user interfaces.

Diversity and genetic structure in natural populations of Hancornia speciosa var. speciosa Gomes in northeastern Brazil

Hancornia speciosa Gomes is a fruit tree native from Brazil that belongs to Apocinaceae family, and is popularly known as Mangabeira. Its fruits are widely consumed raw or processed as fruit jam, juices and ice creams, which have made it a target of intense exploitation. The extractive activities and intense human activity on the environment of natural occurrence of H. speciosa has caused genetic erosion in the species and little is known about the ecology or genetic structure of natural populations. The objective of this research was the evaluation of the genetic diversity and genetic structure of H. speciosa var. speciosa. The genetic variability was assessed using 11 allozyme loci with a sample of 164 individuals distributed in six natural populations located in the States of Pernambuco and Alagoas, Northeastern Brazil. The results showed a high level of genetic diversity within the species (e= 0.36) seeing that the most of the genetic variability of H. speciosa var. speciosa is within its natural populations with low difference among populations ( or = 0.081). The inbreeding values within ( = -0.555) and among populations ( =-0.428) were low showing lacking of endogamy and a surplus of heterozygotes. The estimated gene flow ( m ) was high, ranging from 2.20 to 13.18, indicating to be enough to prevent the effects of genetic drift and genetic differentiation among populations. The multivariate analyses indicated that there is a relationship between genetic and geographical distances, which was confirmed by a spatial pattern analysis using Mantel test (r = 0.3598; p = 0.0920) with 1000 random permutations. The high genetic diversity index in these populations indicates potential for in situ genetic conservation.

Study of variability and genetic structure of European populations of Myotis emarginatus and Myotis capaccinii (Chiroptera, Vespertilionidae)

The Geoffroy’s bat Myotis emarginatus is mainly present in southern, south-eastern and central Europe (Červerný, 1999) and is often recorded from northern Spain (Quetglas, 2002; Flaquer et al., 2004). It has demonstrated the species’ preference for forest. Myotis capaccinii, confined to the Mediterranean (Guille´n, 1999), is classified as ‘vulnerable’ on a global scale (Hutson, Mickleburgh & Racey, 2001). In general, the species preferred calm waters bordered by well-developed riparian vegetation and large (> 5 m) inter-bank distances (Biscardi et al. 2007). In this study we present the first results about population genetic structure of these two species of genus Myotis. We used two methods of sampling: invasive and non-invasive techniques. A total of 323 invasive samples and a total of 107 non-invasive samples were collected and analyzed. For Myotis emarginatus we have individuated for the first time a set of 7 microsatellites, which can work on this species, started from a set developed on Myotis myotis (Castella et al. 2000). We developed also a method for analysis of non-invasive samples, that given a good percentage of positive analyzed samples. The results have highlighted for the species Myotis emarginatus the presence on the European territory of two big groups, discovered by using the microsatellites tracers. On this species, 33 haplotypes of Dloop have been identified, some of them are presented only in some colonies. We identified respectively 33 haplotypes of Dloop and 10 of cytB for Myotis emarginatus and 25 of dloop and 15 of cytB for Myotis capaccinii. Myotis emarginatus’ results, both microsatellites and mtDNA, show that there is a strong genetic flow between different colonies across Europe. The results achieved on Myotis capaccinii are very interesting, in this case either for the microsatellites or the mitochondrial DNA sequences, and it has been highlighted a big difference between different colonies.

Genetic patterns in recovered specimens of Caretta caretta from the Adriatic Sea

The Adriatic Sea is considered a feeding and developmental area for Mediterranean loggerhead turtles, but this area is severely threatened by human impacts. In the Adriatic Sea loggerhead turtles are often found stranded or floating, but they are also recovered as by-catch from fishing activities. Nevertheless, information about population structuring and origin of individuals found in the Adriatic Sea are still limited. Cooperation with fishermen and a good network of voluntary collaborators are essential for understanding their distribution, ecology and for developing conservation strategies in the Adriatic Sea. In this study, a comparative analysis of biometric data and DNA sequence polymorphism of the long fragment of the mitochondrial control region was carried out on ninety-three loggerheads recovered from three feeding areas in the Adriatic Sea: North-western, North-eastern and South Adriatic. Differences in turtles body sizes (e.g. Straight Carapace Length) among the three recovery areas and relationship between SCL and the type of recovery were investigated. The origin of turtles from Mediterranean rookeries and the use of the Adriatic feeding habitats by loggerheads in different life-stages were assessed to understand the migratory pathway of the species. The analysis of biometric data revealed a significant difference in turtle sizes between the Southern and the Northern Adriatic. Moreover, size of captured turtles resulted significantly different from the size of stranded and floating individuals. Actually, neritic sub-adults and adults are more affected by incidental captures than juveniles because of their feeding behavior. The Bayesian mixed-stock analysis showed a strong genetic relationship between the Adriatic aggregates and Mediterranean rookeries, while a low pro¬portion of individuals of Atlantic origin were detected in the Adriatic feeding grounds. The presence of migratory pathways towards the Adriatic Sea due to the surface current system was reinforced by the finding of individuals bearing haplotypes endemic to the nesting populations of Libya, Greece and Israel. A relatively high contribution from Turkey and Cyprus to the Northwest and South Adriatic populations was identified when the three sampled areas were analyzed independently. These results have to be taken in account in a conservative perspective, since coastal hazards, affecting the population of turtles feeding in the Adriatic Sea may also affect the nesting populations of the Eastern Mediterranean with a unique genetic pattern.

Excellent agreement between genetic and hydrogen breath tests for lactase deficiency and the role of extended symptom assessment

Clinical manifestations of lactase (LCT) deficiency include intestinal and extra-intestinal symptoms. Lactose hydrogen breath test (H2-BT) is considered the gold standard to evaluate LCT deficiency (LD). Recently, the single-nucleotide polymorphism C/T(-13910) has been associated with LD. The objectives of the present study were to evaluate the agreement between genetic testing of LCT C/T(-13910) and lactose H2-BT, and the diagnostic value of extended symptom assessment. Of the 201 patients included in the study, 194 (139 females; mean age 38, range 17-79 years, and 55 males, mean age 38, range 18-68 years) patients with clinical suspicion of LD underwent a 3-4 h H2-BT and genetic testing for LCT C/T(-13910). Patients rated five intestinal and four extra-intestinal symptoms during the H2-BT and then at home for the following 48 h. Declaring H2-BT as the gold standard, the CC(-13910) genotype had a sensitivity of 97% and a specificity of 95% with a of 0.9 in diagnosing LCT deficiency. Patients with LD had more intense intestinal symptoms 4 h following the lactose challenge included in the H2-BT. We found no difference in the intensity of extra-intestinal symptoms between patients with and without LD. Symptom assessment yielded differences for intestinal symptoms abdominal pain, bloating, borborygmi and diarrhoea between 120 min and 4 h after oral lactose challenge. Extra-intestinal symptoms (dizziness, headache and myalgia) and extension of symptom assessment up to 48 h did not consistently show different results. In conclusion, genetic testing has an excellent agreement with the standard lactose H2-BT, and it may replace breath testing for the diagnosis of LD. Extended symptom scores and assessment of extra-intestinal symptoms have limited diagnostic value in the evaluation of LD.

A genetic validation study reveals a role of vitamin D metabolism in the response to interferon-alfa-based therapy of chronic hepatitis C

Background To perform a comprehensive study on the relationship between vitamin D metabolism and the response to interferon-α-based therapy of chronic hepatitis C. Methodology/Principal Findings Associations between a functionally relevant polymorphism in the gene encoding the vitamin D 1α-hydroxylase (CYP27B1-1260 rs10877012) and the response to treatment with pegylated interferon-α (PEG-IFN-α) and ribavirin were determined in 701 patients with chronic hepatitis C. In addition, associations between serum concentrations of 25-hydroxyvitamin D3 (25[OH]D3) and treatment outcome were analysed. CYP27B1-1260 rs10877012 was found to be an independent predictor of sustained virologic response (SVR) in patients with poor-response IL28B genotypes (15% difference in SVR for rs10877012 genotype AA vs. CC, p = 0.02, OR = 1.52, 95% CI = 1.061–2.188), but not in patients with favourable IL28B genotype. Patients with chronic hepatitis C showed a high prevalence of vitamin D insufficiency (25[OH]D3<20 ng/mL) during all seasons, but 25(OH)D3 serum levels were not associated with treatment outcome. Conclusions/Significance Our study suggests a role of bioactive vitamin D (1,25[OH]2D3, calcitriol) in the response to treatment of chronic hepatitis C. However, serum concentration of the calcitriol precursor 25(OH)D3 is not a suitable predictor of treatment outcome.

Examining The Genetic, Social-Cognitive, And Neural Correlates Of Autism Spectrum Disorder Behaviors In Typically Developing Chi

Autism spectrum disorders (ASD) are pervasive developmental disorders that affect approximately 1 in 50 children (Blumberg et al., 2013). Due to the social nature of the deficits that characterize the disorders, many have classified them as disorders of social cognition, which is the process that individuals use in order to successfully interact with members of their own species (Frith & Frith, 2007). Previous research has typically neglected the spectrum nature of ASD in favor of a more categorical approach of ¿autistic¿ versus ¿non-autistic,¿ but the spectrum requires a more continuous approach. Thus, the present study sought to examine the genetic, social-cognitive, and neural correlates of ASD-like traits as well as the relationship between these dimensions in typically developing children. Parents and children completed several quantitative measures examining several areas of social-cognitive functioning, including theory of mind and social functioning, restricted/repetitive behaviors and interests, and adaptive/maladaptive functioning. Children were also asked to undergo an EEG and both parents and children contributed a saliva sample that was used to sequence four single nucleotide polymorphisms (SNPs) of the OXTR gene, rs1042778, rs53576, rs2254298, and rs237897. We successfully demonstrated a significant relationship between behavioral measures of social-cognition and differences in face perception via the N170. However, the directionality of these relationships varied based on the behavioral measure and particular N170 difference scores. We also found support for the associations between the G_G allelic combination of rs1042778 and the A_A and A_G allelic combinations of rs2254298 and increased ASD-like behavior with decreased social-cognitive functioning. In contrast, our results contradict previous findings with rs237897 and imply that individuals with the A_A and A_G genotypes are less similar to those with ASD and have higher social cognitive functioning than those with the G_G genotype. In conclusion, we have demonstrated the existence of ASD-like traits in typically developing children and have shown a link between behavioral, genetic, and neural correlates of social-cognition. These findings demonstrate the importance of considering autism as a spectrum disorder and provide support for the move to a more continuous approach to neurodevelopmental disorders.

Do common genetic variants in endotoxin signaling pathway contribute to predisposition to alcoholic liver cirrhosis?

BACKGROUND: Tumor necrosis factor-alpha (TNF-alpha) and interleukin-1beta (IL-1beta), produced by endotoxin-activated Kupffer cells, play a key role in the pathogenesis of alcoholic liver cirrhosis (ALC). Alleles TNFA -238A, IL1B -31T and variant IL1RN*2 of repeat polymorphism in the gene encoding the IL-1 receptor antagonist increase production of TNF-alpha and IL-1beta, respectively. Alleles CD14 -159T, TLR4 c.896G and TLR4 c.1196T modify activation of Kupffer cells by endotoxin. We confirmed the published associations between these common variants and genetic predisposition to ALC by means of a large case-control association study conducted on two Central European populations. METHODS: The study population comprised a Czech sample of 198 ALC patients and 370 controls (MONICA project), and a German sample of 173 ALC patients and 331 controls (KORA-Augsburg), and 109 heavy drinkers without liver disease. RESULTS: Single locus analysis revealed no significant difference between patients and controls in all tested loci. Diplotype [IL1RN 2/ 2; IL1B -31T+] was associated with increased risk of ALC in the pilot study, but not in the validation samples. CONCLUSIONS: Although cytokine mediated immune reactions play a role in the pathogenesis of ALC, hereditary susceptibility caused by variants in the corresponding genes is low in Central European populations.

Expanding the spectrum of mutations in GH1 and GHRHR: genetic screening in a large cohort of patients with congenital isolated growth hormone deficiency

CONTEXT: It is estimated that 3-30% of cases with isolated GH deficiency (IGHD) have a genetic etiology, with a number of mutations being reported in GH1 and GHRHR. The aim of our study was to genetically characterize a cohort of patients with congenital IGHD and analyze their characteristics. PATIENTS AND METHODS: A total of 224 patients (190 pedigrees) with IGHD and a eutopic posterior pituitary were screened for mutations in GH1 and GHRHR. To explore the possibility of an association of GH1 abnormalities with multiple pituitary hormone deficiencies, we have screened 62 patients with either multiple pituitary hormone deficiencies (42 pedigrees), or IGHD with an ectopic posterior pituitary (21 pedigrees). RESULTS: Mutations in GH1 and GHRHR were identified in 41 patients from 21 pedigrees (11.1%), with a higher prevalence in familial cases (38.6%). These included previously described and novel mutations in GH1 (C182X, G120V, R178H, IVS3+4nt, a>t) and GHRHR (W273S, R94L, R162W). Autosomal dominant, type II IGHD was the commonest form (52.4%), followed by type IB (42.8%) and type IA (4.8%). Patients with type II IGHD had highly variable phenotypes. There was no difference in the endocrinology or magnetic resonance imaging appearance between patients with and without mutations, although those with mutations presented with more significant growth failure (height, -4.7 +/- 1.6 SDS vs. -3.4 +/- 1.7 SDS) (P = 0.001). There was no apparent difference between patients with mutations in GH1 and GHRHR. CONCLUSIONS: IGHD patients with severe growth failure and a positive family history should be screened for genetic mutations; the evolving endocrinopathy observed in some of these patients suggests the need for long-term follow-up.

The Individualized Genetic Barrier Predicts Treatment Response in a Large Cohort of HIV-1 Infected Patients

The success of combination antiretroviral therapy is limited by the evolutionary escape dynamics of HIV-1. We used Isotonic Conjunctive Bayesian Networks (I-CBNs), a class of probabilistic graphical models, to describe this process. We employed partial order constraints among viral resistance mutations, which give rise to a limited set of mutational pathways, and we modeled phenotypic drug resistance as monotonically increasing along any escape pathway. Using this model, the individualized genetic barrier (IGB) to each drug is derived as the probability of the virus not acquiring additional mutations that confer resistance. Drug-specific IGBs were combined to obtain the IGB to an entire regimen, which quantifies the virus' genetic potential for developing drug resistance under combination therapy. The IGB was tested as a predictor of therapeutic outcome using between 2,185 and 2,631 treatment change episodes of subtype B infected patients from the Swiss HIV Cohort Study Database, a large observational cohort. Using logistic regression, significant univariate predictors included most of the 18 drugs and single-drug IGBs, the IGB to the entire regimen, the expert rules-based genotypic susceptibility score (GSS), several individual mutations, and the peak viral load before treatment change. In the multivariate analysis, the only genotype-derived variables that remained significantly associated with virological success were GSS and, with 10-fold stronger association, IGB to regimen. When predicting suppression of viral load below 400 cps/ml, IGB outperformed GSS and also improved GSS-containing predictors significantly, but the difference was not significant for suppression below 50 cps/ml. Thus, the IGB to regimen is a novel data-derived predictor of treatment outcome that has potential to improve the interpretation of genotypic drug resistance tests.

DOUBLE-STRAND BREAK REPAIR PATHWAYS IN DNA STRUCTURE-INDUCED GENETIC INSTABILITY

Genetic instability in mammalian cells can occur by many different mechanisms. In the absence of exogenous sources of DNA damage, the DNA structure itself has been implicated in genetic instability. When the canonical B-DNA helix is naturally altered to form a non-canonical DNA structure such as a Z-DNA or H-DNA, this can lead to genetic instability in the form of DNA double-strand breaks (DSBs) (1, 2). Our laboratory found that the stability of these non-B DNA structures was different in mammals versus Escherichia coli (E.coli) bacteria (1, 2). One explanation for the difference between these species may be a result of how DSBs are repaired within each species. Non-homologous end-joining (NHEJ) is primed to repair DSBs in mammalian cells, while bacteria that lack NHEJ (such as E.coli), utilize homologous recombination (HR) to repair DSBs. To investigate the role of the error-prone NHEJ repair pathway in DNA structure-induced genetic instability, E.coli cells were modified to express genes to allow for a functional NHEJ system under different HR backgrounds. The Mycobacterium tuberculosis NHEJ sufficient system is composed of Ku and Ligase D (LigD) (3). These inducible NHEJ components were expressed individually and together in E.coli cells, with or without functional HR (RecA/RecB), and the Z-DNA and H-DNA-induced mutations were characterized. The Z-DNA structure gave rise to higher mutation frequencies compared to the controls, regardless of the DSB repair pathway(s) available; however, the type of mutants produced after repair was greatly dictated on the available DSB repair system, indicated by the shift from 2% large-scale deletions in the total mutant population to 24% large-scale deletions when NHEJ was present (4). This suggests that NHEJ has a role in the large deletions induced by Z-DNA-forming sequences. H-DNA structure, however, did not exhibit an increase in mutagenesis in the newly engineered E.coli environment, suggesting the involvement of other factors in regulating H-DNA formation/stability in bacterial cells. Accurate repair by established DNA DSB repair pathways is essential to maintain the stability of eukaryotic and prokaryotic genomes and our results suggest that an error-prone NHEJ pathway was involved in non-B DNA structure-induced mutagenesis in both prokaryotes and eukaryotes.

Genetic and hormonal factors modulate spreading depression and transient hemiparesis in mouse models of familial hemiplegic migraine type 1.

Familial hemiplegic migraine type 1 (FHM1) is an autosomal dominant subtype of migraine with aura that is associated with hemiparesis. As with other types of migraine, it affects women more frequently than men. FHM1 is caused by mutations in the CACNA1A gene, which encodes the alpha1A subunit of Cav2.1 channels; the R192Q mutation in CACNA1A causes a mild form of FHM1, whereas the S218L mutation causes a severe, often lethal phenotype. Spreading depression (SD), a slowly propagating neuronal and glial cell depolarization that leads to depression of neuronal activity, is the most likely cause of migraine aura. Here, we have shown that transgenic mice expressing R192Q or S218L FHM1 mutations have increased SD frequency and propagation speed; enhanced corticostriatal propagation; and, similar to the human FHM1 phenotype, more severe and prolonged post-SD neurological deficits. The susceptibility to SD and neurological deficits is affected by allele dosage and is higher in S218L than R192Q mutants. Further, female S218L and R192Q mutant mice were more susceptible to SD and neurological deficits than males. This sex difference was abrogated by ovariectomy and senescence and was partially restored by estrogen replacement, implicating ovarian hormones in the observed sex differences in humans with FHM1. These findings demonstrate that genetic and hormonal factors modulate susceptibility to SD and neurological deficits in FHM1 mutant mice, providing a potential mechanism for the phenotypic diversity of human migraine and aura.

Genetic predisposition to prostate cancer

Prostate cancer (PC) is a significant economic and health burden in the U.S. and Europe but its causes are largely unknown. The most significant risk factors (after gender) are age and family history of the disease. A gene with high penetrance but low frequency on chromosome 1q, HPC 1, has been suggested to cause a proportion of the familial aggregation of PC but other more common genes, conferring less risk, are also thought to contribute to disease predisposition. We have pursued a strategy to study both types of genetic risk in PC. To identify high penetrance genes, affected men from thirteen families have been genotyped for genetic linkage analysis at six microsatellite markers spanning 45 cM of 1q24-25. Both LOD score and non-parametric statistics provide no significant support for HPC1 in this genomic region, although 3 of the families did combine to produce a LOD score of 0.9. These families will be included in a genome wide search for other PC predisposition genes as part of a multinational collaboration.^ For study of common genetic factors in PC development, leukocyte DNA samples from an unselected series of 55 patients and 67 controls have been examined for genetic differences in two other candidate genes, the androgen receptor gene, hAR, at Xq11-12, and the vitamin D receptor gene, hVDR, at 12q12-14. hAR was typed for two trinucleotide repeat length polymorphisms, (CAG)$\rm\sb{n}$ and (GGC)$\rm\sb{n},$ encoding polyglutamine and polyglycine tracts, respectively, which have been implicated in PC susceptibility. These data, combined with similarly processed patients and controls from the U.K. show no consistent association of allele length with PC risk. A novel finding, however, has been a significant association between the number of GGC repeats and the length of time between diagnosis and relapse in stage T1-T4 Caucasian patients irrespective of therapy and age of the patient. Of 49 patients who relapsed out of 108 entering the study, those with 16 or fewer GGC repeats had an average relapse-free-period of 101 (+/$-$7.7) months while for those with more than 16 repeats the period averaged 48 (+/$-$2.9) months, a difference of 2.1 fold or 4.4 years.^ The second gene, hVDR, was genotyped at two polymorphisms, a synonymous C/T substitution in exon 9 identified by differential TaqI enzymatic digestion and a variable length polyA tract in the 3$\sp\prime$ UTR. Although these polymorphisms are in strong linkage disequilibrium only the polyA region showed a possible association with PC risk. Men homozygous for alleles with fewer than 18 A's had an increased risk (OR = 3.0, p = 0.0578) compared to controls. This result is opposite to the findings of others and may either indicate off-setting random errors which together balance out to no significant overall effect or reflect more complex genetic and/or environmental associations.^ Overall, this research suggests that single gene familial predisposition may be less prominent in PC than in other cancers and that the characteristics of PC pathology may be useful in identifying the effects of common genetic factors. ^

Genetic diversity of bovine Neospora caninum determined by microsatellite markers

Neospora caninum is one of the most significant parasitic organisms causing bovine abortion worldwide. Despite the economic impact of this infection, relatively little is known about the genetic diversity of this parasite. In this study, using Nc5 and ITS1 nested PCR, N. caninum has been detected in 12 brain samples of aborted fetuses from 298 seropositive dairy cattle collected from four different regions in Tehran, Iran. These specimen (Nc-Iran) were genotyped in multilocus using 9 different microsatellite markers previously described (MS4, MS5, MS6A, MS6B, MS7, MS8, MS10, MS12 and MS21). Microsatellite amplification was completely feasible in 2 samples, semi-completely in 8 samples, and failed in 2 samples. Within the two completely performed allelic profiles of Nc-Iran strains, unique multilocus profiles were obtained for both and novel allelic patterns were found in the MS8 and MS10 microsatellite markers. The Jaccard's similarity index showed significant difference between these two strains and from other standard isolates derived from GenBank such as Nc-Liv, Nc-SweB1, Nc-GER1, KBA1, and KBA2. All samples originating from the same area showed identical allelic numbers and a correlation between the number of repeats and geographic districts was observed.

Minimal residual disease in cancer therapy - Small things make all the difference

Minimal residual disease (MRD) is a major hurdle in the eradication of malignant tumors. Despite the high sensitivity of various cancers to treatment, some residual cancer cells persist and lead to tumor recurrence and treatment failure. Obvious reasons for residual disease include mechanisms of secondary therapy resistance, such as the presence of mutant cells that are insensitive to the drugs, or the presence of cells that become drug resistant due to activation of survival pathways. In addition to such unambiguous resistance modalities, several patients with relapsing tumors do not show refractory disease and respond again when the initial therapy is repeated. These cases cannot be explained by the selection of mutant tumor cells, and the precise mechanisms underlying this clinical drug resistance are ill-defined. In the current review, we put special emphasis on cell-intrinsic and -extrinsic mechanisms that may explain mechanisms of MRD that are independent of secondary therapy resistance. In particular, we show that studying genetically engineered mouse models (GEMMs), which highly resemble the disease in humans, provides a complementary approach to understand MRD. In these animal models, specific mechanisms of secondary resistance can be excluded by targeted genetic modifications. This allows a clear distinction between the selection of cells with stable secondary resistance and mechanisms that result in the survival of residual cells but do not provoke secondary drug resistance. Mechanisms that may explain the latter feature include special biochemical defense properties of cancer stem cells, metabolic peculiarities such as the dependence on autophagy, drug-tolerant persisting cells, intratumoral heterogeneity, secreted factors from the microenvironment, tumor vascularization patterns and immunosurveillance-related factors. We propose in the current review that a common feature of these various mechanisms is cancer cell dormancy. Therefore, dormant cancer cells appear to be an important target in the attempt to eradicate residual cancer cells, and eventually cure patients who repeatedly respond to anticancer therapy but lack complete tumor eradication.