Regulated functional alternative splicing in Drosophila

Improving rigor and reproducibility in chromatin immunoprecipitation assay data analysis workflows with Rocketchip

As genome sequencing technologies advance, the accumulation of sequencing data in public databases necessitates more robust and adaptable data analysis workflows. Here, we present Rocketchip, which aims to offer a solution to this problem by allowing researchers to easily compare and swap out different components of ChIP-seq, CUT&RUN, and CUT&Tag data analysis, thereby facilitating the identification of reliable analysis methodologies. Rocketchip enables researchers to efficiently process large datasets while ensuring reproducibility and allowing for the reanalysis of existing data. By supporting comparative analyses across different datasets and methodologies, Rocketchip contributes to the rigor and reproducibility of scientific findings. Furthermore, Rocketchip serves as a platform for benchmarking algorithms, allowing researchers to identify the most accurate and efficient analytical approaches to be applied to their data. In emphasizing reproducibility and adaptability, Rocketchip represents a significant step towards fostering robust scientific research practices.

High-quality chromosome-level genome assembly of Nicotiana benthamiana

Nicotiana benthamiana is a fundamental model organism in plant research. Recent advancements in genomic sequencing have revealed significant intraspecific genetic variations. This study addresses the pressing need for a precise genome sequence specific to its geographic origin by presenting a comprehensive genome assembly of the N. benthamiana LAB strain from the Republic of Korea (NbKLAB). We compare this assembly with the widely used NbLAB360 strain, shedding light on essential genomic differences between them. The outcome is a high-quality, chromosome-level genome assembly comprising 19 chromosomes, spanning 2,762 Mb, with an N50 of 142.6 Mb. Comparative analyses revealed notable variations, including 46,215 protein-coding genes, with an impressive 99.5% BUSCO completeness score. Furthermore, the NbKLAB assembly substantially improved the QV from 33% for NbLAB360 to 49%. This refined chromosomal genome assembly for N. benthamiana, in conjunction with comparative insights, provides a valuable resource for genomics research and molecular biology. This accomplishment forms a strong foundation for in-depth exploration into the intricacies of plant genetics and genomics, improved precision, and a comparative framework.

Bruno 1/CELF regulates splicing and cytoskeleton dynamics to ensure correct sarcomere assembly in Drosophila flight muscles

Muscles undergo developmental transitions in gene expression and alternative splicing that are necessary to refine sarcomere structure and contractility. CUG-BP and ETR-3-like (CELF) family RNA-binding proteins are important regulators of RNA processing during myogenesis that are misregulated in diseases such as Myotonic Dystrophy Type I (DM1). Here, we report a conserved function for Bruno 1 (Bru1, Arrest), a CELF1/2 family homolog in Drosophila, during early muscle myogenesis. Loss of Bru1 in flight muscles results in disorganization of the actin cytoskeleton leading to aberrant myofiber compaction and defects in pre-myofibril formation. Temporally restricted rescue and RNAi knockdown demonstrate that early cytoskeletal defects interfere with subsequent steps in sarcomere growth and maturation. Early defects are distinct from a later requirement for bru1 to regulate sarcomere assembly dynamics during myofiber maturation. We identify an imbalance in growth in sarcomere length and width during later stages of development as the mechanism driving abnormal radial growth, myofibril fusion, and the formation of hollow myofibrils in bru1 mutant muscle. Molecularly, we characterize a genome-wide transition from immature to mature sarcomere gene isoform expression in flight muscle development that is blocked in bru1 mutants. We further demonstrate that temporally restricted Bru1 rescue can partially alleviate hypercontraction in late pupal and adult stages, but it cannot restore myofiber function or correct structural deficits. Our results reveal the conserved nature of CELF function in regulating cytoskeletal dynamics in muscle development and demonstrate that defective RNA processing due to misexpression of CELF proteins causes wide-reaching structural defects and progressive malfunction of affected muscles that cannot be rescued by late-stage gene replacement.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Transcriptome profile in Drosophila Kc and S2 embryonic cell lines

Drosophila melanogaster cell lines are an important resource for a range of studies spanning genomics, molecular genetics, and cell biology. Amongst these valuable lines are Kc167 (Kc) and Schneider 2 (S2) cells, which were originally isolated in the late 1960s from embryonic sources and have been used extensively to investigate a broad spectrum of biological activities including cell-cell signaling and immune system function. Whole-genome tiling microarray analysis of total RNA from these two cell types was performed as part of the modENCODE project over a decade ago and revealed that they share a number of gene expression features. Here, we expand on these earlier studies by using deep-coverage RNA-sequencing approaches to investigate the transcriptional profile in Kc and S2 cells in detail. Comparison of the transcriptomes reveals that ∼75% of the 13,919 annotated genes are expressed at a detectable level in at least one of the cell lines, with the majority of these genes expressed at high levels in both cell lines. Despite the overall similarity of the transcriptional landscape in the two cell types, 2,588 differentially expressed genes are identified. Many of the genes with the largest fold change are known only by their "CG" designations, indicating that the molecular control of Kc and S2 cell identity may be regulated in part by a cohort of relatively uncharacterized genes. Our data also indicate that both cell lines have distinct hemocyte-like identities, but share active signaling pathways and express a number of genes in the network responsible for dorsal-ventral patterning of the early embryo.

TRA2: The dominant power of alternative splicing in tumors

The dysregulation of alternative splicing (AS) is frequently found in cancer and considered as key markers for cancer progression and therapy. Transformer 2 (TRA2), a nuclear RNA binding protein, consists of transformer 2 alpha homolog (TRA2A) and transformer 2 beta homolog (TRA2B), and plays a role in the regulation of pre-mRNA splicing. Growing evidence has been provided that TRA2A and TRA2B are dysregulated in several types of tumors, and participate in the regulation of proliferation, migration, invasion, and chemotherapy resistance in cancer cells through alteration of AS of cancer-related genes. In this review, we highlight the role of TRA2 in tumorigenesis and metastasis, and discuss potential molecular mechanisms how TRA2 influences tumorigenesis and metastasis via controlling AS of pre-mRNA. We propose that TRA2Ais a novel biomarker and therapeutic target for cancer progression and therapy.

Coordination of Alternative Splicing and Alternative Polyadenylation revealed by Targeted Long-Read Sequencing

Nervous system development is associated with extensive regulation of alternative splicing (AS) and alternative polyadenylation (APA). AS and APA have been extensively studied in isolation, but little is known about how these processes are coordinated. Here, the coordination of cassette exon (CE) splicing and APA in Drosophila was investigated using a targeted long-read sequencing approach we call Pull-a-Long-Seq (PL-Seq). This cost-effective method uses cDNA pulldown and Nanopore sequencing combined with an analysis pipeline to resolve the connectivity of alternative exons to alternative 3' ends. Using PL-Seq, we identified genes that exhibit significant differences in CE splicing depending on connectivity to short versus long 3'UTRs. Genomic long 3'UTR deletion was found to alter upstream CE splicing in short 3'UTR isoforms and ELAV loss differentially affected CE splicing depending on connectivity to alternative 3'UTRs. This work highlights the importance of considering connectivity to alternative 3'UTRs when monitoring AS events.

Cell-specific clock-controlled gene expression program regulates rhythmic fiber cell growth in cotton

BACKGROUND

The epidermis of cotton ovule produces fibers, the most important natural cellulose source for the global textile industry. However, the molecular mechanism of fiber cell growth is still poorly understood.

RESULTS

