BUSCADOR RECOLECTA

Figure S1: GenBank codes for the partial COI sequences of the Aedes albopictus lines used in this study; Figure S2: Mean female fecundity (left) and mean egg fertility (right) in Ae. albopictus ARwPL, ARwPBA, ARwPBN, and BN., Peer reviewed

Additional file 1: Supplementary Fig. 1. A) Box plot showing SLC3A2 gene expression levels in paired normal colon tissue and tumoral CRC tissue, using TNMplot online tool. FC = fold change. P value obtained using Mann Whitney U test. B) Box plot showing SLC3A2 gene expression levels in normal colon tissue and tumoral CRC tissue samples. Data were obtained from the GEPIA2. P value obtained using one-way ANOVA test. C) Box plot showing SLC3A2 gene expression levels in normal colon tissue and tumoral CRC tissue normal, using TNMplot online tool. P value obtained using Mann Whitney U test. D) Box plot showing the CD98hc score analyzed by immunohistochemistry in primary tumor, lymph node and liver metastases. P values were obtained using Student t test (two-sided). E) Immunohistochemical staining of CD98hc in representative samples from the study shown in D, assessed using with anti-CD98hcV509 antibody. Magnification: 40X. F) Box plot showing SLC3A2 gene expression levels in tumor and metastatic tissue of CRC patients. T = tumor, N = normal and M = metastatic. P value obtained using Kruskal–Wallis test. G) Heatmap (top) and boxplot (bottom) representative of the expression of CD98hc in different normal and tumoral tissues, obtained from the TNMplot database. Tissues written in red represent significant differences by the Mann–Whitney test. H) Expression of CD98hc in different normal and tumoral tissues. Data obtained from the GENT2 database., Instituto de Salud Carlos III Consejo Superior de Investigaciones Científicas Junta de Castilla y León ALMOM ACMUMA UCCTA CRIS Cancer Foundation Agencia Estatal de Investigación Consejo Superior de Investigaciones Cientificas (CSIC), Peer reviewed

Additional file 2: Supplementary Fig. 2. A) Internalization of the anti-CD98hcECTO antibody in SW480 cells, analyzed by immunofluorescence. Scale bar = 25 μm. The cells were seeded on coverslips and treated with 10 nM of anti-CD98hcECTO for the times indicated. The images at the right correspond to magnifications of a cell present in the images obtained at 24 h. The colocalization of CD98hc and LAPM1 is show in the merged images. B) Preparation of the antibody–drug conjugate targeting CD98hc. The coupling of DM1 to the anti-CD98hcECTO antibody was evaluated by Western blot using an anti-DM1 antibody. Twenty nanograms of anti-CD98hcECTO-DM1, the nude anti-CD98hcECTO, trastuzumab or T-DM1 were used to detect DM1 (upper panel) and the total amount of protein was evaluated by stain-free blot (lower image). Trastuzumab and T-DM1 were used as a negative and positive controls. C-E) FACS analyses in HT29 (C), cells dissociated from patient PDX BT6224 (D) and human tumoral organoid #48 (E) using anti-CD98hc-DM1 as primary antibody. The red histogram correspond to signals from cells incubated with the secondary antibody alone, whereas the blue histograms represent the fluorescence due to the expression of CD98hc., Instituto de Salud Carlos III Consejo Superior de Investigaciones Científicas Junta de Castilla y León ALMOM ACMUMA UCCTA CRIS Cancer Foundation Agencia Estatal de Investigación Consejo Superior de Investigaciones Cientificas (CSIC), Peer reviewed

Supplementary Fig. 1, Discussion and Tables 1–2., Peer reviewed

