Exploring DNA variant segregation types in pooled genome sequencing enables effective mapping of weeping trait in Malus

OcBSA: An NGS-based bulk segregant analysis tool for outcross populations

Constructing inbred lines for self-incompatible species and species with long generation times is challenging, making the use of F1 outcross/segregating populations the main strategy for genetic studies of such species. However, there is a lack of dedicated algorithms/tools for rapid quantitative trait locus (QTL) mapping using the F1 populations. To this end, we have designed and developed an algorithm/tool called OcBSA specifically for QTL mapping of F1 populations. OcBSA transforms the four-haplotype inheritance problem from the two heterozygous diploid parents of the F1 population into the two-haplotype inheritance problem common in current genetic studies by removing the two haplotypes from the heterozygous parent that do not contribute to phenotype segregation in the F1 population. Testing of OcBSA on 1800 simulated F1 populations demonstrated its advantages over other currently available tools in terms of sensitivity and accuracy. In addition, the broad applicability of OcBSA was validated by QTL mapping using seven reported F1 populations of apple, pear, peach, citrus, grape, tea, and rice. We also used OcBSA to map the QTL for flower color in a newly constructed F1 population of potato generated in this study. The OcBSA mapping result was verified by the insertion or deletion markers to be consistent with a previously reported locus harboring the ANTHOCYANIN 2 gene, which regulates potato flower color. Taken together, these results highlight the power and broad utility of OcBSA for QTL mapping using F1 populations and thus a great potential for functional gene mining in outcrossing species. For ease of use, we have developed both Windows and Linux versions of OcBSA, which are freely available at: https://gitee.com/Bioinformaticslab/OcBSA.

A single amino acid substitution in MdLAZY1A dominantly impairs shoot gravitropism in Malus

Plant architecture is 1 of the most important factors that determines crop yield potential and productivity. In apple (Malus domestica), genetic improvement of tree architecture has been challenging due to a long juvenile phase and growth as complex trees composed of a distinct scion and a rootstock. To better understand the genetic control of apple tree architecture, the dominant weeping growth phenotype was investigated. We report the identification of MdLAZY1A (MD13G1122400) as the genetic determinant underpinning the Weeping (W) locus that largely controls weeping growth in Malus. MdLAZY1A is 1 of the 4 paralogs in apple that are most closely related to AtLAZY1 involved in gravitropism in Arabidopsis (Arabidopsis thaliana). The weeping allele (MdLAZY1A-W) contains a single nucleotide mutation c.584T>C that leads to a leucine to proline (L195P) substitution within a predicted transmembrane domain that colocalizes with Region III, 1 of the 5 conserved regions in LAZY1-like proteins. Subcellular localization revealed that MdLAZY1A localizes to the plasma membrane and nucleus in plant cells. Overexpressing the weeping allele in apple cultivar Royal Gala (RG) with standard growth habit impaired its gravitropic response and altered the growth to weeping-like. Suppressing the standard allele (MdLAZY1A-S) by RNA interference (RNAi) in RG similarly changed the branch growth direction to downward. Overall, the L195P mutation in MdLAZY1A is genetically causal for weeping growth, underscoring not only the crucial roles of residue L195 and Region III in MdLAZY1A-mediated gravitropic response but also a potential DNA base editing target for tree architecture improvement in Malus and other crops.

Quantitative Traits of Interest in Apple Breeding and Their Implications for Selection

Apple breeding is a laborious and long-lasting process that requires qualified resources, land, time, and funds. In this study, more than 5000 F₁ apple hybrids from direct and testcrosses were analyzed. The results revealed how the phenotypic expression of the main quantitative traits of interest assessed in five half-sib families was controlled by the additive genetic effects and by non-additive effects of dominance and epistasis. The statistical number of hybrids required to ensure efficient selection increased exponentially with the number of desirable traits. The minimum number of progenies required to obtain a hybrid with associated quantitative traits of agronomic interest was highly variable. For two independent traits essential in selection (fruit size and quality), but incorporated together in the same hybrid, the statistical number was between about 30 and 300. If three more cumulative traits were added (a large number of fruits per tree, resistance/tolerance to apple scab, and powdery mildew attack), the limits increased to between 1500 and 18,000. The study highlighted the need for new apple varieties due to the narrowing of the genetic diversity of the cultivated species and how the choice of parents used in hybridizations (as well as the objectives pursued in the selection) can increase the efficiency of apple breeding.

The double flower variant of yellowhorn is due to a LINE1 transposon-mediated insertion

