Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Background

The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa.

Results

Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were ‘stress response’ and ‘development’ both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization.

Conclusion

We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-606) contains supplementary material, which is available to authorized users.     

Characterization of the Atg17-Atg29-Atg31 complex specifically required for starvation-induced autophagy in Saccharomyces cerevisiae

Nutrient starvation induces autophagy to degrade cytoplasmic materials in the vacuole/lysosomes. In the yeast, Saccharomyces cerevisiae, Atg17, Atg29, and Atg31/Cis1 are specifically required for autophagosome formation by acting as a scaffold complex essential for pre-autophagosomal structure (PAS) organization. Here, we show that these proteins constitutively form an Atg17-Atg29-Atg31 ternary complex, in which phosphorylated Atg31 is included. Reconstitution analysis of the ternary complex in E. coli indicates that the three proteins are included in equimolar amounts in the complex. The molecular mass of a monomeric Atg17-Atg29-Atg31 complex is calculated at 97 kDa; however, analytical ultracentrifugation shows that the molecular mass of the ternary complex is 198 kDa, suggesting a dimeric complex. We propose that this ternary complex acts as a functional unit for autophagosome formation.     

Two distinct genomic regions, harbouring the period and fruitless genes, affect male courtship song in Drosophila montana

Acoustic signals often have a significant role in pair formation and in species recognition. Determining the genetic basis of signal divergence will help to understand signal evolution by sexual selection and its role in the speciation process. An earlier study investigated quantitative trait locus for male courtship song carrier frequency (FRE) in Drosophila montana using microsatellite markers. We refined this study by adding to the linkage map markers for 10 candidate genes known to affect song production in Drosophila melanogaster. We also extended the analyses to additional song characters (pulse train length (PTL), pulse number (PN), interpulse interval, pulse length (PL) and cycle number (CN)). Our results indicate that loci in two different regions of the genome control distinct features of the courtship song. Pulse train traits (PTL and PN) mapped to the X chromosome, showing significant linkage with the period gene. In contrast, characters related to song pulse properties (PL, CN and carrier FRE) mapped to the region of chromosome 2 near the candidate gene fruitless, identifying these genes as suitable loci for further investigations. In previous studies, the pulse train traits have been found to vary substantially between Drosophila species, and so are potential species recognition signals, while the pulse traits may be more important in intra-specific mate choice.     

Genetics of cuticular hydrocarbon differences between males of the parasitoid wasps Nasonia giraulti and Nasonia vitripennis

Many insects rely on cuticular hydrocarbons (CHCs) as major recognition signals between individuals. Previous research on the genetics of CHCs has focused on Drosophila in which the roles of three desaturases and one elongase were highlighted. Comparable studies in other insect taxa have not been conducted so far. Here, we explore the genetics of CHCs in hybrids of the jewel wasps Nasonia giraulti and Nasonia vitripennis. We analyzed the CHC profiles of pure strain and of F(2) hybrid males using gas chromatography coupled with mass spectrometry and distinguished 54 peaks, of which we identified 52 as straight-chain, monounsaturated, or methyl-branched CHCs. The latter compound class proved to be particularly abundant and diverse in Nasonia. Quantitative trait locus (QTL) analysis suggests fixed genetic differences between the two strains in 42 of the 54 studied traits, making Nasonia a promising genetic model for identifying genes involved in CHC biosynthesis. QTL for methyl-branched CHCs partly clustered in genomic regions with high recombination rate: a possible indication for pleiotropic genes that control their biosynthesis, which is largely unexplored so far. Finally, we identified and mapped genes in the Nasonia genome with high similarity to genes that have been implicated in alkene biosynthesis in Drosophila and discuss those that match in their position with predicted QTL for alkenes.     

Identification of a Novel PNMA-MS1 Gene in Marsupials Suggests the LTR Retrotransposon-Derived PNMA Genes Evolved Differently in Marsupials and Eutherians

Two major gene families derived from Ty3/Gypsy long terminal repeat (LTR) retrotransposons were recently identified in mammals. The sushi-ichi retrotransposon homologue (SIRH) family comprises 12 genes: 11 in eutherians including Peg10 and Peg11/Rtl1 that have essential roles in the eutherian placenta and 1 that is marsupial specific. Fifteen and 12 genes were reported in the second gene family, para-neoplastic antigen MA (PNMA), in humans and mice, respectively, although their biological functions and evolutionary history remain largely unknown. Here, we identified two novel candidate PNMA genes, PNMA-MS1 and -MS2 in marsupials. Like all eutherian-specific PNMA genes, they exhibit the highest homology to a Gypsy12_DR (DR, Danio rerio) Gag protein. PNMA-MS1 is conserved in both Australian and South American marsupial species, the tammar wallaby and grey short-tailed opossum. However, no PNMA-MS1 orthologue was found in eutherians, monotremes or non-mammalian vertebrates. PNMA-MS1 was expressed in the ovary, mammary gland and brain during development and growth in the tammar, suggesting that PNMA-MS1 may have acquired a marsupial-specific function. However, PNMA-MS2 seems to be a pseudogene. The absence of marsupial orthologues of eutherian PNMA genes suggests that the retrotransposition events of the Gypsy12_DR-related retrotransposons that gave rise to the PNMA family occurred after the divergence of marsupials and eutherians.     