Here, we develop an optimized protoplasting method, and integrate single-cell RNA sequencing (scRNA-seq) and single-cell ATAC sequencing (scATAC-seq) to systematically characterize the cells of the outer integument of ovules from wild type and fuzzless/lintless (fl) cotton (Gossypium hirsutum). By jointly analyzing the scRNA-seq data from wildtype and fl, we identify five cell populations including the fiber cell type and construct the development trajectory for fiber lineage cells. Interestingly, by time-course diurnal transcriptomic analysis, we demonstrate that the primary growth of fiber cells is a highly regulated circadian rhythmic process. Moreover, we identify a small peptide GhRALF1 that circadian rhythmically controls fiber growth possibly through oscillating auxin signaling and proton pump activity in the plasma membrane. Combining with scATAC-seq, we further identify two cardinal cis-regulatory elements (CREs, TCP motif, and TCP-like motif) which are bound by the trans factors GhTCP14s to modulate the circadian rhythmic metabolism of mitochondria and protein translation through regulating approximately one third of genes that are highly expressed in fiber cells.

CONCLUSIONS

We uncover a fiber-specific circadian clock-controlled gene expression program in regulating fiber growth. This study unprecedentedly reveals a new route to improve fiber traits by engineering the circadian clock of fiber cells.

Evolution of pathogenicity-associated genes in Rhizoctonia solani AG1-IA by genome duplication and transposon-mediated gene function alterations

BACKGROUND

Rhizoctonia solani is a polyphagous fungal pathogen that causes diseases in crops. The fungal strains are classified into anastomosis groups (AGs); however, genomic complexity, diversification into the AGs and the evolution of pathogenicity-associated genes remain poorly understood.

RESULTS

We report a recent whole-genome duplication and sequential segmental duplications in AG1-IA strains of R. solani. Transposable element (TE) clusters have caused loss of synteny in the duplicated blocks and introduced differential structural alterations in the functional domains of several pathogenicity-associated paralogous gene pairs. We demonstrate that the TE-mediated structural variations in a glycosyl hydrolase domain and a GMC oxidoreductase domain in two paralogous pairs affect the pathogenicity of R. solani. Furthermore, to investigate the association of TEs with the natural selection and evolution of pathogenicity, we sequenced the genomes of forty-two rice field isolates of R. solani AG1-IA. The genomic regions with high population mutation rates and with the lowest nucleotide diversity are enriched with TEs. Genetic diversity analysis predicted the genes that are most likely under diversifying and purifying selections. We present evidence that a smaller variant of a glucosamine phosphate N-acetyltransferase (GNAT) protein, predicted to be under purifying selection, and an LPMP_AA9 domain-containing protein, predicted to be under diversifying selection, are important for the successful pathogenesis of R. solani in rice as well as tomato.

CONCLUSIONS

Our study has unravelled whole-genome duplication, TE-mediated neofunctionalization of genes and evolution of pathogenicity traits in R. solani AG1-IA. The pathogenicity-associated genes identified during the study can serve as novel targets for disease control.

The insect perspective on Z-disc structure and biology

Myofibrils are the intracellular structures formed by actin and myosin filaments. They are paracrystalline contractile cables with unusually well-defined dimensions. The sliding of actin past myosin filaments powers contractions, and the entire system is held in place by a structure called the Z-disc, which anchors the actin filaments. Myosin filaments, in turn, are anchored to another structure called the M-line. Most of the complex architecture of myofibrils can be reduced to studying the Z-disc, and recently, important advances regarding the arrangement and function of Z-discs in insects have been published. On a very small scale, we have detailed protein structure information. At the medium scale, we have cryo-electron microscopy maps, super-resolution microscopy and protein-protein interaction networks, while at the functional scale, phenotypic data are available from precise genetic manipulations. All these data aim to answer how the Z-disc works and how it is assembled. Here, we summarize recent data from insects and explore how it fits into our view of the Z-disc, myofibrils and, ultimately, muscles.

Identification, Phylogenetic and Expression Analyses of the AAAP Gene Family in Liriodendron chinense Reveal Their Putative Functions in Response to Organ and Multiple Abiotic Stresses

In this study, 52 AAAP genes were identified in the L. chinense genome and divided into eight subgroups based on phylogenetic relationships, gene structure, and conserved motif. A total of 48 LcAAAP genes were located on the 14 chromosomes, and the remaining four genes were mapped in the contigs. Multispecies phylogenetic tree and codon usage bias analysis show that the LcAAAP gene family is closer to the AAAP of Amborella trichopoda, indicating that the LcAAAP gene family is relatively primitive in angiosperms. Gene duplication events revealed six pairs of segmental duplications and one pair of tandem duplications, in which many paralogous genes diverged in function before monocotyledonous and dicotyledonous plants differentiation and were strongly purification selected. Gene expression pattern analysis showed that the LcAAAP gene plays a certain role in the development of Liriodendron nectary and somatic embryogenesis. Low temperature, drought, and heat stresses may activate some WRKY/MYB transcription factors to positively regulate the expression of LcAAAP genes to achieve long-distance transport of amino acids in plants to resist the unfavorable external environment. In addition, the GAT and PorT subgroups could involve gamma-aminobutyric acid (GABA) transport under aluminum poisoning. These findings could lay a solid foundation for further study of the biological role of LcAAAP and improvement of the stress resistance of Liriodendron.

U8 variants on the brain: a small nucleolar RNA and human disease

Small nucleolar RNAs (snoRNAs) are non-coding RNAs vital for ribosomal RNA (rRNA) maturation. The U8 snoRNA, encoded by the SNORD118 gene in humans, is an atypical C/D box snoRNA as it promotes rRNA cleavage rather than 2′–O–methylation and is unique to vertebrates. The U8 snoRNA is critical for cleavage events that produce the mature 5.8S and 28S rRNAs of the large ribosomal subunit. Unexpectedly, single nucleotide polymorphisms (SNPs) in the SNORD118 gene were recently found causal to the neurodegenerative disease leukoencephalopathy, brain calcifications, and cysts (LCC; aka Labrune syndrome), but its molecular pathogenesis is unclear. Here, we will review current knowledge on the function of the U8 snoRNA in ribosome biogenesis, and connect it to the preservation of brain function in humans as well as to its dysregulation in inherited white matter disease.

Tissue-specific regulation of translational readthrough tunes functions of the traffic jam transcription factor

Translational readthrough (TR) occurs when the ribosome decodes a stop codon as a sense codon, resulting in two protein isoforms synthesized from the same mRNA. TR has been identified in several eukaryotic organisms; however, its biological significance and mechanism remain unclear. Here, we quantify TR of several candidate genes in Drosophila melanogaster and characterize the regulation of TR in the large Maf transcription factor Traffic jam (Tj). Using CRISPR/Cas9-generated mutant flies, we show that the TR-generated Tj isoform is expressed in a subset of neural cells of the central nervous system and is excluded from the somatic cells of gonads. Control of TR in Tj is critical for preservation of neuronal integrity and maintenance of reproductive health. The tissue-specific distribution of a release factor splice variant, eRF1H, plays a critical role in modulating differential TR of leaky stop codon contexts. Fine-tuning of gene regulatory functions of transcription factors by TR provides a potential mechanism for cell-specific regulation of gene expression.

Genome-wide identification and expression analysis of sucrose nonfermenting-1-related protein kinase (SnRK) genes in Triticum aestivum in response to abiotic stress