As essential organs of reproduction in angiosperms, flowers, and the genetic mechanisms of their development have been well characterized in many plant species but not in the woody tree yellowhorn (Xanthoceras sorbifolium). Here, we focused on the double flower phenotype in yellowhorn, which has high ornamental value. We found a candidate C-class gene, AGAMOUS1 (XsAG1), through bovine serum albumin sequencing and genetics analysis with a Long Interpersed Nuclear Elements 1 (LINE1) transposable element fragment (Xsag1-LINE1-1) inserted into its second intron that caused a loss-of-C-function and therefore the double flower phenotype. In situ hybridization of XsAG1 and analysis of the expression levels of other ABC genes were used to identify differences between single- and double-flower development processes. These findings enrich our understanding of double flower formation in yellowhorn and provide evidence that transposon insertions into genes can reshape plant traits in forest trees.

Gravity sensing and responses in the coordination of the shoot gravitropic setpoint angle

Gravity is one of the fundamental environmental cues that affect plant development. Indeed, the plant architecture in the shoots and roots is modulated by gravity. Stems grow vertically upward, whereas lateral organs, such as the lateral branches in shoots, tend to grow at a specific angle according to a gravity vector known as the gravitropic setpoint angle (GSA). During this process, gravity is sensed in specialised gravity-sensing cells named statocytes, which convert gravity information into biochemical signals, leading to asymmetric auxin distribution and driving asymmetric cell division/expansion in the organs to achieve gravitropism. As a hypothetical offset mechanism against gravitropism to determine the GSA, the anti-gravitropic offset (AGO) has been proposed. According to this concept, the GSA is a balance of two antagonistic growth components, that is gravitropism and the AGO. Although the nature of the AGO has not been clarified, studies have suggested that gravitropism and the AGO share a common gravity-sensing mechanism in statocytes. This review discusses the molecular mechanisms underlying gravitropism as well as the hypothetical AGO in the control of the GSA.

Pedigree reconstruction and genetic analysis of major ornamental characters of ornamental crabapple (Malus spp.) based on paternity analysis

Ornamental crabapple is an important woody ornamental plant in the Northern Hemisphere. Its flowers, fruits, leaves and tree habit are all important ornamental characters. As there has been no research on the selection of superior parents and phenotypic variation, new varieties of ornamental crabapple are mainly selected from open-pollination progeny. In order to explore the transmission rule of ornamental traits between parents and offspring of crabapple, and to provide a basis for the selection of hybrid parents for directional breeding, 14 pairs of SSR markers were used in this study for paternity analysis of 384 offspring from 4 female parents crossed with 91 candidate male parents. And 273 offspring (71.1%) were matched with only the father at a 95% strict confidence level. We reconstructed 7 full-sib families (number of progeny ≥ 10) on the basis of the paternity analysis results. Genetic analysis of characters in the full-sib families revealed that green leaves and white flowers were dominant traits. All the hybrid offspring from the white flower (♀) × non-white flower (♂) cross produced white flowers, while 7.04% produced non-white flowers when both parents had white flowers. The results showed that white flowers might be a dominant qualitative trait in crabapple, while the depth of red was a quantitative trait. The genetic characteristics of green and non-green leaves and the depth of red of the peel were similar to flower color. Compared with the upright and spreading traits, the weeping trait was recessive. Some progeny showed an earlier blooming period, indicating the possibility of breeding for blooming period. Our findings are important for parent screening and improving the breeding efficiency of new varieties in ornamental crabapple hybridization.

Molecular Mechanisms Regulating the Columnar Tree Architecture in Apple

The columnar apple cultivar ‘McIntosh Wijcik’ was discovered as a spontaneous mutant from the top of a ‘McIntosh’ tree in the early 1960s. ‘McIntosh Wijcik’ exhibits the columnar growth phenotype: compact and sturdy growth, short internodes, and very few lateral shoots. Classical genetic analysis revealed that the columnar growth phenotype of ‘McIntosh Wijcik’ is controlled by a single dominant gene, Co. This review focuses on the advances made toward understanding the molecular mechanisms of columnar growth in the last decade. Molecular studies have shown that an 8.2 kb insertion in the intergenic region of the Co locus is responsible for the columnar growth phenotype of ‘McIntosh Wijcik’, implying that the insertion affects the expression patterns of adjacent genes. Among the candidate genes in the Co region, the expression pattern of MdDOX-Co, putatively encoding 2-oxoglutarate-dependent dioxygenase (DOX), was found to vary between columnar and non-columnar apples. Recent studies have found three functions of MdDOX-Co: facilitating bioactive gibberellin deficiency, increasing strigolactone levels, and positively regulating abscisic acid levels. Consequently, changes in these plant hormone levels caused by the ectopic expression of MdDOX-Co in the aerial organs of ‘McIntosh Wijcik’ can lead to dwarf trees with fewer lateral branches. These findings will contribute to the breeding and cultivation of new columnar apple cultivars with improved fruit quality.

