Validation and Comparison of Common Thoracolumbar Injury Classification Treatment Algorithms and a Novel Modification.

BACKGROUND AND OBJECTIVES: The most common thoracolumbar trauma classification systems are the Thoracolumbar Injury Classification and Severity Score (TLICS) and the Thoracolumbar AO Spine Injury Score (TL AOSIS). Predictive accuracy of treatment recommendations is a historical limitation. Our objective was to validate and compare TLICS, TL AOSIS, and a modified TLICS (mTLICS) that awards 2 points for the presence of fractured vertebral body height loss >50% and/or spinal canal stenosis >50% at the fracture site. METHODS: The medical records of adult patients with acute, traumatic thoracolumbar injuries at an urban, Level 1 trauma center were retrospectively reviewed. TLICS, mTLICS, and TL AOSIS scores were calculated for 476 patients using computed tomography, MRI, and the documented neurological examination. Treatment recommendations were compared with treatment received. Standard validity measures were calculated. RESULTS: Treatment recommendations matched actual treatments in 95.6% (455/476) of patients for mTLICS, 91.3% (435/476) for TLICS, and 92.6% (441/476) for TL AOSIS. The differences between the accuracy of mTLICS and TLICS (95.6% vs 91.3%, P < .001) and between mTLICS and TL AOSIS (95.6% vs 91.3%, P = .003) were significant. The sensitivity of mTLICS was higher than that of TLICS (96.3% vs 81.3%, P < .001), and the sensitivity of TL AOSIS was higher than that of TLICS (92.5% vs 81.3%, P < .001). The specificity of mTLICS was equal to that of TLICS (95.3%) and higher than that of TL AOSIS (95.3% vs 92.7%, P = .02). The modifier led to substantial outperformance of mTLICS over TLICS due to 38 patients (20 of whom received surgery) moving from a TLICS score of <4 to a mTLICS score equal to 4. CONCLUSION: All systems performed well. The mTLICS had improved sensitivity and accuracy compared with TLICS and higher accuracy and specificity than TL AOSIS. The sensitivity of TL AOSIS was higher than that of TLICS. Prospective, multi-institutional reliability and validity studies of this mTLICS are needed for adoption.

Novel implant design for comminuted posterior wall acetabular fractures.

BACKGROUND: In recent years, the medical community has witnessed a notable increase in high-energy traumatic injuries, leading to a surge in complex fracture patterns that challenge existing treatment methodologies. Among these, the posterior approach to acetabular fractures stands out for offering direct visualization of the retro-acetabular surface, with current fixation methods relying on 3.5 mm low-profile reconstruction plates and various other implants. Despite the effectiveness of these methods, there is a burgeoning demand for a singular, adaptable implant that not only streamlines the surgical process but also optimizes patient outcomes. METHODS: In an innovative approach to address this need, three-dimensional (3D) models of the posterior acetabular wall were meticulously crafted using AutoCAD® software. The chosen material for the implant was 316L surgical steel for its durability and strength. The design of the implant featured a low-profile mesh structure, which was instrumental in facilitating osteosynthesis. This design allowed for the placement of screws of varying lengths in multiple directions, ensuring the initial reconstruction of the joint in an anatomical position without hindering the placement of the definitive implant. The primary objective was to secure the fixation and stabilization of the fracture by specifically targeting the smaller bone fragments. A comparative analysis was then conducted between this novel plate and a conventional 316L surgical steel, seven-hole, 3.5 mm reconstruction plate through finite element analysis. RESULTS: The comparative analysis unveiled that both plates demonstrated comparable deformation capacities, with no significant differences in load-bearing capabilities observed. This finding suggests that the innovative plate can match the performance of traditional plates used in such surgeries. CONCLUSIONS: The finite element analysis revealed that the newly developed anatomical plate for posterior wall acetabular fractures meets the necessary physical and mechanical criteria for permanent implementation in patients with these fractures. This breakthrough represents a promising advancement that could simplify surgical procedures and potentially elevate patient outcomes. LEVEL OF EVIDENCE II: This study is classified as a Level II, diagnostic study.

Enhanced antimicrobial efficacy of biogenic ZnO nanoparticles through UV-B activation: A novel approach for textile garment.

