Decoding the genome of bloodsucking midge Forcipomyia taiwana (Diptera: Ceratopogonidae): Insights into odorant receptor expansion

Biting midges, notably those within the Ceratopogonidae family, have long been recognized for their epidemiological significance, both as nuisances and vectors for disease transmission in vertebrates. Despite their impact, genomic insights into these insects, particularly beyond the Culicoides genus, remain limited. In this study, we assembled the Forcipomyia taiwana (Shiraki) genome, comprising 113 scaffolds covering 130.4 Mbps-with the longest scaffold reaching 7.6 Mbps and an N50 value of 2.6 Mbps-marking a pivotal advancement in understanding the genetic architecture of ceratopogonid biting midges. Phylogenomic analyses reveal a shared ancestry between F. taiwana and Culicoides sonorensis Wirth & Jones, dating back approximately 124 million years, and highlight a dynamic history of gene family expansions and contractions within the Ceratopogonidae family. Notably, a substantial expansion of the odorant receptor (OR) gene family was observed, which is crucial for the chemosensory capabilities that govern biting midges' interactions with their environment, including host seeking and oviposition behaviors. The distribution of OR genes across the F. taiwana genome displays notable clusters on scaffolds, indicating localized tandem gene duplication events. Additionally, several collinear regions were identified, hinting at segmental duplications, inversions, and translocations, contributing to the olfactory system's evolutionary complexity. Among the 156 ORs identified in F. taiwana, 134 are biting midge-specific ORs, distributed across three distinct clades, each exhibiting unique motif features that distinguish them from the others. Through weighted gene co-expression network analysis, we correlated distinct gene modules with sex and reproductive status, laying the groundwork for future investigations into the interplay between gene expression and adaptive behaviors in F. taiwana. In conclusion, our study not only highlights the unique olfactory repertoire of ceratopogonid biting midges but also sets the stage for future studies into the genetic underpinnings of their unique biological traits and ecological strategies.

Highly efficient capture approach for the identification of diverse inherited retinal disorders

Our study presents a 319-gene panel targeting inherited retinal dystrophy (IRD) genes. Through a multi-center retrospective cohort study, we validated the assay's effectiveness and clinical utility and characterized the mutation spectrum of Taiwanese IRD patients. Between January 2018 and May 2022, 493 patients in 425 unrelated families, all initially suspected of having IRD without prior genetic diagnoses, underwent detailed ophthalmic and physical examinations (with extra-ocular features recorded) and genetic testing with our customized panel. Disease-causing variants were identified by segregation analysis and clinical interpretation, with validation via Sanger sequencing. We achieved a read depth of >200× for 94.2% of the targeted 1.2 Mb region. 68.5% (291/425) of the probands received molecular diagnoses, with 53.9% (229/425) resolved cases. Retinitis pigmentosa (RP) is the most prevalent initial clinical impression (64.2%), and 90.8% of the cohort have the five most prevalent phenotypes (RP, cone-rod syndrome, Usher's syndrome, Leber's congenital amaurosis, Bietti crystalline dystrophy). The most commonly mutated genes of probands that received molecular diagnosis are USH2A (13.7% of the cohort), EYS (11.3%), CYP4V2 (4.8%), ABCA4 (4.5%), RPGR (3.4%), and RP1 (3.1%), collectively accounted for 40.8% of diagnoses. We identify 87 unique unreported variants previously not associated with IRD and refine clinical diagnoses for 21 patients (7.22% of positive cases). We developed a customized gene panel and tested it on the largest Taiwanese cohort, showing that it provides excellent coverage for diverse IRD phenotypes.

Extensive sampling of Saccharomyces cerevisiae in Taiwan reveals ecology and evolution of predomesticated lineages

The ecology and genetic diversity of the model yeast Saccharomyces cerevisiae before human domestication remain poorly understood. Taiwan is regarded as part of this yeast's geographic birthplace, where the most divergent natural lineage was discovered. Here, we extensively sampled the broadleaf forests across this continental island to probe the ancestral species’ diversity. We found that S. cerevisiae is distributed ubiquitously at low abundance in the forests. Whole-genome sequencing of 121 isolates revealed nine distinct lineages that diverged from Asian lineages during the Pleistocene, when a transient continental shelf land bridge connected Taiwan to other major landmasses. Three lineages are endemic to Taiwan and six are widespread in Asia, making this region a focal biodiversity hotspot. Both ancient and recent admixture events were detected between the natural lineages, and a genetic ancestry component associated with isolates from fruits was detected in most admixed isolates. Collectively, Taiwanese isolates harbor genetic diversity comparable to that of the whole Asia continent, and different lineages have coexisted at a fine spatial scale even on the same tree. Patterns of variations within each lineage revealed that S. cerevisiae is highly clonal and predominantly reproduces asexually in nature. We identified different selection patterns shaping the coding sequences of natural lineages and found fewer gene family expansion and contractions that contrast with domesticated lineages. This study establishes that S. cerevisiae has rich natural diversity sheltered from human influences, making it a powerful model system in microbial ecology.

Comparative Genomics of Three Colletotrichum scovillei Strains and Genetic Analysis Revealed Genes Involved in Fungal Growth and Virulence on Chili Pepper