Comparative Association Mapping Reveals Conservation of Major Gene Resistance to White Pine Blister Rust in Southwestern White Pine (Pinus strobiformis) and Limber Pine (P. flexilis)

All native North American white pines are highly susceptible to white pine blister rust (WPBR) caused by Cronartium ribicola. Understanding genomic diversity and molecular mechanisms underlying genetic resistance to WPBR remains one of the great challenges in improvement of white pines. To compare major gene resistance (MGR) present in two species, southwestern white pine (Pinus strobiformis) Cr3 and limber pine (P. flexilis) Cr4, we performed association analyses of Cr3-controlled resistant traits using single nucleotide polymorphism (SNP) assays designed with Cr4-linked polymorphic genes. We found that ∼70% of P. flexilis SNPs were transferable to P. strobiformis. Furthermore, several Cr4-linked SNPs were significantly associated with the Cr3-controlled traits in P. strobiformis families. The most significantly associated SNP (M326511_1126R) almost colocalized with Cr4 on the Pinus consensus linkage group 8, suggesting that Cr3 and Cr4 might be the same R locus, or have localizations very close to each other in the syntenic region of the P. strobiformis and P. flexilis genomes. M326511_1126R was identified as a nonsynonymous SNP, causing amino acid change (Val₃₇₆Ile) in a putative pectin acetylesterase, with coding sequences identical between the two species. Moreover, top Cr3-associated SNPs were further developed as TaqMan genotyping assays, suggesting their usefulness as marker-assisted selection (MAS) tools to distinguish genotypes between quantitative resistance and MGR. This work demonstrates the successful transferability of SNP markers between two closely related white pine species in the hybrid zone, and the possibility for deployment of MAS tools to facilitate long-term WPBR management in P. strobiformis breeding and conservation.

De novo chromosome-level genome of a semi-dwarf cultivar of Prunus persica identifies the aquaporin PpTIP2 as responsible for temperature-sensitive semi-dwarf trait and PpB3-1 for flower type and size

Peach (Prunus persica) is one of the most important fruit crops globally, but its cultivation can be hindered by large tree size. 'Zhongyoutao 14' (CN14) is a temperature-sensitive semi-dwarf (TSSD) cultivar which might be useful as breeding stock. The genome of CN14 was sequenced and assembled de novo using single-molecule real-time sequencing and chromosome conformation capture assembly. A high-quality genome was assembled and annotated, with 228.82 Mb mapped to eight chromosomes. Eighty-six re-sequenced F₁ individuals and 334 previously re-sequenced accessions were used to identify candidate genes controlling TSSD and flower type and size. An aquaporin tonoplast intrinsic protein (PpTIP2) was a strong candidate gene for control of TSSD. Sequence variations in the upstream regulatory region of PpTIP2 correlated with different transcriptional activity at different temperatures. PpB3-1, a candidate gene for flower type (SH) and flower size, contributed to petal development and promoted petal enlargement. The locus of another 12 agronomic traits was identified through genome-wide association study. Most of these loci exhibited consistent and precise association signals, except for flesh texture and flesh adhesion. A 6015-bp insertion in exon 3 and a 26-bp insertion upstream of PpMYB25 were associated with fruit hairless. Along with a 70.5-Kb gap at the F-M locus in CN14, another two new alleles were identified in peach accessions. Our findings will not only promote genomic research and agronomic breeding in peach but also provide a foundation for the peach pan-genome.

De novo chromosome-level genome of a semi-dwarf cultivar of Prunus persica identifies the aquaporin PpTIP2 as responsible for temperature-sensitive semi-dwarf trait and PpB3-1 for flower type and size