Zinc oxide nanoparticles (ZnO NP) are characterized by novel properties which have been attracting the attention of different lines of research due to their wide applicability. Obtaining this nanomaterial is strongly linked to biogenic synthesis methods, which have also been developed in this research, using Coriandrum sativum extract as a reducing agent. ZnO NPs have been properly characterized by techniques to evaluate their morphology by transmission electron microscopy (TEM) and elemental analysis by EDX. The evaluation of the antimicrobial and antifungal effects is linked to the use of a system provided by "locker sanitizer" equipment, which has been designed and built as part of this research, and is intended to treat textile garments by nebulizing the ZnO NP colloid (99.08 µg/mL) + UV-B, water + UV-B, and UV-B only, and also to evaluate the influence of the treatment time for 1, 2 and 3 min. In this sense, it is known that the nanomaterial used shows a better response to UV light because more hydroxyl radicals are produced, leading to a higher reaction rate, which results in greater efficiency in inhibitory processes. The results show that the use of the locker sanitizer is more efficient when using ZnO NP + UV-B light since it achieved 100 % growth inhibition against E. coli, C. albicans, and A. brasiliensis, and >99 % against S. aureus, after 3 min of treatment.

A fast algorithm to factorize high-dimensional tensor product matrices used in genetic models.

Many genetic models (including models for epistatic effects as well as genetic-by-environment) involve covariance structures that are Hadamard products of lower rank matrices. Implementing these models requires factorizing large Hadamard product matrices. The available algorithms for factorization do not scale well for big data, making the use of some of these models not feasible with large sample sizes. Here, based on properties of Hadamard products and (related) Kronecker products, we propose an algorithm that produces an approximate decomposition that is orders of magnitude faster than the standard eigenvalue decomposition. In this article, we describe the algorithm, show how it can be used to factorize large Hadamard product matrices, present benchmarks, and illustrate the use of the method by presenting an analysis of data from the northern testing locations of the G × E project from the Genomes to Fields Initiative (n ∼ 60,000). We implemented the proposed algorithm in the open-source "tensorEVD" R package.

Is there a trend towards semi-extended knee positioning during intramedullary tibial nailing in Latin America? A survey among 990 orthopaedic surgeons.

OBJECTIVE: The aim of this study was to evaluate the Latin American orthopaedic trauma surgeons preference regarding knee positioning and entry portals for IM nailing and identify the reasons of these preferences. METHODS: Using the AO Trauma database, 22.285 surveys were distributed by email to Latin American orthopaedic surgeons. Demographic data and practice patterns, especially regarding knee positioning and approach for tibial nailing, were then evaluated and statistically treated. RESULTS: amongst one thousand five hundred fourteen responses, 990 orthopaedic surgeons (4.4% of response rate) fully responded to the survey. Transpatellar tendon approach (613 / 61.9%,) with the knee in flexion (518 / 52.3%) on a radiolucent table remains the standard practice for intramedullary tibial nailing. Even for proximal and distal tibial nailing, the transpatellar tendon approach (455 / 46%) with the knee in flexion (562 / 56.8%) hold on the most used method. Only 55 (9.36%) orthopaedic surgeons reported that they have migrated to the supra-patellar tibial nailing in the recent years. The main reasons for a low rate of migration were lack of knowledge about the technique and unavailability of specific gigs and cartilage protectors for a safe suprapatellar nailing. CONCLUSION: Even with the potential benefits of the semi-extended knee positioning for tibial nailing, the Latin American orthopaedic community remains using the transpatellar tendon approach with the knee in flexion as the standard technique. Lack of surgical training for suprapatellar and parapatellar approaches with the knee in semi-extension, added by the unavailability of suprapatellar jigs and soft outer protection sleeves contribute to counter the trend towards the semi-extended techniques.

Surgical management of iliac wing fractures: Proposal of a new classification system.

