A comprehensive genomic scan reveals gene dosage balance impacts on quantitative traits in Populus trees

A systems genetic analysis identifies putative mechanisms and candidate genes regulating vessel traits in poplar wood

Wood is the water conducting tissue of tree stems. Like most angiosperm trees, poplar wood contains water-conducting vessel elements whose functional properties affect water transport and growth rates, as well as susceptibility to embolism and hydraulic failure during water stress and drought. Here we used a unique hybrid poplar pedigree carrying genomically characterized chromosomal insertions and deletions to undertake a systems genomics analysis of vessel traits. We assayed gene expression in wood forming tissues from clonal replicates of genotypes covering dosage quantitative trait loci with insertions and deletions, genotypes with extreme vessel trait phenotypes, and control genotypes. A gene co-expression analysis was used to assign genes to modules, which were then used in integrative analyses to identify modules associated with traits, to identify putative molecular and cellular processes associated with each module, and finally to identify candidate genes using multiple criteria including dosage responsiveness. These analyses identified known processes associated with vessel traits including stress response, abscisic acid and cell wall biosynthesis, and in addition identified previously unexplored processes including cell cycle and protein ubiquitination. We discuss our findings relative to component processes contributing to vessel trait variation including signaling, cell cycle, cell expansion, and cell differentiation.

Copy number variation of the restorer Rf4 underlies human selection of three-line hybrid rice breeding

Cytoplasmic male sterility (CMS) lines are important for breeding hybrid crops, and utilization of CMS lines requires strong fertility restorer (Rf) genes. Rf4, a major Rf for Wild-Abortive CMS (CMS-WA), has been cloned in rice. However, the Rf4 evolution and formation of CMS-WA/Rf system remain elusive. Here, we show that the Rf4 locus emerges earlier than the CMS-WA gene WA352 in wild rice, and 69 haplotypes of the Rf4 locus are generated in the Oryza genus through the copy number and sequence variations. Eight of these haplotypes of the Rf4 locus are enriched in modern rice cultivars during natural and human selections, whereas non-functional rf4i is preferentially selected for breeding current CMS-WA lines. We further verify that varieties carrying two-copy Rf4 haplotype have stronger fertility restoration ability and are widely used in three-line hybrid rice breeding. Our findings increase our understanding of CMS/Rf systems and will likely benefit crop breeding.

Unlikely heroes on the long and winding road to potato inbreeding

Conversion of potato from a tetraploid, heterozygous, vegetatively propagated crop to a diploid F1 hybrid, propagated via botanical seed, would constitute a considerable advance for global agriculture, but faces multiple challenges. One such challenge is the difficulty in inbreeding potato, which involves purging deleterious alleles from its genome. This commentary discusses possible reasons for this difficulty and highlights a recent sequence-based effort to classify SNP variation, in potato germplasm, according to its deleterious potential. Tools and strategies connected to this database may facilitate development of F1 hybrids.

Genotypic variation rather than ploidy level determines functional trait expression in a foundation tree species in the presence and absence of environmental stress

BACKGROUND AND AIMS

At the population level, genetic diversity is a key determinant of a tree species' capacity to cope with stress. However, little is known about the relative importance of the different components of genetic diversity for tree stress responses. We compared how two sources of genetic diversity, genotype and cytotype (i.e. differences in ploidy levels), influence growth, phytochemical and physiological traits of Populus tremuloides in the presence and absence of environmental stress.

METHODS

In a series of field studies, we first assessed variation in traits across diploid and triploid aspen genotypes from Utah and Wisconsin under non-stressed conditions. In two follow-up experiments, we exposed diploid and triploid aspen genotypes from Wisconsin to individual and interactive drought stress and defoliation treatments and quantified trait variations under stress.

KEY RESULTS

We found that (1) tree growth and associated traits did not differ significantly between ploidy levels under non-stressed conditions. Instead, variation in tree growth and most other traits was driven by genotypic and population differences. (2) Genotypic differences were critical for explaining variation of most functional traits and their responses to stress. (3) Ploidy level played a subtle role in shaping traits and trait stress responses, as its influence was typically obscured by genotypic differences. (4) As an exception to the third conclusion, we showed that triploid trees expressed 17 % higher foliar defence (tremulacin) levels, 11 % higher photosynthesis levels and 23 % higher rubisco activity under well-watered conditions. Moreover, triploid trees displayed greater drought resilience than diploids as they produced 35 % more new tissue than diploids when recovering from drought stress.

CONCLUSION