Peach (Prunus persica) is one of the most important fruit crops globally, but its cultivation can be hindered by large tree size. 'Zhongyoutao 14' (CN14) is a temperature-sensitive semi-dwarf (TSSD) cultivar which might be useful as breeding stock. The genome of CN14 was sequenced and assembled de novo using single-molecule real-time sequencing and chromosome conformation capture assembly. A high-quality genome was assembled and annotated, with 228.82 Mb mapped to eight chromosomes. Eighty-six re-sequenced F₁ individuals and 334 previously re-sequenced accessions were used to identify candidate genes controlling TSSD and flower type and size. An aquaporin tonoplast intrinsic protein (PpTIP2) was a strong candidate gene for control of TSSD. Sequence variations in the upstream regulatory region of PpTIP2 correlated with different transcriptional activity at different temperatures. PpB3-1, a candidate gene for flower type (SH) and flower size, contributed to petal development and promoted petal enlargement. The locus of another 12 agronomic traits was identified through genome-wide association study. Most of these loci exhibited consistent and precise association signals, except for flesh texture and flesh adhesion. A 6015-bp insertion in exon 3 and a 26-bp insertion upstream of PpMYB25 were associated with fruit hairless. Along with a 70.5-Kb gap at the F-M locus in CN14, another two new alleles were identified in peach accessions. Our findings will not only promote genomic research and agronomic breeding in peach but also provide a foundation for the peach pan-genome.

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

PpMYB36 Encodes a MYB-Type Transcription Factor That Is Involved in Russet Skin Coloration in Pear (Pyrus pyrifolia)

Fruit color is one of the most important external qualities of pear (Pyrus pyrifolia) fruits. However, the mechanisms that control russet skin coloration in pear have not been well characterized. Here, we explored the molecular mechanisms that determine the russet skin trait in pear using the F₁ population derived from a cross between russet skin ('Niitaka') and non-russet skin ('Dangshansu') cultivars. Pigment measurements indicated that the lignin content in the skin of the russet pear fruits was greater than that in the non-russet pear skin. Genetic analysis revealed that the phenotype of the russet skin pear is associated with an allele of the PpRus gene. Using bulked segregant analysis combined with the genome sequencing (BSA-seq), we identified two simple sequence repeat (SSR) marker loci linked with the russet-colored skin trait in pear. Linkage analysis showed that the PpRus locus maps to the scaffold NW_008988489.1: 53297-211921 on chromosome 8 in the pear genome. In the mapped region, the expression level of LOC103929640 was significantly increased in the russet skin pear and showed a correlation with the increase of lignin content during the ripening period. Genotyping results demonstrated that LOC103929640 encoding the transcription factor MYB36 is the causal gene for the russet skin trait in pear. Particularly, a W-box insertion at the PpMYB36 promoter of russet skin pears is essential for PpMYB36-mediated regulation of lignin accumulation and russet coloration in pear. Overall, these results show that PpMYB36 is involved in the regulation of russet skin trait in pear.

Identification of the PmWEEP locus controlling weeping traits in Prunus mume through an integrated genome-wide association study and quantitative trait locus mapping

Weeping Prunus mume (mei) has long been cultivated in East Asia for its specific ornamental value. However, little is known about the regulatory mechanism of the weeping trait in mei, which limits molecular breeding for the improvement of weeping-type cultivars. Here, we quantified the weeping trait in mei using nested phenotyping of 214 accessions and 342 F₁ hybrids. Two major associated loci were identified from the genome-wide association study (GWAS), which was conducted using 3,014,409 single nucleotide polymorphisms (SNPs) derived from resequencing, and 8 QTLs and 55 epistatic loci were identified from QTL mapping using 7,545 specific lengths amplified fragment (SLAF) markers. Notably, an overlapping PmWEEP major QTL was fine mapped within a 0.29 Mb region on chromosome 7 (Pa7), and a core SNP locus closely associated with the weeping trait was screened and validated. Furthermore, a total of 22 genes in the PmWEEP QTL region were expressed in weeping or upright mei based on RNA-seq analysis. Among them, only a novel gene (Pm024213) containing a thioredoxin (Trx) domain was found to be close to the core SNP and specifically expressed in buds and branches of weeping mei. Co-expression analysis of Pm024213 showed that most of the related genes were involved in auxin and lignin biosynthesis. These findings provide insights into the regulatory mechanism of the weeping trait and effective molecular markers for molecular-assisted breeding in Prunus mume.

Comparative Transcriptomics and RNA-Seq-Based Bulked Segregant Analysis Reveals Genomic Basis Underlying Cronartium ribicola vcr2 Virulence

