The karyotype of Aegilops geniculata and its use to identify both addition and substitution lines of wheat
Molecular Cytogenetics volume 12, Article number: 15 (2019)
The annual allotetraploid species Aegilops geniculata harbors a number of traits relevant for wheat improvement. An effective cytogenetic method has yet to be developed to distinguish between each of its 14 chromosomes.
A fluorescence in situ hybridization (FISH) approach was adopted to describe the karyotype of Ae. geniculata. Each of its 14 chromosomes was unequivocally recognized using a cocktail of three probes, namely pTa-713, (AAC)5 and pTa71. FISH karyotyping was then used to detect and characterize selections from an Ae. geniculata × bread wheat wide cross of a chromosome 1Mg disomic addition line and three 4Mg(4B) substitution lines. The identity of the addition line was confirmed by the presence of Glu-M1, detected both using an SDS-PAGE separation of endosperm proteins and by applying a PCR assay directed at the Glu-M1 locus. The status of the substitution lines was validated by genotyping using a wheat single nucleotide polymorphism chip.
FISH karyotyping based on pTa-713, (AAC)5 and pTa71 will be useful for determining the contribution of Ae. geniculata to derivatives of an Ae. geniculata × wheat wide cross. SNP chip-based genotyping is effective for confirming the status of whole chromosome wheat/alien substitution lines.
Aegilops geniculata Roth (syn. Ae. ovata L, 2n = 4x = 28, genome formula UgMg) represents a source of potentially useful genetic variation of relevance to bread wheat improvement. The species is thought to represent a natural allotetraploid between the diploid Ae. umbellulata Zhuk. (U genome carrier) and Ae. comosa Sibth. et Sm. (M genome carrier) . A number of genes conferring resistance to various diseases have been transferred from this species into bread wheat, notably Yr40 (resistance against stripe rust), Lr57 (leaf rust) , Sr53 (stem rust)  and Pm29 (powdery mildew) [4, 5]. Some accessions of Ae. geniculata have displayed high levels of water use efficiency , and the species overall exhibit a higher tolerance to moisture stress than any of the related species Ae. markgrafii (Greuter) K. Hammer, Ae. longissima Schweinf. & Muschl., Ae. searsii Feldman et Kislev ex K. Hammer or Ae. speltoides Tausch . It also carries alleles at the genes encoding endosperm proteins which have been predicted to improve the end-use quality of bread wheat [8, 9].
Genomic in situ hybridization (GISH) has been a very successful technique for discriminating between the chromosomes belonging to the various genomes represented in Triticeae species, while fluorescence in situ hybridization (FISH) tends to be used for identifying individual chromosomes. A FISH-based karyotype of Ae. geniculata has been established, employing the probe combination pSc119.2, Afa family repeats, pAs1 and pTa71 [10,11,12]. However, some segments of the Ae. geniculata genome lack any probe hybridization sites, meaning that FISH karyotyping needs to be supported for the identification of non-intact chromosomes by a GISH-based analysis. For example, the sites of pSc119.2 hybridization are concentrated close to the telomeres of most chromosome arms, making it difficult to differentiate between chromosomes 1Ug, 2Ug, 3Mg, and 4Mg [10, 12]. Here, the objective was to develop a FISH assay able to unequivocally recognize each of the 14 Ae. geniculata chromosomes, and to use this assay to characterize a number of derivatives of an Ae. geniculata × wheat wide cross.
The following taxa were used in these experiments: Ae. umbellulata (2n = 2x = 14, carrier of the U genome) accession AS4, Ae. comosa (2n = 2x = 14, carrier of the M genome) accession PI551068, Ae. geniculata accession AS6, the bread wheat cultivars Yi-yuan 2 (YY2), Chinese Spring and Chuan-mai 41 (CM41), and ten F7 derivatives of a wide cross between Ae. geniculata and wheat (pedigree AS6/YY2//YY2/3/CM41).
GISH and FISH analysis
Cytological preparations were carried out using the methods described by Zhao et al. . For GISH analyses, total genomic DNA was extracted from fresh leaves of Ae. comosa and Ae. geniculata, and labeled with digoxigenin-11-dUTP (Roche Diagnostics GmbH, Mannheim, Germany) via nick translation to use as the probe; non-labeled total genomic DNA of Ae. umbellulata and CM41 were used for blocking. The GISH procedure was based on the protocol described by Hao et al. , with the exception that the concentration of the probe DNA was changed to 0.1 μg/μL and that of blocking DNA to 3.5 μg/μL. FISH experiments were conducted based on the methods given by Zhao et al. . The following probes were essayed: Afa family repeats , pSc119.2, pTa-535 and pTa71 , (AAC)5 , (CTT)5 and pTa-713 . Probes were labeled with either FAM or TAMRA by the TsingKe Biological Technology Company (Chengdu, China). The preparations were stained with DAPI (Vector Laboratories Inc., Burlingame, CA, USA) and the fluorescence signals visualized and captured using an BX-63 epifluorescence microscope equipped with a Photometric SenSys DP70 CCD camera (Olympus, Tokyo, Japan). Raw images were processed using Photoshop v.7.1 (Adobe Systems Inc., San Jose, CA, USA).
Single nucleotide polymorphism (SNP) genotyping
Genomic DNA was extracted from fresh leaves using a plant genomic DNA kit (Tiangen Biotech, Beijing Co. Ltd., Beijing, China). Chip-based genotyping was carried out using the CapitalBio Wheat 55 K SNP array. (www.capitalbio.com); the SNP loci arrayed on this chip represent a sub-set of the Affymetrix® Axiom® Wheat 660 chip, as selected by the Institute of Crop Science, Chinese Academy of Agricultural Sciences (wheat.pw.usda.gov/ggpages/topics/Wheat660_SNP_array_developed_by_CAAS.pdf). The flanking sequences of each locus were used to map each site onto the bread wheat reference sequence (urgi.versailles.inra.fr/download/iwgsc/IWGSC_RefSeq_Assemblies/v1.0/), by imposing a BLASTN E-value threshold of 10− 10, allowing a maximum mismatch of one base. For lines concluded to harbor an Ae. geniculata chromosome substituting for a wheat chromosome, the ratio between the observed and expected number of markers on the wheat chromosome in question (4B) was calculated by considering a series of 10 Mb intervals along the chromosome, applying a sliding window of 10 Mb and a step length of 1 Mb. A graphical representation of these ratios was obtained using the R package ggplot2 v.2.2.1 
High molecular weight glutenin subunit (HMW-GS) analysis
The extraction of protein from single grains and their separation using sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) followed methods given by Yan et al. . After electrophoresis for 150 min (120 V, 20 mA), the gels were stained with Coomassie Brilliant Blue R-250 staining solution for 1 h, then destained in distilled water. The primer pair 5′-ATGGCTAAGCGGYTRGTCCTCTTTG/5′-CTATCACTGGCTRGCGGACAATGG was used to amplify the coding region of the whole set of Glu-1 genes . The methods required for PCR amplification, cloning and sequencing followed those given by Guo et al. .
Development of a PCR assay for Glu-M1x
In a search for informative SNPs, a multiple sequence alignment was carried out of the gene sequences encoding the HMW-GS 1Dx2 (Genbank accession KF466259.1), 1Ax (JQ007589.1), 1Bx13 (EF540764.1), 1Ux (KX375406.1) and 1 Mx (KX375404.1), using DNAMAN v7 software . This permitted the design of a primer pair (5′-CGCCCTCGTGGCTCTCACCC/5′-TTTGCTGCTGGTATTGTCCA) which specifically amplified the encoding sequence of 1 Mx subunit. The amplicon was generated by exposing the reaction to an initial denaturation of 94 °C/5 min, followed by 30 cycles of 94 °C/30 s, 63 °C/30 s, 72 °C/40 s; the PCR was completed with a final extension step of 72 °C/10 min.
The FISH karyotype of ae. Geniculata
When the mitotic chromosomes of Ae. geniculata were probed with labeled total genomic DNA of Ae. comosa and an excess of unlabeled Ae. umbellulata genomic DNA, the M genome chromosomes (labeled green) were readily distinguished from the U genome ones (Fig. 1a). The same preparations were then subjected to FISH using the four probes pSc119.2, Afa family repeats, pTa71 and (CTT)5 following suggestions made elsewhere in the literature [10,11,12, 23, 24]. Most of the pSc119.2 sites were found to lie near the telomeres, while those recognized by the Afa family repeats probe were limited to chromosomes 2Mg, 3Mg and 7Mg (Fig. 1b). The FISH profile of AS6 differed somewhat from previously published profiles of the species: for instance, pTa71 sites were located on seven pairs of chromosomes, of which two were on U genome chromosomes (1Ug, 5Ug) and the other five on the M genome chromosomes 1Mg, 2Mg, 3Mg, 5Mg and 6Mg (Fig. 1b). The (CTT)5 probe detected sites on the chromosomes belonging to both the U and the M genomes (Fig. 1c). A FISH karyotype of AS6 based on probing with pSc119.2, Afa family repeats, pTa71 and (CTT)5 is shown in Fig. 1d. Without the GISH data, it was not easy to differentiate between several of the M and U genome chromosomes: for example, the FISH profiles of chromosomes 3Ug and 4Mg were almost identical to one another (Fig. 1d). For this reason, experiments were carried out to elaborate a more effective set of FISH probes; these led to the choice of the combination pTa-713, (AAC)5 and pTa71. The (ACC)5 sites were largely concentrated around the centromeres and the middle of the chromosome arms, while the pTa-713 sites were distributed across several chromosome arms (Fig. 1e, f). The three probes combination allowed for each of the 14 chromosomes of Ae. geniculata to be discriminated without the need for an accompanying GISH procedure (Fig. 1f).
FISH-based identification of introgression materials
Ten derivatives of a cross between Ae. geniculata and bread wheat were subjected to the newly developed FISH assay. One line had a somatic number of 44, while that of the nine others was 42 chromosomes. A GISH analysis demonstrated that two of the chromosomes present in the 2n = 44 line (Add L-1) had been inherited from Ae. geniculata (Fig. 2a). Applying the FISH procedure confirmed that these chromosomes comprised a pair of 1Mg chromosomes (Fig. 2b, c). Three of the 2n = 42 lines also carried a pair of Ae. geniculata chromosomes (Fig. 2d), and the FISH assay showed that all three (Sub L-1, Sub L-2, and Sub L-3) represented a 4Mg(4B) substitution line (Fig. 2e, f). There was no evidence for the presence of any Ae. geniculata chromatin in any of the other six lines.
The HMW-GS profile of add L-1
Since homeologous group 1 chromosomes harbor the Glu-1 genes which encode HMW-GS, the endosperm protein profile of Add L-1 grain was obtained by SDS-PAGE to identify the presence of the products of the Ae. geniculata homeolog of Glu-1. The Add L-1 and AS6 profiles both included a subunit not represented in either of the wheat parents of Add L-1 or in any of the other sister lines (Fig. 3). When the Glu-1 coding region was PCR-amplified from Ae. geniculata gDNA, the amplicon was found to include a fragment of the same length as that present in the amplicon generated from an Add L-1 template; the fragment was not represented in amplicons produced from euploid wheat (Fig. 4a). When this fragment was cloned and sequenced, both the Add L-1 and AS6 version proved to be a sequence of length 1860 nt and were of identical sequence. This sequence has been deposited in Genbank under accession number MK135469. The sequence differed from that encoding the 1My subunit (KX375405.1) with respect to seven nucleotides, and their predicted products differed for just one residue. The amplified fragment derived from the gene encoding the AS6 1 Mx subunit was likely too similar to one of the wheat fragments (Fig. 4a) to facilitate its isolation, so a PCR assay was designed to target this gene in order to confirm its presence in Add L-1 (Fig. 4b). The conclusion drawn from these assays was that the Ae. geniculata chromosome present in Add L-1 harbored the genes encoding both the x and y subunits of Glu-M1.
SNP genotyping of the 4Mg(4B) substitution lines
The three 4Mg(4B) substitution lines Sub L-1 through L-3, along with their Ae. geniculata parent AS6 and their wheat parents YY2 and CM41 were subjected to SNP genotyping to confirm their FISH-based designation. It was expected that the wheat SNP markers map to 4B (2601) showed a highest ratio to present as missing in the 4Mg(4B) lines. Thus, the missing markers ratio for individual chromosomes were analyzed. As expected, about 60% SNPs (59.0% for Sub L-1; 60.0% for Sub L-2; 60.6% for Sub L-3) for chromosome 4B that covered the whole chromosome (Fig. 5) showed as missing in the three substitution lines and this ratio greatly exceeded for other chromosomes (Table 1).
Of the 51,159 features represented on the SNP chip, 29,537 (57.7%) also hybridized with a sequence(s) present in AS6. Comparability, about 40% 4B SNPs also present in the three 4Mg(4B) substitution lines. It suggested a proportion of the SNP assays also recognized a site on 4Mg. To obtain these SNPs, we filtered 4B SNPs among Sub L-1, Sub L-2, Sub L-3 and AS6. However, ambiguous hybridizing signals may appear using wheat SNP markers to genotype its relative species and introgression lines . Thus, only homozygous SNPs were considered. In all, 240 SNPs covering the whole 4B chromosome presenting in each of the three Sub L lines and AS6, and presenting as same homozygous alleles possibly shared by 4Mg and 4B were obtained (Additional file 1: Table S1).
The ability to identify the extent, location and origin of chromatin introgressed into a crop species genome from a species in its tertiary genepool is important for the successful execution of a chromosome engineering experiment. FISH-based karyotyping based on a panel of multi-copy probes can be an effective means of recognizing individual chromosomes [11, 17, 26, 27]. Here, when a set of seven such probes was applied to derive a FISH-based karyotype of Ae. geniculata, a combination of just three of them, namely pTa-713, (AAC)5 and pTa71, was sufficient to discriminate clearly between each of the species’ 14 chromosomes. The FISH assay could then be applied to detect the origin of Ae. geniculata chromatin introgressed into a number of derivatives of a wide cross between this species and wheat, resulting in the recognition of a chromosome 1Mg disomic addition line and three 4Mg(4B) disomic substitution lines.
High density SNP arrays have been shown to be effective for the genotyping of wheat and certain wheat/alien introgression lines [25, 28]. In species containing genomes closely related to those present in wheat, the expectation is that – depending on the extent of the evolutionary separation involved – a proportion of the SNP assays developed within wheat itself will also detect the presence of alien chromatin. Here, around 60% of the wheat SNP loci recognized a site in the AS6 genome, a consequence of the fairly close relationship between Ae. geniculata and bread wheat [29,30,31]. A set of 240 of the 2601 chromosome 4B SNPs behaved in this manner, highlighting loci on both chromosome 4B and chromosome 4Mg (Additional file 1: Table S1). The conclusion was that a wheat SNP array can be highly informative for identifying the presence of a wheat/alien chromosome substitution line.
The derivatives of the Ae. geniculate × wheat wide cross have inherited their non-nuclear genomes from Ae. geniculata, which has resulted in their exhibiting very late maturity . While this trait has complicated the evaluation of their agronomic performance, it is of some interest that the presence of chromosome 4Mg induces the formation of supernumerary florets, and also introduces a gametocidal mechanism which generates chromosome breakage in gametes which lack the chromosome [33, 34]. The product of Glu-M1 has been associated with a positive effect on the rheological strength of doughs made from flour of grains carrying chromosome 1Mg [8, 9]. Thus, both the 1Mg addition line and the 4Mg(4B) substitution lines can provide the starting material for producing alien translocation lines of interest to the wheat breeding community.
This study has established a FISH protocol able to unequivocally identify each of the 14 chromosomes of Ae. geniculata, avoiding the need to include a parallel GISH procedure. A wheat SNP array was successfully deployed to confirm the cytogenetic status of three independent 4Mg(4B) substitution lines. Both the 1Mg addition line and the 4Mg(4B) substitution lines represent materials of potential utility for wheat improvement.
Kimber G, Sallee PJ, Feiner MM. The interspecific and evolutionary relationships of Triticum ovatum. Genome. 1988;30(2):218–21.
Kuraparthy V, Chhuneja P, Dhaliwal HS, Kaur S, Bowden RL, Gill BS. Characterization and mapping of cryptic introgression from Ae. geniculata with new leaf rust and stripe rust resistance genes Lr57 and Yr40 in wheat. Theor Appl Genet. 2007;114(8):1379–89.
Liu WX, Rouse M, Friebe B, Jin Y, Gill B, Pumphrey MO. Discovery and molecular mapping of a new gene conferring resistance to stem rust, Sr53, derived from Aegilops geniculata, and characterization of spontaneous translocation stocks with reduced alien chromatin. Chromosom Res. 2011;19(5):669–82.
Friebe B, Heun M. C-banding pattern and powdery mildew resistance of Triticum ovatum and four T. aestivum-T. ovatum chromosome addition lines. Theor Appl Genet. 1989;78(3):417–24.
Zeller FJ, Konig L, Hartl L, Mohler V, Hsam SLK. Chromosomal location of genes for resistance to powdery mildew in common wheat (Triticum aestivum L. em Thell.) 7. Gene Pm29 in line Pova. Euphytica. 2002;123(2):187–94.
Zaharieva M, Gaulin E, Havaux M, Acevedo E, Monneveux P. Drought and heat responses in the wild wheat relative Aegilops geniculata Roth: potential interest for wheat improvement. Crop Sci. 2001;41:1321–9.
Pradhan GP, Prasad PVV, Fritz AK, Kirkham MB, Gill BS. Response of Aegilops species to drought stress during reproductive stages of development. Funct Plant Biol. 2011;39(1):51–9.
Grag M, Tsujimoto H, Gupta RK, Kumar A, Kaur N, Kumar R, et al. Chromosome specific substitution lines of Aegilops geniculata alter parameters of bread making quality of wheat. PLoS One. 2016;11(10):e0162350.
Alvarez JB, Guzmán C. Interspecific and intergeneric hybridization as a source of variation for wheat grain quality improvement. Theor Appl Genet. 2018;131(2):225–51.
Molnár I, Cifuentes M, Schneider A, Benavente E, Molnár-Láng M. Association between simple sequence repeat-rich chromosome regions and intergenomic translocation breakpoints in natural populations of allopolyploid wild wheats. Ann Bot. 2011;107(1):65–76.
Kwiatek M, Wiśniewska H, Apolinarska B. Cytogenetic analysis of Aegilops chromosomes, potentially usable in triticale (× Triticosecale Witt.) breeding. J Appl Genet. 2013;54:147–55.
Molnár I, Kubaláková M, Šimková H, Cseh A, Molnár-Láng M, Doležel J. Chromosome isolation by flow sorting in Aegilops umbellulata and Ae. comosa and their allotetraploid hybrids Ae. biuncialis and Ae. geniculata. PLoS One. 2011;6(11):e27708.
Zhao LB, Ning SZ, Yu JJ, Hao M, Zhang LQ, Yuan ZW, et al. Cytological identification of an Aegilops variabilis chromosome carrying stripe rust resistance in wheat. Breeding Sci. 2016;66:522–9.
Hao M, Luo JT, Yang M, Zhang LQ, Yan ZH, Yuan ZW, et al. Comparison of homoeologous chromosome pairing between hybrids of wheat genotypes Chinese spring ph1b and Kaixian-luohanmai with rye. Genome. 2011;54(12):959–64.
Zhao LB, Ning SZ, Yi YJ, Zhang LQ, Yuan ZW, Wang JR, et al. Fluorescence in situ hybridization karyotyping reveals the presence of two distinct genomes in the taxon Aegilops tauschii. BMC Genomics. 2018;19(1):3.
Nagaki K, Tsujimoto H, Isono K, Sasakuma T. Molecular characterization of a tandem repeat, Afa family, and its distribution among Triticeae. Genome. 1995;38(3):479–86.
Tang ZX, Yang ZJ, Fu SL. Oligonucleotides replacing the roles of repetitive sequences pAs1, pSc119.2, pTa-535, pTa71, CCS1, and pAWRC.1 for FISH analysis. J Appl Genet. 2014;55(3):313–8.
Cuadrado A, Schwarzacher T, Jouve N. Identification of different chromatin classes in wheat using in situ hybridization with simple sequence repeat oligonucleotides. Theor Appl Genet. 2000;101(5–6):711–7.
Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer-Verlag; 2016.
Yan ZH, Wan YF, Liu KF, Zheng YL, Wang DW. Identification of a novel HMW glutenin subunit and comparison of its amino acid sequence with those of homologous subunits. Chin Sci Bull. 2002;47(3):222–6.
Guo XH, Wu BH, Hu XG, Bi ZG, Wang ZZ, Liu DC, et al. Molecular characterization of two y-type high molecular weight glutenin subunit alleles 1Ay12* and 1Ay8* from cultivated einkorn wheat (Triticum monococcum ssp. monococcum). Gene. 2013;516(1):1–7.
Wei L, Wang Q, Zhang LL, Ma J, Wang JR, Qi PF, et al. Genetic analyses of Glu - 1S sh in wheat/Aegilops sharonensis hybrid progenies and development of alien HMW-GSs gene-specific markers. Mol Breeding. 2015;35(12):230.
Badaeva ED, Amosova AV, Samatadze TE, Zoshchuk SA, Shostak NG, Chikida NN, et al. Genome differentiation in Aegilops. 4. Evolution of the U-genome cluster. Plant Syst Evol. 2004;246(1/2):45–76.
Landjeva S, Kocheva K, Karceva T, Sepsi A, Molnár I, Schneider A, et al. Molecular cytogenetic identification of a wheat - Aegilops geniculata Roth spontaneous chromosome substitution and its effects on the growth and physiological responses of seedlings to osmotic stress. Plant Breed. 2012;131(1):81–7.
Zhou SH, Zhang JP, Che YH, Liu WH, Lu YQ, Yang XM, et al. Construction of Agropyron Gaertn. Genetic linkage maps using a wheat 660K SNP array reveals a homoeologous relationship with the wheat genome. Plant Biotechnol J. 2017;16:818.
Pedersen C, Langridge P. Identification of the entire chromosome complement of bread wheat by two-colour FISH. Genome. 1997;40(5):589–93.
Komuro S, Endo R, Shikata K, Kato A. Genomic and chromosomal distribution patterns of various repeated DNA sequences in wheat revealed by a fluorescence in situ hybridization procedure. Genome. 2013;56(3):131–7.
Winfield MO, Allen AM, Burridge AJ, Baker GL, Benlow HR, Wikinson PA, et al. High-density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool. Plant Biotechnol J. 2016;14(5):1195–206.
Tiwari VK, Wang SC, Sehgl S, Vrána J, Friebe B, Kubaláková M, et al. SNP discovery for mapping alien introgressions in wheat. BMC Genomics. 2014;15(1):273.
Tiwari VK, Wang SC, Danilova T, Koo D, Vrána j KM, et al. Exploring the tertiary gene pool of bread wheat: sequence assembly and analysis of chromosome 5Mg of Aegilops geniculata. Plant J. 2016;84(4):733–46.
Molnár I, Vrána J, Burešová V, Cápal P, Farkas A, Darkó É, et al. Dissecting the U, M, S and C genomes of wild relatives of bread wheat (Aegilops spp.) into chromosomes and exploring their synteny with wheat. Plant J. 2016;88(3):452–67.
Wu YW, Zhang CL, Liu CG, Ren SX, Zhang Y. Breeding technology of alloplasmic wheat. Sci China Ser C Life Sci. 1998;41(5):449–58.
Friebe BR, Tuleen NA, Gill BS. Development and identification of a complete set of Triticum aestivum-Aegilops geniculata, chromosome addition lines. Genome. 1999;42(3):374–80.
Kwiatek MT, Wiśniewska H, Ślusarkiewicz-Jarzina A, Majka J, Majka M, Belter J, et al. Gametocidal factor transferred from Aegilops geniculata Roth can be adapted for large-scale chromosome manipulations in cereals. Front Plant Sci. 2017;8:409.
We thank the anonymous reviewers of the manuscript for their useful comments.
This research was financially supported by the Chinese government’s National Key Research and Development Program (2016YFD0102000), the Sichuan Provincial Agricultural Department’s Innovative Research Team (wheat-10), and the Sichuan Province’s Science & Technology Department Crops Breeding Project (2016NYZ0030).
Availability of data and materials
The raw SNP data and germplasm are available upon request to Shunzong Ning, Triticeae Research Institute, Sichuan Agricultural University, (email address:email@example.com .).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yi, Y., Zheng, K., Ning, S. et al. The karyotype of Aegilops geniculata and its use to identify both addition and substitution lines of wheat. Mol Cytogenet 12, 15 (2019). https://doi.org/10.1186/s13039-019-0428-2