Colletotrichum scovillei causes anthracnose of chili pepper in many countries. Three strains of this pathogen, Coll-524, Coll-153, and Coll-365, show varied virulence on chili pepper. Among the three strains, Coll-365 showed significant defects in growth and virulence. To decipher the genetic variations among these strains and identify genes contributing to growth and virulence, comparative genomic analysis and gene transformation to show gene function were applied in this study. Compared to Coll-524, Coll-153, and Coll-365 had numerous gene losses including 32 candidate effector genes that are mainly exist in acutatum species complex. A cluster of 14 genes in a 34-kb genomic fragment was lost in Coll-365. Through gene transformation, three genes in the 34-kb fragment were identified to have functions in growth and/or virulence of C. scovillei. CsPLAA encoding a phospholipase A2-activating protein enhanced the growth of Coll-365. A combination of CsPLAA with one transcription factor CsBZTF and one C6 zinc finger domain-containing protein CsCZCP was found to enhance the pathogenicity of Coll-365. Introduction of CsGIP, which encodes a hypothetical protein, into Coll-365 caused a reduction in the germination rate of Coll-365. In conclusion, the highest virulent strain Coll-524 had more genes and encoded more pathogenicity related proteins and transposable elements than the other two strains, which may contribute to the high virulence of Coll-524. In addition, the absence of the 34-kb fragment plays a critical role in the defects of growth and virulence of strain Coll-365.

Whole-Genome Sequence Resource of Calonectria ilicicola, the Casual Pathogen of Soybean Red Crown Rot

Calonectria ilicicola (anamorph: Cylindrocladium parasiticum) is a soilborne plant-pathogenic fungus with a broad host range, and it can cause red crown rot of soybean and Cylindrocladium black rot of peanut, which has become an emerging threat to crop production worldwide. Limited molecular studies have focused on Calonectria ilicicola and one of the possible difficulties is the lack of genomic resources. This study presents the first high quality and near-completed genome of C. ilicicola, using the Oxford Nanopore GridION sequencing platform. A total of 16 contigs were assembled and the genome of C. ilicicola isolate F018 was estimated to have 11 chromosomes. Currently, the C. ilicicola F018 genome represents the most contiguous assembly, which has the lowest contig number and the highest contig N50 among all Calonectria genome resources. Putative protein-coding sequences and secretory proteins were estimated to be 17,308 and 1,930 in the C. ilicicola F018 genome, respectively; and the prediction was close to other plant-pathogenic fungi, such as Fusarium species, within the Nectriaceae family. The availability of this high-quality genome resource is expected to facilitate research on fungal biology and genetics of C. ilicicola and to support advanced understanding of pathogen virulence and disease management.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Whole genome duplication facilitated the evolution of C4 photosynthesis in Gynandropsis gynandra

In higher plants, whole-genome duplication (WGD) is thought to facilitate the evolution of C₄ photosynthesis from C₃ photosynthesis. To understand this issue, we used new and existing leaf-development transcriptomes to construct two coding sequence databases for C₄Gynandropsis gynandra and C₃Tarenaya hassleriana, which shared a WGD before their divergence. We compared duplicated genes in the two species and found that the WGD contributed to four aspects of the evolution of C₄ photosynthesis in G. gynandra. First, G. gynandra has retained the duplicates of ALAAT (alanine aminotransferase) and GOGAT (glutamine oxoglutarate aminotransferase) for nitrogen recycling to establish a photorespiratory CO₂ pump in bundle sheath (BS) cells for increasing photosynthesis efficiency, suggesting that G. gynandra experienced a C₃–C₄ intermediate stage during the C₄ evolution. Second, G. gynandra has retained almost all known vein-development-related paralogous genes derived from the WGD event, likely contributing to the high vein complexity of G. gynandra. Third, the WGD facilitated the evolution of C₄ enzyme genes and their recruitment into the C₄ pathway. Fourth, several genes encoding photosystem I proteins were derived from the WGD and are upregulated in G. gynandra, likely enabling the NADH dehydrogenase-like complex to produce extra ATPs for the C₄ CO₂ concentration mechanism. Thus, the WGD apparently played an enabler role in the evolution of C₄ photosynthesis in G. gynandra. Importantly, an ALAAT duplicate became highly expressed in BS cells in G. gynandra, facilitating nitrogen recycling and transition to the C₄ cycle. This study revealed how WDG may facilitate C₄ photosynthesis evolution.

Chromosomal-level genome assembly of the semi-dwarf rice Taichung Native 1, an initiator of Green Revolution

Here we report the 409.5 Mb chromosome-level assembly of the first bred semi-dwarf rice, the Taichung Native 1 (TN1), which served as the template for the development of the Green Revolution (GR) cultivar IR8 "miracle rice". We sequenced the TN1 genome utilizing multiple platforms and produced PacBio long reads, Illumina paired-end reads, Illumina mate-pair reads and 10x Genomics linked reads. We used a hybrid approach to assemble the 226× coverage of sequences by a combination of de novo and reference-guided approaches. The assembled TN1 genome has an N50 scaffold size of 33.1 Mb with the longest measuring 45.5 Mb. We annotated 37,526 genes, in which 24,102 (64.23%) were assigned Blast2GO annotations. The genome has 4672 or 95.4% complete BUSCOs and a repeat content of 51.52%. We developed our own method of creating a GR pangenome using the orthologous relationships of the proteins of TN1, IR8, MH63 and IR64, identifying 16,999 core orthologue groups of Green Revolution. From the pangenome, we identified a set of shared and unique gene ontology terms for the accessory clusters, characterizing TN1, IR8, MH63 and IR64. This TN1 genome assembly and GR pangenome will be a resource for new genomic discoveries about Green Revolution, and for improving the disease and insect resistances and the yield of rice.

