Syntenic Cell Wall QTLs as Versatile Breeding Tools: Intraspecific Allelic Variability and Predictability of Biomass Quality Loci in Target Plant Species

Syntenic cell wall QTLs (SQTLs) can identify genetic determinants of biomass traits in understudied species based on results from model crops. However, their effective use in plant breeding requires SQTLs to display intraspecific allelic variability and to predict causative loci in other populations/species than the ones used for SQTLs identification. In this study, genome assemblies from different accessions of Arabidopsis, rapeseed, tomato, rice, Brachypodium and maize were used to evaluate the intraspecific variability of SQTLs. In parallel, a genome-wide association study (GWAS) on cell wall quality traits was performed in miscanthus to verify the colocalization between GWAS loci and miscanthus SQTLs. Finally, an analogous approach was applied on a set of switchgrass cell wall QTLs retrieved from the literature. These analyses revealed large SQTLs intraspecific genetic variability, ranging from presence-absence gene variation to SNPs/INDELs and changes in coded proteins. Cell wall genes displaying gene dosage regulation, such as PAL and CAD, displayed presence-absence variation in Brachypodium and rapeseed, while protein INDELs were detected for the Brachypodium homologs of the rice brittle culm-like 8 locus, which may likely impact cell wall quality. Furthermore, SQTLs significantly colocalized with the miscanthus and switchgrass QTLs, with relevant cell wall genes being retained in colocalizing regions. Overall, SQTLs are useful tools to screen germplasm for relevant genes and alleles to improve biomass quality and can increase the efficiency of plant breeding in understudied biomass crops.

A genome wide association study to dissect the genetic architecture of agronomic traits in Andean lupin (Lupinus mutabilis)

Establishing Lupinus mutabilis as a protein and oil crop requires improved varieties adapted to EU climates. The genetic regulation of strategic breeding traits, including plant architecture, growing cycle length and yield, is unknown. This study aimed to identify associations between 16 669 single nucleotide polymorphisms (SNPs) and 9 agronomic traits on a panel of 223 L. mutabilis accessions, grown in four environments, by applying a genome wide association study (GWAS). Seven environment-specific QTLs linked to vegetative yield, plant height, pods number and flowering time, were identified as major effect QTLs, being able to capture 6 to 20% of the phenotypic variation observed in these traits. Furthermore, two QTLs across environments were identified for flowering time on chromosome 8. The genes FAF, GAMYB and LNK, regulating major pathways involved in flowering and growth habit, as well as GA30X1, BIM1, Dr1, HDA15, HAT3, interacting with these pathways in response to hormonal and environmental cues, were prosed as candidate genes. These results are pivotal to accelerate the development of L. mutabilis varieties adapted to European cropping conditions by using marker-assisted selection (MAS), as well as to provide a framework for further functional studies on plant development and phenology in this species.

Genetic dissection of morphological variation in rosette leaves and leafy heads in cabbage (Brassica oleracea var. capitata)

Correlations between morphological traits of cabbage rosette leaves and heads were found. Genome-wide association studies of these traits identified 50 robust quantitative trait loci in multiple years. Half of these loci affect both organs. Cabbage (Brassica oleracea var. capitata) is an economically important vegetable crop cultivated worldwide. Cabbage plants go through four vegetative stages: seedling, rosette, folding and heading. Rosette leaves are the largest leaves of cabbage plants and provide most of the energy needed to produce the leafy head. To understand the relationship and the genetic basis of leaf development and leafy head formation, 308 cabbage accessions were scored for rosette leaf and head traits in three-year field trials. Significant correlations were found between morphological traits of rosette leaves and heads, namely leaf area with the head area, height and width, and leaf width with the head area and head height, when heads were harvested at a fixed number of days after sowing. Fifty robust quantitative trait loci (QTLs) for rosette leaf and head traits distributed over all nine chromosomes were identified with genome-wide association studies. All these 50 loci were identified in multiple years and generally affect multiple traits. Twenty-five of the QTL were associated with both rosette leaf and leafy head traits. We discuss thirteen candidate genes identified in these QTL that are expressed in heading leaves, with an annotation related to auxin and other phytohormones, leaf development, and leaf polarity that likely play a role in leafy head development or rosette leaf expansion.

Diversity and Agronomic Performance of Lupinus mutabilis Germplasm in European and Andean Environments

