Genomic Profiling by ALaP-Seq Reveals Transcriptional Regulation by PML Bodies through DNMT3A Exclusion.

The promyelocytic leukemia (PML) body is a phase-separated nuclear structure physically associated with chromatin, implying its crucial roles in genome functions. However, its role in transcriptional regulation is largely unknown. We developed APEX-mediated chromatin labeling and purification (ALaP) to identify the genomic regions proximal to PML bodies. We found that PML bodies associate with active regulatory regions across the genome and with ∼300 kb of the short arm of the Y chromosome (YS300) in mouse embryonic stem cells. The PML body association with YS300 is essential for the transcriptional activity of the neighboring Y-linked clustered genes. Mechanistically, PML bodies provide specific nuclear spaces that the de novo DNA methyltransferase DNMT3A cannot access, resulting in the steady maintenance of a hypo-methylated state at Y-linked gene promoters. Our study underscores a new mechanism for gene regulation in the 3D nuclear space and provides insights into the functional properties of nuclear structures for genome function.

Overexpression profiling reveals cellular requirements in the context of genetic backgrounds and environments.

Overexpression can help life adapt to stressful environments, making an examination of overexpressed genes valuable for understanding stress tolerance mechanisms. However, a systematic study of genes whose overexpression is functionally adaptive (GOFAs) under stress has yet to be conducted. We developed a new overexpression profiling method and systematically identified GOFAs in Saccharomyces cerevisiae under stress (heat, salt, and oxidative). Our results show that adaptive overexpression compensates for deficiencies and increases fitness under stress, like calcium under salt stress. We also investigated the impact of different genetic backgrounds on GOFAs, which varied among three S. cerevisiae strains reflecting differing calcium and potassium requirements for salt stress tolerance. Our study of a knockout collection also suggested that calcium prevents mitochondrial outbursts under salt stress. Mitochondria-enhancing GOFAs were only adaptive when adequate calcium was available and non-adaptive when calcium was deficient, supporting this idea. Our findings indicate that adaptive overexpression meets the cell's needs for maximizing the organism's adaptive capacity in the given environment and genetic context.

Transcriptomic and functional screening of weapon formation genes implies significance of cell adhesion molecules and female-biased genes in broad-horned flour beetle.

For understanding the evolutionary mechanism of sexually selected exaggerated traits, it is essential to uncover its molecular basis. By using broad-horned flour beetle that has male-specific exaggerated structures (mandibular horn, head horn and gena enlargement), we investigated the transcriptomic and functional characters of sex-biased genes. Comparative transcriptome of male vs. female prepupal heads elucidated 673 sex-biased genes. Counter-intuitively, majority of them were female-biased (584 genes), and GO enrichment analysis showed cell-adhesion molecules were frequently female-biased. This pattern motivated us to hypothesize that female-biased transcripts (i.e. the transcripts diminished in males) may play a role in outgrowth formation. Potentially, female-biased genes may act as suppressors of weapon structure. In order to test the functionality of female-biased genes, we performed RNAi-mediated functional screening for top 20 female-biased genes and 3 genes in the most enriched GO term (cell-cell adhesion, fat1/2/3, fat4 and dachsous). Knockdown of one transcription factor, zinc finger protein 608 (zfp608) resulted in the formation of male-like gena in females, supporting the outgrowth suppression function of this gene. Similarly, knockdown of fat4 induced rudimental, abnormal mandibular horn in female. fat1/2/3RNAi, fat4RNAi and dachsousRNAi males exhibited thick and/or short mandibular horns and legs. These cell adhesion molecules are known to regulate tissue growth direction and known to be involved in the weapon formation in Scarabaeoidea beetles. Functional evidence in phylogenetically distant broad-horned flour beetle suggest that cell adhesion genes are repeatedly deployed in the acquisition of outgrowth. In conclusion, this study clarified the overlooked functions of female-biased genes in weapon development.

GRAS transcription factors regulate cell division planes in moss overriding the default rule.