Prognostic Genetic Biomarkers Based on Oncogenic Signaling Pathways for Outcome Prediction in Patients with Oral Cavity Squamous Cell Carcinoma

Simple Summary

A comprehensive analysis based on mutational signatures and oncogenic signaling pathways to identify a specific subgroup of patients that had a significantly negative impact on both disease-free and overall survival in oral cavity squamous cell carcinoma (OCSCC) from whole exome-sequencing data. This analysis has revealed a variety of biologically relevant candidate target genes. Thirty percent of 165 tumors had multiple targetable alterations in multiple pathways. This suggests the complex interplay and crosstalk of oncogenic signaling pathways play an important role on the outcome of patients with OCSCC, and the candidate genes and pathways identified may include prognostic genetic biomarkers or therapeutic targets for OCSCC.

Abstract

Mutational profiling of patients’ tumors has suggested that the development of oral cavity squamous cell carcinoma (OCSCC) is driven by multiple genes in multiple pathways. This study aimed to examine the association between genomic alterations and clinical outcomes in patients with advanced stages OCSCC to facilitate prognostic stratification. We re-analyzed our previous whole-exome sequencing data from 165 long-term follow-ups of stages III and IV patients with OCSCC. Their frequent mutations were mapped to 10 oncogenic signaling pathways. Clinicopathological risk factors, relapse, and survival were analyzed to identify the genetic factors associated with advanced OCSCC. Frequent genetic alterations included point mutations in TP53, FAT1, NOTCH1, CASP8, CDKN2A, HRAS, PIK3CA, KMT2B (also known as MLL4), and LINC00273; amplified segments in CCND1, EGFR, CTTN, and FGFR1; and lost segments in CDKN2A, ADAM3A, and CFHR1/CFHR4. Comprehensive analysis of genetic alterations revealed that subgroups based on mutational signatures had a significant negative impact on disease-free survival (p = 0.0005) and overall survival (p = 0.0024). Several important signaling pathways were identified to be frequently genetically altered in our cohort. A specific subgroup of patients with alterations in NOTCH, RTK/RAS/MAPK, and TGF-beta pathways that had a significantly negative impact on disease-free survival (p = 0.0009). Thirty percent of samples had multiple targetable mutations in multiple pathways, indicating opportunities for novel therapy.

Molecular asymmetry in the cephalochordate embryo revealed by single-blastomere transcriptome profiling

Studies in various animals have shown that asymmetrically localized maternal transcripts play important roles in axial patterning and cell fate specification in early embryos. However, comprehensive analyses of the maternal transcriptomes with spatial information are scarce and limited to a handful of model organisms. In cephalochordates (amphioxus), an early branching chordate group, maternal transcripts of germline determinants form a compact granule that is inherited by a single blastomere during cleavage stages. Further blastomere separation experiments suggest that other transcripts associated with the granule are likely responsible for organizing the posterior structure in amphioxus; however, the identities of these determinants remain unknown. In this study, we used high-throughput RNA sequencing of separated blastomeres to examine asymmetrically localized transcripts in two-cell and eight-cell stage embryos of the amphioxus Branchiostoma floridae. We identified 111 and 391 differentially enriched transcripts at the 2-cell stage and the 8-cell stage, respectively, and used in situ hybridization to validate the spatial distribution patterns for a subset of these transcripts. The identified transcripts could be categorized into two major groups: (1) vegetal tier/germ granule-enriched and (2) animal tier/anterior-enriched transcripts. Using zebrafish as a surrogate model system, we showed that overexpression of one animal tier/anterior-localized amphioxus transcript, zfp665, causes a dorsalization/anteriorization phenotype in zebrafish embryos by downregulating the expression of the ventral gene, eve1, suggesting a potential function of zfp665 in early axial patterning. Our results provide a global transcriptomic blueprint for early-stage amphioxus embryos. This dataset represents a rich platform to guide future characterization of molecular players in early amphioxus development and to elucidate conservation and divergence of developmental programs during chordate evolution.