The SnRK gene family is a key regulator that plays an important role in plant stress response by phosphorylating the target protein to regulate subsequent signaling pathways. This study was aimed to perform a genome-wide analysis of the SnRK gene family in wheat and the expression profiling of SnRKs in response to abiotic stresses. An in silico analysis identified 174 SnRK genes, which were then categorized into three subgroups (SnRK1/2/3) on the basis of phylogenetic analyses and domain types. The gene intron-exon structure and protein-motif composition of SnRKs were similar within each subgroup but different amongst the groups. Gene duplication and synteny between the wheat and Arabidopsis genomes was also investigated in order to get insight into the evolutionary aspects of the TaSnRK family genes. The result of cis-acting element analysis showed that there were abundant stress- and hormone-related cis-elements in the promoter regions of 129 SnRK genes. Furthermore, quantitative real-time PCR data revealed that heat, salt and drought treatments enhanced TaSnRK2.11 expression, suggesting that it might be a candidate gene for abiotic stress tolerance. We also identified eight microRNAs targeting 16 TaSnRK genes which are playing important role across abiotic stresses and regulation in different pathways. These findings will aid in the functional characterization of TaSnRK genes for further research.

A Candidate RNAi Screen Reveals Diverse RNA-Binding Protein Phenotypes in Drosophila Flight Muscle

The proper regulation of RNA processing is critical for muscle development and the fine-tuning of contractile ability among muscle fiber-types. RNA binding proteins (RBPs) regulate the diverse steps in RNA processing, including alternative splicing, which generates fiber-type specific isoforms of structural proteins that confer contractile sarcomeres with distinct biomechanical properties. Alternative splicing is disrupted in muscle diseases such as myotonic dystrophy and dilated cardiomyopathy and is altered after intense exercise as well as with aging. It is therefore important to understand splicing and RBP function, but currently, only a small fraction of the hundreds of annotated RBPs expressed in muscle have been characterized. Here, we demonstrate the utility of Drosophila as a genetic model system to investigate basic developmental mechanisms of RBP function in myogenesis. We find that RBPs exhibit dynamic temporal and fiber-type specific expression patterns in mRNA-Seq data and display muscle-specific phenotypes. We performed knockdown with 105 RNAi hairpins targeting 35 RBPs and report associated lethality, flight, myofiber and sarcomere defects, including flight muscle phenotypes for Doa, Rm62, mub, mbl, sbr, and clu. Knockdown phenotypes of spliceosome components, as highlighted by phenotypes for A-complex components SF1 and Hrb87F (hnRNPA1), revealed level- and temporal-dependent myofibril defects. We further show that splicing mediated by SF1 and Hrb87F is necessary for Z-disc stability and proper myofibril development, and strong knockdown of either gene results in impaired localization of kettin to the Z-disc. Our results expand the number of RBPs with a described phenotype in muscle and underscore the diversity in myofibril and transcriptomic phenotypes associated with splicing defects. Drosophila is thus a powerful model to gain disease-relevant insight into cellular and molecular phenotypes observed when expression levels of splicing factors, spliceosome components and splicing dynamics are altered.

Genome-wide identification and characterization of bZIP transcription factors and their expression profile under abiotic stresses in Chinese pear (Pyrus bretschneideri)

BACKGROUND

In plants, basic leucine zipper transcription factors (TFs) play important roles in multiple biological processes such as anthesis, fruit growth & development and stress responses. However, systematic investigation and characterization of bZIP-TFs remain unclear in Chinese white pear. Chinese white pear is a fruit crop that has important nutritional and medicinal values.

RESULTS

In this study, 62 bZIP genes were comprehensively identified from Chinese Pear, and 54 genes were distributed among 17 chromosomes. Frequent whole-genome duplication (WGD) and dispersed duplication (DSD) were the major driving forces underlying the bZIP gene family in Chinese white pear. bZIP-TFs are classified into 13 subfamilies according to the phylogenetic tree. Subsequently, purifying selection plays an important role in the evolution process of PbbZIPs. Synteny analysis of bZIP genes revealed that 196 orthologous gene pairs were identified between Pyrus bretschneideri, Fragaria vesca, Prunus mume, and Prunus persica. Moreover, cis-elements that respond to various stresses and hormones were found on the promoter regions of PbbZIP, which were induced by stimuli. Gene structure (intron/exon) and different compositions of motifs revealed that functional divergence among subfamilies. Expression pattern of PbbZIP genes differential expressed under hormonal treatment abscisic acid, salicylic acid, and methyl jasmonate  in pear fruits by real-time qRT-PCR.

CONCLUSIONS

Collectively, a systematic analysis of gene structure, motif composition, subcellular localization, synteny analysis, and calculation of synonymous (Ks) and non-synonymous (Ka) was performed in Chinese white pear. Sixty-two bZIP-TFs in Chinese pear were identified, and their expression profiles were comprehensively analyzed under ABA, SA, and MeJa hormones, which respond to multiple abiotic stresses and fruit growth and development. PbbZIP gene occurred through Whole-genome duplication and dispersed duplication events. These results provide a basic framework for further elucidating the biological function characterizations under multiple developmental stages and abiotic stress responses.

Does the Future of Antibiotics Lie in Secondary Metabolites Produced by Xenorhabdus spp.? A Review

The over-prescription of antibiotics for treatment of infections is primarily to blame for the increase in bacterial resistance. Added to the problem is the slow rate at which novel antibiotics are discovered and the many processes that need to be followed to classify antimicrobials safe for medical use. Xenorhabdus spp. of the family Enterobacteriaceae, mutualistically associated with entomopathogenic nematodes of the genus Steinernema, produce a variety of antibacterial peptides, including bacteriocins, depsipeptides, xenocoumacins and PAX (peptide antimicrobial-Xenorhabdus) peptides, plus additional secondary metabolites with antibacterial and antifungal activity. The secondary metabolites of some strains are active against protozoa and a few have anti-carcinogenic properties. It is thus not surprising that nematodes invaded by a single strain of a Xenorhabdus species are not infected by other microorganisms. In this review, the antimicrobial compounds produced by Xenorhabdus spp. are listed and the gene clusters involved in synthesis of these secondary metabolites are discussed. We also review growth conditions required for increased production of antimicrobial compounds.

Genome-wide identification of the BASS gene family in four Gossypium species and functional characterization of GhBASSs against salt stress

Bile acid sodium symporter (BASS) family proteins encode a class of sodium/solute symporters. Even though the sodium transporting property of BASSs in mammals was well studied, their sodium transportability and functional roles in plant salt tolerance remained largely unknown. Here, BASS family members from 4 cotton species, as well as 30 other species were identified. Then, they were designated as members of BASS1 to BASS5 subfamilies according to their sequence similarity and phylogenetic relationships. There were 8, 11, 16 and 18 putative BASS genes in four cotton species. While whole-genome duplications (WGD) and segmental duplications rendered the expansion of the BASS gene family in cotton, BASS gene losses occurred in the tetraploid cotton during the evolution from diploids to allotetraploids. Concerning functional characterizations, the transcript profiling of GhBASSs revealed that they not only preferred tissue-specific expression but also were differently induced by various stressors and phytohormones. Gene silencing and overexpression experiments showed that GhBASS1 and GhBASS3 positively regulated, whereas GhBASS2, GhBASS4 and GhBASS5 negatively regulated plant salt tolerance. Taken together, BASS family genes have evolved before the divergence from the common ancestor of prokaryotes and eukaryotes, and GhBASSs are plastidial sodium-dependent metabolite co-transporters that can influence plant salt tolerance.

Gene expression and alternative splicing dynamics are perturbed in female head transcriptomes following heterospecific copulation

BACKGROUND

Despite the growing interest in the female side of copulatory interactions, the roles played by differential expression and alternative splicing mechanisms of pre-RNA on tissues outside of the reproductive tract have remained largely unknown. Here we addressed these questions in the context of con- vs heterospecific matings between Drosophila mojavensis and its sister species, D. arizonae. We analyzed transcriptional responses in female heads using an integrated investigation of genome-wide patterns of gene expression, including differential expression (DE), alternative splicing (AS) and intron retention (IR).

RESULTS

