Genetic diversity and population structure of Musa accessions in ex situ conservation

Abstract Background Banana cultivars are mostly derived from hybridization between wild diploid subspecies of Musa acuminata (A genome) and M. balbisiana (B genome), and they exhibit various levels of ploidy and genomic constitution. The Embrapa ex situ Musa collection contains over 220 accessions, of which only a few have been genetically characterized. Knowledge regarding the genetic relationships and diversity between modern cultivars and wild relatives would assist in conservation and breeding strategies. Our objectives were to determine the genomic constitution based on Internal Transcribed Spacer (ITS) regions polymorphism and the ploidy of all accessions by flow cytometry and to investigate the population structure of the collection using Simple Sequence Repeat (SSR) loci as co-dominant markers based on Structure software, not previously performed in Musa. Results From the 221 accessions analyzed by flow cytometry, the correct ploidy was confirmed or established for 212 (95.9%), whereas digestion of the ITS region confirmed the genomic constitution of 209 (94.6%). Neighbor-joining clustering analysis derived from SSR binary data allowed the detection of two major groups, essentially distinguished by the presence or absence of the B genome, while subgroups were formed according to the genomic composition and commercial classification. The co-dominant nature of SSR was explored to analyze the structure of the population based on a Bayesian approach, detecting 21 subpopulations. Most of the subpopulations were in agreement with the clustering analysis. Conclusions The data generated by flow cytometry, ITS and SSR supported the hypothesis about the occurrence of homeologue recombination between A and B genomes, leading to discrepancies in the number of sets or portions from each parental genome. These phenomenons have been largely disregarded in the evolution of banana, as the “single-step domestication” hypothesis had long predominated. These findings will have an impact in future breeding approaches. Structure analysis enabled the efficient detection of ancestry of recently developed tetraploid hybrids by breeding programs, and for some triploids. However, for the main commercial subgroups, Structure appeared to be less efficient to detect the ancestry in diploid groups, possibly due to sampling restrictions. The possibility of inferring the membership among accessions to correct the effects of genetic structure opens possibilities for its use in marker-assisted selection by association mapping.

A note on the use of the generalized odds ratio in meta-analysis of association studies involving bi- and tri-allelic polymorphisms

Abstract Background The generalized odds ratio (GOR) was recently suggested as a genetic model-free measure for association studies. However, its properties were not extensively investigated. We used Monte Carlo simulations to investigate type-I error rates, power and bias in both effect size and between-study variance estimates of meta-analyses using the GOR as a summary effect, and compared these results to those obtained by usual approaches of model specification. We further applied the GOR in a real meta-analysis of three genome-wide association studies in Alzheimer's disease. Findings For bi-allelic polymorphisms, the GOR performs virtually identical to a standard multiplicative model of analysis (e.g. per-allele odds ratio) for variants acting multiplicatively, but augments slightly the power to detect variants with a dominant mode of action, while reducing the probability to detect recessive variants. Although there were differences among the GOR and usual approaches in terms of bias and type-I error rates, both simulation- and real data-based results provided little indication that these differences will be substantial in practice for meta-analyses involving bi-allelic polymorphisms. However, the use of the GOR may be slightly more powerful for the synthesis of data from tri-allelic variants, particularly when susceptibility alleles are less common in the populations (≤10%). This gain in power may depend on knowledge of the direction of the effects. Conclusions For the synthesis of data from bi-allelic variants, the GOR may be regarded as a multiplicative-like model of analysis. The use of the GOR may be slightly more powerful in the tri-allelic case, particularly when susceptibility alleles are less common in the populations.

Functional markers for gene mapping and genetic diversity studies in sugarcane