The introduction of Lupinus mutabilis (Andean lupin) in Europe will provide a new source of protein and oil for plant-based diets and biomass for bio-based products, while contributing to the improvement of marginal soils. This study evaluates for the first time the phenotypic variability of a large panel of L. mutabilis accessions both in their native environment and over two cropping conditions in Europe (winter crop in the Mediterranean region and summer crop in North-Central Europe), paving the way for the selection of accessions adapted to specific environments. The panel of 225 accessions included both germplasm pools from the Andean region and breeding lines from Europe. Notably, we reported higher grain yield in Mediterranean winter-cropping conditions (18 g/plant) than in the native region (9 g/plant). Instead, North European summer-cropping conditions appear more suitable for biomass production (up to 2 kg/plant). The phenotypic evaluation of 16 agronomical traits revealed significant variation in the panel. Principal component analyses pointed out flowering time, yield, and architecture-related traits as the main factors explaining variation between accessions. The Peruvian material stands out among the top-yielding accessions in Europe, characterized by early lines with high grain yield (e.g., LIB065, LIB072, and LIB155). Bolivian and Ecuadorian materials appear more valuable for the selection of genotypes for Andean conditions and for biomass production in Europe. We also observed that flowering time in the different environments is influenced by temperature accumulation. Within the panel, it is possible to identify both early and late genotypes, characterized by different thermal thresholds (600°C-700°C and 1,000-1,200°C GDD, respectively). Indications on top-yielding and early/late accessions, heritability of morpho-physiological traits, and their associations with grain yield are reported and remain largely environmental specific, underlining the importance of selecting useful genetic resources for specific environments. Altogether, these results suggest that the studied panel holds the genetic potential for the adaptation of L. mutabilis to Europe and provide the basis for initiating a breeding program based on exploiting the variation described herein.

Determinants of Common Mental Disorders (CMD) among adolescent girls aged 15-19 years in Indonesia: Analysis of the 2018 National Basic Health Survey Data

Common Mental Disorders (CMD) are distress conditions which manifest themselves with anxiety, somatic, and depressive symptoms. CMD are highly prevalent in Indonesia especially among adolescents. Adolescent girls have a higher risk to develop CMD than boys. This may be related to anaemia, potentially aggravated by early onset of menstruation. This study aimed to determine the association between haemoglobin concentration and other determinants of CMD among adolescent girls in Indonesia. Data of 1,052 adolescent girls aged 15-19 years old from the Indonesian Basic Health Survey 2018 were analysed. CMD was measured using the Self Reporting Questionnaire (SRQ-20). Principal Component Analysis of main determinants was applied and resulting principal components were investigated as risk factors for CMD. The prevalence of CMD among the study population was 16.5%. Anaemia and Age at Menarche (AAM) were not associated with CMD. Three principal components were significantly associated with higher CMD score: 1) higher parental education, better employment of the father, and living in an urban area (β: 0.16, 95%-CI: 0.02; 0.30); 2) higher consumption of salty foods, high-fat foods, and soft drinks (β: 0.23, 95%-CI: 0.05; 0.40); and 3) having asthma, smoking, and a higher haemoglobin concentration (β: 1.74, 95%-CI: 1.59; 1.89). The strongest clustered associates of CMD among adolescent girls in Indonesia were asthma, smoking status, and higher haemoglobin concentration, whereas anaemia and AAM were not associated. Causality of smoking and diet to CMD could not be disentangled in this cross-sectional study. Our findings imply that adolescent girls who have asthma and smoke, as well as those having parents with higher education and secured occupation, are more likely to have mental disorders.

Elucidating the Genetic Architecture of Fiber Quality in Hemp (Cannabis sativa L.) Using a Genome-Wide Association Study

Hemp (Cannabis sativa L.) is a bast-fiber crop with a great potential in the emerging bio-based economy. Yet, hemp breeding for fiber quality is restricted and that is mainly due to the limited knowledge of the genetic architecture of its fiber quality. A panel of 123 hemp accessions, with large phenotypic variability, was used to study the genetic basis of seven cell wall and bast fiber traits relevant to fiber quality. These traits showed large genetic variance components and high values of broad sense heritability in this hemp panel, as concluded from the phenotypic evaluation across three test locations with contrasting environments. The hemp panel was genotyped using restriction site associated DNA sequencing (RAD-seq). Subsequently, a large set (> 600,000) of selected genome-wide single nucleotide polymorphism (SNP) markers was used for a genome-wide association study (GWAS) approach to get insights into quantitative trait loci (QTLs) controlling fiber quality traits. In absence of a complete hemp genome sequence, identification of QTLs was based on the following characteristics: (i) association level to traits, (ii) fraction of explained trait variance, (iii) collinearity between QTLs, and (iv) detection across different environments. Using this approach, 16 QTLs were identified across locations for different fiber quality traits, including contents of glucose, glucuronic acid, mannose, xylose, lignin, and bast fiber content. Among them, six were found across the three environments. The genetic markers composing the QTLs that are common across locations are valuable tools to develop novel genotypes of hemp with improved fiber quality. Underneath the QTLs, 12 candidate genes were identified which are likely to be involved in the biosynthesis and modification of monosaccharides, polysaccharides, and lignin. These candidate genes were suggested to play an important role in determining fiber quality in hemp. This study provides new insights into the genetic architecture of fiber traits, identifies QTLs and candidate genes that form the basis for molecular breeding for high fiber quality hemp cultivars.

