Genetic relatedness does not predict the queen's successors in the primitively eusocial wasp, Ropalidia marginata.

Ropalidia marginata is a social wasp in which colonies consist of a single fertile queen and several sterile workers. If the queen is removed, one of the workers, potential queen (PQ), becomes hyperaggressive and becomes the next queen. The identity of the PQ cannot be predicted in the presence of the queen. The probability of a worker succeeding the queen is uncorrelated with her body size, dominance rank, ovarian or mating status, but imperfectly correlated with her age. Here, we investigate whether genetic relatedness help to predict the queen's successors. We constructed models based on successors being (i) most closely related to the queen, (ii) most closely related to the immediate predecessor queen/PQ, or (iii) having the highest relatedness to the majority of the workers; and (iv) having the highest average relatedness to all the workers. We predicted five successors from each of these models using pair-wise genetic relatedness estimated from polymorphic microsatellite loci. We independently performed serial queen/PQ removal experiments and compared the observed sequence of successors with the predictions from the models. The predictions of none of the models matched the experimental results; on an average 5-6 individuals predicted by the models were bypassed in the experiment. Thus, genetic relatedness is inadequate to predict the queen's successors in this species. We discuss why relatedness sometimes predicts the patterns of altruistic behaviour and sometimes not, and argue that the cost and benefit terms in Hamilton's rule, i.e. ecology, should be vigorously investigated when relatedness does not have adequate explanatory power.

Role of Bmznf-2, a Bombyx mori CCCH zinc finger gene, in masculinisation and differential splicing of Bmtra-2.

Deciphering the regulatory factors involved in Bombyx mori sex determination has been a puzzle, challenging researchers for nearly a century now. The pre-mRNA of B. mori doublesex (Bmdsx), a master regulator gene of sexual differentiation, is differentially spliced, producing Bmdsxm and Bmdsxf transcripts in males and females respectively. The putative proteins encoded by these differential transcripts orchestrate antagonistic functions, which lead to sexual differentiation. A recent study in B. mori illustrated the role of a W-derived fem piRNA in conferring femaleness. In females, the fem piRNA was shown to suppress the activity of a Z-linked CCCH type zinc finger (znf) gene, Masculiniser (masc), which indirectly promotes the Bmdsxm type of splicing. In this study, we report a novel autosomal (Chr 25) CCCH type znf motif encoding gene Bmznf-2 as one of the potential factors in the Bmdsx sex specific differential splicing, and we also provide insights into its role in the alternative splicing of Bmtra2 by using ovary derived BmN cells. Over-expression of Bmznf-2 induced Bmdsxm type of splicing (masculinisation) with a correspondingly reduced expression of Bmdsxf type isoform in BmN cells. Further, the site-directed mutational studies targeting the tandem CCCH znf motifs revealed their indispensability in the observed phenotype of masculinisation. Additionally, the dual luciferase assays in BmN cells using 5' UTR region of the Bmznf-2 strongly implied the existence of a translational repression over this gene. From these findings, we propose Bmznf-2 to be one of the potential factors of masculinisation similar to Masc. From the growing number of Bmdsx splicing regulators, we assume that the sex determination cascade of B. mori is quite intricate in nature; hence, it has to be further investigated for its comprehensive understanding.

Noduler an immune protein augments infection-induced cell proliferation through cross-talking with p38 MAPK.

Noduler, an immune protein that mediates nodule formation by binding to specific bacteria and hemocytes was previously reported in the wild tasar silkworm, Antheraea mylitta. However, the molecular mechanism underlying nodulation in lepidopterans remains unclear. The present study is performed to investigate the functional connection between Noduler with various signalling pathways. It was observed that Noduler is an upstream factor in the phenoloxidase cascade and its knockdown has no direct effect on Toll/Imd pathway inducible genes. Additionally, Noduler was shown to stimulate cell proliferation via activation of p38 mitogen-activated protein kinase (MAPK). Inhibition of p38 in the infected hemocytes cultured in vitro resulted in reduced cell proliferation and melanization. These results suggest that Noduler mediates nodulation via p38/MAPK signalling. This is the first report implicating the p38 MAPK signalling pathway in the nodulation response of insects.

Construction, complete sequence, and annotation of a BAC contig covering the silkworm chorion locus.