Although ploidy level can strongly influence the ecology of tree species, those effects may be relatively small in contrast to the effects of genotypic variation in highly diverse species.

Chromoanagenesis in plants: triggers, mechanisms, and potential impact

Chromoanagenesis is a single catastrophic event that involves, in most cases, localized chromosomal shattering and reorganization, resulting in a dramatically restructured chromosome. First discovered in cancer cells, it has since been observed in various other systems, including plants. In this review, we discuss the origin, characteristics, and potential mechanisms underlying chromoanagenesis in plants. We report that multiple processes, including mutagenesis and genetic engineering, can trigger chromoanagenesis via a variety of mechanisms such as micronucleation, breakage-fusion-bridge (BFB) cycles, or chain-like translocations. The resulting rearranged chromosomes can be preserved during subsequent plant growth, and sometimes inherited to the next generation. Because of their high tolerance to genome restructuring, plants offer a unique system for investigating the evolutionary consequences and potential practical applications of chromoanagenesis.

Is Pollen Production of Birch Controlled by Genetics and Local Conditions?

Intraspecific genetic variation might limit the relevance of environmental factors on plant traits. For example, the interaction between genetics and (a-)biotic factors regulating pollen production are still poorly understood. In this study, we investigated pollen production of 28 birch (Betula pendula Roth) individuals in the years 2019−2021. We sampled catkins of eleven groups of genetically identical trees, which were partially topped, but of the same age and located in a seed plantation in southern Germany characterized by similar microclimatic conditions. Furthermore, we monitored environmental factors such as air temperature, characterized air quality (NO2, NOx and O3), and assessed potential solar radiation. We especially checked for differences between years as well as between and within clones and assessed the synchronicity of years with high/low pollen production. We present a robust mean for the pollen production of Betula pendula (1.66 million pollen grains per catkin). Our findings show temporal (H(2) = 46.29, p < 0.001) and clonal variations (H(4) = 21.44, p < 0.001) in pollen production. We conclude that synchronized high or low pollen production is not utterly site-specific and, in addition, not strictly dependent on genotypes. We suggest that appropriate clone selection based on application (seed plantation, urban planting) might be advantageous and encourage a long-term monitoring.

Variation of Chromosome Composition in a Full-Sib Population Derived From 2x × 3x Interploidy Cross of Populus

Interploidy cross commonly results in complex chromosome number and structural variations. In our previous study, a progeny with segregated ploidy levels was produced by an interploidy cross between diploid female parent Populus tomentosa × Populus bolleana clone TB03 and triploid male parent Populus alba × Populus berolinensis 'Yinzhong'. However, the chromosome compositions of aneuploid genotypes in the progeny were still unclear. In the present study, a microsatellite DNA allele counting-peak ratios (MAC-PR) method was employed to analyze allelic configurations of each genotype to clarify their chromosome compositions, while 45S rDNA fluorescence in situ hybridization (FISH) analysis was used to reveal the mechanism of chromosome number variation. Based on the MAC-PR analysis of 47 polymorphic simple sequence repeat (SSR) markers distributed across all 19 chromosomes of Populus, both chromosomal number and structural variations were detected for the progeny. In the progeny, 26 hypo-triploids, 1 hyper-triploid, 16 hypo-tetraploids, 10 tetraploids, and 5 hyper-tetraploids were found. A total of 13 putative structural variation events (duplications and/or deletions) were detected in 12 genotypes, involved in chromosomes 3, 6, 7, 14, 15, 16, and 18. The 46.2% (six events) structural variation events occurred on chromosome 6, suggesting that there probably is a chromosome breakpoint near the SSR loci of chromosome 6. Based on calculation of the allelic information, the transmission of paternal heterozygosity in the hypo-triploids, hyper-triploid, hypo-tetraploids, tetraploids, and hyper-tetraploids were 0.748, 0.887, 0.830, 0.833, and 0.836, respectively, indicating that the viable pollen gains of the male parent 'Yinzhong' were able to transmit high heterozygosity to progeny. Furthermore, 45S rDNA-FISH analysis showed that specific-chromosome segregation feature during meiosis and chromosome appointment in normal and fused daughter nuclei of telophase II of 'Yinzhong,' which explained that the formation of aneuploids and tetraploids in the progeny could be attributed to imbalanced meiotic chromosomal segregation and division restitution of 'Yinzhong,' The data of chromosomal composition and structural variation of each aneuploid in the full-sib progeny of TB03 × 'Yinzhong' lays a foundation for analyzing mechanisms of trait variation relying on chromosome or gene dosages in Populus.