Our results indicated that early transcriptional responses were largely congruent between con- and heterospecific matings but are substantially perturbed over time. Conspecific matings induced functional pathways related to amino acid balance previously associated with the brain's physiology and female postmating behavior. Heterospecific matings often failed to activate regulation of some of these genes and induced expression of additional genes when compared with those of conspecifically-mated females. These mechanisms showed functional specializations with DE genes mostly linked to pathways of proteolysis and nutrient homeostasis, while AS genes were more related to photoreception and muscle assembly pathways. IR seems to play a more general role in DE regulation during the female postmating response.

CONCLUSIONS

We provide evidence showing that AS genes substantially perturbed by heterospecific matings in female heads evolve at slower evolutionary rates than the genome background. However, DE genes evolve at evolutionary rates similar, or even higher, than those of male reproductive genes, which highlights their potential role in sexual selection and the evolution of reproductive barriers.

Pests, diseases, and aridity have shaped the genome of Corymbia citriodora

Corymbia citriodora is a member of the predominantly Southern Hemisphere Myrtaceae family, which includes the eucalypts (Eucalyptus, Corymbia and Angophora; ~800 species). Corymbia is grown for timber, pulp and paper, and essential oils in Australia, South Africa, Asia, and Brazil, maintaining a high-growth rate under marginal conditions due to drought, poor-quality soil, and biotic stresses. To dissect the genetic basis of these desirable traits, we sequenced and assembled the 408 Mb genome of Corymbia citriodora, anchored into eleven chromosomes. Comparative analysis with Eucalyptus grandis reveals high synteny, although the two diverged approximately 60 million years ago and have different genome sizes (408 vs 641 Mb), with few large intra-chromosomal rearrangements. C. citriodora shares an ancient whole-genome duplication event with E. grandis but has undergone tandem gene family expansions related to terpene biosynthesis, innate pathogen resistance, and leaf wax formation, enabling their successful adaptation to biotic/abiotic stresses and arid conditions of the Australian continent.

Osteocyte transcriptome mapping identifies a molecular landscape controlling skeletal homeostasis and susceptibility to skeletal disease

Osteocytes are master regulators of the skeleton. We mapped the transcriptome of osteocytes from different skeletal sites, across age and sexes in mice to reveal genes and molecular programs that control this complex cellular-network. We define an osteocyte transcriptome signature of 1239 genes that distinguishes osteocytes from other cells. 77% have no previously known role in the skeleton and are enriched for genes regulating neuronal network formation, suggesting this programme is important in osteocyte communication. We evaluated 19 skeletal parameters in 733 knockout mouse lines and reveal 26 osteocyte transcriptome signature genes that control bone structure and function. We showed osteocyte transcriptome signature genes are enriched for human orthologs that cause monogenic skeletal disorders (P = 2.4 × 10-22) and are associated with the polygenic diseases osteoporosis (P = 1.8 × 10-13) and osteoarthritis (P = 1.6 × 10-7). Thus, we reveal the molecular landscape that regulates osteocyte network formation and function and establish the importance of osteocytes in human skeletal disease.

The Corylus mandshurica genome provides insights into the evolution of Betulaceae genomes and hazelnut breeding

Hazelnut is popular for its flavor, and it has also been suggested that hazelnut is beneficial to cardiovascular health because it is rich in oleic acid. Here, we report the first high-quality chromosome-scale genome for the hazelnut species Corylus mandshurica (2n = 22), which has a high concentration of oleic acid in its nuts. The assembled genome is 367.67 Mb in length, and the contig N50 is 14.85 Mb. All contigs were assembled into 11 chromosomes, and 28,409 protein-coding genes were annotated. We reconstructed the evolutionary trajectories of the genomes of Betulaceae species and revealed that the 11 chromosomes of the hazelnut genus were derived from the most ancestral karyotype in Betula pendula, which has 14 protochromosomes, by inferring homology among five Betulaceae genomes. We identified 96 candidate genes involved in oleic acid biosynthesis, and 10 showed rapid evolution or positive selection. These findings will help us to understand the mechanisms of lipid synthesis and storage in hazelnuts. Several gene families related to salicylic acid metabolism and stress responses experienced rapid expansion in this hazelnut species, which may have increased its stress tolerance. The reference genome presented here constitutes a valuable resource for molecular breeding and genetic improvement of the important agronomic properties of hazelnut.

The mitochondrial aspartate/glutamate carrier (AGC or Aralar1) isoforms in D. melanogaster: biochemical characterization, gene structure, and evolutionary analysis

BACKGROUND

In man two mitochondrial aspartate/glutamate carrier (AGC) isoforms, known as aralar and citrin, are required to accomplish several metabolic pathways. In order to fill the existing gap of knowledge in Drosophila melanogaster, we have studied aralar1 gene, orthologue of human AGC-encoding genes in this organism.

METHODS

The blastp algorithm and the "reciprocal best hit" approach have been used to identify the human orthologue of AGCs in Drosophilidae and non-Drosophilidae. Aralar1 proteins have been overexpressed in Escherichia coli and functionally reconstituted into liposomes for transport assays.

RESULTS

The transcriptional organization of aralar1 comprises six isoforms, three constitutively expressed (aralar1-RA, RD and RF), and the remaining three distributed during the development or in different tissues (aralar1-RB, RC and RE). Aralar1-PA and Aralar1-PE, representative of all isoforms, have been biochemically characterized. Recombinant Aralar1-PA and Aralar1-PE proteins share similar efficiency to exchange glutamate against aspartate, and same substrate affinities than the human isoforms. Interestingly, although Aralar1-PA and Aralar1-PE diverge only in their EF-hand 8, they greatly differ in their specific activities and substrate specificity.

CONCLUSIONS

The tight regulation of aralar1 transcripts expression and the high request of aspartate and glutamate during early embryogenesis suggest a crucial role of Aralar1 in this Drosophila developmental stage. Furthermore, biochemical characterization and calcium sensitivity have identified Aralar1-PA and Aralar1-PE as the human aralar and citrin counterparts, respectively.

GENERAL SIGNIFICANCE

The functional characterization of the fruit fly mitochondrial AGC transporter represents a crucial step toward a complete understanding of the metabolic events acting during early embryogenesis.

Comprehensive analysis of NAC transcription factor family uncovers drought and salinity stress response in pearl millet (Pennisetum glaucum)

BACKGROUND

Pearl millet (Pennisetum glaucum) is a cereal crop that possesses the ability to withstand drought, salinity and high temperature stresses. The NAC [NAM (No Apical Meristem), ATAF1 (Arabidopsis thaliana Activation Factor 1), and CUC2 (Cup-shaped Cotyledon)] transcription factor family is one of the largest transcription factor families in plants. NAC family members are known to regulate plant growth and abiotic stress response. Currently, no reports are available on the functions of the NAC family in pearl millet.

RESULTS

Our genome-wide analysis found 151 NAC transcription factor genes (PgNACs) in the pearl millet genome. Thirty-eight and 76 PgNACs were found to be segmental and dispersed duplicated respectively. Phylogenetic analysis divided these NAC transcription factors into 11 groups (A-K). Three PgNACs (- 073, - 29, and - 151) were found to be membrane-associated transcription factors. Seventeen other conserved motifs were found in PgNACs. Based on the similarity of PgNACs to NAC proteins in other species, the functions of PgNACs were predicted. In total, 88 microRNA target sites were predicted in 59 PgNACs. A previously performed transcriptome analysis suggests that the expression of 30 and 42 PgNACs are affected by salinity stress and drought stress, respectively. The expression of 36 randomly selected PgNACs were examined by quantitative reverse transcription-PCR. Many of these genes showed diverse salt- and drought-responsive expression patterns in roots and leaves. These results confirm that PgNACs are potentially involved in regulating abiotic stress tolerance in pearl millet.

