Construction of a high-density genetic map and QTL mapping of leaf traits and plant growth in an interspecific F1 population of Catalpa bungei × Catalpa duclouxii Dode

Genome-Wide Identification and Expression Analysis of the SQUAMOSA Promoter-Binding Protein-like (SPL) Transcription Factor Family in Catalpabungei

As a plant-specific transcription factor, the SPL gene family plays a critical role in plant growth and development. Although the SPL gene family has been identified in diverse plant species, there have been no genome-wide identification or systematic study reports on the SPL gene family in Catalpa bungei. In this study, we identified 19 putative SPL gene family members in the C. bungei genome. According to the phylogenetic relationship, they can be divided into eight groups, and the genes in the same group have a similar gene structure and conserved motifs. Synteny analysis showed that fragment duplication played an important role in the expansion of the CbuSPL gene family. At the same time, CbuSPL genes have cis-acting elements and functions related to light response, hormone response, growth and development, and stress response. Tissue-specific expression and developmental period-specific expression analysis showed that CbuSPL may be involved in flowering initiation and development, flowering transition, and leaf development. In addition, the ectopic expression of CbuSPL4 in Arabidopsis confirmed that it can promote early flowering and induce the expression of related flowering genes. These systematic research results will lay a foundation for further study on the functional analysis of SPL genes in C. bungei.

Genome-wide association study of multiyear dynamic growth traits in hybrid Liriodendron identifies robust genetic loci associated with growth trajectories

The genetic factors underlying growth traits differ over time points or stages. However, most current studies of phenotypes at single time points do not capture all loci or explain the genetic differences underlying growth trajectories. Hybrid Liriodendron exhibits obvious heterosis and is widely cultivated, although its complex genetic mechanism underlying growth traits remains unknown. A genome-wide association study (GWAS) is an effective method for elucidating the genetic architecture by identifying genetic loci underlying complex quantitative traits. In the present study, using a GWAS, we identified robust loci associated with growth trajectories in hybrid Liriodendron populations. We selected 233 hybrid progenies derived from 25 crosses for resequencing, and measured their tree height (H) and diameter at breast height (DBH) for 11 consecutive years; 192 972 high-quality single nucleotide polymorphisms (SNPs) were obtained. The dynamics of the multiyear single-trait GWAS showed that year-specific SNPs predominated, and only five robust SNPs for DBH were identified in at least three different years. Multitrait GWAS analysis with model parameters as latent variables also revealed 62 SNPs for H and 52 for DBH associated with the growth trajectory, displaying different biomass accumulation patterns, among which four SNPs exerted pleiotropic effects. All identified SNPs also exhibited temporal variations in effect sizes and inheritance patterns potentially related to different growth and developmental stages. The haplotypes resulting from these significant SNPs might pyramid favorable loci, benefitting the selection of superior genotypes. The present study provides insights into the genetic architecture of dynamic growth traits and lays a basis for future molecular-assisted breeding.

Association Mapping and Expression Analysis of the Genes Involved in the Wood Formation of Poplar

Xylogenesis is a complex and sequential biosynthetic process controlled by polygenes. Deciphering the genetic architecture of this complex quantitative trait could provide valuable information for increasing wood biomass and improving its properties. Here, we performed genomic resequencing of 64 24-year-old trees (64 hybrids of section Aigeiros and their parents) grown in the same field and conducted full-sib family-based association analyses of two growth and six woody traits using GEMMA as a choice of association model selection. We identified 1342 significantly associated single nucleotide polymorphisms (SNPs), 673 located in the region upstream and downstream of 565 protein-encoding genes. The transcriptional regulation network of secondary cell wall (SCW) biosynthesis was further constructed based on the published data of poplar miRNA, transcriptome, and degradome. These provided a certain scientific basis for the in-depth understanding of the mechanism of poplar timber formation and the molecular-assisted breeding in the future.

High resolution mapping of QTLs for fruit color and firmness in Amrapali/Sensation mango hybrids

Introduction