INTRODUCTION: Isolated fractures of the iliac wing are a rare injury, usually occurring in high-energy trauma, and are associated with other non-musculoskeletal and soft tissue injuries that could compromise the patient's life. Surgical indications are unclear, and there is limited information on the most frequent fracture patterns. METHODS: A descriptive multicenter case series study of isolated fractures of the iliac wing treated surgically in three referral trauma hospitals in Latin America. The different fracture patterns are described, the "iliac ring" concept is proposed, and a classification is made. RESULTS: Twenty-eight patients were included; 24 were male, the median age was 31 years RIQ (24-46), the most frequent trauma mechanism was a traffic accident, and in 14 patients, the ISS >16. The most frequent associated injury was to the appendicular skeleton at another level in 13 patients. In the new classification, according to the number of fragments, 11 patients were classified as type A (1 fragment), ten patients as type B (two fragments), and seven patients as type C (three or more fragments). The most compromised anatomical area was the crest and anterosuperior iliac spines in 26 patients, followed by the fossa and anteroinferior iliac spine in 17 and 8 patients, respectively. DISCUSSION: The patterns of isolated fractures of the iliac wing allow the identification of three types of fractures. Identifying these patterns can help the surgeon decide to perform surgery in these scenarios and choose the fixation technique according to the number and location of the fragments.

Availability and use of resources for emergency fracture care of pelvic trauma associated with haemorrhagic shock in Latin America: A cross-sectional study.

Determining the true availability of resources and understanding the level of training of surgeons involved in the treatment of patients with pelvic fractures and haemorrhagic shock is critical. In the herein study, the availability of technical, technological, and human resources for the care of this injury in Latin America region was analysed, and the preferences of orthopaedic trauma surgeons when performing interventions for the diagnosis and treatment of patients with pelvic trauma and associated haemorrhagic shock was described. A cross sectional web-based survey containing questions on knowledge, attitudes, and practices with respect to imaging resources, emergency pelvic stabilization methods, and interventions used for bleeding control was sent to 948 Latin America orthopaedic trauma surgeons treating pelvic fractures in the emergency department. Differences between regional clusters, level of training, type of hospital, and pelvic surgery volume were assessed. 368 responses were obtained, with 37.5% of respondents reporting formal training in pelvic surgery and 36.0% having available protocol for managing these patients. The most frequently used interventions were the supra-acetabular pelvic external fixator and pelvic packing. Limited hospital and imaging resources are available for the care of patients with pelvic trauma and associated haemorrhagic shock throughout Latin America. In addition, the training of orthopaedic trauma surgeons dealing with this type of injury and the volume of pelvic surgeries per year is heterogeneous. It should be urgently considered to develop management protocols adapted to Latin America according to the availability of resources, as well as to promote training in this severe life-threatening traumatic condition.

Treatment of fracture-related infection in Latin America (FRILA). Proposal for a multicentre regional registry.

Postoperative bone infection is a severe complication in the treatment of fractures. The management of this pathology is challenging, but recent advances have been made to achieve standardization that can help diagnosis and decision-making. However, we are unaware of studies validating these models in Latin America. Therefore, this study aims to collect data from patients with fracture-related infections treated in different institutions in Latin America to create a registry that will assist in future clinical decision-making regarding the diagnostic process and the surgical and medical treatment of these patients.

Leveraging data from the Genomes-to-Fields Initiative to investigate genotype-by-environment interactions in maize in North America.

Genotype-by-environment (G×E) interactions can significantly affect crop performance and stability. Investigating G×E requires extensive data sets with diverse cultivars tested over multiple locations and years. The Genomes-to-Fields (G2F) Initiative has tested maize hybrids in more than 130 year-locations in North America since 2014. Here, we curate and expand this data set by generating environmental covariates (using a crop model) for each of the trials. The resulting data set includes DNA genotypes and environmental data linked to more than 70,000 phenotypic records of grain yield and flowering traits for more than 4000 hybrids. We show how this valuable data set can serve as a benchmark in agricultural modeling and prediction, paving the way for countless G×E investigations in maize. We use multivariate analyses to characterize the data set's genetic and environmental structure, study the association of key environmental factors with traits, and provide benchmarks using genomic prediction models.

Characterization of the influence of spinopelvic parameters on thoracolumbar trauma.

