Search for: All records

Abstract Nanopore signal analysis enables detection of nucleotide modifications from native DNA and RNA sequencing, providing both accurate genetic or transcriptomic and epigenetic information without additional library preparation. At present, only a limited set of modifications can be directly basecalled (for example, 5-methylcytosine), while most others require exploratory methods that often begin with alignment of nanopore signal to a nucleotide reference. We present Uncalled4, a toolkit for nanopore signal alignment, analysis and visualization. Uncalled4 features an efficient banded signal alignment algorithm, BAM signal alignment file format, statistics for comparing signal alignment methods and a reproducible de novo training method fork-mer-based pore models, revealing potential errors in Oxford Nanopore Technologies’ state-of-the-art DNA model. We apply Uncalled4 to RNA 6-methyladenine (m6A) detection in seven human cell lines, identifying 26% more modifications than Nanopolish using m6Anet, including in several genes where m6A has known implications in cancer. Uncalled4 is available open source atgithub.com/skovaka/uncalled4. 

Abstract Plants regenerated from seedling explants (hypocotyls and cotyledons) of the Solanaceae family membersPhysalis grisea(groundcherry),Solanum lycopersicum(tomato), andSolanum prinophyllum(forest nightshade) were used to determine the in vitro culture parameters that contribute to the incidence in polyploidization of tissue culture-derived plants (regenerants) from these species. We examined the possible effects of zeatin concentration in the plant regeneration medium, explant source, and species. Plants were grown to maturity under greenhouse conditions, pollen was collected and germinated. Flow cytometry analysis verified the utility of the pollen germination method for determining differences in ploidy, which was based on the number of pollen tubes produced with one tube representing diploid and two indicating polyploid. As for zeatin concentration, we assessed the effect of our standard method of initiation on medium containing 2 mg/l followed by 1 mg/l 2 weeks after culture initiation in comparison with 0.25, 0.5, and 1 mg/l throughout the culture lifetime. There were no major correlations for zeatin concentration on ploidy status across the species except for plants regenerated fromS. lycopersicumhypocotyl explants where the percentage of polyploid regenerants increased with increasing concentrations. As for species and explant effects,P. griseaplants regenerated from hypocotyl explants had the highest percentage of polyploid plants at 81% compared to 43% and 35% forS. lycopersicumandS. prinophyllum, respectively. From cotyledons, 8% ofS. lycopersicumand 20% ofS. prinophyllumwere polyploid. A comparison withP. griseacould not be made because cotyledon explants do not regenerate on zeatin-containing medium. The results indicated the incidence of polyploidization cannot be generalized for zeatin concentration, however, an influence of explant type and species was observed. Effects of increased ploidy on plant morphology were primarily larger flower and seed size; however, no significant differences were observed in plant or fruit size. 

ABSTRACT Cryptic genetic variants exert minimal or no phenotypic effects alone but have long been hypothesized to form a vast, hidden reservoir of genetic diversity that drives trait evolvability through epistatic interactions. This classical theory has been reinvigorated by pan-genome sequencing, which has revealed pervasive variation within gene families and regulatory networks, including extensive cis-regulatory changes, gene duplication, and divergence between paralogs. Nevertheless, empirical testing of cryptic variation’s capacity to fuel phenotypic diversification has been hindered by intractable genetics, limited allelic diversity, and inadequate phenotypic resolution. Here, guided by natural and engineered cis-regulatory cryptic variants in a recently evolved paralogous gene pair, we identified an additional pair of redundant trans regulators, establishing a regulatory network that controls tomato inflorescence architecture. By combining coding mutations with a cis-regulatory allelic series in populations segregating for all four network genes, we systematically constructed a collection of 216 genotypes spanning the full spectrum of inflorescence complexity and quantified branching in over 27,000 inflorescences. Analysis of the resulting high-resolution genotype-phenotype map revealed a layer of dose-dependent interactions within paralog pairs that enhances branching, culminating in strong, synergistic effects. However, we also uncovered an unexpected layer of antagonism between paralog pairs, where accumulating mutations in one pair progressively diminished the effects of mutations in the other. Our results demonstrate how gene regulatory network architecture and complex dosage effects from paralog diversification converge to shape phenotypic space under a hierarchical model of epistatic interactions. Given the prevalence of paralog evolution in genomes, we propose that paralogous cryptic variation within regulatory networks elicits hierarchies of epistatic interactions, catalyzing bursts of phenotypic change. Keyword:cryptic mutations, paralogs, redundancy, cis-regulatory, tomato, inflorescence, gene regulatory network, modeling, epistasis 