Breeding programs of five-needle pines have documented both major gene resistance (MGR) and quantitative disease resistance (QDR) to Cronartium ribicola (Cri), a non-native, invasive fungal pathogen causing white pine blister rust (WPBR). WPBR is one of the most deadly forest diseases in North America. However, Cri virulent pathotypes have evolved and can successfully infect and kill trees carrying resistance (R) genes, including vcr2 that overcomes MGR conferred by the western white pine (WWP, Pinus monticola) R gene (Cr2). In the absence of a reference genome, the present study generated a vcr2 reference transcriptome, consisting of about 20,000 transcripts with 1,014 being predicted to encode secreted proteins (SPs). Comparative profiling of transcriptomes and secretomes revealed vcr2 was significantly enriched for several gene ontology (GO) terms relating to oxidation-reduction processes and detoxification, suggesting that multiple molecular mechanisms contribute to pathogenicity of the vcr2 pathotype for its overcoming Cr2. RNA-seq-based bulked segregant analysis (BSR-Seq) revealed genome-wide DNA variations, including about 65,617 single nucleotide polymorphism (SNP) loci in 7,749 polymorphic genes shared by vcr2 and avirulent (Avcr2) pathotypes. An examination of the distribution of minor allele frequency (MAF) uncovered a high level of genomic divergence between vcr2 and Avcr2 pathotypes. By integration of extreme-phenotypic genome-wide association (XP-GWAS) analysis and allele frequency directional difference (AFDD) mapping, we identified a set of vcr2-associated SNPs within functional genes, involved in fungal virulence and other molecular functions. These included six SPs that were top candidate effectors with putative activities of reticuline oxidase, proteins with common in several fungal extracellular membrane (CFEM) domain or ferritin-like domain, polysaccharide lyase, rds1p-like stress responsive protein, and two Cri-specific proteins without annotation. Candidate effectors and vcr2-associated genes provide valuable resources for further deciphering molecular mechanisms of virulence and pathogenicity by functional analysis and the subsequent development of diagnostic tools for monitoring the virulence landscape in the WPBR pathosystems.

The roles of the IGT gene family in plant architecture: past, present, and future

Genetic improvement of architectural traits offers tremendous opportunities to dramatically improve crop densities, productivity, and ultimately sustainability. Among these, the orientation, or gravitropic set point angle (GSA), of plant organs is critical to optimize crop profiles, light capture, and nutrient acquisition. Mutant GSA phenotypes have been studied in plants since the 1930's but only recently have the underlying genes been identified. Many of these genes have turned out to fall within the IGT (LAZY1/DRO1/TAC1) family, which initially was not previously recognized due to the lack of sequence conservation of homologous genes across species. Here we discuss recent progress on IGT family genes in various plant species over the past century, review possible functional mechanisms, and provide further analysis of their evolution in land plants and their past and future roles in crop domestication.

Comparative gene expression analysis reveals that multiple mechanisms regulate the weeping trait in Prunus mume

Prunus mume (also known as Mei) is an important ornamental plant that is popular with Asians. The weeping trait in P. mume has attracted the attention of researchers for its high ornamental value. However, the formation of the weeping trait of woody plants is a complex process and the molecular basis of weeping stem development is unclear. Here, the morphological and histochemical characteristics and transcriptome profiles of upright and weeping stems from P. mume were studied. Significant alterations in the histochemical characteristics of upright and weeping stems were observed, and the absence of phloem fibres and less xylem in weeping stems might be responsible for their inability to resist gravity and to grow downward. Transcriptome analysis showed that differentially expressed genes (DEGs) were enriched in phenylpropanoid biosynthesis and phytohormone signal transduction pathways. To investigate the differential responses to hormones, upright and weeping stems were treated with IAA (auxin) and GA₃ (gibberellin A3), respectively, and the results revealed that weeping stems had a weaker IAA response ability and reduced upward bending angles than upright stems. On the contrary, weeping stems had increased upward bending angles than upright stems with GA₃ treatment. Compared to upright stems, interestingly, DEGs associated with diterpenoid biosynthesis and phenylpropanoid biosynthesis were significantly enriched after being treated with IAA, and expression levels of genes associated with phenylpropanoid biosynthesis, ABC transporters, glycosylphosphatidylinositol (GPI)—anchor biosynthesis were altered after being treated with GA₃ in weeping stems. Those results reveal that multiple molecular mechanisms regulate the formation of weeping trait in P. mume, which lays a theoretical foundation for the cultivation of new varieties.

Identification and Characterization of an Early Leaf Senescence Gene ELS1 in Soybean