The silkmoth chorion was studied extensively by F.C. Kafatos' group for almost 40 years. However, the complete structure of the chorion locus was not obtained in the genome sequence of Bombyx mori published in 2008 due to repetitive sequences, resulting in gaps and an incomplete view of the locus. To obtain the complete sequence of the chorion locus, expressed sequence tags (ESTs) derived from follicular epithelium cells were used as probes to screen a bacterial artificial chromosome (BAC) library. Seven BACs were selected to construct a contig which covered the whole chorion locus. By Sanger sequencing, we successfully obtained complete sequences of the chorion locus spanning 871,711 base pairs on chromosome 2, where we annotated 127 chorion genes. The dataset reported here will recruit more researchers to revisit one of the oldest model systems which has been used to study developmentally regulated gene expression. It also provides insights into egg development and fertilization mechanisms and is relevant to applications related to improvements in breeding procedures and transgenesis.

A comprehensive analysis of the chorion locus in silkmoth.

Despite more than 40 years of intense study, essential features of the silkmoth chorion (eggshell) are still not fully understood. To determine the precise structure of the chorion locus, we performed extensive EST analysis, constructed a bacterial artificial chromosome (BAC) contig, and obtained a continuous genomic sequence of 871,711 base pairs. We annotated 127 chorion genes in two segments interrupted by a 164 kb region with 5 non-chorion genes, orthologs of which were on chorion bearing scaffolds in 4 ditrysian families. Detailed transcriptome analysis revealed expression throughout choriogenesis of most chorion genes originally categorized as "middle", and evidence for diverse regulatory mechanisms including cis-elements, alternative splicing and promoter utilization, and antisense RNA. Phylogenetic analysis revealed multigene family associations and faster evolution of early chorion genes and transcriptionally active pseudogenes. Proteomics analysis identified 99 chorion proteins in the eggshell and micropyle localization of 1 early and 6 Hc chorion proteins.

Molecular architecture of silk fibroin of Indian golden silkmoth, Antheraea assama.

The golden silk spun by Indian golden silkmoth Antheraea assama, is regarded for its shimmering golden luster, tenacity and value as biomaterial. This report describes the gene coding for golden silk H-fibroin (AaFhc), its expression, full-length sequence and structurally important motifs discerning the underlying genetic and biochemical factors responsible for its much sought-after properties. The coding region, with biased isocodons, encodes highly repetitious crystalline core, flanked by a pair of 5' and 3' non-repetitious ends. AaFhc mRNA expression is strictly territorial, confined to the posterior silk gland, encoding a protein of size 230 kDa, which makes homodimers making the elementary structural units of the fibrous core of the golden silk. Characteristic polyalanine repeats that make tight ß-sheet crystals alternate with non-polyalanine repeats that make less orderly antiparallel ß-sheets, ß-turns and partial α-helices. Phylogenetic analysis of the conserved N-terminal amorphous motif and the comparative analysis of the crystalline region with other saturniid H-fibroins reveal that AaFhc has longer, numerous and relatively uniform repeat motifs with lower serine content that assume tighter ß-crystals and denser packing, which are speculated to be responsible for its acclaimed properties of higher tensile strength and higher refractive index responsible for golden luster.

Discovery and mapping of genomic regions governing economically important traits of Basmati rice.

BACKGROUND: Basmati rice, originated in the foothills of Himalayas, commands a premium price in the domestic and international markets on account of its unique quality traits. The complex genetic nature of unique traits of Basmati as well as tedious screening methodologies involved in quality testing have been serious constraints to breeding quality Basmati. In the present study, we made an attempt to identify the genomic regions governing unique traits of Basmati rice. RESULTS: A total of 34 Quantitative Trait Loci (QTLs) for 16 economically important traits of Basmati rice were identified employing F(2), F(3) and Recombinant Inbred Line (RIL) mapping populations derived from a cross between Basmati370 (traditional Basmati) and Jaya (semi-dwarf rice). Out of which, 12 QTLs contributing to more than 15 % phenotypic variance were identified and considered as major effect QTLs. Four major effect QTLs coincide with the already known genes viz., sd1, GS3, alk1 and fgr governing plant height, grain size, alkali spreading value and aroma, respectively. For the remaining major QTLs, candidate genes were predicted as auxin response factor for filled grains, soluble starch synthase 3 for chalkiness and VQ domain containing protein for grain breadth and grain weight QTLs, based on the presence of non-synonymous single nucleotide polymorphism (SNPs) that were identified by comparing Basmati genome sequence with that of Nipponbare. CONCLUSIONS: To the best of our knowledge, the current study is the first attempt ever made to carry out genome-wide mapping for the dissection of the genetic basis of economically important traits of Basmati rice. The promising QTLs controlling important traits in Basmati rice, identified in this study, can be used as candidates for future marker-assisted breeding.