CONCLUSION

The pearl millet genome contains 151 NAC transcription factor genes that can be classified into 11 groups. Many of these genes are either upregulated or downregulated by either salinity or drought stress and may therefore contribute to establishing stress tolerance in pearl millet.

The RNA-binding protein Quaking regulates multiciliated and basal cell abundance in the developing lung

RNA-binding proteins (RBPs) form complexes with RNA, changing how the RNA is processed and thereby regulating gene expression. RBPs are important sources of gene regulation during organogenesis, including the development of lungs. The RBP called Quaking (QK) is critical for embryogenesis, yet it has not been studied in the developing lung. Here, we show that QK is widely expressed during rat lung development and into adulthood. The QK isoforms QK5 and QK7 colocalize to the nuclei of nearly all lung cells. QK6 is present in the nuclei and cytoplasm of mesenchymal cells and is only present in the epithelium during branching morphogenesis. QK knockdown in embryonic lung explants caused a greater number of multiciliated cells to appear in the airways, at the expense of basal cells. The mRNA of multiciliated cell genes and the abundance of FOXJ1/SOX2+ cells increased after knockdown, whereas P63/SOX2+ cells decreased. The cytokine IL-6, a known regulator of multiciliated cell differentiation, had increased mRNA levels after QK knockdown, although protein levels remained unchanged. Further studies are necessary to confirm whether QK acts as a blocker for the IL-6-induced differentiation of basal cells into multiciliated cells, and a conditional QK knockout would likely lead to additional discoveries on QK's role during lung development.

Exploring the synthetic biology potential of bacteriophages for engineering non-model bacteria

Non-model bacteria like Pseudomonas putida, Lactococcus lactis and other species have unique and versatile metabolisms, offering unique opportunities for Synthetic Biology (SynBio). However, key genome editing and recombineering tools require optimization and large-scale multiplexing to unlock the full SynBio potential of these bacteria. In addition, the limited availability of a set of characterized, species-specific biological parts hampers the construction of reliable genetic circuitry. Mining of currently available, diverse bacteriophages could complete the SynBio toolbox, as they constitute an unexplored treasure trove for fully adapted metabolic modulators and orthogonally-functioning parts, driven by the longstanding co-evolution between phage and host.

Drosophila melanogaster Mitochondrial Carriers: Similarities and Differences with the Human Carriers

Mitochondrial carriers are a family of structurally related proteins responsible for the exchange of metabolites, cofactors and nucleotides between the cytoplasm and mitochondrial matrix. The in silico analysis of the Drosophila melanogaster genome has highlighted the presence of 48 genes encoding putative mitochondrial carriers, but only 20 have been functionally characterized. Despite most Drosophila mitochondrial carrier genes having human homologs and sharing with them 50% or higher sequence identity, D. melanogaster genes display peculiar differences from their human counterparts: (1) in the fruit fly, many genes encode more transcript isoforms or are duplicated, resulting in the presence of numerous subfamilies in the genome; (2) the expression of the energy-producing genes in D. melanogaster is coordinated from a motif known as Nuclear Respiratory Gene (NRG), a palindromic 8-bp sequence; (3) fruit-fly duplicated genes encoding mitochondrial carriers show a testis-biased expression pattern, probably in order to keep a duplicate copy in the genome. Here, we review the main features, biological activities and role in the metabolism of the D. melanogaster mitochondrial carriers characterized to date, highlighting similarities and differences with their human counterparts. Such knowledge is very important for obtaining an integrated view of mitochondrial function in D. melanogaster metabolism.

Molecular and Vegetative Compatibility Groups Characterization of Aspergillus flavus Isolates from Kenya

The genus Aspergillus contains diverse species and the identification is complicated. Vegetative compatibility groups (VCGs) and molecular mechanisms were deployed to study the species. The study was randomly conducted in four counties in Kenya based on the history of aflatoxicosis and maize cultivation. Thirty-seven Aspergillus flavus isolates from Nandi, Kisumu, Homa Bay and Makueni were characterized to determine their taxonomic status based on their VCGs and genotypes. A phylogenetic analysis of ITS1 and ITS2 sequences of the isolates investigated revealed ITS primers discriminating some of the A. flavus isolates as 100% sequence identity to the RefSeq. Nit mutants' complementation test revealed strong heterokaryon incompatibility between isolates of Nandi region (67%) and Makueni (33%). The trend based on VCGs and molecular findings showed high incidence of toxigenic A. flavus in Makueni, which could be the reason why the region frequently experiences chronic aflatoxicosis incidences over the last few decades as compared to other regions. Interestingly, we have discovered all S and L-morphotypes including the rare S/L-morphotypes, which implies that Kenya is home to all morphotypes of A. flavus. Thus, the analysis provides a deeper understanding of the taxonomic relationship between the A. flavus isolates and could help contextualise the data obtained for each isolate with respect to VCG genetic diversity and genotypes. Determining the primary causal agents of aflatoxin contamination is critical for predicting risk of contamination events and designing and implementing effective management strategies.

Portal of Juglandaceae: A comprehensive platform for Juglandaceae study

Juglandaceae species are plants of great economic value and have been cultivated, domesticated, and utilized by human society for a long time. Their edible, nutrient-rich nuts and tough, durable wood have attracted the attention of botanists and breeders. With the advent of the genomics era, genome sequencing of the Juglandaceae family has been greatly accelerated, and a large amount of data has been generated. In this paper, we introduce the Portal of Juglandaceae (PJU), a tool to bring all these data together. The PJU contains genomes, gene-coding sequences, protein sequences, various types of annotation information, expression data, and miRNA data, which are configured with BLAST, JBrowse, and our self-developed synteny analysis tool. The PJU has a user-friendly and straightforward interface that performs a variety of query tasks with a few simple operations. In the future, we hope that the PJU will serve as a hub for the study of the Juglandaceae family.

Evaluation of Rationally Designed Label-free Stem-loop DNA Probe Opening in the Presence of miR-21 by Circular Dichroism and Fluorescence Techniques

The characteristic features of stem-loop structured probes make them robust tools to detect targets with high sensitivity and selectivity. The basis of the hairpin based sensors operation is a conformational change that occurs upon hybridization of target with stem-loop probe. The design of the stem-loop probe has an important role in target recognition. Therefore, we designed a label-free stem loop probe for targeting miR-21 as a cancer biomarker investigated by web-based tools; its thermodynamic parameters obtained by thermal UV spectroscopy. The efficiency of stem-loop structure opening in the presence of target and non-target sequences was evaluated by fluorescence spectroscopy and circular dichroism spectro-polarimetry. The results showed that the target sequence opens the structure of hairpin efficiently in comparison to non-target sequences. To optimize the stem-loop hybridization to its target, the buffer ionic strength was changed by adding different concentrations of NaCl, KCl and MgCl₂. It was shown that buffering conditions have a significant role in loop structure opening and its optimization, led to an increase in sensitivity detection and have improved LOD from 60 pM to 45 pM.

Contributions of alternative splicing to muscle type development and function

Animals possess a wide variety of muscle types that support different kinds of movements. Different muscles have distinct locations, morphologies and contractile properties, raising the question of how muscle diversity is generated during development. Normal aging processes and muscle disorders differentially affect particular muscle types, thus understanding how muscles normally develop and are maintained provides insight into alterations in disease and senescence. As muscle structure and basic developmental mechanisms are highly conserved, many important insights into disease mechanisms in humans as well as into basic principles of muscle development have come from model organisms such as Drosophila, zebrafish and mouse. While transcriptional regulation has been characterized to play an important role in myogenesis, there is a growing recognition of the contributions of alternative splicing to myogenesis and the refinement of muscle function. Here we review our current understanding of muscle type specific alternative splicing, using examples of isoforms with distinct functions from both vertebrates and Drosophila. Future exploration of the vast potential of alternative splicing to fine-tune muscle development and function will likely uncover novel mechanisms of isoform-specific regulation and a more holistic understanding of muscle development, disease and aging.