Analysis of genomic distributions of SARS-CoV-2 reveals a dominant strain type with strong allelic associations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causal agent of COVID 19, continues to evolve since its first emergence in December 2019. Using the complete sequences of 1,932 SARS-CoV-2 genomes, various clustering analyses consistently identified six types of the strains. Independent of the dendrogram construction, 13 signature variations in the form of single nucleotide variations (SNVs) in protein coding regions and one SNV in the 5' untranslated region (UTR) were identified and provided a direct interpretation for the six types (types I to VI). The six types of the strains and their underlying signature SNVs were validated in two subsequent analyses of 6,228 and 38,248 SARS-CoV-2 genomes which became available later. To date, type VI, characterized by the four signature SNVs C241T (5'UTR), C3037T (nsp3 F924F), C14408T (nsp12 P4715L), and A23403G (Spike D614G), with strong allelic associations, has become the dominant type. Since C241T is in the 5' UTR with uncertain significance and the characteristics can be captured by the other three strongly associated SNVs, we focus on the other three. The increasing frequency of the type VI haplotype 3037T-14408T-23403G in the majority of the submitted samples in various countries suggests a possible fitness gain conferred by the type VI signature SNVs. The fact that strains missing one or two of these signature SNVs fail to persist implies possible interactions among these SNVs. Later SNVs such as G28881A, G28882A, and G28883C have emerged with strong allelic associations, forming new subtypes. This study suggests that SNVs may become an important consideration in SARS-CoV-2 classification and surveillance.

Genomic Analysis of Mycobacterium tuberculosis Isolates and Construction of a Beijing Lineage Reference Genome

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis, kills over 1 million people worldwide annually. Development of drug resistance (DR) in the pathogen is a major challenge for TB control. We conducted whole-genome analysis of seven Taiwan M. tuberculosis isolates: One drug susceptible (DS) and five DR Beijing lineage isolates and one DR Euro-American lineage isolate. Developing a new method for DR mutation identification and applying it to the next-generation sequencing (NGS) data from the 6 Beijing lineage isolates, we identified 13 known and 6 candidate DR mutations and provided experimental support for 4 of them. We assembled the genomes of one DS and two DR Beijing lineage isolates and the Euro-American lineage isolate using NGS data. Moreover, using both PacBio and NGS sequencing data, we obtained a high-quality assembly of an extensive DR Beijing lineage isolate. Comparative analysis of these five newly assembled genomes and two published complete genomes revealed a large number of genetic changes, including gene gains and losses, indels and translocations, suggesting rapid evolution of M. tuberculosis. We found the MazEF toxin–antitoxin system in all the seven isolates studied and several interesting mutations in MazEF proteins. Finally, we used the four assembled Beijing lineage genomes to construct a high-quality Beijing lineage reference genome that is DS and contains all the genes in the four genomes. It contains 212 genes not found in the standard reference H37Rv, which is Euro-American. It is therefore a better reference than H37Rv for the Beijing lineage, the predominant lineage in Asia.

Blast Fungal Genomes Show Frequent Chromosomal Changes, Gene Gains and Losses, and Effector Gene Turnover

Pyricularia is a fungal genus comprising several pathogenic species causing the blast disease in monocots. Pyricularia oryzae, the best-known species, infects rice, wheat, finger millet, and other crops. As past comparative and population genomics studies mainly focused on isolates of P. oryzae, the genomes of the other Pyricularia species have not been well explored. In this study, we obtained a chromosomal-level genome assembly of the finger millet isolate P. oryzae MZ5-1-6 and also highly contiguous assemblies of Pyricularia sp. LS, P. grisea, and P. pennisetigena. The differences in the genomic content of repetitive DNA sequences could largely explain the variation in genome size among these new genomes. Moreover, we found extensive gene gains and losses and structural changes among Pyricularia genomes, including a large interchromosomal translocation. We searched for homologs of known blast effectors across fungal taxa and found that most avirulence effectors are specific to Pyricularia, whereas many other effectors share homologs with distant fungal taxa. In particular, we discovered a novel effector family with metalloprotease activity, distinct from the well-known AVR-Pita family. We predicted 751 gene families containing putative effectors in 7 Pyricularia genomes and found that 60 of them showed differential expression in the P. oryzae MZ5-1-6 transcriptomes obtained under experimental conditions mimicking the pathogen infection process. In summary, this study increased our understanding of the structural, functional, and evolutionary genomics of the blast pathogen and identified new potential effector genes, providing useful data for developing crops with durable resistance.

The SnRK1-eIFiso4G1 signaling relay regulates the translation of specific mRNAs in Arabidopsis under submergence

Cellular responses to oxygen deprivation are essential for survival during energy crises in plants and animals. Hypoxia caused by submergence results in reprogramming of translation dynamic in plants, but the molecular mechanisms are not well understood. Here we show that Arabidopsis Snf1-related protein kinase 1 (SnRK1) phosphorylates the translation initiation factor eIFiso4G to regulate translation dynamic under submergence. In Arabidopsis, there are two eIFiso4G genes, eIFiso4G1 and eIFiso4G2, which belong to the eIF4G family. Both eIFiso4Gs were phosphorylated by SnRK1 under submergence. Interestingly, the eIFiso4G1 knockout mutant, but not the eIFiso4G2 mutant, became more sensitive to submergence, implying that eIFiso4G1 is involved in regulating submergence tolerance in Arabidopsis. Comparison of RNA sequences in the polysome fraction and the RNAs immunoprecipitated by eIFiso4G1 from Col-0 and the SnRK1 and eIFiso4G1 mutants revealed that lack of eIFiso4G1 phosphorylation disrupts the translation of specific mRNAs under submergence. Taken together, our findings suggest that the SnRK1-eIFiso4G1 relay controls the translation of an array of genes under hypoxia, including core hypoxia response genes and genes related to stress response and biosynthetic process.