De novo assembly and transcriptome analysis of the Mediterranean fruit fly Ceratitis capitata early embryos.

The agricultural pest Ceratitis capitata, also known as the Mediterranean fruit fly or Medfly, belongs to the Tephritidae family, which includes a large number of other damaging pest species. The Medfly has been the first non-drosophilid fly species which has been genetically transformed paving the way for designing genetic-based pest control strategies. Furthermore, it is an experimentally tractable model, in which transient and transgene-mediated RNAi have been successfully used. We applied Illumina sequencing to total RNA preparations of 8-10 hours old embryos of C. capitata, This developmental window corresponds to the blastoderm cellularization stage. In summary, we assembled 42,614 transcripts which cluster in 26,319 unique transcripts of which 11,045 correspond to protein coding genes; we identified several hundreds of long ncRNAs; we found an enrichment of transcripts encoding RNA binding proteins among the highly expressed transcripts, such as CcTRA-2, known to be necessary to establish and, most likely, to maintain female sex of C. capitata. Our study is the first de novo assembly performed for Ceratitis capitata based on Illumina NGS technology during embryogenesis and it adds novel data to the previously published C. capitata EST databases. We expect that it will be useful for a variety of applications such as gene cloning and phylogenetic analyses, as well as to advance genetic research and biotechnological applications in the Medfly and other related Tephritidae.

Positive selection drives faster-Z evolution in silkmoths.

Genes linked to X or Z chromosomes, which are hemizygous in the heterogametic sex, are predicted to evolve at different rates than those on autosomes. This "faster-X effect" can arise either as a consequence of hemizygosity, which leads to more efficient selection for recessive beneficial mutations in the heterogametic sex, or as a consequence of reduced effective population size of the hemizygous chromosome, which leads to increased fixation of weakly deleterious mutations due to genetic drift. Empirical results to date suggest that, while the overall pattern across taxa is complicated, systems with male heterogamy show a faster-X effect attributable to more efficient selection, whereas the faster-Z effect in female-heterogametic taxa is attributable to increased drift. To test the generality of the faster-Z pattern seen in birds and snakes, we sequenced the genome of the lepidopteran silkmoth Bombyx huttoni. We show that silkmoths experience faster-Z evolution, but unlike in birds and snakes, the faster-Z effect appears to be attributable to more efficient positive selection. These results suggest that female heterogamy alone is unlikely to explain the reduced efficacy of selection on vertebrate Z chromosomes. It is likely that many factors, including differences in overall effective population size, influence Z chromosome evolution.

Nodulation: an unexplored cellular defense mechanism in insects.

Nodulation is a highly conserved process that involves aggregation of cells around microorganisms, leading to their entrapment with the help of cellular milieu. In insects upon infection, the humoral and cellular arms of the innate immune system orchestrate recognition of pathogens facilitating effector responses through various signaling pathways. Existing data suggests a wide range of immune functions for multiple pattern recognition molecules but their role in nodulation is not known. Hence, an in-depth knowledge of components implicated in the signaling pathways across diverse species is crucial for understanding their evolutionary conservation. Here, we attempted to consolidate available information on the nodulation response in insects and made an analogy with other known systems.

Computational analyses of protein coded by rice (Oryza sativa japonica) cDNA (GI: 32984786) indicate lectin like Ca(2+) binding properties for Eicosapenta Peptide Repeats (EPRs).