One Hundred Years of Gene Balance: How Stoichiometric Issues Affect Gene Expression, Genome Evolution, and Quantitative Traits

A century ago experiments with the flowering plant Datura stramonium and the fruit fly Drosophila melanogaster revealed that adding an extra chromosome to a karyotype was much more detrimental than adding a whole set of chromosomes. This phenomenon was referred to as gene balance and has been recapitulated across eukaryotic species. Here, we retrace some developments in this field. Molecular studies suggest that the basis of balance involves stoichiometric relationships of multi-component interactions. This concept has implication for the mechanisms controlling gene expression, genome evolution, sex chromosome evolution/dosage compensation, speciation mechanisms, and the underlying genetics of quantitative traits.

LD-CNV: rapid and simple discovery of chromosomal translocations using linkage disequilibrium between copy number variable loci

Large-scale structural variations, such as chromosomal translocations, can have profound effects on fitness and phenotype, but are difficult to identify and characterize. Here, we describe a simple and effective method aimed at identifying translocations using only the dosage of sequence reads mapped on the reference genome. We binned reads on genomic segments sized according to sequencing coverage and identified instances when copy number segregated in populations. For each dosage-polymorphic 1 Mb bin, we tested independence, effectively an apparent linkage disequilibrium (LD), with other variable bins. In nine potato (Solanum tuberosum) dihaploid families translocations affecting pericentromeric regions were common and in two cases were due to genomic misassembly. In two populations, we found evidence for translocation affecting euchromatic arms. In cv. PI 310467, a nonreciprocal translocation between chromosomes (chr.) 7 and 8 resulted in a 5-3 copy number change affecting several Mb at the respective chromosome tips. In cv. "Alca Tarma," the terminal arm of chr. 4 translocated to the tip of chr. 1. Using oligonucleotide-based fluorescent in situ hybridization painting probes (oligo-FISH), we tested and confirmed the predicted arrangement in PI 310467. In 192 natural accessions of Arabidopsis thaliana, dosage haplotypes tended to vary continuously and resulted in higher noise, while apparent LD between pericentromeric regions suggested the effect of repeats. This method, LD-CNV, should be useful in species where translocations are suspected because it tests linkage without the need for genotyping.

Chromoanagenesis from radiation-induced genome damage in Populus

Chromoanagenesis is a genomic catastrophe that results in chromosomal shattering and reassembly. These extreme single chromosome events were first identified in cancer, and have since been observed in other systems, but have so far only been formally documented in plants in the context of haploid induction crosses. The frequency, origins, consequences, and evolutionary impact of such major chromosomal remodeling in other situations remain obscure. Here, we demonstrate the occurrence of chromoanagenesis in poplar (Populus sp.) trees produced from gamma-irradiated pollen. Specifically, in this population of siblings carrying indel mutations, two individuals exhibited highly frequent copy number variation (CNV) clustered on a single chromosome, one of the hallmarks of chromoanagenesis. Using short-read sequencing, we confirmed the presence of clustered segmental rearrangement. Independently, we identified and validated novel DNA junctions and confirmed that they were clustered and corresponded to these rearrangements. Our reconstruction of the novel sequences suggests that the chromosomal segments have reorganized randomly to produce a novel rearranged chromosome but that two different mechanisms might be at play. Our results indicate that gamma irradiation can trigger chromoanagenesis, suggesting that this may also occur when natural or induced mutagens cause DNA breaks. We further demonstrate that such events can be tolerated in poplar, and even replicated clonally, providing an attractive system for more in-depth investigations of their consequences.

Plant breeders often use radiation treatment to produce variation, with the goal of identifying new varieties with superior traits. We studied a population of poplar trees produced by gamma irradiation of pollen, and asked what kind of DNA changes were associated with this variation. We found many changes, most often in the form of added (insertions) or removed (deletions) pieces of DNA. We also found two lines with much more drastic changes. In those lines, we observed massive reorganization. We characterized these two lines in detail and found that catastrophic pulverization and random reassembly only occurred on a single chromosome. Looking closely at how the pieces were put back together suggest that the rearrangements in these two lines may have resulted from two slightly different mechanisms. This type of rearrangement is commonly observed in human cancer cells, but has rarely been observed in plants. We demonstrated here that they can be induced by gamma irradiation, indicating this type of event might be more widespread than we expected. Characterizing such genome restructuring instances helps to understand how genome instability can remodel chromosomes and affect genome function.

Genetic Regulation of Vessel Morphology in Populus