Genetic Variability of Morphological, Flowering, and Biomass Quality Traits in Hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is a bast-fiber crop well-known for the great potential to produce sustainable fibers. Nevertheless, hemp fiber quality is a complex trait, and little is known about the phenotypic variability and heritability of fiber quality traits in hemp. The aim of this study is to gain insights into the variability in fiber quality within the hemp germplasm and to estimate the genetic components, environmental components, and genotype-by-environment (G×E) interactions on fiber quality traits in hemp. To investigate these parameters, a panel of 123 hemp accessions was phenotyped for 28 traits relevant to fiber quality at three locations in Europe, corresponding to climates of northern, central, and southern Europe. In general, hemp cultivated in northern latitudes showed a larger plant vigor while earlier flowering was characteristic of plants cultivated in southern latitudes. Extensive variability between accessions was observed for all traits. Most cell wall components (contents of monosaccharides derived from cellulose and hemicellulose; and lignin content), bast fiber content, and flowering traits revealed large genetic components with low G×E interactions and high broad-sense heritability values, making these traits suitable to maximize the genetic gains of fiber quality. In contrast, contents of pectin-related monosaccharides, most agronomic traits, and several fiber traits (fineness and decortication efficiency) showed low genetic components with large G×E interactions affecting the rankings across locations. These results suggest that pectin, agronomic traits, and fiber traits are unsuitable targets in breeding programs of hemp, as their large G×E interactions might lead to unexpected phenotypes in untested locations. Furthermore, all environmental effects on the 28 traits were statistically significant, suggesting a strong adaptive behavior of fiber quality in hemp to specific environments. The high variability in fiber quality observed in the hemp panel, the broad range in heritability, and adaptability among all traits prescribe positive prospects for the development of new hemp cultivars of excellent fiber quality.

Genetic Architecture of Flowering Time and Sex Determination in Hemp (Cannabis sativa L.): A Genome-Wide Association Study

Flowering time and sex determination in hemp (Cannabis sativa L.) strongly influence fiber quality and seed production of this crop. The control of these traits is paramount for the breeding of new cultivars. Yet, little is known about the genetics underlying such complex traits and a better understanding requires in depth knowledge of the molecular mechanisms responsible for these traits. In this report, the genetic architecture of flowering time and sex determination in hemp was studied using a Genome-Wide Association Studies (GWAS) approach. Association studies were performed on a panel of 123 hemp accessions, tested in three contrasting environments, using a set of 600 K SNP markers. Altogether, eight QTLs were identified across environments; six for flowering time traits and two for sex determination. These QTLs covered genomic regions with 33 transcripts predicted to be involved in flowering and sex determination as well as a microRNA, miR156. Genes related to perception and transduction of light and transcription factors well-known to regulate flowering were identified in QTLs for flowering time traits. Transcription factors and genes involved in regulating the balance of phytohormones, specially auxins and gibberellic acid, were identified in QTLs for sex determination. Sex determination QTLs were associated with the development of male flowers in female plants and thus with the stability of sex determination in monecious plants. The present study elucidates relevant knowledge on the genetic mechanisms of flowering and sex determination traits in hemp, and provides new tools for hemp breeding.

Population structure and genome-wide association analysis for frost tolerance in oat using continuous SNP array signal intensity ratios

Infinium SNP data analysed as continuous intensity ratios enabled associating genotypic and phenotypic data from heterogeneous oat samples, showing that association mapping for frost tolerance is a feasible option. Oat is sensitive to freezing temperatures, which restricts the cultivation of fall-sown or winter oats to regions with milder winters. Fall-sown oats have a longer growth cycle, mature earlier, and have a higher productivity than spring-sown oats, therefore improving frost tolerance is an important goal in oat breeding. Our aim was to test the effectiveness of a Genome-Wide Association Study (GWAS) for mapping QTLs related to frost tolerance, using an approach that tolerates continuously distributed signals from SNPs in bulked samples from heterogeneous accessions. A collection of 138 European oat accessions, including landraces, old and modern varieties from 27 countries was genotyped using the Infinium 6K SNP array. The SNP data were analyzed as continuous intensity ratios, rather than converting them into discrete values by genotype calling. PCA and Ward's clustering of genetic similarities revealed the presence of two main groups of accessions, which roughly corresponded to Continental Europe and Mediterranean/Atlantic Europe, although a total of eight subgroups can be distinguished. The accessions were phenotyped for frost tolerance under controlled conditions by measuring fluorescence quantum yield of photosystem II after a freezing stress. GWAS were performed by a linear mixed model approach, comparing different corrections for population structure. All models detected three robust QTLs, two of which co-mapped with QTLs identified earlier in bi-parental mapping populations. The approach used in the present work shows that SNP array data of heterogeneous hexaploid oat samples can be successfully used to determine genetic similarities and to map associations to quantitative phenotypic traits.