Identification of the soybean HyPRP family and specific gene response to Asian soybean rust disease

Soybean [Glycine max (L.) Merril], one of the most important crop species in the world, is very susceptible to abiotic and biotic stress. Soybean plants have developed a variety of molecular mechanisms that help them survive stressful conditions. Hybrid proline-rich proteins (HyPRPs) constitute a family of cell-wall proteins with a variable N-terminal domain and conserved C-terminal domain that is phylogenetically related to non-specific lipid transfer proteins. Members of the HyPRP family are involved in basic cellular processes and their expression and activity are modulated by environmental factors. In this study, microarray analysis and real time RT-qPCR were used to identify putative HyPRP genes in the soybean genome and to assess their expression in different plant tissues. Some of the genes were also analyzed by time-course real time RT-qPCR in response to infection by Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease. Our findings indicate that the time of induction of a defense pathway is crucial in triggering the soybean resistance response to P. pachyrhizi. This is the first study to identify the soybean HyPRP group B family and to analyze disease-responsive GmHyPRP during infection by P. pachyrhizi.     

The Transient Inactivation of the Master Cell Cycle Phosphatase Cdc14 Causes Genomic Instability in Diploid Cells of Saccharomyces cerevisiae

Genomic instability is a common feature found in cancer cells . Accordingly, many tumor suppressor genes identified in familiar cancer syndromes are involved in the maintenance of the stability of the genome during every cell division and are commonly referred to as caretakers. Inactivating mutations and epigenetic silencing of caretakers are thought to be the most important mechanisms that explain cancer-related genome instability. However, little is known of whether transient inactivation of caretaker proteins could trigger genome instability and, if so, what types of instability would occur. In this work, we show that a brief and reversible inactivation, during just one cell cycle, of the key phosphatase Cdc14 in the model organism Saccharomyces cerevisiae is enough to result in diploid cells with multiple gross chromosomal rearrangements and changes in ploidy. Interestingly, we observed that such transient loss yields a characteristic fingerprint whereby trisomies are often found in small-sized chromosomes, and gross chromosome rearrangements, often associated with concomitant loss of heterozygosity, are detected mainly on the ribosomal DNA-bearing chromosome XII. Taking into account the key role of Cdc14 in preventing anaphase bridges, resetting replication origins, and controlling spindle dynamics in a well-defined window within anaphase, we speculate that the transient loss of Cdc14 activity causes cells to go through a single mitotic catastrophe with irreversible consequences for the genome stability of the progeny.     

Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome

Despite growing evidence of rapid evolution in protein coding genes, the contribution of positive selection to intra- and interspecific differences in protein coding regions of the genome is unclear. We attempted to see if genes coding for secreted proteins and genes with narrow expression, specifically those preferentially expressed in the mammary gland, have diverged at a faster rate between domestic cattle (Bos taurus) and humans (Homo sapiens) than other genes and whether positive selection is responsible. Using a large data set, we identified groups of genes based on secretion and expression patterns and compared them for the rate of nonsynonymous (dN) and synonymous (dS) substitutions per site and the number of radical (Dr) and conservative (Dc) amino acid substitutions. We found evidence of rapid evolution in genes with narrow expression, especially for those expressed in the liver and mammary gland and for genes coding for secreted proteins. We compared common human polymorphism data with human-cattle divergence and found that genes with high evolutionary rates in human-cattle divergence also had a large number of common human polymorphisms. This argues against positive selection causing rapid divergence in these groups of genes. In most cases dN/dS ratios were lower in human-cattle divergence than in common human polymorphism presumably due to differences in the effectiveness of purifying selection between long-term divergence and short-term polymorphism.     

Overexpression of a LAM Domain Containing RNA-Binding Protein LARP1c Induces Precocious Leaf Senescence in Arabidopsis