Early leaf senescence phenotype in soybean could be helpful to shorten the maturation period and prevent green stem disorder. From a high-density mutation library, we identified two early leaf senescence soybean mutant lines, els1-1 (early leaf senescence 1) and els1-2. The chlorophyll contents of both els1-1 and els1-2 were low in pre-senescent leaves. They degraded rapidly in senescent leaves, revealing that ELS1 is involved in chlorophyll biosynthesis during leaf development and chlorophyll degradation during leaf senescence. The causal mutations in els1 were identified by next-generation sequencing-based bulked segregant analysis. ELS1 encodes the ortholog of the Arabidopsis CaaX-like protease BCM1, which is localized in chloroplasts. Soybean ELS1 was highly expressed in green tissue, especially in mature leaves. The accumulation of photosystem I core proteins and light-harvesting proteins in els1 was low even in pre-senescent leaves, and their degradation was accelerated during leaf senescence. These results suggest that soybean ELS1 is involved in both chlorophyll synthesis and degradation, consistent with the findings in Arabidopsis BCM1. The gene els1, characterized by early leaf senescence and subsequent early maturation, does not affect the flowering time. Hence, the early leaf senescence trait regulated by els1 helps shorten the harvesting period because of early maturation characteristics. The els1-1 allele with weakly impaired function of ELS1 has only a small effect on agricultural traits and could contribute to practical breeding.

Weeping candidate genes screened using comparative transcriptomic analysis of weeping and upright progeny in an F1 population of Prunus mume

Weeping is a specific plant architecture with high ornamental value. Despite the considerable importance of the weeping habit to landscaping applications and knowledge of plant architecture biology, little is known regarding the underlying molecular mechanisms. In this study, growth and phytohormone content were analyzed among the progeny of different branch types in an F1 mapping population of Prunus mume with varying plant architecture. Bulked segregant RNA sequencing was conducted to compare differences among progeny at a transcriptional level. The weeping habit appears to be a complex process regulated by a series of metabolic pathways, with photosynthesis and flavonoid biosynthesis highly enriched in differentially expressed genes between weeping and upright progeny. Based on functional annotation and homologous analyses, we identified 30 candidate genes related to weeping that merit further analysis, including 10 genes related to IAA and GA₃ biosynthesis, together with 6 genes related to secondary branch growth. The results of this study will facilitate further studies of the molecular mechanisms underlying the weeping habit in P. mume.

Identification of two QTLs associated with high fruit acidity in apple using pooled genome sequencing analysis

Acidity is a critical component determining apple fruit quality. Previous studies reported two major acidity quantitative trait loci (QTLs) on linkage groups (LGs) 16 (Ma) and 8 (Ma3), respectively, and their homozygous genotypes mama and ma3ma3 usually confer low titratable acidity (TA) (<3.0 mg ml-1) to apple fruit. However, apples of genotypes Ma- (MaMa and Mama) or Ma3- (Ma3Ma3 and Ma3ma3) frequently show an acidity range spanning both regular (TA 3.0-10.0 mg ml-1) and high (TA > 10 mg ml-1) acidity levels. To date, the genetic control for high-acidity apples remains essentially unknown. In order to map QTLs associated with high acidity, two genomic DNA pools, one for high acidity and the other for regular acidity, were created in an interspecific F1 population Royal Gala (Malus domestica) × PI 613988 (M. sieversii) of 191 fruit-bearing progenies. By Illumina paired-end sequencing of the high and regular acidity pools, 1,261,640 single-nucleotide variants (SNVs) commonly present in both pools were detected. Using allele frequency directional difference and density (AFDDD) mapping approach, one region on chromosome 4 and another on chromosome 6 were identified to be putatively associated with high acidity, and were named Ma6 and Ma4, respectively. Trait association analysis of DNA markers independently developed from the Ma6 and Ma4 regions confirmed the mapping of Ma6 and Ma4. In the background of MaMa, 20.6% of acidity variation could be explained by Ma6, 28.5% by Ma4, and 50.7% by the combination of both. The effects of Ma6 and Ma4 in the background of Mama were also significant, but lower. These findings provide important genetic insight into high acidity in apple.

Genetic background mutations drive neural circuit hyperconnectivity in a fragile X syndrome model

BACKGROUND

Neural circuits are initially assembled during development when neurons synapse with potential partners and later refined as appropriate connections stabilize into mature synapses while inappropriate contacts are eliminated. Disruptions to this synaptogenic process impair connectivity optimization and can cause neurodevelopmental disorders. Intellectual disability (ID) and autism spectrum disorder (ASD) are often characterized by synaptic overgrowth, with the maintenance of immature or inappropriate synapses. Such synaptogenic defects can occur through mutation of a single gene, such as fragile X mental retardation protein (FMRP) loss causing the neurodevelopmental disorder fragile X syndrome (FXS). FXS represents the leading heritable cause of ID and ASD, but many other genes that play roles in ID and ASD have yet to be identified.

RESULTS