Degradome Sequencing in Plants

Degradome sequencing provides large amounts of data regarding RNA degradation. The degradome library construction described here is modified from the 5'-rapid amplification of cDNA ends (5'-RACE), and each degradome cDNA is sequenced by next-generation sequencing (NGS). Degradome profiles provide information confirming miRNA-mediated cleavage of target genes and allow the identification of novel targets. Furthermore, degradome sequencing provides additional information for the study of RNA processing, such as information regarding RNA-binding proteins. In this chapter, we describe a detailed optimized protocol for constructing a degradome library with high yield and quality, along with NGS and data mining procedures. We hope that the degradome approach will be applied to further studies of non-model organisms.

Elevated auxin biosynthesis and transport underlie high vein density in C4 leaves

High vein density, a distinctive trait of C4 leaves, is central to both C3-to-C4 evolution and conversion of C3 to C4-like crops. We tested the hypothesis that high vein density in C4 leaves is due to elevated auxin biosynthesis and transport in developing leaves. Up-regulation of genes in auxin biosynthesis pathways and higher auxin content were found in developing C4 leaves compared with developing C3 leaves. The same observation held for maize foliar (C4) and husk (C3) leaf primordia. Moreover, auxin content and vein density were increased in loss-of-function mutants of Arabidopsis MYC2, a suppressor of auxin biosynthesis. Treatment with an auxin biosynthesis inhibitor or an auxin transport inhibitor led to much fewer veins in new leaves. Finally, both Arabidopsis thaliana auxin efflux transporter pin1 and influx transporter lax2 mutants showed reduced vein numbers. Thus, development of high leaf vein density requires elevated auxin biosynthesis and transport.

Identification and evolutionary analysis of long non-coding RNAs in zebra finch

Background

Long non-coding RNAs (lncRNAs) are important in various biological processes, but very few studies on lncRNA have been conducted in birds. To identify IncRNAs expressed during feather development, we analyzed single-stranded RNA-seq (ssRNA-seq) data from the anterior and posterior dorsal regions during zebra finch (Taeniopygia guttata) embryonic development. Using published transcriptomic data, we further analyzed the evolutionary conservation of IncRNAs in birds and amniotes.

Results

A total of 1,081 lncRNAs, including 965 intergenic lncRNAs (lincRNAs), 59 intronic lncRNAs, and 57 antisense lncRNAs (lncNATs), were identified using our newly developed pipeline. These avian IncRNAs share similar characteristics with lncRNAs in mammals, such as shorter transcript length, lower exon number, lower average expression level and less sequence conservation than mRNAs. However, the proportion of lncRNAs overlapping with transposable elements in birds is much lower than that in mammals. We predicted the functions of IncRNAs based on the enriched functions of co-expressed protein-coding genes. Clusters of lncRNAs associated with natal down development were identified. The sequences and expression levels of candidate lncRNAs that shared conserved sequences among birds were validated by qPCR in both zebra finch and chicken. Finally, we identified three highly conserved lncRNAs that may be associated with natal down development.

Conclusions

Our study provides the first systematical identification of avian lncRNAs using ssRNA-seq analysis and offers a resource of embryonically expressed lncRNAs in zebra finch. We also predicted the biological function of identified lncRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3506-z) contains supplementary material, which is available to authorized users.

Assembling the Setaria italica L. Beauv. genome into nine chromosomes and insights into regions affecting growth and drought tolerance

The diploid C4 plant foxtail millet (Setaria italica L. Beauv.) is an important crop in many parts of Africa and Asia for the vast consumption of its grain and ability to grow in harsh environments, but remains understudied in terms of complete genomic architecture. To date, there have been only two genome assembly and annotation efforts with neither assembly reaching over 86% of the estimated genome size. We have combined de novo assembly with custom reference-guided improvements on a popular cultivar of foxtail millet and have achieved a genome assembly of 477 Mbp in length, which represents over 97% of the estimated 490 Mbp. The assembly anchors over 98% of the predicted genes to the nine assembled nuclear chromosomes and contains more functional annotation gene models than previous assemblies. Our annotation has identified a large number of unique gene ontology terms related to metabolic activities, a region of chromosome 9 with several growth factor proteins, and regions syntenic with pearl millet or maize genomic regions that have been previously shown to affect growth. The new assembly and annotation for this important species can be used for detailed investigation and future innovations in growth for millet and other grains.

Regulatory Differences in Natal Down Development between Altricial Zebra Finch and Precocial Chicken

Birds can be classified into altricial and precocial. The hatchlings of altricial birds are almost naked, whereas those of precocial birds are covered with natal down. This regulatory divergence is thought to reflect environmental adaptation, but the molecular basis of the divergence is unclear. To address this issue, we chose the altricial zebra finch and the precocial chicken as the model animals. We noted that zebra finch hatchlings show natal down growth suppressed anterior dorsal (AD) skin but partially down-covered posterior dorsal (PD) skin. Comparing the transcriptomes of AD and PD skins, we found that the feather growth promoter SHH (sonic hedgehog) was expressed higher in PD skin than in AD skin. Moreover, the data suggested that the FGF (fibroblast growth factor)/Mitogen-activated protein kinase (MAPK) signaling pathway is involved in natal down growth suppression and that FGF16 is a candidate upstream signaling suppressor. Ectopic expression of FGF16 on chicken leg skin showed downregulation of SHH, upregulation of the feather growth suppressor FGF10, and suppression of feather bud elongation, similar to the phenotype found in zebra finch embryonic AD skin. Therefore, we propose that FGF16-related signals suppress natal down elongation and cause the naked AD skin in zebra finch. Our study provides insights into the regulatory divergence in natal down formation between precocial and altricial birds.