Abstract Background The database of sugarcane expressed sequence tags (EST) offers a great opportunity for developing molecular markers that are directly associated with important agronomic traits. The development of new EST-SSR markers represents an important tool for genetic analysis. In sugarcane breeding programs, functional markers can be used to accelerate the process and select important agronomic traits, especially in the mapping of quantitative traits loci (QTL) and plant resistant pathogens or qualitative resistance loci (QRL). The aim of this work was to develop new simple sequence repeat (SSR) markers in sugarcane using the sugarcane expressed sequence tag (SUCEST database). Findings A total of 365 EST-SSR molecular markers with trinucleotide motifs were developed and evaluated in a collection of 18 genotypes of sugarcane (15 varieties and 3 species). In total, 287 of the EST-SSRs markers amplified fragments of the expected size and were polymorphic in the analyzed sugarcane varieties. The number of alleles ranged from 2-18, with an average of 6 alleles per locus, while polymorphism information content values ranged from 0.21-0.92, with an average of 0.69. The discrimination power was high for the majority of the EST-SSRs, with an average value of 0.80. Among the markers characterized in this study some have particular interest, those that are related to bacterial defense responses, generation of precursor metabolites and energy and those involved in carbohydrate metabolic process. Conclusions These EST-SSR markers presented in this work can be efficiently used for genetic mapping studies of segregating sugarcane populations. The high Polymorphism Information Content (PIC) and Discriminant Power (DP) presented facilitate the QTL identification and marker-assisted selection due the association with functional regions of the genome became an important tool for the sugarcane breeding program.

Genetic and environmental influences on blood pressure and physical activity: a study of nuclear families from Muzambinho, Brazil

Blood pressure (BP) and physical activity (PA) levels are inversely associated. Since genetic factors account for the observed variation in each of these traits, it is possible that part of their association may be related to common genetic and/or environmental influences. Thus, this study was designed to estimate the genetic and environmental correlations of BP and PA phenotypes in nuclear families from Muzambinho, Brazil. Families including 236 offspring (6 to 24 years) and their 82 fathers and 122 mothers (24 to 65 years) were evaluated. BP was measured, and total PA (TPA) was assessed by an interview (commuting, occupational, leisure time, and school time PA). Quantitative genetic modeling was used to estimate maximal heritability (h²), and genetic and environmental correlations. Heritability was significant for all phenotypes (systolic BP: h² = 0.37 ± 0.10, P < 0.05; diastolic BP: h² = 0.39 ± 0.09, P < 0.05; TPA: h² = 0.24 ± 0.09, P < 0.05). Significant genetic (r g) and environmental (r e) correlations were detected between systolic and diastolic BP (r g = 0.67 ± 0.12 and r e = 0.48 ± 0.08, P < 0.05). Genetic correlations between BP and TPA were not significant, while a tendency to an environmental cross-trait correlation was found between diastolic BP and TPA (r e = -0.18 ± 0.09, P = 0.057). In conclusion, BP and PA are under genetic influences. Systolic and diastolic BP share common genes and environmental influences. Diastolic BP and TPA are probably under similar environmental influences.

Association between celiac disease and Crohn's disease - a challenge to the coloproctologist

Over the past few years, many studies on the association between celiac disease and inflammatory bowel disease have been reported. The genetic origin of this association has prompted research that searches for a common link for the concomitant manifestation of these pathologies. Clinical studies aim not only to demonstrate this relation, but also to establish the epidemiological frequencies among affected individuals and their relatives as compared to the general population. The similar clinical symptoms, difficulties, diagnoses, and therapeutics are still a challenge, since this association is unknown to most coloproctologists, thereby culminating in treatments and surgical procedures with no benefits for the patient.

Association between IL8 haplotypes and pathogen levels in chronic periodontitis

Chronic periodontitis (CP) is considered to be a multifactorial disease influenced by microbial and genetic factors. The aim of the present study was to investigate whether the genetic susceptibility to CP in individuals with the IL8 ATC/TTC haplotype is associated with subgingival levels of periodontopathogens. Sixty-five individuals, grouped according to the presence (n=28) or absence (n=37) of the IL8 haplotype, were evaluated. After clinical periodontal evaluation, each group was subdivided according to the presence (CP) or absence (H) of periodontitis. Four subgingival samples were obtained from CP and two samples per subject from H patients. The levels and proportions of Porphyromonas gingivalis, Tannerella forsythia, and Treponema denticola were analyzed using quantitative real-time polymerase chain reaction (q-PCR). No differences were found in the proportion of periodontopathogenic bacteria between groups with the presence or absence of the IL8 haplotype. However, in the CP groups, the levels of periodontopathogens were significantly higher in the individuals withou the IL8 haplotype than in the individuals with the IL8 haplotype. These results suggest that periodontal destruction may occur in patients who are considered to be genetically susceptible to CP with a lower microbial challenge because of the presence of the IL8 ATC/TTC haplotype than in patients without this haplotype.

