Skip to main content

Genetic mapping of novel modifiers for ApcMin induced intestinal polyps’ development using the genetic architecture power of the collaborative cross mice

Abstract

Background

Familial adenomatous polyposis is an inherited genetic disease, characterized by colorectal polyps. It is caused by inactivating mutations in the Adenomatous polyposis coli (Apc) gene. Mice carrying a nonsense mutation in the Apc gene at R850, which is designated ApcMin/+ (Multiple intestinal neoplasia), develop intestinal adenomas. Several genetic modifier loci of Min (Mom) were previously mapped, but so far, most of the underlying genes have not been identified. To identify novel modifier loci associated with ApcMin/+, we performed quantitative trait loci (QTL) analysis for polyp development using 49 F1 crosses between different Collaborative Cross (CC) lines and C57BL/6 J-ApcMin/+mice. The CC population is a genetic reference panel of recombinant inbred lines, each line independently descended from eight genetically diverse founder strains. C57BL/6 J-ApcMin/+ males were mated with females from 49 CC lines. F1 offspring were terminated at 23 weeks and polyp counts from three sub-regions (SB1–3) of small intestinal and colon were recorded.

Results

The number of polyps in all these sub-regions and colon varied significantly between the different CC lines. At 95% genome-wide significance, we mapped nine novel QTL for variation in polyp number, with distinct QTL associated with each intestinal sub-region. QTL confidence intervals varied in width between 2.63–17.79 Mb. We extracted all genes in the mapped QTL at 90 and 95% CI levels using the BioInfoMiner online platform to extract, significantly enriched pathways and key linker genes, that act as regulatory and orchestrators of the phenotypic landscape associated with the ApcMin/+ mutation.

Conclusions

Genomic structure of the CC lines has allowed us to identify novel modifiers and confirmed some of the previously mapped modifiers. Key genes involved mainly in metabolic and immunological processes were identified. Future steps in this analysis will be to identify regulatory elements – and possible epistatic effects – located in the mapped QTL.

Background

Colorectal cancer (CRC) is a complex genetic disease, with many genes influencing the expression of the disease [1]. Mutations in the Adenomatous polyposis coli (Apc) gene are relevant for > 80% of sporadic colorectal adenomas and inherited mutations in the Apc gene cause Familial Adenomatous Polyposis (FAP) syndrome [2]. However, the presence of a mutation in the Apc gene alone cannot explain the wide range of different clinical features of CRC, such as number/size/specific location and onset of polyp development. Environmental factors may contribute to these phenotypic differences, as do modify genes that modulate and regulate the expression and severity of the cancer development [3].

Experimental mouse models of cancer are ideal for examining the effects of genetic modifiers. Modifiers include loci that act, epistatically with known susceptibility loci, (i.e. a mutation in the Apc gene in CRC). Epistasis is difficult to detect in human genome wide association studies (GWAS), due to the very large sample size required. However, it is straightforward approach to engineer mutant mice in which a known susceptibility locus is altered to increase the risk of disease. By crossing the mutant into a population of mice with different genetic backgrounds of naturally occurring variation, it is theoretically possible to unearth the modifier loci.

Three decades ago, a mouse model for intestinal and colorectal cancer research was introduced by Moser et al. 1990; it was created by mutagenesis in germline of C57BL/6 (B6) mice strain and called Min (Multiple intestinal Neoplasia). Mice progeny from this mutated germline suffered from progressive anemia and had visible polyps in large and small intestine. This mice model has allowed further research of intestinal tumorgenesis [4].

Thus far, several genetic modifiers of Min, called Moms have been identified in various mice models, containing mutant versions of ApcMin/+ [5,6,7,8]. The genomic confidence intervals of most of the reported Moms, with few exceptions, were large, which limits the identification of candidate genes underlying these quantitative trait loci (QTL). So far, only two genes underlying two of these Moms were cloned, Pla2g2a for Mom1 and Atp5a1 for Mom2, although their clinical significance is still not clear [7, 8].

In this study, we used a mouse panel formed by crossing ApcMin/+mice with Collaborative Cross (CC) mice [9, 10], in order to map novel Moms. Nowadays, the completed CC comprises a set of ~ 70 Recombinant Inbred Lines (RILs) that were created by full reciprocal matings between 8 different mice strains (the CC founders). These 8 founder strains are genetically diverse, including 5 common laboratory strains: A/J, C57BL/6 J, 129S1/SvImJ, NOD/LtJ, NZO/HiLtJ, and 3 wild-derived strains: CAST/Ei, PWK/PhJ, and WSB/EiJ [11].

The advantages of using CC F1 hybrids for modifier mapping include the numerous genetic variants segregating in the population (there are over 50 million SNPs) [12, 13] e.g. only ~ 4.4 million SNPs segregate between the founders of the BXD panel of RILs [14], and the relatively high level of recombination events compared to two-parent mouse RILs. The three wild-derived founders of the CC represent different subspecies, M.m castaneus, M.m musculus and M.m. domesticus, and contribute many novel sequence variants, not segregating among classical laboratory strains descended from M.m. domesticus [13,14,15]. Many QTLs mapped in CC mice involve allelic contrasts between the wild-derived and laboratory strains [16, 17]. Previous simulation of QTL mapping in CC mice has shown that confidence intervals are typically shorter than 1 Mb [18], and our recent results from variety of studies have shown that it was possible to map the QTL even within less than 1 MB genomic intervals [16, 17].

Methods

Generation of CC- B/6-min mice

In total, 957 F1 mice were produced by a cross of females from 49 CC lines to C57B/6 J-ApcMin/+ males and after PCR analysis for Min genotype, 402 F1 CC-C57BL/6-ApcMin/+ (CC-B/6-ApcMin/+) mice were identified and included in the study for further assessment and analysis. Table 1 shows the list of all the used 49 CC lines and number of mice used from each line. The CC mouse lines were developed and maintained at conventional environmental conditions at the small animal facility of Tel-Aviv University (TAU) and were between generations of G10 to G28 of inbreeding by full-sib mating as, fully described, earlier [11]. The C57BL/6 J- ApcMin/+ mouse line was purchased from the Jackson Laboratory (Bar Harbor, Maine, USA). All experimental mice and protocols were approved by the Institutional Animal Care and Use Committee (IACUC) of Tel-Aviv University (TAU), approval numbers: M-08-075; M-12-024, which adheres to the Israeli guidelines that follow NIH/USA animal care and use protocols.