Leaf senescence is the final stage of leaf life history, and it can be regulated by multiple internal and external cues. La-related proteins (LARPs), which contain a well-conserved La motif (LAM) domain and normally a canonical RNA recognition motif (RRM) or noncanonical RRM-like motif, are widely present in eukaryotes. Six LARP genes (LARP1a-1c and LARP6a-6c) are present in Arabidopsis, but their biological functions have not been studied previously. In this study, we investigated the biological roles of LARP1c from the LARP1 family. Constitutive or inducible overexpression of LARP1c caused premature leaf senescence. Expression levels of several senescence-associated genes and defense-related genes were elevated upon overexpression of LARP1c. The LARP1c null mutant 1c-1 impaired ABA-, SA-, and MeJA-induced leaf senescence in detached leaves. Gene expression profiles of LARP1c showed age-dependent expression in rosette leaves. Taken together, our results suggest LARP1c is involved in regulation of leaf senescence.     

Reduced mycorrhizal colonization (rmc) tomato mutant lacks expression of SymRK signaling pathway genes

Comparison of the expression of 13 genes involved in arbuscular mycorrhizal (AM) symbiosis was performed in a wild type tomato (Solanum lycopersicum cv 76R) and its reduced mycorrhizal colonization mutant rmc in response to colonization with Glomus fasiculatum. Four defense-related genes were induced to a similar extent in the mutant and wild type AM colonized plants, indicating a systemic response to AM colonization. Genes related to nutrient exchange between the symbiont partners showed higher expression in the AM roots of wild type plants than the mutant plants, which correlated with their arbuscular frequency. A symbiosis receptor kinase that is involved in both nodulation and AM symbiosis was not expressed in the rmc mutant. The fact that some colonization was observed in rmc was suggestive of the existence of an alternate colonization signaling pathway for AM symbiosis in this mutant.     

Defining Specificity Determinants of cGMP Mediated Gustatory Sensory Transduction in Caenorhabditis elegans

Cyclic guanosine monophosphate (cGMP) is a key secondary messenger used in signal transduction in various types of sensory neurons. The importance of cGMP in the ASE gustatory receptor neurons of the nematode Caenorhabditis elegans was deduced by the observation that multiple receptor-type guanylyl cyclases (rGCs), encoded by the gcy genes, and two presently known cyclic nucleotide-gated ion channel subunits, encoded by the tax-2 and tax-4 genes, are essential for ASE-mediated gustatory behavior. We describe here specific mechanistic features of cGMP-mediated signal transduction in the ASE neurons. First, we assess the specificity of the sensory functions of individual rGC proteins. We have previously shown that multiple rGC proteins are expressed in a left/right asymmetric manner in the functionally lateralized ASE neurons and are required to sense distinct salt cues. Through domain swap experiments among three different rGC proteins, we show here that the specificity of individual rGC proteins lies in their extracellular domains and not in their intracellular, signal-transducing domains. Furthermore, we find that rGC proteins are also sufficient to confer salt sensory responses to other neurons. Both findings support the hypothesis that rGC proteins are salt receptor proteins. Second, we identify a novel, likely downstream effector of the rGC proteins in gustatory signal transduction, a previously uncharacterized cyclic nucleotide-gated (CNG) ion channel, encoded by the che-6 locus. che-6 mutants show defects in gustatory sensory transduction that are similar to defects observed in animals lacking the tax-2 and tax-4 CNG channels. In contrast, thermosensory signal transduction, which also requires tax-2 and tax-4, does not require che-6, but requires another CNG, cng-3. We propose that CHE-6 may form together with two other CNG subunits, TAX-2 and TAX-4, a gustatory neuron-specific heteromeric CNG channel complex.     

Identification of Genes That Promote or Inhibit Olfactory Memory Formation in Drosophila

Genetic screens in Drosophila melanogaster and other organisms have been pursued to filter the genome for genetic functions important for memory formation. Such screens have employed primarily chemical or transposon-mediated mutagenesis and have identified numerous mutants including classical memory mutants, dunce and rutabaga. Here, we report the results of a large screen using panneuronal RNAi expression to identify additional genes critical for memory formation. We identified >500 genes that compromise memory when inhibited (low hits), either by disrupting the development and normal function of the adult animal or by participating in the neurophysiological mechanisms underlying memory formation. We also identified >40 genes that enhance memory when inhibited (high hits). The dunce gene was identified as one of the low hits and further experiments were performed to map the effects of the dunce RNAi to the α/β and γ mushroom body neurons. Additional behavioral experiments suggest that dunce knockdown in the mushroom body neurons impairs memory without significantly affecting acquisition. We also characterized one high hit, sickie, to show that RNAi knockdown of this gene enhances memory through effects in dopaminergic neurons without apparent effects on acquisition. These studies further our understanding of two genes involved in memory formation, provide a valuable list of genes that impair memory that may be important for understanding the neurophysiology of memory or neurodevelopmental disorders, and offer a new resource of memory suppressor genes that will aid in understanding restraint mechanisms employed by the brain to optimize resources.