Plant cells are surrounded by a cell wall and do not migrate, which makes the regulation of cell division orientation crucial for development. Regulatory mechanisms controlling cell division orientation may have contributed to the evolution of body organization in land plants. The GRAS family of transcription factors was transferred horizontally from soil bacteria to an algal common ancestor of land plants. SHORTROOT (SHR) and SCARECROW (SCR) genes in this family regulate formative periclinal cell divisions in the roots of flowering plants, but their roles in nonflowering plants and their evolution have not been studied in relation to body organization. Here, we show that SHR cell autonomously inhibits formative periclinal cell divisions indispensable for leaf vein formation in the moss Physcomitrium patens, and SHR expression is positively and negatively regulated by SCR and the GRAS member LATERAL SUPPRESSOR, respectively. While precursor cells of a leaf vein lacking SHR usually follow the geometry rule of dividing along the division plane with the minimum surface area, SHR overrides this rule and forces cells to divide nonpericlinally. Together, these results imply that these bacterially derived GRAS transcription factors were involved in the establishment of the genetic regulatory networks modulating cell division orientation in the common ancestor of land plants and were later adapted to function in flowering plant and moss lineages for their specific body organizations.

Odd-Paired is Involved in Morphological Divergence of Snail-Feeding Beetles.

Body shape and size diversity and their evolutionary rates correlate with species richness at the macroevolutionary scale. However, the molecular genetic mechanisms underlying the morphological diversification across related species are poorly understood. In beetles, which account for one-fourth of the known species, adaptation to different trophic niches through morphological diversification appears to have contributed to species radiation. Here, we explored the key genes for the morphological divergence of the slender to stout body shape related to divergent feeding methods on large to small snails within the genus Carabus. We show that the zinc-finger transcription factor encoded by odd-paired (opa) controls morphological variation in the snail-feeding ground beetle Carabus blaptoides. Specifically, opa was identified as the gene underlying the slender to stout morphological difference between subspecies through genetic mapping and functional analysis via gene knockdown. Further analyses revealed that changes in opa cis-regulatory sequences likely contributed to the differences in body shape and size between C. blaptoides subspecies. Among opa cis-regulatory sequences, single nucleotide polymorphisms on the transcription factor binding sites may be associated with the morphological differences between C. blaptoides subspecies. opa was highly conserved in a wide range of taxa, especially in beetles. Therefore, opa may play an important role in adaptive morphological divergence in beetles.

Newt Hoxa13 has an essential and predominant role in digit formation during development and regeneration.

The 5'Hox genes play crucial roles in limb development and specify regions in the proximal-distal axis of limbs. However, there is no direct genetic evidence that Hox genes are essential for limb development in non-mammalian tetrapods or for limb regeneration. Here, we produced single to quadruple Hox13 paralog mutants using the CRISPR/Cas9 system in newts (Pleurodeles waltl), which have strong regenerative capacities, and also produced germline mutants. We show that Hox13 genes are essential for digit formation in development, as in mice. In addition, Hoxa13 has a predominant role in digit formation, unlike in mice. The predominance is probably due to the restricted expression pattern of Hoxd13 in limb buds and the strong dependence of Hoxd13 expression on Hoxa13. Finally, we demonstrate that Hox13 genes are also necessary for digit formation in limb regeneration. Our findings reveal that the general function of Hox13 genes is conserved between limb development and regeneration, and across taxa. The predominance of Hoxa13 function both in newt limbs and fish fins, but not in mouse limbs, suggests a potential contribution of Hoxa13 function in fin-to-limb transition.

Genomic and transcriptomic analyses of the subterranean termite Reticulitermes speratus: Gene duplication facilitates social evolution.

Termites are model social organisms characterized by a polyphenic caste system. Subterranean termites (Rhinotermitidae) are ecologically and economically important species, including acting as destructive pests. Rhinotermitidae occupies an important evolutionary position within the clade representing a transitional taxon between the higher (Termitidae) and lower (other families) termites. Here, we report the genome, transcriptome, and methylome of the Japanese subterranean termite Reticulitermes speratus Our analyses highlight the significance of gene duplication in social evolution in this termite. Gene duplication associated with caste-biased gene expression was prevalent in the R. speratus genome. The duplicated genes comprised diverse categories related to social functions, including lipocalins (chemical communication), cellulases (wood digestion and social interaction), lysozymes (social immunity), geranylgeranyl diphosphate synthase (social defense), and a novel class of termite lineage-specific genes with unknown functions. Paralogous genes were often observed in tandem in the genome, but their expression patterns were highly variable, exhibiting caste biases. Some of the assayed duplicated genes were expressed in caste-specific organs, such as the accessory glands of the queen ovary and the frontal glands of soldier heads. We propose that gene duplication facilitates social evolution through regulatory diversification, leading to caste-biased expression and subfunctionalization and/or neofunctionalization conferring caste-specialized functions.

PECT1, a rate-limiting enzyme in phosphatidylethanolamine biosynthesis, is involved in the regulation of stomatal movement in Arabidopsis.