OBJECTIVE: The current Roussouly classification identifies four groups of "normal" sagittal spine morphology, which has greatly expanded the understanding of normal heterogeneity of the spine. While there has been extensive characterization of the influence of spinopelvic parameters on outcomes after degenerative spine surgery, the influence of spinopelvic parameters on thoracolumbar trauma has yet to be described. The goal of this study was to determine if spinopelvic parameters and global spine morphology influence fracture location, fracture morphology, and rate of neurological deficit in the setting of thoracolumbar trauma. METHODS: Of 2896 patients reviewed in the authors' institutional spine database between January 2014 and April 2020 with an ICD-9/10 diagnosis of thoracolumbar trauma, 514 met the inclusion criteria of acute thoracolumbar fracture on CT and visible femoral heads on sagittal CT. Pelvic incidence (PI) was calculated on sagittal CT. Demographic and clinical data including age, sex, BMI, smoking status, concomitant cervical fracture, mechanism of injury, major fracture location, neurological deficit, AO Spine thoracolumbar injury classification, and management type (operative vs nonoperative) were collected. Patients were stratified into high-PI (≥ 50°) and low-PI (< 50°) groups. RESULTS: Patients with high PI had a lower incidence of fractures in the lower lumbar spine (below L2) compared with patients with low PI (16% vs 8%, p < 0.01). The last lordotic vertebrae were observed between T10 and L4, and of fractures that occurred at these levels, 75% were at the last lordotic vertebrae. Fall from height was the most common cause of neurological deficit, accounting for 47%. Of the patients presenting with a fall from height, AO Spine type B distraction injuries were more common in the high-PI group (41% vs 18%, p = 0.01). Similarly, within the same subgroup, AO Spine type A compression injuries were more common in the low-PI group (73% vs 53%, p = 0.01). CONCLUSIONS: Spinopelvic parameters and sagittal balance influence the location and morphology of thoracolumbar fractures. Fractures of the thoracolumbar junction are strongly associated with the inflection point, which is defined by sagittal alignment. While the importance of considering sagittal balance is known for decision-making in degenerative spinal pathology, further studies are required to determine if spinopelvic parameters and sagittal balance should play a role in the decision-making for management of thoracolumbar fractures.

Genomes to Fields 2022 Maize genotype by Environment Prediction Competition.

OBJECTIVES: The Genomes to Fields (G2F) 2022 Maize Genotype by Environment (GxE) Prediction Competition aimed to develop models for predicting grain yield for the 2022 Maize GxE project field trials, leveraging the datasets previously generated by this project and other publicly available data. DATA DESCRIPTION: This resource used data from the Maize GxE project within the G2F Initiative [1]. The dataset included phenotypic and genotypic data of the hybrids evaluated in 45 locations from 2014 to 2022. Also, soil, weather, environmental covariates data and metadata information for all environments (combination of year and location). Competitors also had access to ReadMe files which described all the files provided. The Maize GxE is a collaborative project and all the data generated becomes publicly available [2]. The dataset used in the 2022 Prediction Competition was curated and lightly filtered for quality and to ensure naming uniformity across years.

Simultaneous ipsilateral floating hip and knee: the double floating extremity-a systematic review and proposal of a treatment algorithm.

PURPOSE: To systematically review the currently available existing evidence related to the presentation and management of simultaneous floating hip and knee injuries to identify injury characteristics, treatment strategies, and complications. METHODS: Data sources: Relevant articles were identified by searching Medline, PubMed, and Google Scholar databases with no language restrictions. Manual searches of other relevant databases (SciELO and grey literature databases) and reference lists of primary articles found from initial searches were also conducted. STUDY SELECTION: All types of study designs published from January 1st, 2000 to October 1st, 2022 involving skeletally mature patients with simultaneous floating hip and knee injuries were included. DATA EXTRACTION: Basic information and specific injury-related information were collected. RESULTS: Eight case reports were included. No study adequately reported the case with sufficient detail to allow other investigators to make inferences, nor was the result properly calculated, nor was the follow-up considered adequate for adequate functional assessment to occur in 80% of the studies. CONCLUSION: The exact treatment strategy and the follow-up time are not uniform across the included studies; therefore, they are not sufficient to adequately recommend surgical approach, timing of fixation, and fixation method. Our findings warrant the need for better documentation and reporting information about the mode of treatment of simultaneous floating hip and knee injuries.

Mitotic DNA damage promotes chromokinesin-mediated missegregation of polar chromosomes in cancer cells.