In a Drosophila FXS disease model, one dfmr150M null mutant stock exhibits previously unreported axonal overgrowths at developmental and mature stages in the giant fiber (GF) escape circuit. These excess axon projections contain both chemical and electrical synapse markers, indicating mixed synaptic connections. Extensive analyses show these supernumerary synapses connect known GF circuit neurons, rather than new, inappropriate partners, indicating hyperconnectivity within the circuit. Despite the striking similarities to well-characterized FXS synaptic defects, this new GF circuit hyperconnectivity phenotype is driven by genetic background mutations in this dfmr150M stock. Similar GF circuit synaptic overgrowth is not observed in independent dfmr1 null alleles. Bulked segregant analysis (BSA) was combined with whole genome sequencing (WGS) to identify the quantitative trait loci (QTL) linked to neural circuit hyperconnectivity. The results reveal 8 QTL associated with inappropriate synapse formation and maintenance in the dfmr150M mutant background.

CONCLUSIONS

Synaptogenesis is a complex, precisely orchestrated neurodevelopmental process with a large cohort of gene products coordinating the connectivity, synaptic strength, and excitatory/inhibitory balance between neuronal partners. This work identifies a number of genetic regions that contain mutations disrupting proper synaptogenesis within a particularly well-mapped neural circuit. These QTL regions contain potential new genes involved in synapse formation and refinement. Given the similarity of the synaptic overgrowth phenotype to known ID and ASD inherited conditions, identifying these genes should increase our understanding of these devastating neurodevelopmental disease states.

Exploring DNA Variant Segregation Types Enables Mapping Loci for Recessive Phenotypic Suppression of Columnar Growth in Apple

Columnar apples trees, originated from a bud mutation 'Wijcik McIntosh,' develop a simple canopy and set fruit on spurs. These characteristics make them an important genetic resource for improvement of tree architecture. Genetic studies have uncovered that columnar growth habit is a dominant trait and is caused by a retroposon insertion that induces the expression of the neighboring gene Co encoding a 2OG-Fe(II) oxygenase. Here we report the genetic mapping of two loci of recessive suppressors (genes) c2 (on Chr10) and c3 (on Chr9) that are linked to repression of the retroposon-induced Co gene expression and associated columnar phenotype in 275 F₁ seedlings, which were developed from a reciprocal cross between two columnar selections heterozygous at the Co locus. The mapping was accomplished by sequencing a genomic pool comprising 18 columnar seedlings and another pool of 16 standard seedlings that also carry the retroposon insertion, and by exploring DNA variants of segregation types that are informative for mapping recessive traits in apple. The informative segregation types include <hk × hk>, <lm × ll>, <nn × np>, <lm × mm>, and <pp × np>, where each letter denotes one of the four DNA bases and the letters in bold represent variants in relation to the reference genome. The alleles in each first and third positions are assumed in linkage with the recessive suppressors' allele in the two parents, respectively. Using RNA-seq analysis, we further revealed that the Co gene together with the differentially expressed genes under loci c2 and c3 formed a co-expression gene-network module associated with growth habit, in which 12 MapMan Bins were enriched.

Identification of Differentially Up-regulated Genes in Apple with White Rot Disease

Fuji, a major apple cultivar in Korea, is susceptible to white rot. Apple white rot disease appears on the stem and fruit; the development of which deteriorates fruit quality, resulting in decreases in farmers' income. Thus, it is necessary to characterize molecular markers related to apple white rot resistance. In this study, we screened for differentially expressed genes between uninfected apple fruits and those infected with Botryosphaeria dothidea, the fungal pathogen that causes white rot. Antimicrobial tests suggest that a gene expression involved in the synthesis of the substance inhibiting the growth of B. dothidea in apples was induced by pathogen infection. We identified seven transcripts induced by the infection. The seven transcripts were homologous to genes encoding a flavonoid glucosyltransferase, a metallothionein-like protein, a senescence-induced protein, a chitinase, a wound-induced protein, and proteins of unknown function. These genes have functions related to responses to environmental stresses, including pathogen infections. Our results can be useful for the development of molecular markers for early detection of the disease or for use in breeding white rotresistant cultivars.

Screening of Applicable SSR Molecular Markers Linked to Creeping Trait in Crape Myrtle