Mango (Mangifera indica L.), acclaimed as the 'king of fruits' in the tropical world, has historical, religious, and economic values. It is grown commercially in more than 100 countries, and fresh mango world trade accounts for ~3,200 million US dollars for the year 2020. Mango is widely cultivated in sub-tropical and tropical regions of the world, with India, China, and Thailand being the top three producers. Mango fruit is adored for its taste, color, flavor, and aroma. Fruit color and firmness are important fruit quality traits for consumer acceptance, but their genetics is poorly understood.

Methods

For mapping of fruit color and firmness, mango varieties Amrapali and Sensation, having contrasting fruit quality traits, were crossed for the development of a mapping population. Ninety-two bi-parental progenies obtained from this cross were used for the construction of a high-density linkage map and identification of QTLs. Genotyping was carried out using an 80K SNP chip array.

Results and discussion

Initially, we constructed two high-density linkage maps based on the segregation of female and male parents. A female map with 3,213 SNPs and male map with 1,781 SNPs were distributed on 20 linkages groups covering map lengths of 2,844.39 and 2,684.22cM, respectively. Finally, the integrated map was constructed comprised of 4,361 SNP markers distributed on 20 linkage groups, which consisted of the chromosome haploid number in Mangifera indica (n =20). The integrated genetic map covered the entire genome of Mangifera indica cv. Dashehari, with a total genetic distance of 2,982.75 cM and an average distance between markers of 0.68 cM. The length of LGs varied from 85.78 to 218.28 cM, with a mean size of 149.14 cM. Phenotyping for fruit color and firmness traits was done for two consecutive seasons. We identified important consistent QTLs for 12 out of 20 traits, with integrated genetic linkages having significant LOD scores in at least one season. Important consistent QTLs for fruit peel color are located at Chr 3 and 18, and firmness on Chr 11 and 20. The QTLs mapped in this study would be useful in the marker-assisted breeding of mango for improved efficiency.

Genetic architecture of leaf morphology revealed by integrated trait module in Catalpa bungei

Leaves are crucial for maintaining plant growth and development via photosynthesis, and their function is simultaneously regulated by a suite of phenotypic traits. Although much is known about the genetic architecture of individual leaf traits, unraveling the genetic basis of complex leaf morphology remains a challenge. Based on the functional correlation and coordination of multi-traits, we divided 15 leaf morphological traits into three modules, comprising size (area, length, width, and perimeter), shape (leaf lobes, aspect ratio, circularity, rectangularity, and the relevant ratios), and color (red, green, and blue) for an ornamental tree species, Catalpa bungei. A total of 189 significant single-nucleotide polymorphisms were identified in the leaves of C. bungei: 35, 82, and 76 in the size, shape, and color modules, respectively. Four quantitative trait loci were common between the size and shape modules, which were closely related according to phenotype correlation, genetic mapping, and mRNA analysis. The color module was independent of them. Synergistic changes in the aspect ratio, leaf lobe, and circularity suggest that these traits could be the core indicators of the leaf shape module. The LAS and SRK genes, associated with leaf lobe and circularity, were found to function in plant defense mechanisms and the growth of leaves. The associations between the SRK and CRK2 genes and the leaf lobe and circularity traits were further verified by RT-qPCR. Our findings demonstrate the importance of integrating multi-trait modules to characterize leaf morphology and facilitate a holistic understanding of the genetic architecture of intraspecific leaf morphology diversity.

Multiple dynamic models reveal the genetic architecture for growth in height of Catalpa bungei in the field

Growth in height (GH) is a critical determinant for tree survival and development in forests and can be depicted using logistic growth curves. Our understanding of the genetic mechanism underlying dynamic GH, however, is limited, particularly under field conditions. We applied two mapping models (Funmap and FVTmap) to find quantitative trait loci responsible for dynamic GH and two epistatic models (2HiGWAS and 1HiGWAS) to detect epistasis in Catalpa bungei grown in the field. We identified 13 co-located quantitative trait loci influencing the growth curve by Funmap and three heterochronic parameters (the timing of the inflection point, maximum acceleration and maximum deceleration) by FVTmap. The combined use of FVTmap and Funmap reduced the number of candidate genes by >70%. We detected 76 significant epistatic interactions, amongst which a key gene, COMT14, co-located by three models (but not 1HiGWAS) interacted with three other genes, implying that a novel network of protein interaction centered on COMT14 may control the dynamic GH of C. bungei. These findings provide new insights into the genetic mechanisms underlying the dynamic growth in tree height in natural environments and emphasize the necessity of incorporating multiple dynamic models for screening more reliable candidate genes.