Genomic and transcriptomic insights into molecular basis of sexually dimorphic nuptial spines in Leptobrachium leishanense

Sexually dimorphic (SD) traits are important in sexual selection and species survival, yet the molecular basis remains elusive, especially in amphibians where SD traits have evolved repeatedly. We focus on the Leishan moustache toad (Leptobrachium leishanense), in which males develop nuptial spines on their maxillary skin. Here we report a 3.5 Gb genome assembly with a contig N50 of 1.93 Mb. We find a specific expansion of the intermediate filament gene family including numerous keratin genes. Within these genes, a cluster of duplicated hair keratin genes exhibits male-biased and maxillary skin-specific expression, suggesting a role in developing nuptial spines. We identify a module of coexpressed genes significantly associated with spine formation. In addition, we find several hormones likely to be involved in regulating spine development. This study not only presents a high-quality anuran genome but also provides a reference for studying skin-derived SD traits in amphibians.

The basis of sexual dimorphism in non-model species may be elusive, in part due to a lack of genomic and molecular resources. Here, Li et al. report a high-quality anuran genome and reveal candidate genes and pathways associated with shaping sexually dimorphic nuptial spines in a moustache toad.

The Combined Strategy for iron uptake is not exclusive to domesticated rice (Oryza sativa)

Iron (Fe) is an essential micronutrient that is frequently inaccessible to plants. Rice (Oryza sativa L.) plants employ the Combined Strategy for Fe uptake, which is composed by all features of Strategy II, common to all Poaceae species, and some features of Strategy I, common to non-Poaceae species. To understand the evolution of Fe uptake mechanisms, we analyzed the root transcriptomic response to Fe deficiency in O. sativa and its wild progenitor O. rufipogon. We identified 622 and 2,017 differentially expressed genes in O. sativa and O. rufipogon, respectively. Among the genes up-regulated in both species, we found Fe transporters associated with Strategy I, such as IRT1, IRT2 and NRAMP1; and genes associated with Strategy II, such as YSL15 and IRO2. In order to evaluate the conservation of these Strategies among other Poaceae, we identified the orthologs of these genes in nine species from the Oryza genus, maize and sorghum, and evaluated their expression profile in response to low Fe condition. Our results indicate that the Combined Strategy is not specific to O. sativa as previously proposed, but also present in species of the Oryza genus closely related to domesticated rice, and originated around the same time the AA genome lineage within Oryza diversified. Therefore, adaptation to Fe2+ acquisition via IRT1 in flooded soils precedes O. sativa domestication.

Drosophila Sister-of-Sex-lethal reinforces a male-specific gene expression pattern by controlling Sex-lethal alternative splicing

In Drosophila, female development is governed by a single RNA-binding protein, Sex-lethal (Sxl), that controls the expression of key factors involved in dosage compensation, germline homeostasis and the establishment of female morphology and behaviour. Sxl expression in female flies is maintained by an auto-regulatory, positive feedback loop with Sxl controlling splicing of its own mRNA. Until now, it remained unclear how males prevent accidental triggering of the Sxl expression cascade and protect themselves against runaway protein production. Here, we identify the protein Sister-of-Sex-lethal (Ssx) as an inhibitor of Sxl auto-regulatory splicing. Sxl and Ssx have a comparable RNA-binding specificity and compete for binding to RNA regulatory elements present in the Sxl transcript. In cultured Drosophila cells, Sxl-induced changes to alternative splicing can be reverted by the expression of Ssx. Moreover, in adult male flies ablation of the ssx gene results in a low level of productive Sxl mRNA splicing and Sxl protein production in isolated, clonal cell populations. In sum, this demonstrates that Ssx safeguards male animals against Sxl protein production to reinforce a stable, male-specific gene expression pattern.

Novel insights into PARPs in gene expression: regulation of RNA metabolism

Poly(ADP-ribosyl)ation (PARylation) is an important post-translational modification in which an ADP-ribose group is transferred to the target protein by poly(ADP-riboses) polymerases (PARPs). Since the discovery of poly-ADP-ribose (PAR) 50 years ago, its roles in cellular processes have been extensively explored. Although research initially focused on the functions of PAR and PARPs in DNA damage detection and repair, our understanding of the roles of PARPs in various nuclear and cytoplasmic processes, particularly in gene expression, has increased significantly. In this review, we discuss the current advances in understanding the roles of PARylation with a particular emphasis in gene expression through RNA biogenesis and processing. In addition to updating PARP's significance in transcriptional regulation, we specifically focus on how PARPs and PARylation affect gene expression, especially inflammation-related genes, at the post-transcriptional levels by modulating RNA processing and degrading. Increasing evidence suggests that PARP inhibition is a promising treatment for inflammation-related diseases besides conventional chemotherapy for cancer.

Evolutionary characteristics of intergenic transcribed regions indicate rare novel genes and widespread noisy transcription in the Poaceae

Extensive transcriptional activity occurring in intergenic regions of genomes has raised the question whether intergenic transcription represents the activity of novel genes or noisy expression. To address this, we evaluated cross-species and post-duplication sequence and expression conservation of intergenic transcribed regions (ITRs) in four Poaceae species. Among 43,301 ITRs across the four species, 34,460 (80%) are species-specific. ITRs found across species tend to be more divergent in expression and have more recent duplicates compared to annotated genes. To assess if ITRs are functional (under selection), machine learning models were established in Oryza sativa (rice) that could accurately distinguish between phenotype genes and pseudogenes (area under curve-receiver operating characteristic = 0.94). Based on the models, 584 (8%) and 4391 (61%) rice ITRs are classified as likely functional and nonfunctional with high confidence, respectively. ITRs with conserved expression and ancient retained duplicates, features that were not part of the model, are frequently classified as likely-functional, suggesting these characteristics could serve as pragmatic rules of thumb for identifying candidate sequences likely to be under selection. This study also provides a framework to identify novel genes using comparative transcriptomic data to improve genome annotation that is fundamental for connecting genotype to phenotype in crop and model systems.

Genome-wide identification and expression profiling of the C2H2-type zinc finger protein genes in the silkworm Bombyx mori

Cys2-His2 zinc finger (C2H2-ZF) proteins comprise the largest class of putative eukaryotic transcription factors. The zinc finger motif array is highly divergent, indicating that most proteins will have distinctive binding sites and perform different functions. However, the binding sites and functions of the majority of C2H2-ZF proteins remain unknown. In this study, we identified 327 C2H2-ZF protein genes in the silkworm, 290 in the monarch butterfly, 243 in the fruit fly, 107 in elegans, 673 in mouse, and 1,082 in human. The C2H2-ZF protein genes of the silkworm were classified into three main grouping clades according to a phylogenetic classification, and 312 of these genes could be mapped onto 27 chromosomes. Most silkworm C2H2-ZF protein genes exhibited specific expression in larval tissues. Furthermore, several C2H2-ZF protein genes had sex-specific expression during metamorphosis. In addition, we found that some C2H2-ZF protein genes are involved in metamorphosis and female reproduction by using expression clustering and gene annotation analysis. Among them, five genes were selected, BGIBMGA002091 (CTCF), BGIBMGA006492 (fru), BGIBMGA006230 (wor), BGIBMGA004640 (lola), and BIGBMGA004569, for quantitative real-time PCR analysis from larvae to adult ovaries. The results showed that the five genes had different expression patterns in ovaries, among which BGIBMGA002091 (CTCF) gene expression level was the highest, and its expression level increased rapidly in late pupae and adult stages. These findings provide a basis for further investigation of the functions of C2H2-ZF protein genes in the silkworm, and the results offer clues for further research into the development of metamorphosis and female reproduction in the silkworm.