Analysis of genetic lineages and their correlation with virulence genes in enterococcus faecalis clinical isolates from root canal and systemic infections

Introduction: Enterococcus faecalis is a member of the mammalian gastrointestinal microbiota but has been considered a leading cause of hospital-acquired infections. In the oral cavity, it is commonly detected from root canals of teeth with failed endodontic treatment. However, little is known about the virulence and genetic relatedness among E. faecalis isolates from different clinical sources. This study compared the presence of enterococcal virulence factors among root canal strains and clinical isolates from hospitalized patients to identify virulent clusters of E. faecalis. Methods: Multilocus sequence typing analysis was used to determine genetic lineages of 40 E. faecalis clinical isolates from different sources. Virulence clusters were determined by evaluating capsule (cps) locus polymorphisms, pathogenicity island gene content, and antibiotic resistance genes by polymerase chain reaction. Results: The clinical isolates from hospitalized patients formed a phylogenetically separate group and were mostly grouped in the clonal complex 2, which is a known virulent cluster of E. faecalis that has caused infection outbreaks globally. The clonal complex 2 group comprised capsule-producing strains harboring multiple antibiotic resistance and pathogenicity island genes. On the other hand, the endodontic isolates were more diverse and harbored few virulence and antibiotic resistance genes. In particular, although more closely related to isolates from hospitalized patients, capsuleproducing E. faecalis strains from root canals did not carry more virulence/antibiotic genes than other endodontic isolates. Conclusions: E. faecalis isolates from endodontic infections have a genetic and virulence profile different from pathogenic clusters of hospitalized patients’ isolates, which is most likely due to niche specialization conferred mainly by variable regions in the genome.

The humankind genome: from genetic diversity to the origin of human diseases

Genome-wide association studies have failed to establish common variant risk for the majority of common human diseases. The underlying reasons for this failure are explained by recent studies of resequencing and comparison of over 1200 human genomes and 10 000 exomes, together with the delineation of DNA methylation patterns (epigenome) and full characterization of coding and noncoding RNAs (transcriptome) being transcribed. These studies have provided the most comprehensive catalogues of functional elements and genetic variants that are now available for global integrative analysis and experimental validation in prospective cohort studies. With these datasets, researchers will have unparalleled opportunities for the alignment, mining, and testing of hypotheses for the roles of specific genetic variants, including copy number variations, single nucleotide polymorphisms, and indels as the cause of specific phenotypes and diseases. Through the use of next-generation sequencing technologies for genotyping and standardized ontological annotation to systematically analyze the effects of genomic variation on humans and model organism phenotypes, we will be able to find candidate genes and new clues for disease’s etiology and treatment. This article describes essential concepts in genetics and genomic technologies as well as the emerging computational framework to comprehensively search websites and platforms available for the analysis and interpretation of genomic data.

Exome analysis of HIV patients submitted to dendritic cells therapeutic vaccine reveals an association of CNOT1 gene with response to the treatment

INTRODUCTION: With the aim of searching genetic factors associated with the response to an immune treatment based on autologous monocyte-derived dendritic cells pulsed with autologous inactivated HIV, we performed exome analysis by screening more than 240,000 putative functional exonic variants in 18 HIV-positive Brazilian patients that underwent the immune treatment. METHODS: Exome analysis has been performed using the ILLUMINA Infinium HumanExome BeadChip. zCall algorithm allowed us to recall rare variants. Quality control and SNP-centred analysis were done with GenABEL R package. An in-house implementation of the Wang method permitted gene-centred analysis. RESULTS: CCR4-NOT transcription complex, subunit 1 (CNOT1) gene (16q21), showed the strongest association with the modification of the response to the therapeutic vaccine (p=0.00075). CNOT1 SNP rs7188697 A/G was significantly associated with DC treatment response. The presence of a G allele indicated poor response to the therapeutic vaccine (p=0.0031; OR=33.00; CI=1.74-624.66), and the SNP behaved in a dominant model (A/A vs. A/G+G/G p=0.0009; OR=107.66; 95% CI=3.85-3013.31), being the A/G genotype present only in weak/transient responders, conferring susceptibility to poor response to the immune treatment. DISCUSSION: CNOT1 is known to be involved in the control of mRNA deadenylation and mRNA decay. Moreover, CNOT1 has been recently described as being involved in the regulation of inflammatory processes mediated by tristetraprolin (TTP). The TTP-CCR4-NOT complex (CNOT1 in the CCR4-NOT complex is the binding site for TTP) has been reported as interfering with HIV replication, through post-transcriptional control. Therefore, we can hypothesize that genetic variation occurring in the CNOT1 gene could impair the TTP-CCR4-NOT complex, thus interfering with HIV replication and/or host immune response. CONCLUSIONS: Being aware that our findings are exclusive to the 18 patients studied with a need for replication, and that the genetic variant of CNOT1 gene, localized at intron 3, has no known functional effect, we propose a novel potential candidate locus for the modulation of the response to the immune treatment, and open a discussion on the necessity to consider the host genome as another potential variant to be evaluated when designing an immune therapy study