Transcriptomic analyses of regenerating adult feathers in chicken

Background

Feathers have diverse forms with hierarchical branching patterns and are an excellent model for studying the development and evolution of morphological traits. The complex structure of feathers allows for various types of morphological changes to occur. The genetic basis of the structural differences between different parts of a feather and between different types of feather is a fundamental question in the study of feather diversity, yet there is only limited relevant information for gene expression during feather development.

Results

We conducted transcriptomic analysis of five zones of feather morphologies from two feather types at different times during their regeneration after plucking. The expression profiles of genes associated with the development of feather structure were examined. We compared the gene expression patterns in different types of feathers and different portions of a feather and identified morphotype-specific gene expression patterns. Many candidate genes were identified for growth control, morphogenesis, or the differentiation of specific structures of different feather types.

Conclusion

This study laid the ground work for studying the evolutionary origin and diversification of feathers as abundant data were produced for the study of feather morphogenesis. It significantly increased our understanding of the complex molecular and cellular events in feather development processes and provided a foundation for future studies on the development of other skin appendages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1966-6) contains supplementary material, which is available to authorized users.

Expression Divergence of Chemosensory Genes between Drosophila sechellia and Its Sibling Species and Its Implications for Host Shift

Drosophila sechellia relies exclusively on the fruits of Morinda citrifolia, which are toxic to most insects, including its sibling species Drosophila melanogaster and Drosophila simulans. Although several odorant binding protein (Obp) genes and olfactory receptor (Or) genes have been suggested to be associated with the D. sechellia host shift, a broad view of how chemosensory genes have contributed to this shift is still lacking. We therefore studied the transcriptomes of antennae, the main organ responsible for detecting food resource and oviposition, of D. sechellia and its two sibling species. We wanted to know whether gene expression, particularly chemosensory genes, has diverged between D. sechellia and its two sibling species. Using a very stringent definition of differential gene expression, we found a higher percentage of chemosensory genes differentially expressed in the D. sechellia lineage (7.8%) than in the D. simulans lineage (5.4%); for upregulated chemosensory genes, the percentages were 8.8% in D. sechellia and 5.2% in D. simulans. Interestingly, Obp50a exhibited the highest upregulation, an approximately 100-fold increase, and Or85c--previously reported to be a larva-specific gene--showed approximately 20-fold upregulation in D. sechellia. Furthermore, Ir84a (ionotropic receptor 84a), which has been proposed to be associated with male courtship behavior, was significantly upregulated in D. sechellia. We also found expression divergence in most of the chemosensory gene families between D. sechellia and the two sibling species. Our observations suggest that the host shift of D. sechellia was associated with the enrichment of differentially expressed, particularly upregulated, chemosensory genes.

Genomic organization, transcriptomic analysis, and functional characterization of avian α- and β-keratins in diverse feather forms

Feathers are hallmark avian integument appendages, although they were also present on theropods. They are composed of flexible corneous materials made of α- and β-keratins, but their genomic organization and their functional roles in feathers have not been well studied. First, we made an exhaustive search of α- and β-keratin genes in the new chicken genome assembly (Galgal4). Then, using transcriptomic analysis, we studied α- and β-keratin gene expression patterns in five types of feather epidermis. The expression patterns of β-keratin genes were different in different feather types, whereas those of α-keratin genes were less variable. In addition, we obtained extensive α- and β-keratin mRNA in situ hybridization data, showing that α-keratins and β-keratins are preferentially expressed in different parts of the feather components. Together, our data suggest that feather morphological and structural diversity can largely be attributed to differential combinations of α- and β-keratin genes in different intrafeather regions and/or feather types from different body parts. The expression profiles provide new insights into the evolutionary origin and diversification of feathers. Finally, functional analysis using mutant chicken keratin forms based on those found in the human α-keratin mutation database led to abnormal phenotypes. This demonstrates that the chicken can be a convenient model for studying the molecular biology of human keratin-based diseases.

Genome-wide patterns of genetic variation in two domestic chickens