Table 1 Summary of number of all the used male and female mice of the 49 different lines of the Collaborativ Cross mouse population. # shows the 49 lines; TAU CC lines, shows the TAU designation i.e. Ilxxxx; JAX CCxxx shows the current international CC designation available at JAX laboratory; Male, shows the number of used male mice per line; Female, shows the number of used females per line

All experimental mice were weaned at age of 3 weeks old, housed separately by sex, maximum five mice per cage, with standard rodents’ chow diet (TD.2018SC, Teklad Global, Harlan Inc., Madison, WI, USA, containing % Kcal from Fat 18%, Protein 24%, and Carbohydrates 58%) and water ad libitum. All animals housed in TAU animal facility at conventional open environment conditions, in clean polycarbonate cages with stainless metal covers, and bedded with wood shavings. A Light: dark cycles of 12:12 h, and constant room temperature of 220c (±2). Due to genetic variations between the CC lines, breeding rate, number and sex of litters in each cycle might vary.

Genotyping of CC-B/6-min mice

At 4 weeks old, 0.5 cm tail biopsies were collected from CC-X B/6- ApcMin/+ mice and DNA extracted by NaOH boiling protocol [19]. Mice were genotyped by Polymerase chain reaction (PCR) for the ApcMin/+ mutant allele, using the primers: MAPC-min (TTCTGAGAAAGACAGAAGTTA), MAPC-15 (TTCCACTTTGGCATAAGG), and MAPC-9 (GCCATCCCTTCACGT). For Apc wild type alleles, we used the primers MAPC-15 and MAPC-9, while for the mutant allele we used MAPC-min and MAPC-15 primers [20]. For later identification each mouse was labeled with ear clipping.

Intestinal preparations for polyps count

At the terminal point of the experiment (when mice were 23 weeks old), 402 mice (215 males and 187 females), from 49 CC-B/6-ApcMin/+ lines (n = 1–18 mice per line) were sacrificed by CO2 protocol. Subsequently, small intestines and colons were extracted and washed with Phosphate Buffered Saline (PBS). The small intestines were divided into three segments (SB1-proximal, SB2-middle, and SB3-distal), and the colon was kept as a whole and spread over 3 mm paper. The intestines were fixed in 10% Neutral Buffered Formalin (NBF) overnight and stained by 0.02% methylene blue. The samples were then examined by binocular. The counts and sizes (< 1 mm, 1-2 mm, 2-3 mm, > 3 mm) of polyps in each of the four intestinal sub-regions were recorded as described in Rudling et al. 2006 [21].

Data analysis

Initial statistical analyses were performed using a statistical software package SPSS version 19. One-way Analysis of variance (ANOVA) was performed to test the significance levels of variations in total polyp counts between the different CC-B/6-Min crosses.

CC lines genotype data

High molecular genomic DNA of the CC lines were initially genotyped with the mouse diversity array (MDA), which consists of 620,000 SNPs [22] and re-genotyped by mouse universal genotype array (MUGA-7500 markers) and eventually with MegaMuga (77,800 markers) SNP arrays to confirm their genotype status [12]. The genotype database used in this study is, publically available at: http://mtweb.cs.ucl.ac.uk/mus/www/preCC/CC-2018/LIFTOVER/CONDENSED/.

Data analysis was performed using the statistical software R (R Development Core Team 2009), including the R package HAPPY.HBREM [23].

Reconstruction of CC ancestral genome mosaics