Abstract Pan-genomics and genome-editing technologies are revolutionizing breeding of global crops^1,2. A transformative opportunity lies in exchanging genotype-to-phenotype knowledge between major crops (that is, those cultivated globally) and indigenous crops (that is, those locally cultivated within a circumscribed area)^3–5to enhance our food system. However, species-specific genetic variants and their interactions with desirable natural or engineered mutations pose barriers to achieving predictable phenotypic effects, even between related crops^6,7. Here, by establishing a pan-genome of the crop-rich genusSolanum⁸and integrating functional genomics and pan-genetics, we show that gene duplication and subsequent paralogue diversification are major obstacles to genotype-to-phenotype predictability. Despite broad conservation of gene macrosynteny among chromosome-scale references for 22 species, including 13 indigenous crops, thousands of gene duplications, particularly within key domestication gene families, exhibited dynamic trajectories in sequence, expression and function. By augmenting our pan-genome with African eggplant cultivars⁹and applying quantitative genetics and genome editing, we dissected an intricate history of paralogue evolution affecting fruit size. The loss of a redundant paralogue of the classical fruit size regulatorCLAVATA3(CLV3)^10,11was compensated by a lineage-specific tandem duplication. Subsequent pseudogenization of the derived copy, followed by a large cultivar-specific deletion, created a single fusedCLV3allele that modulates fruit organ number alongside an enzymatic gene controlling the same trait. Our findings demonstrate that paralogue diversifications over short timescales are underexplored contingencies in trait evolvability. Exposing and navigating these contingencies is crucial for translating genotype-to-phenotype relationships across species. 

Societal Impact StatementGroundcherry (Physalis grisea) is a plant species grown for its flavorful fruit. The fruit drops from the plant, hence the common name groundcherry. This makes harvest cumbersome and puts the fruit at risk for carrying soil‐borne pathogens, therefore making them unsellable. Furthermore, insects often damage the plants, reducing yield. Advances in gene editing offer promise for addressing these issues and aiding home gardeners and farmers. Improvement will expand access to this nutritious fruit, rich in potassium, vitamin C, and antioxidants. Additionally, studies of its biology could serve as a model for improving other fruiting plants, particularly underutilized species. SummaryP. griseais an underutilized, semidomesticated fruit crop with rising agronomic value. Several resources have been developed for its use in fundamental biological research, including a plant transformation system and a high‐quality reference genome. Already,P. griseahas been used as a model to investigate biological phenomena including inflated calyx syndrome and gene compensation.P. griseahas also been used to demonstrate the potential of fast‐tracking domestication trait improvement through approaches such as CRISPR/Cas9 gene editing. This work has led to thePhysalisImprovement Project, which relies on reverse genetics to understand the mechanisms that underlie fruit abscission and plant–herbivore interactions to guide approaches for improvement of undesirable characteristics. CRISPR/Cas9 gene editing has been used to targetP. griseagenes that are suspected to act in fruit abscission, particularly orthologs of those that are reported in tomato abscission zone development. A similar approach is being taken to targetP. griseagenes involved in the withanolide biosynthetic pathway to determine the impact of withanolides on plant–herbivore interactions. Results from these research projects will lead to a greater understanding of important biological processes and will also generate knowledge needed to develop cultivars with reduced fruit drop and increased resistance to insect herbivory. 