Domestic chickens are excellent models for investigating the genetic basis of phenotypic diversity, as numerous phenotypic changes in physiology, morphology, and behavior in chickens have been artificially selected. Genomic study is required to study genome-wide patterns of DNA variation for dissecting the genetic basis of phenotypic traits. We sequenced the genomes of the Silkie and the Taiwanese native chicken L2 at ∼23- and 25-fold average coverage depth, respectively, using Illumina sequencing. The reads were mapped onto the chicken reference genome (including 5.1% Ns) to 92.32% genome coverage for the two breeds. Using a stringent filter, we identified ∼7.6 million single-nucleotide polymorphisms (SNPs) and 8,839 copy number variations (CNVs) in the mapped regions; 42% of the SNPs have not found in other chickens before. Among the 68,906 SNPs annotated in the chicken sequence assembly, 27,852 were nonsynonymous SNPs located in 13,537 genes. We also identified hundreds of shared and divergent structural and copy number variants in intronic and intergenic regions and in coding regions in the two breeds. Functional enrichments of identified genetic variants were discussed. Radical nsSNP-containing immunity genes were enriched in the QTL regions associated with some economic traits for both breeds. Moreover, genetic changes involved in selective sweeps were detected. From the selective sweeps identified in our two breeds, several genes associated with growth, appetite, and metabolic regulation were identified. Our study provides a framework for genetic and genomic research of domestic chickens and facilitates the domestic chicken as an avian model for genomic, biomedical, and evolutionary studies.

Transcriptional profiling of adult Drosophila antennae by high-throughput sequencing

Background

Antennae of fruit flies are the major organs responsible for detecting environmental volatiles, e.g., egg-laying substrates. An adult antenna contains many sensilla full of olfactory sensory neurons, where olfactory receptor (Or) genes are expressed. Each sensory neuron only expresses up to three receptors, making it difficult to estimate expression levels by conventional methods. In this study, we applied Illumina RNA sequencing (RNA-seq) to study the expression levels of Or and other genes in fly antennae.

Results

RNA from approximately 1,200 pairs of adult antennae from each sex of Drosophila melanogaster was used to obtain the antennal transcriptome of each sex. We detected approximately 12,000 genes expressed in antennae of either sex. The most highly expressed genes included pheromone-binding genes, transmembrane transporter genes, and sensory reception genes. Among the 61 annotated Or genes, we observed 53 and 54 genes (approximately 90%) expressed (fragments per kilobase of exon per million fragments mapped (FPKM) > 0.05) in male and female antennae, respectively; approximately 25 genes were expressed with FPKM > 15. Compared to previous studies, which extracted RNA from the whole body or head and used microarrays, antenna-specific transcriptomes obtained by RNA-seq provided more reliable estimates of gene expression levels and revealed many lowly expressed genes. Ninty-one genes, including one odorant-binding protein (Obp) gene and four Or genes, were differentially expressed between male and female antennae. These sexually biased genes were enriched on the X chromosome and showed enrichment in different gene ontology categories for male and female flies. The present and previous data together suggest that a gene family with putative immune response functions is related to pheromone detection and involved in the courtship behavior of male flies.

Conclusions

Tissue-specific RNA-seq is powerful for detecting lowly expressed genes. Our study provides new insight into the expression of olfactory-related genes in Drosophila antennae.

Inheritance of gene expression level and selective constraints on trans- and cis-regulatory changes in yeast

Gene expression evolution can be caused by changes in cis- or trans-regulatory elements or both. As cis and trans regulation operate through different molecular mechanisms, cis and trans mutations may show different inheritance patterns and may be subjected to different selective constraints. To investigate these issues, we obtained and analyzed gene expression data from two Saccharomyces cerevisiae strains and their hybrid, using high-throughput sequencing. Our data indicate that compared with other types of genes, those with antagonistic cis-trans interactions are more likely to exhibit over- or underdominant inheritance of expression level. Moreover, in accordance with previous studies, genes with trans variants tend to have a dominant inheritance pattern, whereas cis variants are enriched for additive inheritance. In addition, cis regulatory differences contribute more to expression differences between species than within species, whereas trans regulatory differences show a stronger association between divergence and polymorphism. Our data indicate that in the trans component of gene expression differences genes subjected to weaker selective constraints tend to have an excess of polymorphism over divergence compared with those subjected to stronger selective constraints. In contrast, in the cis component, this difference between genes under stronger and weaker selective constraint is mostly absent. To explain these observations, we propose that purifying selection more strongly shapes trans changes than cis changes and that positive selection may have significantly contributed to cis regulatory divergence.

Characterizing regulatory and functional differentiation between maize mesophyll and bundle sheath cells by transcriptomic analysis

To study the regulatory and functional differentiation between the mesophyll (M) and bundle sheath (BS) cells of maize (Zea mays), we isolated large quantities of highly homogeneous M and BS cells from newly matured second leaves for transcriptome profiling by RNA sequencing. A total of 52,421 annotated genes with at least one read were found in the two transcriptomes. Defining a gene with more than one read per kilobase per million mapped reads as expressed, we identified 18,482 expressed genes; 14,972 were expressed in M cells, including 53 M-enriched transcription factor (TF) genes, whereas 17,269 were expressed in BS cells, including 214 BS-enriched TF genes. Interestingly, many TF gene families show a conspicuous BS preference in expression. Pathway analyses reveal differentiation between the two cell types in various functional categories, with the M cells playing more important roles in light reaction, protein synthesis and folding, tetrapyrrole synthesis, and RNA binding, while the BS cells specialize in transport, signaling, protein degradation and posttranslational modification, major carbon, hydrogen, and oxygen metabolism, cell division and organization, and development. Genes coding for several transporters involved in the shuttle of C(4) metabolites and BS cell wall development have been identified, to our knowledge, for the first time. This comprehensive data set will be useful for studying M/BS differentiation in regulation and function.