During secondary growth, forest trees can modify the anatomy of the wood produced by the vascular cambium in response to environmental conditions. Notably, the trees of the model angiosperm genus, Populus, reduce the risk of cavitation and hydraulic failure under water stress by producing water-conducting vessel elements with narrow lumens, which are more numerous and more interconnected with each other. Here, we determined the genetic architecture of vessel traits affecting hydraulic physiology and resilience to water stress. Vessel traits were measured for clonally replicated genotypes of a unique Populus deltoides x nigra population carrying genomically defined insertions and deletions that create gene dosage variation. We found significant phenotypic variation for all traits measured (mean vessel diameter, height-corrected mean vessel diameter, vessel frequency, height-corrected vessel frequency, vessel grouping index, and mean vessel circularity), and that all traits were under genetic control and showed moderate heritability values, ranging from 0.32 to 0.53. Whole-genome scans of correlations between gene dosage and phenotypic traits identified quantitative trait loci for tree height, mean vessel diameter, height-corrected mean vessel diameter, height-corrected vessel frequency, and vessel grouping index. Our results demonstrate that vessel traits affecting hydraulic physiology are under genetic control, and both pleiotropic and trait-specific quantitative trait loci are found for these traits.

Genomic imbalance determines positive and negative modulation of gene expression in diploid maize

Genomic imbalance caused by changing the dosage of individual chromosomes (aneuploidy) has a more detrimental effect than varying the dosage of complete sets of chromosomes (ploidy). We examined the impact of both increased and decreased dosage of 15 distal and 1 interstitial chromosomal regions via RNA-seq of maize (Zea mays) mature leaf tissue to reveal new aspects of genomic imbalance. The results indicate that significant changes in gene expression in aneuploids occur both on the varied chromosome (cis) and the remainder of the genome (trans), with a wider spread of modulation compared with the whole-ploidy series of haploid to tetraploid. In general, cis genes in aneuploids range from a gene-dosage effect to dosage compensation, whereas for trans genes the most common effect is an inverse correlation in that expression is modulated toward the opposite direction of the varied chromosomal dosage, although positive modulations also occur. Furthermore, this analysis revealed the existence of increased and decreased effects in which the expression of many genes under genome imbalance are modulated toward the same direction regardless of increased or decreased chromosomal dosage, which is predicted from kinetic considerations of multicomponent molecular interactions. The findings provide novel insights into understanding mechanistic aspects of gene regulation.

Predominantly inverse modulation of gene expression in genomically unbalanced disomic haploid maize

The phenotypic consequences of the addition or subtraction of part of a chromosome is more severe than changing the dosage of the whole genome. By crossing diploid trisomies to a haploid inducer, we identified 17 distal segmental haploid disomies that cover ∼80% of the maize genome. Disomic haploids provide a level of genomic imbalance that is not ordinarily achievable in multicellular eukaryotes, allowing the impact to be stronger and more easily studied. Transcriptome size estimates revealed that a few disomies inversely modulate most of the transcriptome. Based on RNA sequencing, the expression levels of genes located on the varied chromosome arms (cis) in disomies ranged from being proportional to chromosomal dosage (dosage effect) to showing dosage compensation with no expression change with dosage. For genes not located on the varied chromosome arm (trans), an obvious trans-acting effect can be observed, with the majority showing a decreased modulation (inverse effect). The extent of dosage compensation of varied cis genes correlates with the extent of trans inverse effects across the 17 genomic regions studied. The results also have implications for the role of stoichiometry in gene expression, the control of quantitative traits, and the evolution of dosage-sensitive genes.

A systems genetics approach to deciphering the effect of dosage variation on leaf morphology in Populus

Gene copy number variation is frequent in plant genomes of various species, but the impact of such gene dosage variation on morphological traits is poorly understood. We used a large population of Populus carrying genomically characterized insertions and deletions across the genome to systematically assay the effect of gene dosage variation on a suite of leaf morphology traits. A systems genetics approach was used to integrate insertion and deletion locations, leaf morphology phenotypes, gene expression, and transcriptional network data, to provide an overview of how gene dosage influences morphology. Dosage-sensitive genomic regions were identified that influenced individual or pleiotropic morphological traits. We also identified cis-expression quantitative trait loci (QTL) within these dosage QTL regions, a subset of which modulated trans-expression QTL as well. Integration of data types within a gene co-expression framework identified co-expressed gene modules that are dosage sensitive, enriched for dosage expression QTL, and associated with morphological traits. Functional description of these modules linked dosage-sensitive morphological variation to specific cellular processes, as well as candidate regulatory genes. Together, these results show that gene dosage variation can influence morphological variation through complex changes in gene expression, and suggest that frequently occurring gene dosage variation has the potential to likewise influence quantitative traits in nature.