We removed SNPs with heterozygous or missing genotypes in the 8 CC founders, or were not in common between the arrays, leaving 170,935 SNPs. The SNPs were mapped onto build 37 of the mouse genome. We reconstructed the genome mosaic of each CC line in terms of the eight CC founders using a hidden Markov Model HAPPY ([23] across the genotypes to compute probabilities of descent from founders, setting the generation parameter to g = 7. To allow for genotyping error, we configured the HMM to allow a small probability of 0.001 that any founder was consistent with any SNP allele. The HAPPY HMM computed a descent probability distribution for each of the 170 k SNP intervals, which we reduced to 8533 intervals by averaging the matrices in groups of n = 20 consecutive SNPs. This reduction reduced further the effects of genotyping error and made analyses faster. Mean heterozygosity was computed across each window of 20 SNPs.

The locus-specific fraction of CC lines carrying each of the founders was estimated by summing the HMM posterior probabilities at each interval across all lines. Genome wide thresholds for significance were computed by permuting the identities of the founders separately within each line, then recomputing the locus-specific fractions and recording the genome wide maximum and minimum fractions in the permuted data. This process was repeated 200 times to estimate the upper and lower thresholds exceeded in 10% of permutations.

QTL analysis

The genome of each CC line is a mosaic of the inbred founders, which we reconstructed using a hidden Markov model implemented in the HAPPY R package across the genotypes to compute probabilities of descent from the founders [13, 23]. The presence of a QTL at a given locus was tested using the probabilities of descent from each founder calculated through HAPPY and testing for association between the founder haplotype at each locus and the median polyp count within each CC line, using multiple linear regression. Sex was included as a covariate. QTL effect sizes were estimated as the proportion of the log-likelihood explained by the locus effects at the QTL. Genome-wide significance was estimated by permutation, where the CC line labels were permuted between the phenotypes. Permutation-based false discovery rate (FDR) was calculated for a given P-value threshold, following the formula: (expected number of false discoveries)/ (number of observed discoveries).

Testing sequence variation segregating between the CC founders

Except for a small number of de-novo mutations arising during breeding, all sequence variants segregating in the CC should also segregate in the CC founders. Therefore we use the merge analysis methodology [24] to test which variants under a QTL peak were compatible with the pattern of action at the QTL. A variant with A alleles inside the locus L merges the 8 CC founders into A < 8 groups according to whether they share the same allele at the variant (A = 2 in the case of SNPs). This merging is characterized by an 8xA merge matrix Msa defined to be 1 when strain s carries allele a, and 0 otherwise. The effect of this merging is tested by comparing the fit of the QTL model above with one in which the Nx8 matrix XLis is replaced by the NxA matrix Zia = Σs XLis Msa. We use the Perlegen SNP database to test sequence variants globally and the Sanger SNP database for individual genes. This approach was, successfully applied in our previous studies [16, 17, 24].

Estimation of QTL confidence intervals

The confidence intervals of the QTL were estimated through simulation of a QTL with a similar logP and strain effects in the neighborhood (5 Mb) of the observed QTL peak, using a similar approach as presented in our previous studies [16, 17] to take into account local patterns of linkage disequilibrium. Briefly, accurate estimates of QTL mapping resolution should take into account local patterns of linkage disequilibrium. We devised a method that preserved the genotypes of the data, whilst simulating survival times caused by a QTL in the neighborhood (5 Mb) of the observed QTL peak, and with a similar logP to that observed. We first extracted the parameter estimates \( {\hat{\beta}}_s \) and residuals \( {\hat{r}}_i \) of the fitted polyp counts model at the QTL peak. Let \( {\hat{t}}_i \) be a random permutation of \( {\hat{r}}_i \). Then in a marker interval K within 5 Mb of the QTL peak L we simulated a set of survival times ZiK caused by a QTL at K by substituting the parameter estimates and permuted residuals:

$$ {Z}_{iK}=t\hat{\mkern6mu} i\kern0.5em \exp \left(\upmu \hat{\mkern6mu} +\Sigma s\ {X}_{Kis}\upbeta \hat{\mkern6mu} s\right) $$

We then rescanned the region and found the interval with the highest logP. We simulated 1000 QTLs at each interval K and estimated the p% CI from interval containing p% of the simulated local maxima.

Founder effects

Except for a small number of de-novo mutations arising during breeding, all sequence variants segregating in the CC lines should also segregate in the CC founders. The founder strain trait effects at each QTL were shown relatively to WSB/EiJ, using a similar approach as presented in our previous studies [16, 17]. Briefly, except for a small number of de-novo mutations arising during breeding, all sequence variants segregating in the CC should also segregate in the CC founders. Therefore, we use the merge analysis methodology [24] to test which variants under a QTL peak were compatible with the pattern of action at the QTL. A variant with A alleles inside the locus L merges the 8 CC founders into A < 8 groups according to whether they share the same allele at the variant (A = 2 in the case of SNPs). This merging is characterized by an 8xA merge matrix Msa defined to be 1 when strain s carries allele a, and 0 otherwise. The effect of this merging is tested by comparing the fit of the QTL model above with one in which the Nx8 matrix XLis is replaced by the NxA matrix Zia = Σs XLis Msa. We use the Perlegen SNP database (http://mouse.perlegen.com/mouse/download.html) to test sequence variants globally and the Sanger mouse genomes database (http://www.sanger.ac.uk/resources/mouse/genomes/) for individual genes.

Within the QTLs we classified the sequence variants according to the genome annotation as repetitive, intergenic, upstream, downstream, UTR, intronic or coding. We then classified variants according to whether their merge logP was greater or less than the corresponding haplotype-based logP. The enrichment of variants with high logP values within each category was computed.

List of suggested candidate genes

We used the SNP tools package in R, and the MGI database (http://www.informatics.jax.org) to find all the genes in the 95% confidence interval for each QTL. We focused on protein-coding genes in these regions, but also non-coding RNA genes, such as miRNA loci. Also, if the 3′ UTR or the 5′ UTR of a gene were inside the interval then we included the gene in our list. We used these candidate gene lists as an input for BioInfoMiner.

Functional analysis with BioInfoMiner

We performed functional pathway analysis using BioInfoMiner [25]. BioInfoMiner (https://bioinfominer.com) performs statistical and network analysis on biological hierarchical vocabularies to detect and rank significantly enriched processes and the underlying hub genes involved in these processes. For our analysis, we used Gene Ontology (GO) [26], Reactome [27] and MGI Mammalian Phenotype (MGI) [28]. The BioInfoMiner algorithm maps the genes in the supplied gene list to a semantic network created from ontological data, corrected through AI-inspired semantic network pruning and clustering and then prioritizes the genes based on the topological properties of the thus corrected network. This analysis prioritized genes with central functional and regulatory roles in enriched processes, underlying the studied phenotype. The correction for potential semantic inconsistencies on the selected ontological scheme and bias mitigation regarding the different depth of the branches of the semantic tree, as a result of differences in knowledge representation for distinct scientific concepts, was performed by restoring the order of the resolution of annotation of each gene with its ancestral ontological terms.

Results

Polyp counts

We mapped QTL modifiers of ApcMin/+ based on polyp counts in the small intestine and colon, in 49 CC-B/6-ApcMin/+ lines at 23 weeks old (n = 402 mice) (see Fig. 1). The overall population mean of total polyp counts was 32.48 ± 1.36 polyps, ranging widely from 9 polyps (IL1286) to 88 polyps (IL2288). Based on one-way ANOVA, significant variation (p < 0.01) was found between different 49 CC-B/6-ApcMin/+ lines in their total counted polyps. Polyp counts were approximately normally distributed, suggesting the intervention of numerous genetic and environmental factors in this trait. The mean of polyp number for the parental line B/6-ApcMin/+(n = 5) (first column Fig. 1) was 64.25 ± 6.65 polyps. The majority of CC-B/6-ApcMin/+ lines (30/49, 61%) had lower polyp counts compared to the B/6-ApcMin/+ parental line suggesting that resistant alleles for intestinal tumorigenesis segregate among the CC lines. We also investigated if different segments of the intestine exhibited differential polyp distribution and different genetic architectures. The small intestine was subdivided into 3 sections (small intestine proximal-SB1, middle-SB2, and distal-SB3), and the colon was treated separately. Overall polyps were distributed approximately equally between 3 segments of the small intestine: SB1 with 8.12 ± 0.45 polyps (25%), SB2 with 9.25 ± 0.53 (28.48%), SB3 with 9.37 ± 0.48 (28.85%) and the colon was with 5.7 ± 0.19 (17.47%).

Fig. 1
figure 1

Polyp’s count (±SE) inSB1, SB2, SB3, and Colon of CC-F1 mice crosses at the age of 23 wks. Old (terminal time point). Data analysis of 49 CC-B/6-Min F1 hybrid lines (n = 1–18 mice/line) and C57BL/6-Minstrain (4 mice). The Y-axis represents the number of polyps; X-axis represents different APC-min F1 hybrid lines. First column represents C57BL/6 carrying the ApcMin/+ mutation (first column) and mean of the CC-B/6-Min population. Data analyzed by One-way Analysis of Variation (ANOVA), *p-value< 0.05

QTL analysis

QTL analysis using HAPPY [16, 17, 23] was performed for polyp count traits sub-divided into three parts of the small intestinal tract (SB1, SB2, and SB3) and colon, for the 402 mice in 49 CC-B/6- ApcMin/+ F1 crosses, including males and females. Nine significant QTLs at the genome-wide significance threshold of 90% (i.e. where < 10% of permutations had a genome-wide maximum exceeding an observed QTL score) were detected (Table 2). Five of these QTLs were significant at the more stringent 95% level of genome-wide significance.

Table 2 Genomic location of the significant Quantitative Trait Loci (QTL) at 90 and 95% genome wide significant thresholds associated with polyp counts in SB1, SB2, SB3, Colon and total polyps in the entire intestines (SB1–3 and colon) regions of different CC lines. QTL associated with polyp counts detected on different chromosomal regions. Experiment-wide thresholds of significance at *P% of 50, 90 and 95% levels are presented for each trait, accordingly

In the proximal section of the small intestine, SB1, (Fig. 2A), a significant QTL (95%) was found on chromosome 3, peak at 13.839 Mb, logP =4.43, designated Mom19. Another significant QTL (90%) was found on chromosome 12, peak at 111.37 Mb, logP = 3.71, designated Mom20. For SB2, (Fig. 2B), a significant QTL (95%) was found on chromosome 10, peak at 18.805 Mb, logP =4.11, designated Mom21. Additionally, two well-separated significant QTLs (95%) for SB2 were found on chromosome 16, peak at 53.51 Mb (Mom22) and 73.216 Mb (Mom23), logP > 4. For SB3, (Fig. 2C), two significant QTLs (95%) were found on chromosome 6 and chromosome 12, peak at 146.203 Mb (Mom24) and 113.449 Mb (Mom25) respectively, logP > 4.2. Further, two QTLs (90%) were found on chromosome 9, peak at 37.55 Mb, logP = 3.9, on chromosome 10 same location as Mom21. For polyp’s count in colon, Fig. 2D, a solo significant QTL (95%) was found on chromosome 6, peak at 35.91 Mb, logP = 4.19, designated Mom27. For total polyp counts, Fig. 2E, a significant QTL (95%) was mapped to same locations of Mom20, Mom22 and Mom23.

Fig. 2
figure 2

Genome scans for counted polyps with significant QTLs. A Polyps counted in SB1. B Polyps counted in SB2. C Polyps counted in SB3. D Polyps counted in the colon. E Total counted polyps. The X-axis represents the 19 mouse chromosomes and the position of mapped QTLs. Y-axis represents the logP of the test of association between locus and polyps count

In summary, nine distinct and novel QTLs at 90 and 95% genome-wide significant thresholds levels (GWSL). These QTLs were designated as modifiers of Min gene (Mom) numbers 19–27, respectively, presented in Table 2 the location of the peak, the interval and number of genes within each of QTLs.

Finally, we lowered the genome-wide significant thresholds to the 50% level (i.e. where there was a probability of 50% that a QTL exceeding the threshold was a false positive), to identify potential candidate QTLs, which might be genome-wide significant in a larger and more powerful mapping populations. Table 3 summarizes all the mapped QTL at 50% threshold level, their chromosomal locations, 90 and 95% confidence intervals and number of genes identified within these intervals. The analysis has detected four QTL associated with SB1 trait, where two of them were mapped at different positions on chromosome 7 (peaks at 87.273 and 103,732 Mb), one QTL mapped at chromosome 6 (peak at 145.956 Mb), and one on chromosome 14 (peak at 73.794 Mb).

Table 3 Chromosomal locations of the mapped Quantitative Trait Loci (QTL) at 50% genome wide significant threshold level and found to be associated with polyp counts in SB1, SB2, SB3, Colon and total polyps in the entire intestines (SB1–3 and colon) regions of different CC lines. QTL associated with polyp counts detected on different chromosomal regions. Experiment- confidence intervals (CI) of significance at *P% of 90 and 95% levels are presented for each trait, accordingly

Two QTL were detected with SB2 trait and mapped on chromosomes 6 and 12 at positions of 146.284 and 113.669 Mb, respectively. Six QTL were detected with SB3 trait, while two of them were mapped at different positions on chromosome 16 (peaks at 53.511 and 73.271 Mb), a single QTL was mapped on chromosomes 3, 4, 7 and 18 at positions 12.197, 140.994, 145.941 and 40.491, respectively. One QTL was detected with polyps in the colon and mapped at chromosome 1 its peak was located at 34.282 Mb. Finally, six QTL were detected with total polyp’s trait, and mapped at chromosomes 3, 4, 6, 10, 14 and 18, and its peaks were located at 11.832, 126.165, 146.233, 19.133, 76.654 and 10.631, respectively. These 90 and 95% confidence intervals of the identified QTL were ranged between 10 and 20 Mb, and number of genes identified within these intervals were ranged between 150 to 670.

Founder effects

The effects of each founder haplotype on polyp counts for the mapped QTLs were evaluated as deviation relative to the WSB/EiJ parental strain, which was arbitrary assigned the baseline zero effect. All the data presented in Table 4. For Mom19 there were slight positive effects on poly counts for CAST/EiJ, NZO/HILtJ, 129S1/SvImJ strains and minimal negative effects for A/J, C57BL/6 J, NOD/LtJ, and PWK/PhJ. For Mom20 all founder strains have positive effects except A/J and C57BL/6 J. For Mom21 all the founder strains contributed a positive effect on polyps count (i.e. this QTL involved a contrast between WSB/eiJ vs the rest). For Mom22and Mom23 all the founder strains except PWK/PhJ contributed positive effects. For Mom24 positive effects were seen in all the founder strains, except CAST/EiJ and 129S1/SvImJ. For Mom25 only A/J and C57BL/6 J strain had a minor negative effect on polyps count. For Mom26 and Mom27 C57BL/6 J, CAST/EiJ and 129S1/SvImJ had a negative effect.

Table 4 The estimated strain effects on polyp count for the 8 CC founder strains for each of the mapped Quantitative Trait Loci (QTL) Mom19-Mom27, which were mapped at 90 and 95% genome wide significant thresholds levels. Effects are shown as deviations relative to WSB/EiJ, which is arbitrarily assigned the trait effect

Merge analysis

The haplotype QTL analysis was then refined by merge analysis in order to identify SNPs within each QTL whose strain distribution patterns among the founder strains were consistent with the patterns of action at the QTL. Results are presented in Fig. 3. In every plot two vertical lines that delineate the location of the mapped QTL. In three plots of Fig. 3 of SB2 Chr16, SB3 Chr6 and SB3 Chr9, the locations of SNPs were mapped at the same interval as the mapped QTLs, Mom23, Mom24, Mom26. For SB1, we did not find any SNPs with logP> 4 within the mapped QTLs, suggesting that the effect was not driven by a single biallelic variant but instead was haplotype-based. For SB2 significant SNPs were found on chromosome 16, same location as Mom23. For SB3 significant SNPs were found on chromosome 9 same location as Mom26, on chromosomes 10 and 12 significant SNPs were found but outside of mapped QTLs. It is interesting to note that significant SNPs on chromosome 10 were mapped to the same location as in the previously mapped Mom17 [29]. For colon polyps count, unfortunately, we did not find significant SNPs within the mapped QTLs. For total polys count, significant SNPs were found on chromosome 12 but outside of mapped QTL interval, on chromosome 16 significant SNPs were found on same locations as Moms22 and 23.

Fig. 3
figure 3

Merge analysis of sequence variants around the mapped QTL, for susceptibility to polyp development in the CC lines. The X-axis represents genomic location, the Y-axis represents the logP of the test of association between locus and polyps count. The continuous red line is the genome scan of different mapped QTL in Fig. 2. The black dots are the results of merge analysis tests of sequence variants segregating in the 8 founders of the CC lines

Candidate genes

Our results have shown that based on the mapped QTL analysis, we have 1862 unique genes in the intervals of all the mapped QTL. When the Bioinfominer analysis was performed, while searching for genes based on functional pathway analysis, which prioritized genes with central functional and regulatory roles in enriched processes, underlying the studied phenotype, we identified 411 unique genes related to the studied traits (supplement 1).

Merge analysis identified candidate functional SNPs that may play a role in each QTL, some of which were outside the haplotype-based QTL regions. We used the MGI database to select genes nearest these SNPs. Being outside of the original QTL intervals, some of the are not in the original candidate gene lists based on the BioInfoMiner analysis. We found a total of 123 candidate genes, presented in Table 5, but only two of them are topic relevant. These two genes are mapped to a SNP on chromosome 10 which is associated with SB3 polyp counts. The first is colon tumor susceptibility 9 (Scc9) locus, which was previously mapped in BALB/c × CcS crosses [30], and the second gene is angiogenesis by VEGF QTL 1(Angvq1), which was mapped by using BXD strains [31].

Table 5 Candidate gene list identified within the mapped QTL, Mom19–27, based on merge analysis and using MGI data base

Finally, when we combined these results with the merge analysis, we shortened the list to 123 genes (Table 5), only, while some are overlapping between the three approaches (i.e. QTL analysis at 90 and 95% genome-wide significance threshold (1862 genes), BioInfoMiner (411 genes), and Merge analysis (123 genes). These approaches have shown the power of identifying of candidate genes, which may lead to future plans for further studies with these genes.

Discussion

Colorectal cancer is a complex disease, with many genes modifying the expression of the phenotype. The presence of mutations in the Apc gene alone cannot explain the wide range of different clinical features observed. It is well documented that modifying genes (host genes that modulate and regulate the expression and severity of the cancer development) have a crucial role on tumorgenesis [32]. The modifier genes of CRC have been studied in human GWAS and numerous loci found [33,34,35,36]. The combined effects of these mapped and identified alleles are currently too small to explain the bulk of heritable disease risk [33,34,35,36]. These suggest that genetic influence towards cancer susceptibility cannot be unraveled solely using approaches designed to identify the main effects of individual alleles in human populations.

Experimental mouse models are ideal for examining the effects of genetic modifiers. By crossing the mutant into a population of mice with different genetic backgrounds of naturally occurring variation, it is possible to map modifier loci. Several earlier studies [5, 29, 37, 38] on mouse models have been performed that mapped 18 Moms. The most tightly mapped of these QTLs (width 7.4 Mb) was found for Mom7 [37], while the rest are between 16 to 53 Mb.

Here, we present a genetic analysis of intestinal polyp counts in 49 F1 CC-B/6-ApcMin/+crosses to search for Moms. We observed wide heritable variations in polyp counts between the 49 crosses, in accordance with our previous study that showed mice with different genetic backgrounds vary in their progression of intestinal polyp development [39]. Polyp counts in F1 CC-B/6-ApcMin/+crosses differ from those in the parental line (B/6-ApcMin/+). This suggests that the CC population contains modifiers that either suppressor enhance ApcMin/+ mutation, which might be caused by the high genetic diversity of the three wild-derived strains. This variation enabled us to identify new Mom QTL and improve the resolution of previously mapped modifiers. By using founder effect analysis, we found the haplotypic effects of founder CC stains varied between QTL and could have positive or negative effects.

The density of polyps varies across the gastrointestinal tract, suggesting tumorigenesis in different parts of the intestine is controlled by different genes and with distinct genetic architectures. Each part of the intestinal tract has specific physiological functions, with different gene expression profiles, pH and microbiota [40]. In our study we counted polyps in each of the proximal (SB1), middle (SB2) or distal (SB3) parts of small intestine and colon. Most polyps were found within the small intestine, in accordance to a previous study that showed ApcMin/+ mice usually develop polyps in the small intestine, unlike Familial adenomatous polyposis (FAP) patients [41]. Within the small intestine, there was no preferred location for polyp development.

We mapped nine distinct and novel Mom QTLs at the 90% genome-wide significance threshold, and at least an additional 16 more potential QTL at 50% threshold, but will not be fully discussed here. These QTLs are designated as modifiers of Min gene (Mom) numbers 19–27, respectively. We found different Moms to be responsible for polyp development in different parts of the intestinal tract. The wild-derived CC founder strains contain genetic variations absent from standard laboratory mouse strain (SLMS), explaining why we were able to map novel QTLs. However, in this study we were able to map some QTL associated with SLMS, as well.

Most of the Moms mapped in this study are distinct from those identified previously [29]. However, we mapped Mom27 for colon polyp development on chromosome 6:34.72–38.33 Mb, which overlap with a previously mapped Mom12 with wider range of 6:17.3–50.8 Mb [38]. Additionally, Mom26 which mapped on chromosome 9:32.56–42.56 Mb partly overlaps Colon Cancer loci susceptibility 4 (Ccs4) mapped to 36.84–49.23 Mb [42]. We mapped Mom21 to 10:8.90–28.47 Mb which does not overlap Mom17 (10:69–89 Mb) [8], so these probably represent distinct loci. The rest of our genome-wide significant QTLs are novel. Some of the 16 potential QTLs at 50% genome-wide significance [Table 3] overlap with previously mapped Moms i.e. Mom1 [5] and Mom2 [45] on chromosomes 4 and 18, respectively. |.

We have identified candidate genes underlying these traits. In this study we adopted three approaches for identifying and suggesting candidate genes (i.e. QTL analysis at 90 and 95% genome-wide significance threshold (1862 genes), BioInfoMiner (411 genes), and Merge analysis (123 genes), including list of genes identified in the genomic intervals of the mapped QTL, based on mouse genome data base, Bioinfominer analysis using Gene Ontology (GO) [26], Reactome [27] and MGI Mammalian Phenotype (MGI) data, so to identify prioritized genes with central functional and regulatory roles in enriched processes, underlying the studied phenotype, and finally based on the merge analysis. Indeed, each approach has suggested different number of genes, while the lowest number was obtained by merge analysis (123 genes).

This report and our previous study [16, 17, 39] demonstrate the utility of the CC lines in the analysis of complex traits in mouse models of human disease. This is study showed the power of using CC mice to dissect the genetic response to intestinal cancer development and the first to use of the CC F1 cross design for modifier mapping. Even a modest number of lines [16, 17] are useful with sufficient replication (3–5 mice) within each line. Similar recent studies were reported by using this CCXMutant F1 approach for defining the genetic mechanisms of host susceptibility to melanoma [43, 44].

The genomic intervals of the mapped QTL in this study were small enough to suggest candidate genes, although further confirmation work is required, including knockout or knockdown of specific candidate genes analysis. Many of these candidates are involved in innate and adaptive immune responses.

Conclusions

Variation in polyp development is heritable and controlled, to an appreciable extent, by genetic factors segregating in the CC population which is therefore well-suited for identifying novel modifier genes associated with ApcMin/+mutation. The expected findings from this study may be used for early prediction of potential intestine cancer development in host carrying susceptible genetic factors, thus can be applied for better control and sufficient application therapy tools and approaches.

Availability of data and materials

Availability of Data:

All phenotype, genotype, and bioinformatics analysis programs should be freely following DOIs:

phenoPolyps.txt https://0-doi-org.brum.beds.ac.uk/10.5522/04/12790100

R codes https://0-doi-org.brum.beds.ac.uk/10.5522/04/12790160

Genotypes https://0-doi-org.brum.beds.ac.uk/10.5522/04/12790187

Abbreviations

CEO:

Chief executive officer

SME:

Small medium enterprise

UK:

United of Kingdom

APC:

Adenomatous polyposis coli gene

Min:

Multiple intestinal neoplasia

APC Min/+ :

Mice carrying a nonsense mutation in the Apc gene at R850, which is designated ApcMin/

CRC:

Colorectal cancer

MOM:

modifier loci of Min

FAP:

Familial Adenomatous Polyposis

Mb:

Megabases

CC:

Collaborative Cross

GRP:

Genetic Reference Panel

HFD:

High-Fat Diet

RIL:

Recombinant Inbred Lines

GWAS:

Genome wide association studies

QTL:

Quantitative Trait Loci

TAU:

Tel-Aviv University

SE:

Standard Error

GO:

Gene Ontology

MGI:

Mouse Genome Informatics

CI:

Confidence Interval

BXD:

C57BL/6 JXDBA2/J mice

M.m :

Mus musculus

ANOVA:

Analysis of variance

SB:

Small intestinal tract

GWSL:

genome-wide significant thresholds levels

Go:

Gene Ontology

References

  1. Weigl K, Chang-Claude J, Knebel P, Hsu L, Hoffmeister M, Brenner H. Strongly enhanced colorectal cancer risk stratification by combining family history and genetic risk score. Clin Epidemiol. 2018;10:143–52. https://0-doi-org.brum.beds.ac.uk/10.2147/CLEP.S145636.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Fodde R. The APC gene in colorectal cancer. Eur J Cancer. 2002;38(7):867–71. https://0-doi-org.brum.beds.ac.uk/10.1016/S0959-8049(02)00040-0.

    Article  CAS  PubMed  Google Scholar 

  3. Wang F, Wong SC, Chan LW, Cho WC, Yip SP, Yung BY. Multiple regression analysis of mRNA-miRNA associations in colorectal cancer pathway. Biomed Res Int. 2014;2014:676724.

    PubMed  PubMed Central  Google Scholar 

  4. Moser AP, Pitot HC, Dove FW. A dominant mutation that predisposes to multiple intestinal neoplasia in the mouse. Science. 1990;247(4940):322–4. https://0-doi-org.brum.beds.ac.uk/10.1126/science.2296722.

    Article  CAS  PubMed  Google Scholar 

  5. Dietrich WF, Lander ES, Smith JS, Moser AR, Gould KA, Luongo C. Genetic identification of mom-1 a major modifier locus affecting min induced intestinal neoplasia in the mouse. Cell. 1993;75(4):631–9. https://0-doi-org.brum.beds.ac.uk/10.1016/0092-8674(93)90484-8.

    Article  CAS  PubMed  Google Scholar 

  6. Gould KA, Dove WF. Analysis of the Mom1 modifier of intestinal neoplasia in mice. Exp Lung Res. 1998;24(4):437–53. https://0-doi-org.brum.beds.ac.uk/10.3109/01902149809087379.

    Article  CAS  PubMed  Google Scholar 

  7. MacPhee M, Chepenik KP, Liddell RA, Nelson KK, Siracusa LD, Buchberg AM. The secretory phospholipase A2 gene is a candidate for the Mom1 locus, a major modifier of ApcMin-induced intestinal neoplasia. Cell. 1995;81(6):957–66. https://0-doi-org.brum.beds.ac.uk/10.1016/0092-8674(95)90015-2.

    Article  CAS  PubMed  Google Scholar 

  8. Prami C, Savelyeva L, Le Pastier D, Siracusa LD, Buchberg AM, Schwab M, et al. Human homologue of a candidate for the Mom1 locus, the secretory type II phospholipase A2 (PLA2S-II), maps to 1p35-36.1/D1S19. Cancer Res. 1995;55:5504–6.

    Google Scholar 

  9. The Complex Trait Consortium. The collaborative cross, a community resource for the genetic analysis of complex traits. Nat Genet. 2004;36(11):1133–7. https://0-doi-org.brum.beds.ac.uk/10.1038/ng1104-1133.

    Article  CAS  Google Scholar 

  10. Complex Trait Community website: http://www.complextrait.org.

  11. Chesler EJ, Miller DM, Branstetter LR, Galloway LD, Jackson BL. The collaborative cross at oak Ridge National Laboratory: developing a powerful resource for systems genetics. Mamm Genome. 2008;19(6):382–9. https://0-doi-org.brum.beds.ac.uk/10.1007/s00335-008-9135-8.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Collaborative Cross Consortium. The genome architecture of the collaborative cross mouse genetic reference population. Genetics. 2011;190:389–401.

    Article  Google Scholar 

  13. Keane TM, Goodstadt L, Danecek P, White MA, Wong K, Yalcin B, et al. Mouse genomic variation and its effect on phenotypes and gene regulation. Nature. 2011;77(7364):289–94.

    Article  Google Scholar 

  14. Roberts A, Pardo-Manuel de Villena F, Wang W, McMillan L, Threadgill DW. The polymorphism architecture of mouse genetic resources elucidated using genome-wide resequencing data: implications for QTL discovery and systems genetics. Mamm Genome. 2007;18(6–7):473–81. https://0-doi-org.brum.beds.ac.uk/10.1007/s00335-007-9045-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Beck JA, Lloyd S, Hafezparast M, Lennon-Pierce M, Eppig JT, Festing MFW, et al. Genealogies of mouse inbred strains. Nat Genet. 2000;24(1):23–5. https://0-doi-org.brum.beds.ac.uk/10.1038/71641.

    Article  CAS  PubMed  Google Scholar 

  16. Durrant C, Tayem H, Yalcin B, Cleak J, Goodstadt L, Pardo-Manuel de Villena F. Collaborative cross mice and their power to map host susceptibility to Aspergillus fumigates infection. Genome Res. 2011;21(8):1239–48. https://0-doi-org.brum.beds.ac.uk/10.1101/gr.118786.110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Vered K, Durrant C, Mott R, Iraqi FA. Susceptibility to klebsiella pneumonaie infection in collaborative cross mice is a complex trait controlled by at least three loci acting at different time points. BMC Genomics. 2014;15(1):865. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-15-865.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Valdar W, Flint J, Mott R. Simulating the collaborative cross: power of quantitative trait loci detection and mapping resolution in large sets of recombinant inbred strains of mice. Genetics. 2006;172(3):1783–97. https://0-doi-org.brum.beds.ac.uk/10.1534/genetics.104.039313.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Truett GE, Heeger P, Mynatt RL, Truett AA, Walker JA, Warman ML. Preparation of PCR-quality mouse genomic DNA with hot sodium hydroxide and TRIS (HotSHOT). Biotechniques. 2000;29(1):52–4. https://0-doi-org.brum.beds.ac.uk/10.2144/00291bm09.

    Article  CAS  PubMed  Google Scholar 

  20. Ju J, Hong J, Zhou JN, Pan Z, Bose M, Liao J. Inhibition of intestinal tumorigenesis in ApcMin/+ mice by (−)-epigallocatechin-3-gallate, the major catechin in green tea. Cancer Res. 2005;65(22):10623–31. https://0-doi-org.brum.beds.ac.uk/10.1158/0008-5472.CAN-05-1949.

    Article  CAS  PubMed  Google Scholar 

  21. Rudling R, Hassan AB, Kitau J, Mandir N, Goodlad RA. A simple device to rapidly prepare whole mounts of murine intestine. Cell Prolif. 2006;39(5):415–20. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1365-2184.2006.00391.x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Yang H, Ding Y, Hutchins LN, Szatkiewicz J, Bell TA, Paigen B. A customized and versatile high-density genotyping array for the mouse. Nat Methods. 2009;6(9):663–6. https://0-doi-org.brum.beds.ac.uk/10.1038/nmeth.1359.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Mott R, Talbot CJ, Turri MG, Collins AC, Flint J. A method for fine mapping quantitative trait loci in outbred animal stocks. Proc Natl Acad Sci. 2000;97(23):12649–54. https://0-doi-org.brum.beds.ac.uk/10.1073/pnas.230304397.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Yalcin B, Flint J, Mott R. Using progenitor strain information to identify quantitative trait nucleotides in outbred mice. Genetics. 2005;171(2):673–81. https://0-doi-org.brum.beds.ac.uk/10.1534/genetics.104.028902.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Koutsandreas T, Binenbaum I, Pilalis E, Valavanis I, Papadodima O, Chatziioannou A. Analyzing and visualizing genomic complexity for the derivation of the emergent molecular networks. Int J Monit Surveill Technol Res. 2016;4(2):30–49. https://0-doi-org.brum.beds.ac.uk/10.4018/IJMSTR.2016040103.

    Article  Google Scholar 

  26. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9. https://0-doi-org.brum.beds.ac.uk/10.1038/75556.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Fabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 2018;46(D1):D649–55. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkx1132.

    Article  CAS  PubMed  Google Scholar 

  28. Smith CL, Eppig JT. The mammalian phenotype ontology as a unifying standard for experimental and high-throughput phenotyping data. Mamm Genome. 2012;23(9–10):653–68. https://0-doi-org.brum.beds.ac.uk/10.1007/s00335-012-9421-3.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Nnadi SC, Watson R, Innocent J, Gonye GE, Buchberg AM, Siracusa LD. Identification of five novel modifier loci of Apc (min) harbored in the BXH14 recombinant inbred strain. Carcinogenesis. 2012;33(8):1589–97. https://0-doi-org.brum.beds.ac.uk/10.1093/carcin/bgs185.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. van Wezel T, Ruivenkamp CAL, Stassen APM, et al. Four New Colon Cancer Susceptibility Loci, Scc6 to Scc9 in the Mouse. Cancer Res. 1999;59:4216–8.

    PubMed  Google Scholar 

  31. Nakai K, Rogers MS, Baba T, Funakoshi T, Birsner AE, Luyindula DS, et al. Genetic loci that control the size of laser-induced choroidal neovascularization. FASEB J. 2009;23(7):2235–43. https://0-doi-org.brum.beds.ac.uk/10.1096/fj.08-124321 Epub 2009 Feb 23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kinzler WK, Vogelstein B. Lessons from hereditary colorectal Cancer. Cell. 1996;87(2):159–70. https://0-doi-org.brum.beds.ac.uk/10.1016/S0092-8674(00)81333-1.

    Article  CAS  PubMed  Google Scholar 

  33. Houlston RS, Webb E, Broderick P, Pittman AM, Di Bernardo MC, Lubbe S, et al. Meta-analysis of genome-wide association data identifies four new susceptibility loci for colorectal cancer. Nat Genet. 2008;40(12):1426–35.

    Article  CAS  PubMed  Google Scholar 

  34. Tenesa A, Farrington SM, Prendergast JGD, Porteous ME, Walker M, Haq N, et al. Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21. Nat Genet. 2008;40(5):631–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Tomlinson IPM, Webb E, Carvajal-Carmona L, Broderick P, Howarth K, Pittman AM, et al. A genome-wide association study identifies colorectal cancer susceptibility loci on chromosomes 10p14 and 8q23.3. Nat Genet. 40:623–30.

  36. Tomlinson I, Webb E, Carvajal-Carmona L, Broderick P, Kemp Z, Spain S, et al. A genome-wide association scan of tag SNPs identifies a susceptibility variant for colorectal cancer at 8q24.21. Nat Genet. 39:984–8.

  37. Kwong LN, Shedlovsky A, Biehl BS, Clipson L, Pasch CA, Dove WF. Identification of Mom7, a novel modifier of ApcMin/+ on mouse chromosome 18. Genetics. 2007;176(2):1237–44. https://0-doi-org.brum.beds.ac.uk/10.1534/genetics.107.071217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Crist RC, Roth JJ, Lisanti MP, Siracusa LD, Buchberg AM. Identification of Mom12 and Mom13, two novel modifier loci of Apc (min) -mediated intestinal tumorigenesis. Cell Cycle. 2011;10(7):1092–9. https://0-doi-org.brum.beds.ac.uk/10.4161/cc.10.7.15089.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Dorman A, Baerr D, Tomlinson I, Mott R, Iraqi FA. Genetic analysis of intestinal polyp development in collaborative cross mice carrying the ApcMin/+ mutation. BMC Genet. 2016;17(1):46. https://0-doi-org.brum.beds.ac.uk/10.1186/s12863-016-0349-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Gu S, Chen D, Zhang JN, Lv X, Wang K, Duan LP, et al. Bacterial community mapping of the mouse gastrointestinal tract. PLoS One. 2013;8(10):e47957.

    Google Scholar 

  41. Lawrence NK, William FD. APC and its modifiers in colon cancer. APC Proteins. Adv Exp Med Biol. 2009;656:85–106.

    Article  Google Scholar 

  42. Van Der Kraak L, Meunier C, Turbide C, et al. A two-locus system controls susceptibility to colitis-associated colon cancer in mice. Oncotarget. 2010;1(6):436–46. https://0-doi-org.brum.beds.ac.uk/10.18632/oncotarget.101009.

    Article  Google Scholar 

  43. Ferguson B, Ram R, Handoko HY, Mukhopadhyay P, Muller HK, Soyer HP, et al. Melanoma susceptibility as a complex trait: genetic variation controls all stages of tumor progression. Oncogene. 2015;34(22):2879–86. https://0-doi-org.brum.beds.ac.uk/10.1038/onc.2014.227.

    Article  CAS  PubMed  Google Scholar 

  44. Ferguson B, Handoko HY, Mukhopadhyay P, Chitsazan A, Balmer L, Morah Walker GJ. Different genetic mechanisms mediate spontaneous versus UVR-induced malignant melanoma. Elife. 2019;8:E42424. https://0-doi-org.brum.beds.ac.uk/10.7554/eLife.42424.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

This work was supported by the Wellcome Trust grants 090532/Z/09/Z, 085906/Z/08/Z, 083573/Z/07/Z, and 075491/Z/04, Cancer research counsel-UK, and the Cancer biology research center (CBRC) of Tel Aviv University, Core funding by Tel-Aviv University (TAU), Israeli Science foundation (ISF) grants 429/09, 961/15 and 1085/18, the Israeli Cancer research fund (ICRF). Ilona Binenbaum’s PhD thesis is supported by a scholarship from the state scholarship Foundation in Greece (IKY) (operational program “Human Resources Development – Education and Lifelong Learning”, partnership agreement (PA) 2014–2020)

Author information

Authors and Affiliations

Authors

Contributions

This work was part of AD Master thesis. AD performed the experiments on the mice and collected the data, and drafted the manuscript. IB, TK, EP and AC were involved with QTL and bioinformatics analysis, AD and IB performed the analysis, and HJA helped drafting the manuscript. IT, RM and FAI conceived of and designed the study, helped with collection, analysis and interpretation of the data, drafted and finalized the manuscript. All authors approved the final version for publication.

Corresponding author

Correspondence to Fuad A. Iraqi.

Ethics declarations

Competing of interests

AC is founder and CEO of e-NIOS Applications SME. The other authors declare no potential conflict of interest with respect to financial or Non-financial competing interests, the authorship and/or publication of this article.

Ethics approval and consent to participate

The Institutional Animal Care approved all experiments in this study and Use Committee (IACUC) at TAU, which adheres to the Israeli guidelines, which follow the NIH/USA animal care and use protocols (approved experiment number M-08-075; M-12-024).

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplement 1.

Identified 411 unique genes related to the studied traits.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dorman, A., Binenbaum, I., Abu-Toamih Atamni, H.J. et al. Genetic mapping of novel modifiers for ApcMin induced intestinal polyps’ development using the genetic architecture power of the collaborative cross mice. BMC Genomics 22, 566 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-021-07890-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-021-07890-x

Keywords