Additional file 3: Supplementary Fig. 3. A) Dose–response analyses of the effect of anti-CD98hc-DM1 on the proliferation of parental and CD98hc CRISPR #5 and #11 HT29 cells. Cells were treated with anti-CD98hc-DM1 at the indicated doses for four days. Results are shown as the mean ± SD of quadruplicates of an experiment repeated three times. B) Expression of CD98hc in normal human fibroblasts (NHF) and immortalized keratinocytes (HaCaT), compared to CRC cell lines. Cell extracts (20 µg) were used to identify CD98hc by Western blot with the anti-CD98hcV509 antibody. Calnexin was used as a loading control. C) Dose–response analyses of the anti-CD98hc-DM1 ADC on NHF and HaCaT, compared to HT29 cells. Cells were treated with the ADC for four days at the indicated doses. The data are plotted as the percentage of MTT metabolization with respect to control. Results are shown as the mean ± SD of quadruplicates of an experiment repeated two times. D) Evaluation of the effect of anti-CD98hc-DM1 (10 nM, 48 h) on the distribution of the different cell cycle phases in HT29 and HCT116 cell lines. E) Immunofluorecescence analyses of the action of anti-CD98hc-DM1 on spindle assembly and organization on HCT116 cells treated with CD98hc-DM1 (10 nM) for 48 h. β-Tubulin (green), DAPI (blue). Scale bars = 7.5 µm. F). Detection of giant multinucleated cells or altered nuclear structures after anti-CD98hc-DM1 treatment. HCT116 and HT29 cells were treated with 10 nM anti-CD98hcECTO-DM1 for 72 h, fixed and stained for nucleoporin p62 (red) and DNA (blue). Scale bar = 10 μm. The arrows indicate giant multinucleated cells., Instituto de Salud Carlos III Consejo Superior de Investigaciones Científicas Junta de Castilla y León ALMOM ACMUMA UCCTA CRIS Cancer Foundation Agencia Estatal de Investigación Consejo Superior de Investigaciones Cientificas (CSIC), Peer reviewed

Source Data Fig. 1: Raw data for resistivity measurements. Source Data Fig. 3: Data for scaled magneto-anisotropy. Source Data Fig. 4: (a) S1 to S4, data for the temperature dependence of resistivity anisotropy. (b) Theory, data for theoretically predicted anisotropy versus scaled temperature., Peer reviewed