Creeping plants have unique ornamental value because they have more branches and flowers and the creeping trait is rare in crape myrtle (Lagerstroemia indica L.). In this study, the first filial generation (F1) population was derived from Lagerstroemia fauriei Koehne (standard) and L. indica “Creole” (creeping) and the backcross1 (BC1) population was derived from the backcross of F1 individual S82 (creeping) and L. fauriei. The segregation of the creeping trait was analyzed for 174 seedlings of the BC1 population to examine the linkage relationship between simple sequence repeat (SSR) molecular markers and the creeping trait. Creeping genes were screened using bulked segregant analysis combined with 322 SSR primers, which were detected with good polymorphism. The results show that two SSR markers (S364 and LYS12) were detected, with genetic distances of 23.49 centimorgan (cM) and 25.86 cM from the loci controlling the plant opening angle trait and the branching angle trait, respectively. The accuracy rate for phenotypic verification using S364 and LYS12 was 76.51% and 74.14%, respectively. Our results provide basic information for the molecular marker-assisted selective breeding and cloning of the creeping gene to improve architecture diversity in the breeding of crape myrtle.

Branching out: new insights into the genetic regulation of shoot architecture in trees

Directional growth in all plants involves both phototropic and gravitropic responses. Accordingly, mechanisms controlling shoot architecture throughout the plant kingdom are likely similar. However, as forms vary between species due in part to gene copy number and functional divergence, some aspects of how plants predetermine and regulate architecture can differ. This is especially true when comparing annual herbaceous species (e.g. model plants) to woody perennials such as trees. In the past decade, inexpensive genomic sequencing and technological advances enabled gene discovery and functional analyses in trees. This led to the identification of genes associated with tree shoot architecture control. Here, we present recent discoveries on the regulation of shoot architectures for which causative genes have been identified, including dwarf, weeping, columnar, and pillar growth habits. We also discuss potential applications of these findings.

Apple fruit acidity is genetically diversified by natural variations in three hierarchical epistatic genes: MdSAUR37, MdPP2CH and MdALMTII

Many efforts have been made to map quantitative trait loci (QTLs) to facilitate practical marker-assisted selection (MAS) in plants. In the present study, using MapQTL and BSA-seq (bulk segregant analysis using next generation sequencing) with two independent pedigree-based populations, we identified four major genome-wide QTLs responsible for apple fruit acidity. Candidate genes were screened in major QTL regions, and three functional gene markers, including a non-synonymous A/G single-nucleotide polymorphism (SNP) in the coding region of MdPP2CH, a 36-bp insertion in the promoter of MdSAUR37 and a previously reported SNP in MdALMTII, were validated to influence the malate content of apple fruits. In addition, MdPP2CH inactivated three vacuolar H⁺ -ATPases (MdVHA-A3, MdVHA-B2 and MdVHA-D2) and one aluminium-activated malate transporter (MdALMTII) via dephosphorylation and negatively influenced fruit malate accumulation. The dephosphotase activity of MdPP2CH was suppressed by MdSAUR37, which implied a higher hierarchy of genetic interaction. Therefore, the MdSAUR37/MdPP2CH/MdALMTII chain cascaded hierarchical epistatic genetic effects to precisely determine apple fruit malate content. An A/G SNP (-1010) on the MdMYB44 promoter region from a major QTL (qtl08.1) was closely associated with fruit malate content. The predicted phenotype values (PPVs) were estimated using the tentative genotype values of the gene markers, and the PPVs were significantly correlated with the observed phenotype values. Our findings provide an insight into plant genome-based selection in apples and will aid in conducting research to understand the fundamental physiological basis of quantitative genetics.

Loss of a highly conserved sterile alpha motif domain gene (WEEP) results in pendulous branch growth in peach trees

Plant shoots typically grow upward in opposition to the pull of gravity. However, exceptions exist throughout the plant kingdom. Most conspicuous are trees with weeping or pendulous branches. While such trees have long been cultivated and appreciated for their ornamental value, the molecular basis behind the weeping habit is not known. Here, we characterized a weeping tree phenotype in Prunus persica (peach) and identified the underlying genetic mutation using a genomic sequencing approach. Weeping peach tree shoots exhibited a downward elliptical growth pattern and did not exhibit an upward bending in response to 90° reorientation. The causative allele was found to be an uncharacterized gene, Ppa013325, having a 1.8-Kb deletion spanning the 5' end. This gene, dubbed WEEP, was predominantly expressed in phloem tissues and encodes a highly conserved 129-amino acid protein containing a sterile alpha motif (SAM) domain. Silencing WEEP in the related tree species Prunus domestica (plum) resulted in more outward, downward, and wandering shoot orientations compared to standard trees, supporting a role for WEEP in directing lateral shoot growth in trees. This previously unknown regulator of branch orientation, which may also be a regulator of gravity perception or response, provides insights into our understanding of how tree branches grow in opposition to gravity and could serve as a critical target for manipulating tree architecture for improved tree shape in agricultural and horticulture applications.