A High-Density Genetic Map Enables Genome Synteny and QTL Mapping of Vegetative Growth and Leaf Traits in Gardenia

The gardenia is a traditional medicinal horticultural plant in China, but its molecular genetic research has been largely hysteretic. Here, we constructed an F₁ population with 200 true hybrid individuals. Using the genotyping-by-sequencing method, a high-density sex-average genetic map was generated that contained 4,249 SNPs with a total length of 1956.28 cM and an average genetic distance of 0.46 cM. We developed 17 SNP-based Kompetitive Allele-Specific PCR markers and found that 15 SNPs were successfully genotyped, of which 13 single-nucleotide polymorphism genotypings of 96 F₁ individuals showed genotypes consistent with GBS-mined genotypes. A genomic collinearity analysis between gardenia and the Rubiaceae species Coffea arabica, Coffea canephora and Ophiorrhiza pumila showed the relativity strong conservation of LG11 with NC_039,919.1, HG974438.1 and Bliw01000011.1, respectively. Lastly, a quantitative trait loci analysis at three phenotyping time points (2019, 2020, and 2021) yielded 18 QTLs for growth-related traits and 31 QTLs for leaf-related traits, of which qBSBN7-1, qCD8 and qLNP2-1 could be repeatably detected. Five QTL regions (qCD8 and qSBD8, qBSBN7 and qSI7, qCD4-1 and qLLLS4, qLNP10 and qSLWS10-2, qSBD10 and qLLLS10) with potential pleiotropic effects were also observed. This study provides novel insight into molecular genetic research and could be helpful for further gene cloning and marker-assisted selection for early growth and development traits in the gardenia.

Efficient Transformation of Catalpa bungei Shows Crystal Genes Conferring Resistance to the Shoot Borer Omphisa plagialis

Although Catalpa bungei is a forest plant with considerable economic and ornamental value in China, its wood and decorative qualities are constrained by insect pests such as the shoot borer Omphisa plagialis (Lepidoptera). Overexpressing insect resistance genes such as crystal genes to develop an insect-resistant variety of C. bungei is an environmental and ecological approach. However, genotype limitations and low regeneration rates of embryogenic calli (EC) inhibit the development of transformation and the insect-resistant gene expression system in C. bungei. Here, we first established embryogenic callus induction and regeneration systems of five genotypes using mature seed and stem segment explants; the highest induction and regeneration rates of EC were 39.89 and 100%, respectively. Next, an efficient and stable Agrobacterium-mediated genetic transformation system was developed from EC and its positive frequency was up to 92.31%. Finally, using the transformation system, 15 and 22 transgenic C. bungei lines that expressed Cry2A and Cry9Aa-like were generated, respectively. These transgenic lines that exhibited significantly higher resistance to O. plagialis in the laboratory and field have great promise for meeting the challenge of future pest management under changing climatic conditions. Additionally, this efficient, fast, and stable transformation system could be a potential tool for gene function analysis and forest tree genetic improvement.

A Bivariate Mapping Model Identifies Major Covariation QTLs for Biomass Allocation Between Leaf and Stem Growth of Catalpa bungei