Eicosapenta peptide repeats (EPRs) occur exclusively in flowering plant genomes and exhibit very high amino acid residue conservation across occurrence. DNA and amino acid sequence searches yielded no indications about the function due to absence of similarity to known sequences. Tertiary structure of an EPR protein coded by rice (Oryza sativa japonica) cDNA (GI: 32984786) was determined based on ab initio methodology in order to draw clues on functional significance of EPRs. The resultant structure comprised of seven α-helices and thirteen anti-parallel ß-sheets. Surface-mapping of conserved residues onto the structure deduced that (i) regions equivalent to ß α4- the primary function of EPR protein could be Ca(2+) binding, and (iii) the putative EPR Ca(2+) binding domain is structurally similar to calcium-binding domains of plant lectins. Additionally, the phylogenetic analysis showed an evolving taxa-specific distribution of EPR proteins observed in some GNA-like lectins.

Large scale full-length cDNA sequencing reveals a unique genomic landscape in a lepidopteran model insect, Bombyx mori.

The establishment of a complete genomic sequence of silkworm, the model species of Lepidoptera, laid a foundation for its functional genomics. A more complete annotation of the genome will benefit functional and comparative studies and accelerate extensive industrial applications for this insect. To realize these goals, we embarked upon a large-scale full-length cDNA collection from 21 full-length cDNA libraries derived from 14 tissues of the domesticated silkworm and performed full sequencing by primer walking for 11,104 full-length cDNAs. The large average intron size was 1904 bp, resulting from a high accumulation of transposons. Using gene models predicted by GLEAN and published mRNAs, we identified 16,823 gene loci on the silkworm genome assembly. Orthology analysis of 153 species, including 11 insects, revealed that among three Lepidoptera including Monarch and Heliconius butterflies, the 403 largest silkworm-specific genes were composed mainly of protective immunity, hormone-related, and characteristic structural proteins. Analysis of testis-/ovary-specific genes revealed distinctive features of sexual dimorphism, including depletion of ovary-specific genes on the Z chromosome in contrast to an enrichment of testis-specific genes. More than 40% of genes expressed in specific tissues mapped in tissue-specific chromosomal clusters. The newly obtained FL-cDNA sequences enabled us to annotate the genome of this lepidopteran model insect more accurately, enhancing genomic and functional studies of Lepidoptera and comparative analyses with other insect orders, and yielding new insights into the evolution and organization of lepidopteran-specific genes.

Engineering silkworms for resistance to baculovirus through multigene RNA interference.

Bombyx mori nucleopolyhedrovirus (BmNPV) that infects the silkworm, B. mori, accounts for >50% of silk cocoon crop losses globally. We speculated that simultaneous targeting of several BmNPV essential genes in transgenic silkworm would elicit a stable defense against the virus. We introduced into the silkworm germline the vectors carrying short sequences of four essential BmNPV genes in tandem, either in sense or antisense or in inverted-repeat arrangement. The transgenic silkworms carrying the inverted repeat-containing transgene showed stable protection against high doses of baculovirus infection. Further, the antiviral trait was incorporated to a commercially productive silkworm strain highly susceptible to BmNPV. This led to combining the high-yielding cocoon and silk traits of the parental commercial strain and a very high level of refractoriness (>75% survival rate as compared to <15% in nontransgenic lines) to baculovirus infection conferred by the transgene. We also observed impaired infectivity of the occlusion bodies derived from the transgenic lines as compared to the wild-type ones. Currently, large-scale exploitation of these transgenic lines is underway to bring about economic transformation of sericulture.

Pyrosequencing-based transcriptome analysis of the asian rice gall midge reveals differential response during compatible and incompatible interaction.

The Asian rice gall midge (Orseolia oryzae) is a major pest responsible for immense loss in rice productivity. Currently, very little knowledge exists with regard to this insect at the molecular level. The present study was initiated with the aim of developing molecular resources as well as identifying alterations at the transcriptome level in the gall midge maggots that are in a compatible (SH) or in an incompatible interaction (RH) with their rice host. Roche 454 pyrosequencing strategy was used to develop both transcriptomics and genomics resources that led to the identification of 79,028 and 85,395 EST sequences from gall midge biotype 4 (GMB4) maggots feeding on a susceptible and resistant rice variety, TN1 (SH) and Suraksha (RH), respectively. Comparative transcriptome analysis of the maggots in SH and RH revealed over-representation of transcripts from proteolysis and protein phosphorylation in maggots from RH. In contrast, over-representation of transcripts for translation, regulation of transcription and transcripts involved in electron transport chain were observed in maggots from SH. This investigation, besides unveiling various mechanisms underlying insect-plant interactions, will also lead to a better understanding of strategies adopted by insects in general, and the Asian rice gall midge in particular, to overcome host defense.