DNA damage response (DDR) during interphase involves active signaling and repair to ensure genomic stability. However, how mitotic cells respond to DNA damage remains poorly understood. Supported by correlative live-/fixed-cell microscopy, it was found that mitotic cells exposed to several cancer chemotherapy compounds acquire and signal DNA damage, regardless of how they interact with DNA. In-depth analysis upon DNA damage during mitosis revealed a spindle assembly checkpoint (SAC)-dependent, but ataxia telangiectasia mutated-independent, mitotic delay. This delay was due to the presence of misaligned chromosomes that ultimately satisfy the SAC and missegregate, leading to micronuclei formation. Mechanistically, it is shown that mitotic DNA damage causes missegregation of polar chromosomes due to the action of arm-ejection forces by chromokinesins. Importantly, with the exception of DNA damage induced by etoposide-a topoisomerase II inhibitor-this outcome was independent of a general effect on kinetochore microtubule stability. Colony formation assays in pan-cancer cell line models revealed that mitotic DNA damage causes distinct cytotoxic effects, depending on the nature and extent of the damage. Overall, these findings unveil and raise awareness that therapeutic DNA damage regimens may contribute to genomic instability through a surprising link with chromokinesin-mediated missegregation of polar chromosomes in cancer cells.

Micronuclei from misaligned chromosomes that satisfy the spindle assembly checkpoint in cancer cells.

Chromosome alignment to the spindle equator is a hallmark of mitosis thought to promote chromosome segregation fidelity in metazoans. Yet chromosome alignment is only indirectly supervised by the spindle assembly checkpoint (SAC) as a byproduct of chromosome bi-orientation, and the consequences of defective chromosome alignment remain unclear. Here, we investigated how human cells respond to chromosome alignment defects of distinct molecular nature by following the fate of live HeLa cells after RNAi-mediated depletion of 125 proteins previously implicated in chromosome alignment. We confirmed chromosome alignment defects upon depletion of 108/125 proteins. Surprisingly, in all confirmed cases, depleted cells frequently entered anaphase after a delay with misaligned chromosomes. Using depletion of prototype proteins resulting in defective chromosome alignment, we show that misaligned chromosomes often satisfy the SAC and directly missegregate without lagging behind in anaphase. In-depth analysis of specific molecular perturbations that prevent proper kinetochore-microtubule attachments revealed that misaligned chromosomes that missegregate frequently result in micronuclei. Higher-resolution live-cell imaging indicated that, contrary to most anaphase lagging chromosomes that correct and reintegrate the main nuclei, misaligned chromosomes are a strong predictor of micronuclei formation in a cancer cell model of chromosomal instability, but not in non-transformed near-diploid cells. We provide evidence supporting that intrinsic differences in kinetochore-microtubule attachment stability on misaligned chromosomes account for this distinct outcome. Thus, misaligned chromosomes that satisfy the SAC may represent a previously overlooked mechanism driving chromosomal/genomic instability during cancer cell division, and we unveil genetic conditions predisposing for these events.

Sparse kernel models provide optimization of training set design for genomic prediction in multiyear wheat breeding data.

The success of genomic selection (GS) in breeding schemes relies on its ability to provide accurate predictions of unobserved lines at early stages. Multigeneration data provides opportunities to increase the training data size and thus, the likelihood of extracting useful information from ancestors to improve prediction accuracy. The genomic best linear unbiased predictions (GBLUPs) are performed by borrowing information through kinship relationships between individuals. Multigeneration data usually becomes heterogeneous with complex family relationship patterns that are increasingly entangled with each generation. Under these conditions, historical data may not be optimal for model training as the accuracy could be compromised. The sparse selection index (SSI) is a method for training set (TRN) optimization, in which training individuals provide predictions to some but not all predicted subjects. We added an additional trimming process to the original SSI (trimmed SSI) to remove less important training individuals for prediction. Using a large multigeneration (8 yr) wheat (Triticum aestivum L.) grain yield dataset (n = 68,836), we found increases in accuracy as more years are included in the TRN, with improvements of ∼0.05 in the GBLUP accuracy when using 5 yr of historical data relative to when using only 1 yr. The SSI method showed a small gain over the GBLUP accuracy but with an important reduction on the TRN size. These reduced TRNs were formed with a similar number of subjects from each training generation. Our results suggest that the SSI provides a more stable ranking of genotypes than the GBLUP as the TRN becomes larger.

Plantas medicinales utilizadas en la región de Teziutlán, Puebla, México / Medicinal plants used in the region of Teziutlán, Puebla, México