A highly efficient β-glucosidase from the buffalo rumen fungus Neocallimastix patriciarum W5

BACKGROUND

Cellulose, which is the most abundant renewable biomass on earth, is a potential bio-resource of alternative energy. The hydrolysis of plant polysaccharides is catalyzed by microbial cellulases, including endo-β-1,4-glucanases, cellobiohydrolases, cellodextrinases, and β-glucosidases. Converting cellobiose by β-glucosidases is the key factor for reducing cellobiose inhibition and enhancing the efficiency of cellulolytic enzymes for cellulosic ethanol production.

RESULTS

In this study, a cDNA encoding β-glucosidase was isolated from the buffalo rumen fungus Neocallimastix patriciarum W5 and is named NpaBGS. It has a length of 2,331 bp with an open reading frame coding for a protein of 776 amino acid residues, corresponding to a theoretical molecular mass of 85.1 kDa and isoelectric point of 4.4. Two GH3 catalytic domains were found at the N and C terminals of NpaBGS by sequence analysis. The cDNA was expressed in Pichia pastoris and after protein purification, the enzyme displayed a specific activity of 34.5 U/mg against cellobiose as the substrate. Enzymatic assays showed that NpaBGS was active on short cello-oligosaccharides from various substrates. A weak activity in carboxymethyl cellulose (CMC) digestion indicated that the enzyme might also have the function of an endoglucanase. The optimal activity was detected at 40°C and pH 5 ~ 6, showing that the enzyme prefers a weak acid condition. Moreover, its activity could be enhanced at 50°C by adding Mg2+ or Mn2+ ions. Interestingly, in simultaneous saccharification and fermentation (SSF) experiments using Saccharomyces cerevisiae BY4741 or Kluyveromyces marxianus KY3 as the fermentation yeast, NpaBGS showed advantages in cell growth, glucose production, and ethanol production over the commercial enzyme Novo 188. Moreover, we showed that the KY3 strain engineered with the NpaNGS gene can utilize 2 % dry napiergrass as the sole carbon source to produce 3.32 mg/ml ethanol when Celluclast 1.5 L was added to the SSF system.

CONCLUSION

Our characterizations of the novel β-glucosidase NpaBGS revealed that it has a preference of weak acidity for optimal yeast fermentation and an optimal temperature of ~40°C. Since NpaBGS performs better than Novo 188 under the living conditions of fermentation yeasts, it has the potential to be a suitable enzyme for SSF.

Transcriptomes of mouse olfactory epithelium reveal sexual differences in odorant detection

To sense numerous odorants and chemicals, animals have evolved a large number of olfactory receptor genes (Olfrs) in their genome. In particular, the house mouse has ∼1,100 genes in the Olfr gene family. This makes the mouse a good model organism to study Olfr genes and olfaction-related genes. To date, whether male and female mice possess the same ability in detecting environmental odorants is still unknown. Using the next generation sequencing technology (paired-end mRNA-seq), we detected 1,088 expressed Olfr genes in both male and female olfactory epithelium. We found that not only Olfr genes but also odorant-binding protein (Obp) genes have evolved rapidly in the mouse lineage. Interestingly, Olfr genes tend to express at a higher level in males than in females, whereas the Obp genes clustered on the X chromosome show the opposite trend. These observations may imply a more efficient odorant-transporting system in females, whereas a more active Olfr gene expressing system in males. In addition, we detected the expression of two genes encoding major urinary proteins, which have been proposed to bind and transport pheromones or act as pheromones in mouse urine. This observation suggests a role of main olfactory system (MOS) in pheromone detection, contrary to the view that only accessory olfactory system (AOS) is involved in pheromone detection. This study suggests the sexual differences in detecting environmental odorants in MOS and demonstrates that mRNA-seq provides a powerful tool for detecting genes with low expression levels and with high sequence similarities.

Natural selection on cis and trans regulation in yeasts

Gene expression is regulated both by cis elements, which are DNA segments closely linked to the genes they regulate, and by trans factors, which are usually proteins capable of diffusing to unlinked genes. Understanding the patterns and sources of regulatory variation is crucial for understanding phenotypic and genome evolution. Here, we measure genome-wide allele-specific expression by deep sequencing to investigate the patterns of cis and trans expression variation between two strains of Saccharomyces cerevisiae. We propose a statistical modeling framework based on the binomial distribution that simultaneously addresses normalization of read counts derived from different parents and estimating the cis and trans expression variation parameters. We find that expression polymorphism in yeast is common for both cis and trans, though trans variation is more common. Constraint in expression evolution is correlated with other hallmarks of constraint, including gene essentiality, number of protein interaction partners, and constraint in amino acid substitution, indicating that both cis and trans polymorphism are clearly under purifying selection, though trans variation appears to be more sensitive to selective constraint. Comparing interspecific expression divergence between S. cerevisiae and S. paradoxus to our intraspecific variation suggests a significant departure from a neutral model of molecular evolution. A further examination of correlation between polymorphism and divergence within each category suggests that cis divergence is more frequently mediated by positive Darwinian selection than is trans divergence.