Genetic and environmental factors associated with the risk of cognitive decline and dementia

AD is the most common age related neurodegenerative disease in the industrialized world. Clinically AD is defined as a progressing decline of cognitive functions. Neuropathologically, AD is characterized by the aggregation of b-amyloid (Ab) peptide in the form of extracellular senile plaques, and hyperphosphorlylated tau protein in the form of intracellular neurofibrillary tangles. These neuropathological hallmarks are often accompanied by abundant microvascular damage and pronounced inflammation of the affected brain regions. In this thesis we investigated several aspects of AD focusing on the genetic aspect. We confirmed that Alpha 1 antichymotrypsin (ACT), an acute phase protein, was associated to AD subjects, being plasma levels higher in AD cases than controls. In addition, in a GWA study we demonstrated that two different gene, Clusterin and CR1 were strongly associated to AD. A single gene association not explain such a complex disease like AD. The goal should be to created a network of genetic, phenotypic and clinical data associated to AD. We used a new algorithm, the ANNs, aimed to map variables and search for connectivity among variables. We found specific variables associated to AD like cholesterol levels, the presence of variation in HMGCR enzyme and the age. Other factors such as the BMI, the amount of HDL and blood folate levels were also associated with AD. Pathogen infections, above all viral infections, have been previously associated to AD. The hypothesis suggests that virus and in particular herpes virus could enter the brain when an individual becomes older, perhaps because of a decline in the immune system. Our new hypothesis is that the presence of SNPs in our GWA gene study results in a genetic signature that might affect individual brain susceptibility to infection by herpes virus family during aging.

Influence of genetic variants and post transcriptional factors on imatinib transporters as determinants of the pharmacological response

Introduction – Although imatinib (IM) is a recognized gold standard in chronic myeloid leukemia (CML) therapy, resistance has emerged in a significant proportion of patients. Aim – The aim of this study was: (1) to investigate the role of genetic variants in genes encoding for IM transporters, as candidate of IM responsiveness and (2) to test the influence of miRNAs on IM response, focusing on efflux transporters. Methods – As a first step, a panel of polymorphisms (SNPs) was genotyped in a subgroup population of 189 patients enrolled in the Tyrosine Kinase Inhibitor Optimization and Selectivity (TOPS) trial. The association with cytogenetic response and molecular response (MR) was assessed for each SNP. As a second step, an in vitro IM-resistant model (K-562 CML cell line) was established. miRNAs profiles were analyzed using Taqman arrays and in silico search was performed for miRNAs deregulated after IM treatment. mRNA and protein expression were quantified using TaqMan realtime PCR and Western blotting, respectively. Results – (1) Among Caucasian patients, ABCB1 rs60023214 significantly correlated with complete MR (P = 0.005). Concerning SNPs combination in IM uptake transporters, the associations with treatment outcomes were statistically significant for both major and complete MR (P = 0.005 and P = 0.01, respectively). (2) ABCB1 protein was not expressed under any conditions of treatment, differently from ABCG2. Two deregulated miRNAs, namely miR-212 and miR-328, were identified to be inversely correlated with ABCG2 (r2= 0.57; p=0.03 and r2=0.47; p=0.06, respectively). Experiments of loss and gain of function confirmed the functional influence of these miRNAs on ABCG2. Conclusion – The multiple candidate gene approach identified single and combination of SNPs that can be proposed as predictor of IM response. The in vitro study suggested that IM resistance could be mediated by miRNA-dependent mechanism. Further studies are needed to validate these preliminary findings.