[Description of methods used for collection/generation of data] 12870_2023_4431_MOESM1_ESM.xlsx. Orthologue/homoeologue sequences obtained from EnsemblPlants after BLASTn. 12870_2023_4431_MOESM2_ESM.xlsx. Genomic sequences of genes TraesCS7D02G094000, TraesCS4A02G397900, TRITD4Av1G231840 and TRITD4Av1G231510 were retrieved from Ensembl Plants and aligned using the multiple sequence alignment tool ClustalW. Primers were designed in the conserved 5’ and 3’ regions by using the NCBI Primer-Blast tool. SNP markers were designed following a Tetra-Primer ARMS (amplification refractory mutation system) strategy for SNPs detection. Primer1 web service (http://primer1.soton.ac.uk/primer1.html) was used for primer design. 12870_2023_4431_MOESM4_ESM.xlsx. Sequence editing, alignment and assembly were performed with SeqMan Pro Lasergene Software v17 (DNAStar, WI, US). The identity of the clones as GDSL esterase-lipase- like sequences was confirmed by BLASTn at NCBI. The coding sequences were predicted based on the exon-intron structure of TraesCS7D02G094000, and open reading frames were searched by using ORFfinder at NCBI (https://www.ncbi.nlm.nih.gov/orffinder). Signal peptide and cellular location of the expected proteins were predicted by SignalP 6.0 software 12870_2023_4431_MOESM5_ESM.xlsx. The genetic map was constructed using Joinmap(R) 5.0 [Kyazma(R), The Netherlands]. 12870_2023_4431_MOESM6_ESM.xlsx. Carotenoid content and profile was determined according to Rodríguez-Suárez C, Requena-Ramírez MD, Hornero-Méndez D, Atienza SG (2022). Chapter 4. The breeder’s tool-box for enhancing the content of esterified carotenoids in wheat: From extraction and profiling of carotenoids to marker-assisted selection of candidate genes. En: Carotenoids: Carotenoid and apocarotenoid biosynthesis, metabolic engineering and synthetic biology. pp. 99-125. Editor. Eleanore T. Wurtzel. Book series: Methods in Enzymology. Editorial: Academic Press. Elsevier. Hardcover ISBN: 9780323913539., This research was financed by project PID2021-122152NB-I00 funded by MCIN/AEI/10.13039/501100011033/ and by ERDF “ERDF A way of making Europe”. M.D.R.-R. was supported by PRE2018-084037 funded by MCIN/AEI/10.13039/501100011033 and ESF “ESF investing in your future”. DH-M is member of the Spanish Carotenoid Network (CaRed), grant RED2022-134577-T. SA, CR-S and MDR-R are members of CeReS Network, grant RED2022-134922-T. Both networks are funded by MCIN/AEI/10.13039/501100011033., 12870_2023_4431_MOESM1_ESM.xlsx 12870_2023_4431_MOESM2_ESM.xlsx 12870_2023_4431_MOESM3_ESM.xlsx 12870_2023_4431_MOESM4_ESM.xlsx 12870_2023_4431_MOESM5_ESM.xlsx 12870_2023_4431_MOESM6_ESM.xlsx 12870_2023_4431_MOESM1_ESM.xlsx contains the results of the Search for XAT candidate genes in durum wheat genomes using XAT-7D (TraesCS7D02G094000) from common wheat as gene model. 12870_2023_4431_MOESM2_ESM.xlsx contains the list of primers designed in this work. 12870_2023_4431_MOESM3_ESM.xlsx contains the list of landraces showing positive amplification for XAT candidate gene. 12870_2023_4431_MOESM4_ESM.pdf contains the alignment of XAT-7A1 and XAT-7D proteins. The predicted active sites are shown with black triangles. Signal peptide in positions 1–22 is highlighted. 12870_2023_4431_MOESM5_ESM.xlsx contains the genetic map for the F2 population derived from the cross BGE047535 × 'Athoris' using DArTSeq markers. 12870_2023_4431_MOESM6_ESM.xlsx contains the carotenoid content and profile of the F2 population derived from the cross BGE047535 x Athoris., Peer reviewed

-allR.oggu: md5:5c51eed351afe17c4c062353ae1f8c54 -CVS1_strain_glue1.mph: md5:9edd31d718dd390f47ff65d116c34deb -CVS1_strain_mormal.mph: md5:573f3bdaf751bd3c1728ad43d2c81dce -CVS_membrane.mph: md5:0f44822ed62478ac8d193df77ef33966 -Data for NP.opju: md5:ea7b0b52b4d77c405215adb211fc4faa -Raw-data.zip: md5:302f15aa75ce00fdc3598d0a2acd8036 -Springs.mph: md5:8538155fac19b2db8c4b723d633861be, Data deposite for the manuscript entitled "Correlated order at the tipping point in the kagome metal CsV3Sb5"., Peer reviewed

20 figures, 2 tables.-- TEM, size distribution and elemental analysis of CNC: TEM images of the type II CNC are presented in Figure S1 and the distribution of diameters and heights of the CNC are in Figure S2. It can be observed the characteristic shape and size of this type of CNC, with a mean diameter 28 ± 13 nm and a length of 58 ± 15 nm. Type II CNCs were analyzed to determine their composition by elemental analysis and the results are presented in Table S1. The most characteristic result is the relative high mass percentage of sulfur of 3.45%, which is clearly higher than the value in type I CNCs prepared by the same approach. The elemental composition (C, H, N and S) of the CNCs was determined using a LECO 628 elemental analyzer (Velp Scientifica). The elemental analysis was performed in triplicate to ensure reproducibility, and the average values were reported. Commercial TiO2 paste (TiO2-P): Film fabrication: The commercial paste (TiO2-P) was dried in an oven at 120 ºC overnight to remove organic solvents and get a material ready for subsequent solid stated characterization such as XRD and TGA measurements. For the preparation of the photoelectrodes using the commercial paste, an optimized screen-printing procedure was employed following the instructions provided by the supplier. The paste was applied to cover a 1 cm2 surface area of the FTO substrates. Subsequently, the electrode was thermally sintered in an oven with the following temperature profile: 5 minutes at 325 ºC, 5 minutes at 375 ºC, 5 minutes at 450 ºC, and 15 minutes at 500 ºC, under an air atmosphere, according to the instructions from the provider. Prior to the photoelectrochemical (PEC) evaluation, the TiO2 film was activated by treating it at 500 ºC in air for 30 minutes. Preparation of TiO2(NH4OH) and solid material: TiO2(NH4OH) dispersions were prepared by mixing 50 mg of anatase powder, 300 µL of commercial aqueous ammonia (30%) and 19.7 mL of ultrapure water. Then, the resulting mixture was homogenized in an ultrasonic bath for 1 h. Dispersions were freeze-dried to obtain powder materials for further characterization. UV-Vis of the employed materials: Figure S5a shows the transmittance curves of various aqueous dispersions, including freshly prepared TiO2-NPs, TiO2(NH4OH), CNC materials, and the TiO2(NH4OH) material after 24 hours of dispersion preparation. At the selected wavelength of 360 nm, the TiO2- NPs dispersion exhibits the highest transmittance, indicating an unstable system with most of the material settling down. Conversely, the TiO2(NH4OH) material shows lower transmittance values, which moderately increases after 24 hours, evidencing successful dispersion of the TiO2-NPs when ammonia was used. Figure S5b shows the variation of the transmittance for the TiO2(NH4OH), TiO2-CNC and TiO2-CNC(NH4OH) materials with time. The results demonstrate the effective role of ammonia in facilitating the dispersion of TiO2-NPs, although CNC exhibited superior stabilization efficiency. Notably, the combination of ammonia and CNC yields the most stable aqueous dispersions, as evidenced by consistent transmittance values even after 24 hours. TEM of the TiO2(NH4OH) material: TEM images in Figure S6 reveal the presence of TiO2 aggregates, typically smaller than 100 nm. It seems that the addition of NH4OH during the preparation of the TiO2-NPs dispersion effectively disrupt the formation of large aggregates. Thermogravimetric analysis of CNC: All thermogravimetric analysis (TGA) of powder materials were carried out under air atmosphere in a Libra F1 (Netzsch) thermobalance using a ramp of 10 ºC/min. Figure S7 confirms that CNC have been completely eliminated, as there is no residual mass after the experiment. Thermogravimetric analysis of the TiO2-CNC material: As the CNC residue is 0% at 800 ºC and approximately of 5% at 450 ºC, it can be deduced that CNC have been completely removed from the TiO2 matrix during the sintering step (450 ºC, 2 hours). The TiO2-CNC material reveals a 50% of residual mass, which confirms that all CNC were removed at a 1:1 TiO2 : CNC ratio (Figure S8). Differences in surface morphology between films prepared from commercial TiO2- P and TiO2-CNC(NH4OH): Figure S9 shows the results from the profilometry measurements of TiO2-P and TiO2- CNC(NH4OH) films. Both films displayed remarkable differences in surface characteristics. While TiO2-P film shows a smooth surface, the film obtained from the TiO2-CNC(NH4OH) dispersion exhibits a significant higher roughness. This disparity is clearly a consequence of the respective fabrication methodologies. In the case of the TiO2-P, screen printing was employed, allowing the particles to accommodate under the gentle pressure of the printing blade, resulting in a roughness of 110 nm (Figure S9a). By contrast, the spray coating process immobilized the TiO2-NPs upon contact with the hot substrate, leading to rapid droplet evaporation and the formation of a film with a pronounced roughness of 1300 nm (Figure S9b). Effect of CNC in the surface morphology. Comparison between films prepared from commercial TiO2-P and TiO2(NH4OH): To gain more insight into the effect of film processing, namely screen printing and spray coating, a comparison between TiO2 photoanodes of the same thickness (~3.5 µm) prepared by both techniques is herein shown. The screen-printed TiO2 photoanode was prepared with the GreatCell® paste, whereas the spray-coated TiO2 photoanode was fabricated from the TiO2(NH4OH) dispersions in order to discard the CNC effect upon sintering, as described in the main article. In terms of surface morphology, it is of great interest the study of such property according to the followed fabrication procedure. As commented before, the screen-printed films display low Rq values (Figure S10a) whereas the TiO2(NH4OH) ones show a higher roughness (Figure S10b). These differences directly arise from the film fabrication method followed. Furthermore, the addition of CNC clearly influences the morphology of the film, with the bare TiO2(NH4OH) film showing lower roughness (900 nm, Figure S10b) compared to the TiO2-CNC(NH4OH) one (1300 nm, Figure S9b). Gas physisorption of the TiO2 materials (N2 isotherms): To further explore the effect of the CNC on the macroporous structure of the TiO2-CNC(NH4OH) material used as photoanode, physisorption measurements were conducted on both bare TiO2-NPs and TiO2-CNC(NH4OH) powder materials after sintering. N2 adsorption−desorption at -196 °C (Quantachrome Autosorb-6B Instrument) was measured after sample degassing (250 °C, 4 h) to characterize the porous texture and the equivalent Brunauer−Emmett−Teller (BET) specific surface area (SBET). Figure S11 shows the N2 isotherms of the employed materials. Both materials exhibit type II isotherms according to the IUPAC classification, typical of non-porous solids.The isotherms (Figure S11) reveal an initial increase at low relative pressure values. When it comes to the intermediate region of the isotherms, it is important to note that the TiO2 CNC hybrid has a slightly increased adsorption due to the removal of the biopolymer during the thermal treatment. This probably refers to an enlarged separation between the solid TiO2 particles. A narrow hysteresis loop appears at very high relative pressures, around p/p0 = 0.9, which is commonly ascribed to the capillary condensation taking place within the interstitial pores between TiO2 particles. The BET specific surface area of the TiO2-CNC(NH4OH) powder material is marginally higher (42 m2/g) to the observed for the bare TiO2 nanoparticles (35 m2 /g), evidencingnot too much influence of the CNC in the final internal porosity. Pore size distribution (DFT and BJH methods): Pore size distribution has been calculated for both samples from their N2 adsorption isotherms using the density functional theory (DFT) (Figure S12) and the Barrett-JoynerHalenda (BJH) method (Figure S13), which uses the Kelvin model of pore filling. X-Ray diffraction: X-ray diffraction results of the employed materials, namely CNC (type-II), TiO2-NPs (anatase), TiO2-CNC, and the commercial TiO2 paste (GreatCell®) are shown in Figure S14. The diffractogram of CNC (type-II) is in agreement with literature. The synthesized TiO2-NPs exhibits the characteristic profile of anatase NPs, evidencing a comparable crystal phase composition and crystallite size to the commercial TiO2 paste. The average crystallite size of the employed TiO2-NPs is 25 nm,3 whereas the commercial TiO2 paste is composed of particles distributed in two sizes: 20 and 300 nm. Notably, the XRD analysis of the prepared TiO2-CNC(NH4OH) hybrid shows a combined pattern from TiO2 and CNC, not showing any additional peaks. Thermogravimetric analysis of the commercial TiO2 paste: The concentration of TiO2 in the commercial paste is somewhat higher (62%, Figure S15) to that of the TiO2-CNC. This mass loss is ascribed to the removal of alkylated celluloses from the paste. Scanning electron microscopy of photoanodes from commercial TiO2-P based: Figure S16 shows SEM images of the film obtained from TiO2 commercial paste. The TiO2-P based photoelectrode film displays a smooth surface, both before and after air sintering, despite having a similar content of cellulose derivatives (62 wt.% of TiO2) to our TiO2-CNC(NH4OH) hybrid (50 wt.%) (Figure S8). Photoelectrochemical characterization of the films prepared from TiO2(NH4OH) and commercial TiO2-P: Figure S17 displays the CV profiles of the TiO2(NH4OH) and TiO2-P films, both under dark and illumination. In the absence of light, the photoelectrodes exhibit the characteristic reversible redox behavior of TiO2 electrodes. A cathodic current is observed at more negative potentials, and a nearly symmetric positive current during the backward scan. The voltammograms exhibited an accumulation region at approximately -0.8 V, and a depletion region at higher potentials, -0.4 V. Under illumination conditions, both materials show similar photocurrent values (~ 46 µA·cm2). This observation is further supported by transient photocurrent measurements (Figure S18), which reveal slight differences between the two electrodes (from 32 µA·cm-2 to 43 µA·cm-2), highlighting the significant improvement of the TiO2 electrodes when using CNC and ammonia for their fabrication. 5-hour photocurrent measurements: Aiming to study the stability of the TiO2-CNC(NH4OH) and TiO2-P photoanodes, 5-hour experiments at a constant potential (0 V vs. Ag/AgCl) were performed (Figure S19). A photocurrent decay is observed in both cases, mainly caused by the blocking of the TiO2 active sites due to parasitic redox processes (see references 33, 34 and 35). Nevertheless, the TiO2-CNC(NH4OH) retained a higher PEC performance, even after 5 hours under operation conditions due to its specific morphology that arises from the CNC processing. Electrochemical impedance spectroscopy (EIS): Equivalent circuit: EIS spectra were analyzed according to the model circuit shown in Figure S20. Resistances (RS, Rct and Rsc) and constant phase elements, i.e. non-ideal capacitors with phase angle <90, (CPEdl and CPEsc) are calculated as model fitting parameters to the experimental data.-- Under a Creative Commons license BY-NC 3.0., TEM, size distribution and elemental analysis of CNC. Commercial TiO2 paste (TiO2-P): Film fabrication. Preparation of TiO2(NH4OH) and solid material. DLS and ζ-potential measurements. UV-Vis of the employed materials. TEM of the TiO2(NH4OH) material. Thermogravimetric analysis of CNC. Thermogravimetric analysis of the TiO2-CNC material. Differences in surface morphology between films prepared from commercial TiO2-P and TiO2-CNC(NH4OH). Effect of CNC in the surface morphology. Comparison between films prepared from commercial TiO2-P and TiO2(NH4OH). Physisorption of the TiO2 materials (N2 isotherms). Pore size distribution (DFT and BJH methods). X-Ray diffraction. Thermogravimetric analysis of the commercial TiO2 paste. SEM of photoanodes from commercial TiO2-P. Photoelectrochemical characterization of films prepared from TiO2(NH4OH) and commercial TiO2-P. 5-hour photocurrent measurements. Electrochemical impedance spectroscopy (EIS). Equivalent circuit. References., Financial support from Spanish MCIN/AEI/10.13039/501100011033 and “ERDF A way of making Europe” under project grants PID2022-139671OB-I00 and PID2020-120439-RA-I00, as well as by the Gobierno de Aragón (DGA) under projects T03_23R and E47_23R (Grupos de Investigación Reconocidos) is acknowledged. V. C. is thankful for his PhD contract funded by DGA (Ref. CUS/581/2020)., Peer reviewed

Additional file 4: Supplementary Fig. 4. A) Kaplan–Meier survival curve of mice from the experiment performed in Fig. 6A. The Kaplan–Meier survival plot was created using a tumor volume threshold of 1,000 mm3. P values were calculated using one-sided log-rank tests. B) Effect of the anti-CD98hc ADC on the weight of mice xenografted with HT29 cells. Data are plotted as mean ± SD of six mice/group. C) Analysis of the antitumoral effect of naked anti-CD98hc and DM1 on tumor growth in nude mice implanted with HT29 cells. Arrows indicate days of administration of anti-CD98hc (15 mg/Kg) or DM1 (0.14 mg/Kg). Data are plotted as mean tumor volumes ± SEM. P values were calculated using Student’s t test (two-sided). D) Kaplan–Meier survival curve of mice from the experiment of panel C. The Kaplan–Meier survival plot was created using a tumor volume threshold of 650 mm3. P values were calculated using one-sided log-rank tests. E) Expression levels or phosphorylation of proteins involved in cell cycle and apoptosis in the tumors of the experiment performed in Fig. 6A. Tumor samples were obtained on day 21 after initiation of treatments (seven days after the last treatment). Tissue extracts of the tumors were used to analyze the levels of expression of pH3, PARP, pH2AX, pCDK1 and cleaved Caspase 3. Stain free blot was analyzed to verify equal loading. F) Quantitation of the levels of DM1 (data shown in Fig. 6B), pH3, PARP, pH2AX, pCDK1 and cleaved Caspase 3 of the experiments shown in panel E. The graphs represent the mean intensity (arbitrary units) ± SD of the different proteins present in control (C) or treated (anti-CD98hc-DM1) mice groups. P values were calculated using Student’s t test (two-sided). G) Immunohistochemical detection of CD98hc in PDX BT6224 (P2M1: passage 2, mouse #1) using the anti-CD98hcV509 antibody. H) Kaplan–Meier survival curve of mice from the experiment performed in Fig. 6E. The Kaplan–Meier survival plot was created using a tumor volume threshold of 800 mm3. P values were calculated using one-sided log-rank tests. I) Effect of the anti-CD98hc ADC on the weight of mice xenografted with PDX BT6224. Data are plotted as mean ± SD of four mice/group. J) Quantitation of the levels of DM1 of the experiments performed in Fig. 6F. The graphs represent the mean intensity (arbitrary units) ± SD of IgG Heavy and Light chains coupled to DM1 present in control (C) or treated (anti-CD98hc-DM1) groups. P values were calculated using Student’s t test (two-sided)., Instituto de Salud Carlos III. Consejo Superior de Investigaciones Científicas (España). Junta de Castilla y León. Fundación CRIS contra el Cáncer., Peer reviewed