Genetic variation within native populations of endemic silkmoth Antheraea assamensis (Helfer) from Northeast India indicates need for in situ conservation.

A. assamensis is a phytophagous Lepidoptera from Northeast India reared on host trees of Lauraceae family for its characteristic cocoon silk. Source of these cocoons are domesticated farm stocks that crash frequently and/or wild insect populations that provide new cultures. The need to reduce dependence on wild populations for cocoons necessitates assessment of genetic diversity in cultivated and wild populations. Molecular markers based on PCR of Inter-simple sequence repeats (ISSR) and simple sequence repeats (SSR) were used with four populations of wild insects and eleven populations of cultivated insects. Wild populations had high genetic diversity estimates (H(i) = 0.25; H(S) = 0.28; H(E) = 0.42) and at least one population contained private alleles. Both marker systems indicated that genetic variability within populations examined was significantly high. Among cultivated populations, insects of the Upper Assam region (H(i) = 0.19; H(S) = 0.18; H(E) = 0) were genetically distinct (F(ST) = 0.38 with both marker systems) from insects of Lower Assam (H(i) =0.24; H(S) =0.25; H(E) = 0.3). Sequencing of polymorphic amplicons suggested transposition as a mechanism for maintaining genomic diversity. Implications for conservation of native populations in the wild and preserving in-farm diversity are discussed.

Genetic diversity and population structure of Indian golden silkmoth (Antheraea assama).

BACKGROUND: The Indian golden saturniid silkmoth (Antheraea assama), popularly known as muga silkmoth, is a semi-domesticated silk producing insect confined to a narrow habitat range of the northeastern region of India. Owing to the prevailing socio-political problems, the muga silkworm habitats in the northeastern region have not been accessible hampering the phylogeography studies of this rare silkmoth. Recently, we have been successful in our attempt to collect muga cocoon samples, although to a limited extent, from their natural habitats. Out of 87 microsatellite markers developed previously for A. assama, 13 informative markers were employed to genotype 97 individuals from six populations and analyzed their population structure and genetic variation. METHODOLOGY/PRINCIPAL FINDINGS: We observed highly significant genetic diversity in one of the populations (WWS-1, a population derived from West Garo Hills region of Meghalaya state). Further analysis with and without WWS-1 population revealed that dramatic genetic differentiation (global F(ST) = 0.301) was due to high genetic diversity contributed by WWS-1 population. Analysis of the remaining five populations (excluding WWS-1) showed a marked reduction in the number of alleles at all the employed loci. Structure analysis showed the presence of only two clusters: one formed by WWS-1 population and the other included the remaining five populations, inferring that there is no significant genetic diversity within and between these five populations, and suggesting that these five populations are probably derived from a single population. Patterns of recent population bottlenecks were not evident in any of the six populations studied. CONCLUSIONS/SIGNIFICANCE: A. assama inhabiting the WWS-1 region revealed very high genetic diversity, and was genetically divergent from the five populations studied. The efforts should be continued to identify and study such populations from this region as well as other muga silkworm habitats. The information generated will be very useful in conservation of dwindling muga culture in Northeast India.

Role of sequence encoded κB DNA geometry in gene regulation by Dorsal.

Many proteins of the Rel family can act as both transcriptional activators and repressors. However, mechanism that discerns the 'activator/repressor' functions of Rel-proteins such as Dorsal (Drosophila homologue of mammalian NFκB) is not understood. Using genomic, biophysical and biochemical approaches, we demonstrate that the underlying principle of this functional specificity lies in the 'sequence-encoded structure' of the κB-DNA. We show that Dorsal-binding motifs exist in distinct activator and repressor conformations. Molecular dynamics of DNA-Dorsal complexes revealed that repressor κB-motifs typically have A-tract and flexible conformation that facilitates interaction with co-repressors. Deformable structure of repressor motifs, is due to changes in the hydrogen bonding in A:T pair in the 'A-tract' core. The sixth nucleotide in the nonameric κB-motif, 'A' (A(6)) in the repressor motifs and 'T' (T(6)) in the activator motifs, is critical to confer this functional specificity as A(6) → T(6) mutation transformed flexible repressor conformation into a rigid activator conformation. These results highlight that 'sequence encoded κB DNA-geometry' regulates gene expression by exerting allosteric effect on binding of Rel proteins which in turn regulates interaction with co-regulators. Further, we identified and characterized putative repressor motifs in Dl-target genes, which can potentially aid in functional annotation of Dorsal gene regulatory network.