Analysis of differential gene expression and alternative splicing is significantly influenced by choice of reference genome

RNA-seq analysis has enabled the evaluation of transcriptional changes in many species including nonmodel organisms. However, in most species only a single reference genome is available and RNA-seq reads from highly divergent varieties are typically aligned to this reference. Here, we quantify the impacts of the choice of mapping genome in rice where three high-quality reference genomes are available. We aligned RNA-seq data from a popular productive rice variety to three different reference genomes and found that the identification of differentially expressed genes differed depending on which reference genome was used for mapping. Furthermore, the ability to detect differentially used transcript isoforms was profoundly affected by the choice of reference genome: Only 30% of the differentially used splicing features were detected when reads were mapped to the more commonly used, but more distantly related reference genome. This demonstrated that gene expression and splicing analysis varies considerably depending on the mapping reference genome, and that analysis of individuals that are distantly related to an available reference genome may be improved by acquisition of new genomic reference material. We observed that these differences in transcriptome analysis are, in part, due to the presence of single nucleotide polymorphisms between the sequenced individual and each respective reference genome, as well as annotation differences between the reference genomes that exist even between syntenic orthologs. We conclude that even between two closely related genomes of similar quality, using the reference genome that is most closely related to the species being sampled significantly improves transcriptome analysis.

Conserved functions of RNA-binding proteins in muscle

Animals require different types of muscle for survival, for example for circulation, motility, reproduction and digestion. Much emphasis in the muscle field has been placed on understanding how transcriptional regulation generates diverse types of muscle during development. Recent work indicates that alternative splicing and RNA regulation are as critical to muscle development, and altered function of RNA-binding proteins causes muscle disease. Although hundreds of genes predicted to bind RNA are expressed in muscles, many fewer have been functionally characterized. We present a cross-species view summarizing what is known about RNA-binding protein function in muscle, from worms and flies to zebrafish, mice and humans. In particular, we focus on alternative splicing regulated by the CELF, MBNL and RBFOX families of proteins. We discuss the systemic nature of diseases associated with loss of RNA-binding proteins in muscle, focusing on mis-regulation of CELF and MBNL in myotonic dystrophy. These examples illustrate the conservation of RNA-binding protein function and the marked utility of genetic model systems in understanding mechanisms of RNA regulation.

Comparative ploidy response to experimental hydrogen peroxide exposure in Atlantic salmon (Salmo salar)

While research into the growth, survival, nutrition and, more recently, disease susceptibility of triploid Atlantic salmon has expanded, there remains an overall lack of studies assessing the response of triploids to chemical treatments. It is essential that the response of triploids to disease treatments be characterised to validate their suitability for commercial production. This study aimed to investigate and compare the stress and immune responses of triploid and diploid Atlantic salmon following an experimental treatment with hydrogen peroxide (H₂O₂). A dose response test was first undertaken to determine a suitable test dose for both diploid and triploid Atlantic salmon. Following this, diploids and triploids were exposed to H₂O₂ (1800 ppm) for 20 min, as per commercial practices, after which blood glucose and lactate, and plasma cortisol and lysozyme were measured, along with the expression of oxidative stress and immune-related genes. In the first 6 h post-exposure to H₂O₂, comparable mortalities occurred in both diploid and triploid Atlantic salmon. Cortisol, glucose and lactate were not significantly influenced by ploidy suggesting that, physiologically, triploid Atlantic salmon are able to cope with the stress associated with H₂O₂ exposure as well as their diploid counterparts. Exposure to H₂O₂ significantly elevated the expression of cat and sod2 in diploid livers and gr, il1β and crp/sap1b in diploid gills, while it significantly decreased the expression of saa5 and crp/sap1a in diploid gills. In triploids, the expression levels of cat, hsp70, sod1, saa5, crp/sap1a and crp/sap1b in liver was significantly higher in fish exposed to H₂O₂ compared to control fish. The expression of gr, sod1 and il1β in triploid gills was also elevated in response to H₂O₂ exposure. This study represents the first experimental evidence of the effects of H₂O₂ exposure on triploid Atlantic salmon and continues to support their application into commercial production.

Specification of Drosophila neuropeptidergic neurons by the splicing component brr2

During embryonic development, a number of genetic cues act to generate neuronal diversity. While intrinsic transcriptional cascades are well-known to control neuronal sub-type cell fate, the target cells can also provide critical input to specific neuronal cell fates. Such signals, denoted retrograde signals, are known to provide critical survival cues for neurons, but have also been found to trigger terminal differentiation of neurons. One salient example of such target-derived instructive signals pertains to the specification of the Drosophila FMRFamide neuropeptide neurons, the Tv4 neurons of the ventral nerve cord. Tv4 neurons receive a BMP signal from their target cells, which acts as the final trigger to activate the FMRFa gene. A recent FMRFa-eGFP genetic screen identified several genes involved in Tv4 specification, two of which encode components of the U5 subunit of the spliceosome: brr2 (l(3)72Ab) and Prp8. In this study, we focus on the role of RNA processing during target-derived signaling. We found that brr2 and Prp8 play crucial roles in controlling the expression of the FMRFa neuropeptide specifically in six neurons of the VNC (Tv4 neurons). Detailed analysis of brr2 revealed that this control is executed by two independent mechanisms, both of which are required for the activation of the BMP retrograde signaling pathway in Tv4 neurons: (1) Proper axonal pathfinding to the target tissue in order to receive the BMP ligand. (2) Proper RNA splicing of two genes in the BMP pathway: the thickveins (tkv) gene, encoding a BMP receptor subunit, and the Medea gene, encoding a co-Smad. These results reveal involvement of specific RNA processing in diversifying neuronal identity within the central nervous system.

The nervous system displays daunting cellular diversity, largely generated through complex regulatory input operating on stem cells and their neural lineages during development. Most of the reported mechanisms acting to generate neural diversity pertain to transcriptional regulation. In contrast, little is known regarding the post-transcriptional mechanisms involved. Here, we use a specific group of neurons, Apterous neurons, in the ventral nerve cord of Drosophila melanogaster as our model, to analyze the function of two essential components of the spliceosome; Brr2 and Prp8. Apterous neurons require a BMP retrograde signal for terminal differentiation, and we find that brr2 and Prp8 play crucial roles during this process. brr2 is critical for two independent events; axon pathfinding and BMP signaling, both of which are required for the activation of the retrograde signaling pathway necessary for Apterous neurons. These results identify a post-transcriptional mechanism as key for specifying neuronal identity, by ensuring the execution of a retrograde signal.

Successive duplication-divergence mechanisms at the RCO locus contributed to leaf shape diversity in the Brassicaceae

Gene duplication is a major driver for the increase of biological complexity. The divergence of newly duplicated paralogs may allow novel functions to evolve, while maintaining the ancestral one. Alternatively, partitioning the ancestral function among paralogs may allow parts of that role to follow independent evolutionary trajectories. We studied the REDUCED COMPLEXITY (RCO) locus, which contains three paralogs that have evolved through two independent events of gene duplication, and which underlies repeated events of leaf shape evolution within the Brassicaceae. In particular, we took advantage of the presence of three potentially functional paralogs in Capsella to investigate the extent of functional divergence among them. We demonstrate that the RCO copies control growth in different areas of the leaf. Consequently, the copies that are retained active in the different Brassicaceae lineages contribute to define the leaf dissection pattern. Our results further illustrate how successive gene duplication events and subsequent functional divergence can increase trait evolvability by providing independent evolutionary trajectories to specialized functions that have an additive effect on a given trait.