Evolutionary network genomics of wood formation in a phylogenetic survey of angiosperm forest trees

Wood formation was present in early angiosperms, but has been highly modified through evolution to generate the anatomical diversity seen in extant angiosperm lineages. In this project, we modeled changes in gene coexpression relationships associated with the evolution of wood formation in a phylogenetic survey of 13 angiosperm tree species. Gravitropic stimulation was used as an experimental treatment to alter wood formation and also perturb gene expression. Gene transcript abundances were determined using RNA sequencing of developing wood tissues from upright trees, and from the top (tension wood) and bottom (opposite wood) tissues of gravistimulated trees. A network-based approach was employed to align gene coexpression networks across species based on orthologous relationships. A large-scale, multilayer network was modeled that identified both lineage-specific gene coexpression modules and modules conserved across multiple species. Functional annotation and analysis of modules identified specific regulatory processes associated with conserved modules, including regulation of hormones, protein phosphorylation, meristem development and epigenetic processes. Our results provide novel insights into the evolution and development of wood formation, and demonstrate the ability to identify biological processes and genes important for the evolution of a foundational trait in nonmodel, undomesticated forest trees.

Magnitude of modulation of gene expression in aneuploid maize depends on the extent of genomic imbalance

Aneuploidy has profound effects on an organism, typically more so than polyploidy, and the basis of this contrast is not fully understood. A dosage series of the maize long arm of chromosome 1 (1L) was used to compare relative global gene expression in different types and degrees of aneuploidy to gain insights into how the magnitude of genomic imbalance as well as hypoploidy affects global gene expression. While previously available methods require a selective examination of specific genes, RNA sequencing provides a whole-genome view of gene expression in aneuploids. Most studies of global aneuploidy effects have concentrated on individual types of aneuploids because multiple dose aneuploidies of the same genomic region are difficult to produce in most model genetic organisms. The genetic toolkit of maize allows the examination of multiple ploidies and 1-4 doses of chromosome arms. Thus, a detailed examination of expression changes both on the varied chromosome arms and elsewhere in the genome is possible, in both hypoploids and hyperploids, compared with euploid controls. Previous studies observed the inverse trans effect, in which genes not varied in DNA dosage were expressed in a negative relationship to the varied chromosomal region. This response was also the major type of changes found globally in this study. Many genes varied in dosage showed proportional expression changes, though some were seen to be partly or fully dosage compensated. It was also found that the effects of aneuploidy were progressive, with more severe aneuploids producing effects of greater magnitude.

Transcriptome Profiling Provides Insight into the Genes in Carotenoid Biosynthesis during the Mesocarp and Seed Developmental Stages of Avocado (Persea americana)

Avocado (Persea americana Mill.) is an economically important crop because of its high nutritional value. However, the absence of a sequenced avocado reference genome has hindered investigations of secondary metabolism. For next-generation high-throughput transcriptome sequencing, we obtained 365,615,152 and 348,623,402 clean reads as well as 109.13 and 104.10 Gb of sequencing data for avocado mesocarp and seed, respectively, during five developmental stages. High-quality reads were assembled into 100,837 unigenes with an average length of 847.40 bp (N50 = 1725 bp). Additionally, 16,903 differentially expressed genes (DEGs) were detected, 17 of which were related to carotenoid biosynthesis. The expression levels of most of these 17 DEGs were higher in the mesocarp than in the seed during five developmental stages. In this study, the avocado mesocarp and seed transcriptome were also sequenced using single-molecule long-read sequencing to acquired 25.79 and 17.67 Gb clean data, respectively. We identified 233,014 and 238,219 consensus isoforms in avocado mesocarp and seed, respectively. Furthermore, 104 and 59 isoforms were found to correspond to the putative 11 carotenoid biosynthetic-related genes in the avocado mesocarp and seed, respectively. The isoform numbers of 10 out of the putative 11 genes involved in the carotenoid biosynthetic pathway were higher in the mesocarp than those in the seed. Besides, alpha- and beta-carotene contents in the avocado mesocarp and seed during five developmental stages were also measured, and they were higher in the mesocarp than in the seed, which validated the results of transcriptome profiling. Gene expression changes and the associated variations in gene dosage could influence carotenoid biosynthesis. These results will help to further elucidate carotenoid biosynthesis in avocado.