An Arabidopsis mutant displaying impaired stomatal responses to CO2 , cdi4, was isolated by a leaf thermal imaging screening. The mutated gene PECT1 encodes CTP:phosphorylethanolamine cytidylyltransferase. The cdi4 exhibited a decrease in phosphatidylethanolamine levels and a defect in light-induced stomatal opening as well as low-CO2 -induced stomatal opening. We created RNAi lines in which PECT1 was specifically repressed in guard cells. These lines are impaired in their stomatal responses to low-CO2 concentrations or light. Fungal toxin fusicoccin (FC) promotes stomatal opening by activating plasma membrane H+ -ATPases in guard cells via phosphorylation. Arabidopsis H+ -ATPase1 (AHA1) has been reported to be highly expressed in guard cells, and its activation by FC induces stomatal opening. The cdi4 and PECT1 RNAi lines displayed a reduced stomatal opening response to FC. However, similar to in the wild-type, cdi4 maintained normal levels of phosphorylation and activation of the stomatal H+ -ATPases after FC treatment. Furthermore, the cdi4 displayed normal localization of GFP-AHA1 fusion protein and normal levels of AHA1 transcripts. Based on these results, we discuss how PECT1 could regulate CO2 - and light-induced stomatal movements in guard cells in a manner that is independent and downstream of the activation of H+ -ATPases. [Correction added on 15 May 2023, after first online publication: The third sentence is revised in this version.].

Slow decrease in temperature produces readthrough transcripts in mammalian hibernation.

Accumulating evidence suggests that various cellular stresses interfere with the end processing of mRNA synthesis and lead to the production of abnormally long transcripts, known as readthrough transcripts (RTTs), which extend beyond the termination sites. Small mammalian hibernators repeatedly enter a state referred to as deep torpor (DT), where the metabolic rate, respiration rate, and core body temperature become extremely low, which produces various types of cellular stresses and therefore induces RTTs. However, the types of stresses and processes around the DT that cause RTTs are unclear. In the present study, we showed that RTTs are produced from different gene loci in the livers of Syrian hamsters under DT and summer-like conditions. Moreover, in vitro analysis using hamster primary hepatocytes revealed that DT-specific RTTs are induced by a slow decline in temperature, as seen in body temperature in the entrance phase of DT, but not by rapid cold treatment or hypoxia. In addition, it was observed that RTTs were not elongated under a significantly cold temperature (4 °C). These results indicate that DT-specific RTTs are produced during the entrance phase of torpor by a slow decrease in body temperature.

Experimental evidence for the correlation between RNA structural fluctuations and the frequency of beneficial mutations.

RNA has been used as a model molecule to understand the adaptive evolution process owing to the simple relationship between the structure (i.e., phenotype) and sequence (i.e., genotype). RNA usually forms multiple substructures with similar thermodynamic stabilities, called structural fluctuations. Ancel and Fontana theoretically proposed that structural fluctuation is directly related to the ease of change in structures by mutations and thus works as a source of adaptive evolution; however, experimental verification is limited. Here, we analyzed 76 RNA genotypes that appeared in our previous in vitro evolution to examine whether (i) RNA fluctuation decreases as adaptive evolution proceeds and (ii) RNAs that have larger fluctuations tend to have higher frequencies of beneficial mutations. We first computationally estimated the structural fluctuations of all RNAs and observed that they tended to decrease as their fitness increased. We next measured the frequency of beneficial mutations for 10 RNA genotypes and observed that the total number of beneficial mutations was correlated with the size of the structural fluctuations. These results consistently support the idea that the structural fluctuation of RNA, at least those evaluated in our study, works as a source of adaptive evolution.

W chromosome sequences of two bombycid moths provide an insight into the origin of Fem.

Fem is a W-linked gene that encodes a piRNA precursor, and its product, Fem piRNA, is a master factor of female determination in Bombyx mori. Fem has low similarity to any known sequences, and the origin of Fem remains unclear. So far, two hypotheses have been proposed for the origin of Fem: The first hypothesis is that Fem is an allele of Masc, which assumes that the W chromosome was originally a homologous chromosome of the Z chromosome. The second hypothesis is that Fem arose by the transposition of Masc to the W chromosome. To explore the origin of Fem, we determined the W chromosome sequences of B. mori and, as a comparison, a closely relative bombycid species of Trilocha varians with a Fem-independent sex determination system. To our surprise, although the sequences of W and Z chromosomes show no homology to each other, a few pairs of homologues are shared by W and Z chromosomes, indicating the W chromosome of both species originated from Z chromosome. In addition, the W chromosome of T. varians lacks Fem, while the W chromosome of B. mori has over 100 copies of Fem. The high-quality assembly of the W chromosome of B. mori arose the third hypothesis about the origin of Fem: Fem is a chimeric sequence of multiple transposons. More than half of one transcriptional unit of Fem shows a significant homology to RTE-BovB. Moreover, the Fem piRNA-producing region could correspond to the boundary of the two transposons, gypsy and satellite DNA.