Biomass allocation plays a critical role in plant morphological formation and phenotypic plasticity, which greatly impact plant adaptability and competitiveness. While empirical studies on plant biomass allocation have focused on molecular biology and ecology approaches, detailed insight into the genetic basis of biomass allocation between leaf and stem growth is still lacking. Herein, we constructed a bivariate mapping model to identify covariation QTLs governing carbon (C) allocation between the leaves and stem as well as the covariation of traits within and between organs in a full-sib mapping population of C. bungei. A total of 123 covQTLs were detected for 23 trait pairs, including six leaf traits (leaf length, width, area, perimeter, length/width ratio and petiole length) and five stem traits (height, diameter at breast height, wood density, stemwood volume and stemwood biomass). The candidate genes were further identified in tissue-specific gene expression data, which provided insights into the genetic architecture underlying C allocation for traits or organs. The key QTLs related to growth and biomass allocation, which included UVH1, CLPT2, GAD/SPL, COG1 and MTERF4, were characterised and verified via gene function annotation and expression profiling. The integration of a bivariate Quantitative trait locus mapping model and gene expression profiling will enable the elucidation of genetic architecture underlying biomass allocation and covariation growth, in turn providing a theoretical basis for forest molecular marker-assisted breeding with specific C allocation strategies for adaptation to heterogeneous environments.

InDel marker development and QTL analysis of agronomic traits in mung bean [Vigna radiate (L.) Wilczek]

The stem color of young mung bean is a very useful tool in germplasm identification. Flowering time and plant height (PH) are known to be strongly correlated with crop adaption and yield. However, few studies have focused on elucidating the genetic mechanisms that regulate these five particular traits: young stem color (YSC), days to first flowering (DFF), days to maturity (DM), PH, and nodes on the main stem (NMS). In this study, a genetic linkage map for the F2 population was constructed using 129 InDel markers that were developed based on the sequence variations between parents. A total of 14 QTLs related to YSC, DFF, DM, PH, and NMS were detected. These QTLs were distributed on six chromosomes (1, 3, 4, 6, 7, and 10), which individually accounted for 1.32 to 90.07% of the total phenotypic variation. Using a short and high-density linkage map for the F3 population, six of the seven QTLs which clustered at two intervals on chromosomes 3 and 10 were detected again. Further analysis found that four QTLs between InDel markers R3-15 and R3-19 controlled DFF, DM, PH, and NMS, and each QTL accounted for a large percent of the total phenotypic variation. Analysis of two separated F2:3 lines also found that the phenotype was highly corresponded to its genotype which was between R3-15 and R3-19. Phenotype and genotype analysis for 30 mung bean accessions showed that the major effect QTL qDFF3 was a key regulator for DFF. Using a map-based cloning method, the major effect QTL qYSC4 for YSC was mapped in a 347 Kb interval on chromosome 4. Candidate gene analysis showed that sequence variations and expression level differences existed in the predicted candidate gene between the parents. These results provide a theoretical basis for cloning these QTLs and marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01233-0.

Construction of a High-Density Genetic Map and Identification of Leaf Trait-Related QTLs in Chinese Bayberry (Myrica rubra)

Chinese bayberry (Myrica rubra) is an economically important fruit tree that is grown in southern China. Owing to its over 10-year seedling period, the crossbreeding of bayberry is challenging. The characteristics of plant leaves are among the primary factors that control plant architecture and potential yields, making the analysis of leaf trait-related genetic factors crucial to the hybrid breeding of any plant. In the present study, molecular markers associated with leaf traits were identified via a whole-genome re-sequencing approach, and a genetic map was thereby constructed. In total, this effort yielded 902.11 Gb of raw data that led to the identification of 2,242,353 single nucleotide polymorphisms (SNPs) in 140 F1 individuals and parents (Myrica rubra cv. Biqizhong × Myrica rubra cv. 2012LXRM). The final genetic map ultimately incorporated 31,431 SNPs in eight linkage groups, spanning 1,351.85 cM. This map was then used to assemble and update previous scaffold genomic data at the chromosomal level. The genome size of M. rubra was thereby established to be 275.37 Mb, with 94.98% of sequences being assembled into eight pseudo-chromosomes. Additionally, 18 quantitative trait loci (QTLs) associated with nine leaf and growth-related traits were identified. Two QTL clusters were detected (the LG3 and LG5 clusters). Functional annotations further suggested two chlorophyll content-related candidate genes being identified in the LG5 cluster. Overall, this is the first study on the QTL mapping and identification of loci responsible for the regulation of leaf traits in M. rubra, offering an invaluable scientific for future marker-assisted selection breeding and candidate gene analyses.