Study of Candidate Genes for production and carcass traits in Italian Heavy Pigs: identification of a new DNA Markers, Expression Studies and Association Analysis using Different Experimental Design

Heavy pig breeding in Italy is mainly oriented for the production of high quality processed products. Of particular importance is the dry cured ham production, which is strictly regulated and requires specific carcass characteristics correlated with green leg characteristics. Furthermore, as pigs are slaughtered at about 160 kg live weight, the Italian pig breeding sector faces severe problems of production efficiency that are related to all biological aspects linked to growth, feed conversion, fat deposition and so on. It is well known that production and carcass traits are in part genetically determined. Therefore, as a first step to understand genetic basis of traits that could have a direct or indirect impact on dry cured ham production, a candidate gene approach can be used to identify DNA markers associated with parameters of economic importance. In this thesis, we investigated three candidate genes for carcass and production traits (TRIB3, PCSK1, MUC4) in pig breeds used for dry cured ham production, using different experimental approaches in order to find molecular markers associated with these parameters.

Comparative DNA sequence and methylation analyses of orthologous genes in humans and non-human primates: Genetic and epigenetic footprints of evolution

Welche genetische Unterschiede machen uns verschieden von unseren nächsten Verwandten, den Schimpansen, und andererseits so ähnlich zu den Schimpansen? Was wir untersuchen und auch verstehen wollen, ist die komplexe Beziehung zwischen den multiplen genetischen und epigenetischen Unterschieden, deren Interaktion mit diversen Umwelt- und Kulturfaktoren in den beobachteten phänotypischen Unterschieden resultieren. Um aufzuklären, ob chromosomale Rearrangements zur Divergenz zwischen Mensch und Schimpanse beigetragen haben und welche selektiven Kräfte ihre Evolution geprägt haben, habe ich die kodierenden Sequenzen von 2 Mb umfassenden, die perizentrischen Inversionsbruchpunkte flankierenden Regionen auf den Chromosomen 1, 4, 5, 9, 12, 17 und 18 untersucht. Als Kontrolle dienten dabei 4 Mb umfassende kollineare Regionen auf den rearrangierten Chromosomen, welche mindestens 10 Mb von den Bruchpunktregionen entfernt lagen. Dabei konnte ich in den Bruchpunkten flankierenden Regionen im Vergleich zu den Kontrollregionen keine höhere Proteinevolutionsrate feststellen. Meine Ergebnisse unterstützen nicht die chromosomale Speziationshypothese für Mensch und Schimpanse, da der Anteil der positiv selektierten Gene (5,1% in den Bruchpunkten flankierenden Regionen und 7% in den Kontrollregionen) in beiden Regionen ähnlich war. Durch den Vergleich der Anzahl der positiv und negativ selektierten Gene per Chromosom konnte ich feststellen, dass Chromosom 9 die meisten und Chromosom 5 die wenigsten positiv selektierten Gene in den Bruchpunkt flankierenden Regionen und Kontrollregionen enthalten. Die Anzahl der negativ selektierten Gene (68) war dabei viel höher als die Anzahl der positiv selektierten Gene (17). Eine bioinformatische Analyse von publizierten Microarray-Expressionsdaten (Affymetrix Chip U95 und U133v2) ergab 31 Gene, die zwischen Mensch und Schimpanse differentiell exprimiert sind. Durch Untersuchung des dN/dS-Verhältnisses dieser 31 Gene konnte ich 7 Gene als negativ selektiert und nur 1 Gen als positiv selektiert identifizieren. Dieser Befund steht im Einklang mit dem Konzept, dass Genexpressionslevel unter stabilisierender Selektion evolvieren. Die meisten positiv selektierten Gene spielen überdies eine Rolle bei der Fortpflanzung. Viele dieser Speziesunterschiede resultieren eher aus Änderungen in der Genregulation als aus strukturellen Änderungen der Genprodukte. Man nimmt an, dass die meisten Unterschiede in der Genregulation sich auf transkriptioneller Ebene manifestieren. Im Rahmen dieser Arbeit wurden die Unterschiede in der DNA-Methylierung zwischen Mensch und Schimpanse untersucht. Dazu wurden die Methylierungsmuster der Promotor-CpG-Inseln von 12 Genen im Cortex von Menschen und Schimpansen mittels klassischer Bisulfit-Sequenzierung und Bisulfit-Pyrosequenzierung analysiert. Die Kandidatengene wurden wegen ihrer differentiellen Expressionsmuster zwischen Mensch und Schimpanse sowie wegen Ihrer Assoziation mit menschlichen Krankheiten oder dem genomischen Imprinting ausgewählt. Mit Ausnahme einiger individueller Positionen zeigte die Mehrzahl der analysierten Gene keine hohe intra- oder interspezifische Variation der DNA-Methylierung zwischen den beiden Spezies. Nur bei einem Gen, CCRK, waren deutliche intraspezifische und interspezifische Unterschiede im Grad der DNA-Methylierung festzustellen. Die differentiell methylierten CpG-Positionen lagen innerhalb eines repetitiven Alu-Sg1-Elements. Die Untersuchung des CCRK-Gens liefert eine umfassende Analyse der intra- und interspezifischen Variabilität der DNA-Methylierung einer Alu-Insertion in eine regulatorische Region. Die beobachteten Speziesunterschiede deuten darauf hin, dass die Methylierungsmuster des CCRK-Gens wahrscheinlich in Adaption an spezifische Anforderungen zur Feinabstimmung der CCRK-Regulation unter positiver Selektion evolvieren. Der Promotor des CCRK-Gens ist anfällig für epigenetische Modifikationen durch DNA-Methylierung, welche zu komplexen Transkriptionsmustern führen können. Durch ihre genomische Mobilität, ihren hohen CpG-Anteil und ihren Einfluss auf die Genexpression sind Alu-Insertionen exzellente Kandidaten für die Förderung von Veränderungen während der Entwicklungsregulation von Primatengenen. Der Vergleich der intra- und interspezifischen Methylierung von spezifischen Alu-Insertionen in anderen Genen und Geweben stellt eine erfolgversprechende Strategie dar.