Monoterpene glucosides in Eustoma grandiflorum roots promote hyphal branching in arbuscular mycorrhizal fungi.

Host plant-derived strigolactones trigger hyphal branching in arbuscular mycorrhizal (AM) fungi, initiating a symbiotic interaction between land plants and AM fungi. However, our previous studies revealed that gibberellin-treated lisianthus (Eustoma grandiflorum, Gentianaceae) activates rhizospheric hyphal branching in AM fungi using unidentified molecules other than strigolactones. In this study, we analyzed independent transcriptomic data of E. grandiflorum and found that the biosynthesis of gentiopicroside (GPS) and swertiamarin (SWM), characteristic monoterpene glucosides in Gentianaceae, was upregulated in gibberellin-treated E. grandiflorum roots. Moreover, these metabolites considerably promoted hyphal branching in the Glomeraceae AM fungi Rhizophagus irregularis and Rhizophagus clarus. GPS treatment also enhanced R. irregularis colonization of the monocotyledonous crop chive (Allium schoenoprasum). Interestingly, these metabolites did not provoke the germination of the root parasitic plant common broomrape (Orobanche minor). Altogether, our study unveiled the role of GPS and SWM in activating the symbiotic relationship between AM fungi and E. grandiflorum.

Autoactivation of mycorrhizal symbiosis signaling through gibberellin deactivation in orchid seed germination.

Orchids parasitically depend on external nutrients from mycorrhizal fungi for seed germination. Previous findings suggest that orchids utilize a genetic system of mutualistic arbuscular mycorrhizal (AM) symbiosis, in which the plant hormone gibberellin (GA) negatively affects fungal colonization and development, to establish parasitic symbiosis. Although GA generally promotes seed germination in photosynthetic plants, previous studies have reported low sensitivity of GA in seed germination of mycoheterotrophic orchids where mycorrhizal symbiosis occurs concurrently. To elucidate the connecting mechanisms of orchid seed germination and mycorrhizal symbiosis at the molecular level, we investigated the effect of GA on a hyacinth orchid (Bletilla striata) seed germination and mycorrhizal symbiosis using asymbiotic and symbiotic germination methods. Additionally, we compared the transcriptome profiles between asymbiotically and symbiotically germinated seeds. Exogenous GA negatively affected seed germination and fungal colonization, and endogenous bioactive GA was actively converted to the inactive form during seed germination. Transcriptome analysis showed that B. striata shared many of the induced genes between asymbiotically and symbiotically germinated seeds, including GA metabolism- and signaling-related genes and AM-specific marker homologs. Our study suggests that orchids have evolved in a manner that they do not use bioactive GA as a positive regulator of seed germination and instead autoactivate the mycorrhizal symbiosis pathway through GA inactivation to accept the fungal partner immediately during seed germination.

Dual oxidase enables insect gut symbiosis by mediating respiratory network formation.

Most animals harbor a gut microbiota that consists of potentially pathogenic, commensal, and mutualistic microorganisms. Dual oxidase (Duox) is a well described enzyme involved in gut mucosal immunity by the production of reactive oxygen species (ROS) that antagonizes pathogenic bacteria and maintains gut homeostasis in insects. However, despite its nonspecific harmful activity on microorganisms, little is known about the role of Duox in the maintenance of mutualistic gut symbionts. Here we show that, in the bean bug Riptortus pedestris, Duox-dependent ROS did not directly contribute to epithelial immunity in the midgut in response to its mutualistic gut symbiont, Burkholderia insecticola Instead, we found that the expression of Duox is tracheae-specific and its down-regulation by RNAi results in the loss of dityrosine cross-links in the tracheal protein matrix and a collapse of the respiratory system. We further demonstrated that the establishment of symbiosis is a strong oxygen sink triggering the formation of an extensive network of tracheae enveloping the midgut symbiotic organ as well as other organs, and that tracheal breakdown by Duox RNAi provokes a disruption of the gut symbiosis. Down-regulation of the hypoxia-responsive transcription factor Sima or the regulators of tracheae formation Trachealess and Branchless produces similar phenotypes. Thus, in addition to known roles in immunity and in the formation of dityrosine networks in diverse extracellular matrices, Duox is also a crucial enzyme for tracheal integrity, which is crucial to sustain mutualistic symbionts and gut homeostasis. We expect that this is a conserved function in insects.