Abstract Stem cell homeostasis is pivotal for continuous and programmed formation of organs in plants. The precise control of meristem proliferation is mediated by the evolutionarily conserved signaling that encompasses complex interactions among multiple peptide ligands and their receptor-like kinases. Here, we identified compensation mechanisms involving the CLAVATA1 (CLV1) receptor and its paralogs, BARELY ANY MERISTEMs (BAMs), for stem cell proliferation in two Solanaceae species, tomato and groundcherry. Genetic analyses of higher-order mutants deficient in multiple receptor genes, generated via CRISPR-Cas9 genome editing, reveal that tomato SlBAM1 and SlBAM2 compensate for slclv1 mutations. Unlike the compensatory responses between orthologous receptors observed in Arabidopsis, tomato slclv1 mutations do not trigger transcriptional upregulation of four SlBAM genes. The compensation mechanisms within receptors are also conserved in groundcherry, and critical amino acid residues of the receptors associated with the physical interaction with peptide ligands are highly conserved in Solanaceae plants. Our findings demonstrate that the evolutionary conservation of both compensation mechanisms and critical coding sequences between receptor-like kinases provides a strong buffering capacity during stem cell homeostasis in tomato and groundcherry. 

Abstract An enduring question in evolutionary biology concerns the degree to which episodes of convergent trait evolution depend on the same genetic programs, particularly over long timescales. Here we genetically dissected repeated origins and losses of prickles, sharp epidermal projections, that convergently evolved in numerous plant lineages. Mutations in a cytokinin hormone biosynthetic gene caused at least 16 independent losses of prickles in eggplants and wild relatives in the genusSolanum. Strikingly, homologs promote prickle formation across angiosperms that collectively diverged over 150 million years ago. By developing newSolanumgenetic systems, we leveraged this discovery to eliminate prickles in a wild species and an indigenously foraged berry. Our findings implicate a shared hormone-activation genetic program underlying evolutionarily widespread and recurrent instances of plant morphological innovation. 

Abstract Advancing crop genomics requires efficient genetic systems enabled by high-quality personalized genome assemblies. Here, we introduce RagTag, a toolset for automating assembly scaffolding and patching, and we establish chromosome-scale reference genomes for the widely used tomato genotype M82 along with Sweet-100, a new rapid-cycling genotype that we developed to accelerate functional genomics and genome editing in tomato. This work outlines strategies to rapidly expand genetic systems and genomic resources in other plant species. 

Developmental transitions require precise temporal and spatial control of gene expression. In plants, such regulation is critical for flower formation, which involves the progressive maturation of stem cell populations within shoot meristems to floral meristems, followed by rapid sequential differentiation into floral organs. Across plant taxa, these transitions are orchestrated by the F-box transcriptional cofactor geneUNUSUAL FLORAL ORGANS(UFO). The conserved and pleiotropic functions ofUFOoffer a useful framework for investigating how evolutionary processes have shaped the intricatecis-regulation of key developmental genes. By pinpointing a conserved promoter sequence in an accessible chromatin region of the tomato ortholog ofUFO, we engineered in vivo a series ofcis-regulatory alleles that caused both loss- and gain-of-function floral defects. These mutant phenotypes were linked to disruptions in predicted transcription factor binding sites for known transcriptional activators and repressors. Allelic combinations revealed dosage-dependent interactions between opposing alleles, influencing the penetrance and expressivity of gain-of-function phenotypes. These phenotypic differences support that robustness in tomato flower development requires precise temporal control ofUFOexpression dosage. Bridging our analysis toArabidopsis, we found that although homologous sequences to the tomato regulatory region are dispersed within theUFOpromoter, they maintain similar control over floral development. However, phenotypes from disrupting these sequences differ due to the differing expression patterns ofUFO. Our study underscores the complexcis-regulatory control of dynamic developmental genes and demonstrates that critical short stretches of regulatory sequences that recruit both activating and repressing machinery are conserved to maintain developmental robustness. 

The wide array of currently available genomes displays a wonderful diversity in size, composition, and structure and is quickly expanding thanks to several global biodiversity genomics initiatives. However, sequencing of genomes, even with the latest technologies, can still be challenging for both technical (e.g., small physical size, contaminated samples, or access to appropriate sequencing platforms) and biological reasons (e.g., germline-restricted DNA, variable ploidy levels, sex chromosomes, or very large genomes). In recent years,k-mer-based techniques have become popular to overcome some of these challenges. They are based on the simple process of dividing the analyzed sequences (e.g., raw reads or genomes) into a set of subsequences of lengthk, calledk-mers, and then analyzing the frequency or sequences of thosek-mers. Analyses based onk-mers allow for a rapid and intuitive assessment of complex sequencing data sets. Here, we provide a comprehensive review to the theoretical properties and practical applications ofk-mers in biodiversity genomics with a special focus on genome modeling.