Isolation and characterization of microsatellite loci in the Asian rice gall midge (Orseolia oryzae) (Diptera: Cecidomyiidae).

Microsatellite loci were isolated from the genomic DNA of the Asian rice gall midge, Orseolia oryzae (Wood-Mason) using a hybridization capture approach. A total of 90 non-redundant primer pairs, representing unique loci, were designed. These simple sequence repeat (SSR) markers represented di (72%), tri (15.3%), and complex repeats (12.7%). Three biotypes of gall midge (20 individuals for each biotype) were screened using these SSRs. The results revealed that 15 loci were hyper variable and showed polymorphism among different biotypes of this pest. The number of alleles ranged from two to 11 and expected heterozygosity was above 0.5. Inheritance studies with three markers (observed to be polymorphic between sexes) revealed sex linked inheritance of two SSRs (Oosat55 and Oosat59) and autosomal inheritance of one marker (Oosat43). These markers will prove to be a useful tool to devise strategies for integrated pest management and in the study of biotype evolution in this important rice pest.

Novel female-specific splice form of dsx in the silkworm, Bombyx mori.

The Bombyx mori doublesex (Bmdsx), a homologue of doublesex of Drosophila, is the bottom most gene of the sex determination cascade. Bmdsx plays a very crucial role in somatic sexual development. Its pre-mRNA sex-specifically splices to generate two splice variants; one encodes female-specific and the other encodes male-specific polypeptides which differ only at their C-termini. The open reading frame of Bmdsx consists of 5 exons, of which exons 3 and 4 are female-specific and are skipped in males. In the present study, we have identified a third splice form of the Bmdsx which is specific only to females and differs from the previously reported Bmdsxf isoform by the presence of 15 bp sequence. This new female splice form is generated as a result of alternative 5' splice site selection in the third exon adding additional 15 bp sequence in exon 3 which results in alteration of the reading frame leading to incorporation of an early stop codon. Thus the protein encoded by this splice form is 20 aa shorter than the known BmDsxF. Initial results obtained from the study of dsx homologues in Saturniid silkmoths suggest that both the female-specific Dsx proteins are essential for female sexual differentiation. It remains to be seen whether female-specific multiple splice forms of dsx are characteristic feature of only silkmoths or widespread among lepidopterans. The findings that sex determination mechanism is unique in lepidopterans offer an opportunity to develop genetic sexing methods in beneficial as well as economically destructive lepidopteran pests.

Dynamic repositioning of dorsal to two different kappaB motifs controls its autoregulation during immune response in Drosophila.

Autoregulation is one of the mechanisms of imparting feedback control on gene expression. Positive autoregulatory feedback results in induction of a gene, and negative feedback leads to its suppression. Here, we report an interesting mechanism of autoregulation operating on Drosophila Rel gene dorsal that can activate as well as repress its expression. Using biochemical and genetic approaches, we show that upon immune challenge Dorsal regulates its activation as well as repression by dynamically binding to two different kappaB motifs, kappaB(I) (intronic kappaB) and kappaB(P) (promoter kappaB), present in the dorsal gene. Although the kappaB(I) motif functions as an enhancer, the kappaB(P) motif acts as a transcriptional repressor. Interestingly, Dorsal binding to these two motifs is dynamic; immediately upon immune challenge, Dorsal binds to the kappaB(I) leading to auto-activation, whereas at the terminal phase of the immune response, it is removed from the kappaB(I) and repositioned at the kappaB(P), resulting in its repression. Furthermore, we show that repression of Dorsal as well as its binding to the kappaB(P) depends on the transcription factor AP1. Depletion of AP1 by RNA interference resulted in constitutive expression of Dorsal. In conclusion, this study suggests that during acute phase response dorsal is regulated by following two subcircuits: (i) Dl-kappaB(I) for activation and (ii) Dl-AP1-kappaB(P) for repression. These two subcircuits are temporally delineated and bring about overall regulation of dorsal during immune response. These results suggest the presence of a previously unknown mechanism of Dorsal autoregulation in immune-challenged Drosophila.