Molecular strategies for genetic diversity analysis and development of markers linked to resistance traits in apple

The goal of many plant scientists’ research is to explain natural phenotypic variation in term of simple changes in DNA sequence. DNA-based molecular markers are extensively used for the construction of genome-wide molecular maps and to perform genetic analysis for simple and complex traits. The PhD thesis was divided into two main research lines according to the different approaches adopted. The first research line is to analyze the genetic diversity in an Italian apple germplasm collection for the identification of markers tightly linked to targeted genes by an association genetic method. This made it possible to identify synomym and homonym accessions and triploids. The fruit red skin color trait has been used to test the reliability of the genetic approaches in this species. The second line is related to the development of molecular markers closely linked to the Rvi13 and Rvi5 scab resistance genes, previously mapped on apple’s chromosome 10 and 17 respectively by using the traditional linkage mapping method. Both region have been fine-mapped with various type of markers that could be used for marker-assisted selection in future breeding programs and to isolate the two resistance genes.

Genetic basis of variation for root traits and response to heat stress in durum wheat

Durum wheat is the second most important wheat species worldwide and the most important crop in several Mediterranean countries including Italy. Durum wheat is primarily grown under rainfed conditions where episodes of drought and heat stress are major factors limiting grain yield. The research presented in this thesis aimed at the identification of traits and genes that underlie root system architecture (RSA) and tolerance to heat stress in durum wheat, in order to eventually contribute to the genetic improvement of this species. In the first two experiments we aimed at the identification of QTLs for root trait architecture at the seedling level by studying a bi-parental population of 176 recombinant inbred lines (from the cross Meridiano x Claudio) and a collection of 183 durum elite accessions. Forty-eight novel QTLs for RSA traits were identified in each of the two experiments, by means of linkage- and association mapping-based QTL analysis, respectively. Important QTLs controlling the angle of root growth in the seedling were identified. In a third experiment, we investigated the phenotypic variation of root anatomical traits by means of microscope-based analysis of root cross sections in 10 elite durum cultivars. The results showed the presence of sizeable genetic variation in aerenchyma-related traits, prompting for additional studies aimed at mapping the QTLs governing such variation and to test the role of aerenchyma in the adaptive response to abiotic stresses. In the fourth experiment, an association mapping experiment for cell membrane stability at the seedling stage (as a proxy trait for heat tolerance) was carried out by means of association mapping. A total of 34 QTLs (including five major ones), were detected. Our study provides information on QTLs for root architecture and heat tolerance which could potentially be considered in durum wheat breeding programs.