Evolution and functionalization of vitellogenin genes in the termite Reticulitermes speratus.

Eusociality has been commonly observed in distinct animal lineages. The reproductive division of labor is a particular feature, achieved by the coordination between fertile and sterile castes within the same nest. The sociogenomic approach in social hymenopteran insects indicates that vitellogenin (Vg) has undergone neo-functionalization in sterile castes. Here, to know whether Vgs have distinct roles in nonreproductive castes in termites, we investigated the unique characteristics of Vgs in the rhinotermitid termite Reticulitermes speratus. The four Vgs were identified from R. speratus (RsVg1-4), and RsVg3 sequences were newly identified using the RACE method. Molecular phylogenetic analysis supported the monophyly of the four termite Vgs. Moreover, the termites Vg1-3 and Vg4 were positioned in two different clades. The  dN/dS ratios indicated that the branch leading to the common ancestor of termite Vg4 was under weak purifying selection. Expression analyses among castes (reproductives, workers, and soldiers) and females (nymphs, winged alates, and queens) showed that RsVg1-3 was highly expressed in fertile queens. In contrast, RsVg4 was highly expressed in workers and female nonreproductives (nymphs and winged adults). Localization of RsVg4 messenger RNA was confirmed in the fat body of worker heads and abdomens. These results suggest that Vg genes are functionalized after gene duplication during termite eusocial transition and that Vg4 is involved in nonreproductive roles in termites.

Gene expression profiles of chemosensory genes of termite soldier and worker antennae.

Termites have an elaborate social system that involves cooperation and division of labour among colony members. Although this social system is regulated by chemical signals produced in the colony, it remains unclear how these signals are perceived by other members. Signal transduction is well known to be triggered by the reception of odorant molecules by some binding proteins in the antennae, after which, a signal is transmitted to chemosensory receptors. However, there is insufficient information on the role of chemosensory genes involved in signal transduction in termites. Here, we identified the genes involved in chemosensory reception in the termite Reticulitermes speratus and performed a genome-wide comparative transcriptome analysis of worker and soldier antennae. First, we identified 31 odorant-binding proteins (OBPs), and three chemosensory protein A (CheA) from the genome data. Thereafter, we performed RNA sequencing to compare the expression levels of OBPs, CheAs, and previously identified chemosensory receptor genes between worker and soldier antennae. There were no receptor genes with significant differences in expression between castes. However, the expression levels of three non-receptor odorant-detection/binding proteins (OBP, CheA, and Sensory neuron membrane protein) were significantly different between castes. Real-time qPCR (RT-qPCR) analysis using antennae and other head parts confirmed that these genes were highly expressed in soldier antennae. Finally, independent RT-qPCR analysis showed that the expression patterns of these genes were altered in soldiers from different social contexts. Present results suggest that gene expression levels of some non-receptors are affected by both castes and behavioural interactions among colony members in termites.

Intestinal expression patterns of transcription factors and markers for interstitial cells in the larval zebrafish.

For the digestion of food, it is important for the gut to be differentiated regionally and to have proper motor control. However, the number of transcription factors that regulate its development is still limited. Meanwhile, the interstitial cells of the gastrointestinal (GI) tract are necessary for intestinal motility in addition to the enteric nervous system. There are anoctamine1 (Ano1)-positive and platelet-derived growth factor receptor α (Pdgfra)-positive interstitial cells in mammal, but Pdgfra-positive cells have not been reported in the zebrafish. To identify new transcription factors involved in GI tract development, we used RNA sequencing comparing between larval and adult gut. We isolated 40 transcription factors that were more highly expressed in the larval gut. We demonstrated expression patterns of the 13 genes, 7 of which were newly found to be expressed in the zebrafish larval gut. Six of the 13 genes encode nuclear receptors. The osr2 is expressed in the anterior part, while foxP4 in its distal part. Also, we reported the expression pattern of pdgfra for the first time in the larval zebrafish gut. Our data provide fundamental knowledge for studying vertebrate gut regionalization and motility by live imaging using zebrafish.

Behavioral and brain- transcriptomic synchronization between the two opponents of a fighting pair of the fish Betta splendens.