The objective of this work was to collect information on the curative use of plants in the municipality of Teziutlán, Puebla through semi-structured interviews. Thus, 78 plants used for medicinal purposes were identified, of which 40 are native to Mexico and 38 introduced; The value of use (UV) of each one and the Informant's Consensus Factor (FCI) of 10 categories of diseases were calculated. The five most frequently used plants are Ruta chalepensis L., Rosmarinus officinalis, Arnica montana, Loeselia mexicana (Lam.) Brandegee and Sambucus cerulea var., Neomexicana, which underwent a chemical and pharmacological review. On the other hand, the most frequent preparations are infusion and decoction, using mainly leaves (49.34%) and flowers (19.51%). It is concluded that the inhabitants of Teziutlán have a fairly homogeneous ethnomedical knowledge, setting the standard for research on its pharmacological properties.

El objetivo del presente trabajo fue recopilar información sobre el uso curativo de las plantas del municipio de Teziutlán, Puebla por medio de entrevistas semiestructuradas. Así, se identificaron 78 plantas utilizadas con fines medicinales, de las cuales 40 son propias de México y 38 introducidas; se calculó el valor de uso (UV) de cada una y el Factor de Consenso del Informante (FCI) de 10 categorías de padecimientos. Las cinco plantas utilizadas con mayor frecuencia son Ruta chalepensis L., Rosmarinus officinalis, Arnica montana, Loeselia mexicana (Lam.) Brandegee y Sambucus cerulea var., neomexicana, a las que se les hizo una revisión química y farmacológica. Por otro lado, las preparaciones más frecuentes son infusión y decocción, utilizando principalmente hojas (49.34 %) y flores (19.51 %). Se concluye que los habitantes de Teziutlán cuentan con un conocimiento etnomédico bastante homogéneo, dando la pauta para investigaciones sobre sus propiedades farmacológicas.

A Comparative Study of Post-Quantum Cryptosystems for Internet-of-Things Applications.

The existence of quantum computers and Shor's algorithm poses an imminent threat to classical public-key cryptosystems. These cryptosystems are currently used for the exchange of keys between servers and clients over the Internet. The Internet of Things (IoT) is the next step in the evolution of the Internet, and it involves the connection of millions of low-powered and resource-constrained devices to the network. Because quantum computers are becoming more capable, the creation of a new cryptographic standard that cannot be compromised by them is indispensable. There are several current proposals of quantum-resistant or post-quantum algorithms that are being considered for future standards. Given that the IoT is increasing in popularity, and given its resource-constrained nature, it is worth adapting those new standards to IoT devices. In this work, we study some post-quantum cryptosystems that could be suitable for IoT devices, adapting them to work with current cryptography and communication software, and conduct a performance measurement on them, obtaining guidelines for selecting the best for different applications in resource-constrained hardware. Our results show that many of these algorithms can be efficiently executed in current IoT hardware, providing adequate protection from the attacks that quantum computers will eventually be capable of.

Genomic selection for morphological and yield-related traits using genome-wide SNPs in oil palm.

Oil palm is the most important oil crop worldwide. Colombia is the fourth largest producer, primarily relying on production from interspecific hybrids, derived from crosses between Elaeis oleifera and Elaeis guineensis (OxG). However, conventional breeding can take up to 20 years to generate a new variety. Therefore, reducing the breeding cycle while improving the genetic gain for complex traits is desirable. Genomic selection (GS) is an approach with the potential to achieve this goal. In this study, we evaluated 431 F1 interspecific hybrids (OxG) and 444 backcrosses (BC1) for morphological and yield-related traits. Genomic predictions were performed with the G-BLUP model using three different population datasets for training the model: the same population (TRN1), the other population (TRN2), and both populations (TRN1+2). Higher multi-family prediction accuracies were obtained for foliar area (0.3 in OxG) and trunk height (0.47 in BC1) when the model was trained with TRN1. Single-family prediction accuracies were lower in the OxG compared to BC1 families for traits such as trunk diameter, trunk height, bunch number, and yield using TRN1. Conversely, lower prediction accuracies were obtained for most traits when the model was trained using TRN2 (< 0.1). Multi-trait models showed a substantial increase of the predictions for traits such as yield (0.22 for OxG and 0.44 for BC1), because of the genetic correlations between traits. The results herein highlighted the potential of GS for parental selection in OxG and BC1 populations, but further studies are required to improve the models to select individuals by their genetic value. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01341-5.