Statistical epistasis between candidate gene alleles for complex tuber traits in an association mapping population of tetraploid potato

Association mapping using DNA-based markers is a novel tool in plant genetics for the analysis of complex traits. Potato tuber yield, starch content, starch yield and chip color are complex traits of agronomic relevance, for which carbohydrate metabolism plays an important role. At the functional level, the genes and biochemical pathways involved in carbohydrate metabolism are among the best studied in plants. Quantitative traits such as tuber starch and sugar content are therefore models for association genetics in potato based on candidate genes. In an association mapping experiment conducted with a population of 243 tetraploid potato varieties and breeding clones, we previously identified associations between individual candidate gene alleles and tuber starch content, starch yield and chip quality. In the present paper, we tested 190 DNA markers at 36 loci scored in the same association mapping population for pairwise statistical epistatic interactions. Fifty marker pairs were associated mainly with tuber starch content and/or starch yield, at a cut-off value of q ≤ 0.20 for the experiment-wide false discovery rate (FDR). Thirteen marker pairs had an FDR of q ≤ 0.10. Alleles at loci encoding ribulose-bisphosphate carboxylase/oxygenase activase (Rca), sucrose phosphate synthase (Sps) and vacuolar invertase (Pain1) were most frequently involved in statistical epistatic interactions. The largest effect on tuber starch content and starch yield was observed for the paired alleles Pain1-8c and Rca-1a, explaining 9 and 10% of the total variance, respectively. The combination of these two alleles increased the means of tuber starch content and starch yield. Biological models to explain the observed statistical epistatic interactions are discussed.

Natural DNA variation at candidate loci is associated with potato chip color, tuber starch content, yield and starch yield

Complex characters of plants such as starch and sugar content of seeds, fruits, tubers and roots are controlled by multiple genetic and environmental factors. Understanding their molecular basis will facilitate diagnosis and combination of superior alleles in crop improvement programs ("precision breeding"). Association genetics based on candidate genes is one approach toward this goal. Tetraploid potato varieties and breeding clones related by descent were evaluated for 2 years for chip quality before and after cold storage, tuber starch content, yield and starch yield. Chip quality is inversely correlated with tuber sugar content. A total of 36 loci on 11 potato chromosomes were evaluated for natural DNA variation in 243 individuals. These loci included microsatellites and genes coding for enzymes that function in carbohydrate metabolism or transport (candidate loci). The markers were used to analyze population structure and were tested for association with the tuber quality traits. Highly significant and robust associations of markers with 1-4 traits were identified. Most frequent were associations with chip quality and tuber starch content. Alleles increasing tuber starch content improved chip quality and vice versa. With two exceptions, the most significant and robust associations (q < 0.01) were observed with DNA variants in genes encoding enzymes that function in starch and sugar metabolism or transport. Comparing linkage and linkage disequilibrium between loci provided evidence for the existence of large haplotype blocks in the breeding materials analyzed.

Association mapping of leaf traits, flowering time, and phytate content in Brassica rapa

Association mapping was used to investigate the genetic basis of variation within Brassica rapa, which is an important vegetable and oil crop. We analyzed the variation of phytate and phosphate levels in seeds and leaves and additional developmental and morphological traits in a set of diverse B. rapa accessions and tested association of these traits with AFLP markers. The analysis of population structure revealed four subgroups in the population. Trait values differed between these subgroups, thus defining associations between population structure and trait values, even for traits such as phytate and phosphate levels. Marker-trait associations were investigated both with and without taking population structure into account. One hundred and seventy markers were found to be associated with the observed traits without correction for population structure. Association analysis with correction for population structure led to the identification of 27 markers, 6 of which had known map positions; 3 of these were confirmed in additional QTL mapping studies.

Analysis of multivariate extreme intakes of food chemicals

A recently published multivariate Extreme Value Theory (EVT) model is applied to the estimation of population risks associated with dietary intake of pesticides. The objective is to quantify the acute risk of pesticide intake above a threshold and relate it to the consumption of specific primary food products. As an example daily intakes of a pesticide from three foods are considered. The method models and extrapolates simultaneous intakes of pesticide, and estimates probability of exceeding unobserved large intakes. Multivariate analysis was helpful in identifying whether the avoidance of certain food combinations would reduce the likelihood of exceeding a threshold. We argue that the presented method can be an important contribution to exposure assessment studies.