Conspecific male animals fight for resources such as food and mating opportunities but typically stop fighting after assessing their relative fighting abilities to avoid serious injuries. Physiologically, how the fighting behavior is controlled remains unknown. Using the fighting fish Betta splendens, we studied behavioral and brain-transcriptomic changes during the fight between the two opponents. At the behavioral level, surface-breathing, and biting/striking occurred only during intervals between mouth-locking. Eventually, the behaviors of the two opponents became synchronized, with each pair showing a unique behavioral pattern. At the physiological level, we examined the expression patterns of 23,306 brain transcripts using RNA-sequencing data from brains of fighting pairs after a 20-min (D20) and a 60-min (D60) fight. The two opponents in each D60 fighting pair showed a strong gene expression correlation, whereas those in D20 fighting pairs showed a weak correlation. Moreover, each fighting pair in the D60 group showed pair-specific gene expression patterns in a grade of membership analysis (GoM) and were grouped as a pair in the heatmap clustering. The observed pair-specific individualization in brain-transcriptomic synchronization (PIBS) suggested that this synchronization provides a physiological basis for the behavioral synchronization. An analysis using the synchronized genes in fighting pairs of the D60 group found genes enriched for ion transport, synaptic function, and learning and memory. Brain-transcriptomic synchronization could be a general phenomenon and may provide a new cornerstone with which to investigate coordinating and sustaining social interactions between two interacting partners of vertebrates.

Potential contribution of intrinsic developmental stability toward body plan conservation.

BACKGROUND: Despite the morphological diversity of animals, their basic anatomical patterns-the body plans in each animal phylum-have remained highly conserved over hundreds of millions of evolutionary years. This is attributed to conservation of the body plan-establishing developmental period (the phylotypic period) in each lineage. However, the evolutionary mechanism behind this phylotypic period conservation remains under debate. A variety of hypotheses based on the concept of modern synthesis have been proposed, such as negative selection in the phylotypic period through its vulnerability to embryonic lethality. Here we tested a new hypothesis that the phylotypic period is developmentally stable; it has less potential to produce phenotypic variations than the other stages, and this has most likely led to the evolutionary conservation of body plans. RESULTS: By analyzing the embryos of inbred Japanese medaka embryos raised under the same laboratory conditions and measuring the whole embryonic transcriptome as a phenotype, we found that the phylotypic period has greater developmental stability than other stages. Comparison of phenotypic differences between two wild medaka populations indicated that the phylotypic period and its genes in this period remained less variational, even after environmental and mutational modifications accumulated during intraspecies evolution. Genes with stable expression levels were enriched with those involved in cell-cell signalling and morphological specification such as Wnt and Hox, implying possible involvement in body plan development of these genes. CONCLUSIONS: This study demonstrated the correspondence between the developmental stage with low potential to produce phenotypic variations and that with low diversity in micro- and macroevolution, namely the phylotypic period. Whereas modern synthesis explains evolution as a process of shaping of phenotypic variations caused by mutations, our results highlight the possibility that phenotypic variations are readily limited by the intrinsic nature of organisms, namely developmental stability, thus biasing evolutionary outcomes.

Arabidopsis Glucan Synthase-Like1 (GSL1) Is Required for Tolerance to Low-Calcium Conditions and Exhibits a Function Comparable to GSL10.

Crops that exhibit symptoms of calcium (Ca) deficiency constitute a major agricultural problem. Molecular breeding of resistant cultivars is a promising method for overcoming this problem. However, the involved genes must first be identified. Here, we show that the glucan synthase-like (GSL) 1 gene is essential for low-Ca tolerance in Arabidopsis thaliana. GSL1 is homologous to GSL10, which we previously showed was essential for low-Ca tolerance. Under low-Ca conditions, gsl1 mutants exhibit reduced growth and the onset of necrosis in new leaves. These symptoms are typical of Ca-deficient crops. A grafting experiment suggested that the shoot genotype, but not the root genotype, was important for the suppression of shoot necrosis. The ectopic accumulation of callose under low-Ca conditions was significantly reduced in gsl1 mutants compared with wild-type plants. Because the corresponding single-mutant phenotypes are similar, we investigated the interaction between GSL1 and GSL10 by testing the gsl1 gsl10 double mutant for sensitivity to low-Ca conditions. The double mutant exhibited a more severe phenotype than did the single mutants, indicating that the effects of GSL1 and GSL10 on low-Ca tolerance are additive. Because GSL genes are highly conserved within the plant kingdom, the GSL loci may be useful for breeding low-Ca tolerant crops.