In silico identification and experimental validation of amino acid motifs required for the Rho-of-plants GTPase-mediated activation of receptor-like cytoplasmic kinases

Several amino acid motifs required for Rop-dependent activity were found to form a common surface on RLCKVI_A kinases. This indicates a unique mechanism for Rho-type GTPase-mediated kinase activation in plants. Rho-of-plants (Rop) G-proteins are implicated in the regulation of various cellular processes, including cell growth, cell polarity, hormonal and pathogen responses. Our knowledge about the signalling pathways downstream of Rops is continuously increasing. However, there are still substantial gaps in this knowledge. One reason for this is that these pathways are considerably different from those described for yeast and/or animal Rho-type GTPases. Among others, plants lack all Rho/Rac/Cdc42-activated kinase families. Only a small group of plant-specific receptor-like cytoplasmic kinases (RLCK VI_A) has been shown to exhibit Rop-binding-dependent in vitro activity. These kinases do not carry any known GTPase-binding motifs. Based on the sequence comparison of the Rop-activated RLCK VI_A and the closely related but constitutively active RLCK VI_B kinases, several distinguishing amino acid residues/motifs were identified. All but one of these were found to be required for the Rop-mediated regulation of the in vitro activity of two RLCK VI_A kinases. Structural modelling indicated that these motifs might form a common Rop-binding surface. Based on in silico data mining, kinases that have the identified Rop-binding motifs are present in Embryophyta but not in unicellular green algae. It can, therefore, be supposed that Rops recruited these plant-specific kinases for signalling at an early stage of land plant evolution.

Drosophila Sister-of-Sex-lethal is a repressor of translation

The RNA-binding protein Sex-lethal (Sxl) is an important post-transcriptional regulator of sex determination and dosage compensation in female Drosophila To prevent the assembly of the MSL dosage compensation complex in female flies, Sxl acts as a repressor of male-specific lethal-2 (msl-2) mRNA translation. It uses two distinct and mutually reinforcing blocks to translation that operate on the 5' and 3' untranslated regions (UTRs) of msl-2 mRNA, respectively. While 5' UTR-mediated translational control involves an upstream open reading frame, 3' UTR-mediated regulation strictly requires the co-repressor protein Upstream of N-ras (Unr), which is recruited to the transcript by Sxl. We have identified the protein Sister-of-Sex-lethal (Ssx) as a novel repressor of translation with Sxl-like activity. Both proteins have a comparable RNA-binding specificity and can associate with uracil-rich RNA regulatory elements present in msl-2 mRNA. Moreover, both repress translation when bound to the 5' UTR of msl-2 However, Ssx is inactive in 3' UTR-mediated regulation, as it cannot engage the co-repressor protein Unr. The difference in activity maps to the first RNA-recognition motif (RRM) of Ssx. Conversion of three amino acids within this domain into their Sxl counterpart results in a gain of function and repression via the 3' UTR, allowing detailed insights into the evolutionary origin of the two proteins and into the molecular requirements of an important translation regulatory pathway.

Inferring synteny between genome assemblies: a systematic evaluation

BACKGROUND

Genome assemblies across all domains of life are being produced routinely. Initial analysis of a new genome usually includes annotation and comparative genomics. Synteny provides a framework in which conservation of homologous genes and gene order is identified between genomes of different species. The availability of human and mouse genomes paved the way for algorithm development in large-scale synteny mapping, which eventually became an integral part of comparative genomics. Synteny analysis is regularly performed on assembled sequences that are fragmented, neglecting the fact that most methods were developed using complete genomes. It is unknown to what extent draft assemblies lead to errors in such analysis.

RESULTS

We fragmented genome assemblies of model nematodes to various extents and conducted synteny identification and downstream analysis. We first show that synteny between species can be underestimated up to 40% and find disagreements between popular tools that infer synteny blocks. This inconsistency and further demonstration of erroneous gene ontology enrichment tests raise questions about the robustness of previous synteny analysis when gold standard genome sequences remain limited. In addition, assembly scaffolding using a reference guided approach with a closely related species may result in chimeric scaffolds with inflated assembly metrics if a true evolutionary relationship was overlooked. Annotation quality, however, has minimal effect on synteny if the assembled genome is highly contiguous.

CONCLUSIONS

Our results show that a minimum N50 of 1 Mb is required for robust downstream synteny analysis, which emphasizes the importance of gold standard genomes to the science community, and should be achieved given the current progress in sequencing technology.

Genetic Architecture of Natural Variation Underlying Adult Foraging Behavior That Is Essential for Survival of Drosophila melanogaster

Foraging behavior is critical for the fitness of individuals. However, the genetic basis of variation in foraging behavior and the evolutionary forces underlying such natural variation have rarely been investigated. We developed a systematic approach to assay the variation in survival rate in a foraging environment for adult flies derived from a wild Drosophila melanogaster population. Despite being such an essential trait, there is substantial variation of foraging behavior among D. melanogaster strains. Importantly, we provided the first evaluation of the potential caveats of using inbred Drosophila strains to perform genome-wide association studies on life-history traits, and concluded that inbreeding depression is unlikely a major contributor for the observed large variation in adult foraging behavior. We found that adult foraging behavior has a strong genetic component and, unlike larval foraging behavior, depends on multiple loci. Identified candidate genes are enriched in those with high expression in adult heads and, demonstrated by expression knock down assay, are involved in maintaining normal functions of the nervous system. Our study not only identified candidate genes for foraging behavior that is relevant to individual fitness, but also shed light on the initial stage underlying the evolution of the behavior.

Genome Editing: CRISPR-Cas9

In the present chapter, we present the protocols and guidelines to facilitate implementation of CRISPR-Cas9 technology in fungi where few or no genetic tools are in place. Hence, we firstly explain how to identify dominant markers for genetic transformation. Secondly, we provide a guide for construction of Cas9/sgRNA episomal expression vectors. Thirdly, we present how to mutagenize reporter genes to explore the efficiency of CRISPR-Cas9 in the relevant fungus and to ease subsequent CRISPR-mediated genetic engineering. Lastly, we describe how to make CRISPR-mediated marker-dependent and marker-free gene targeting.

The Y Chromosome Modulates Splicing and Sex-Biased Intron Retention Rates in Drosophila

The Drosophila Y chromosome is a 40-Mb segment of mostly repetitive DNA; it harbors a handful of protein-coding genes and a disproportionate amount of satellite repeats, transposable elements, and multicopy DNA arrays. Intron retention (IR) is a type of alternative splicing (AS) event by which one or more introns remain within the mature transcript. IR recently emerged as a deliberate cellular mechanism to modulate gene expression levels and has been implicated in multiple biological processes. However, the extent of sex differences in IR and the contribution of the Y chromosome to the modulation of AS and IR rates has not been addressed. Here we showed pervasive IR in the fruit fly Drosophila melanogaster with thousands of novel IR events, hundreds of which displayed extensive sex bias. The data also revealed an unsuspected role for the Y chromosome in the modulation of AS and IR. The majority of sex-biased IR events introduced premature termination codons and the magnitude of sex bias was associated with gene expression differences between the sexes. Surprisingly, an extra Y chromosome in males (X^YY genotype) or the presence of a Y chromosome in females (X^XY genotype) significantly modulated IR and recapitulated natural differences in IR between the sexes. Our results highlight the significance of sex-biased IR in tuning sex differences and the role of the Y chromosome as a source of variable IR rates between the sexes. Modulation of splicing and IR rates across the genome represent new and unexpected outcomes of the Drosophila Y chromosome.