Mapping QTLs underpin nutrition components in aromatic rice germplasm.

As rice is an important staple food globally, research for development and enhancement of its nutritional value it is an imperative task. Identification of nutrient enriched rice germplasm and exploiting them for breeding programme is the easiest way to develop better quality rice. In this study, we analyzed 113 aromatic rice germplasm in order to identify quantitative trait loci (QTL) underpinning nutrition components and determined by measuring the normal frequency distribution for Fe, Zn, amylose, and protein content in those rice germplasm. Comparatively, the germplasm Radhuni pagal, Kalobakri, Thakurbhog (26.6 ppm) and Hatisail exhibited the highest mean values for Fe (16.9 ppm), Zn (34.1 ppm), amylose (26.6 ppm) and protein content (11.0 ppm), respectively. Moreover, a significant linear relationship (R2 = 0.693) was observed between Fe and Zn contents. Cluster analysis based on Mahalanobis D2 distances revealed four major clusters of 113 rice germplasm, with cluster III containing a maximum 37 germplasm and a maximum inter-cluster distance between clusters III and IV. The 45 polymorphic SSRs and four trait associations exhibited eight significant quantitative trait loci (QTL) located on eight different chromosomes using composite interval mapping (CIM). The highly significant QTL (variance 7.89%, LOD 2.02) for protein content (QTL.pro.1) was observed on chromosome 1 at 94.9cM position. Also, four QTLs for amylose content were observed with the highly significant QTL.amy.8 located on chromosome 8 exhibiting 7.2% variance with LOD 1.83. Only one QTL (QTL.Fe.9) for Fe content was located on chromosome 9 (LOD 1.24), and two (QTL.Zn.4 and QTL.Zn.5) for Zn on chromosome 4 (LOD 1.71) and 5 (LOD 1.18), respectively. Overall, germplasm from clusters III and IV might offer higher heterotic response with the identified QTLs playing a significant role in any rice biofortification breeding program and released with development of new varieties.

Variability Assessment of Aromatic Rice Germplasm by Pheno-Genomic traits and Population Structure Analysis.

While the pleasant scent of aromatic rice is making it more popular, with demand for aromatic rice expected to rise in future, varieties of this have low yield potential. Genetic diversity and population structure of aromatic germplasm provide valuable information for yield improvement which has potential market value and farm profit. Here, we show diversity and population structure of 113 rice germplasm based on phenotypic and genotypic traits. Phenotypic traits showed that considerable variation existed across the germplasm. Based on Shannon-Weaver index, the most variable phenotypic trait was lemma-palea color. Detecting 140 alleles, 11 were unique and suitable as a germplasm diagnostic tool. Phylogenetic cluster analysis using genotypic traits classified germplasm into three major groups. Moreover, model-based population structure analysis divided all germplasm into three groups, confirmed by principal component and neighbors joining tree analyses. An analysis of molecular variance (AMOVA) and pairwise FST test showed significant differentiation among all population pairs, ranging from 0.023 to 0.068, suggesting that all three groups differed. Significant correlation coefficient was detected between phenotypic and genotypic traits which could be valuable to select further improvement of germplasm. Findings from this study have the potential for future use in aromatic rice molecular breeding programs.

Diversity and population structure of red rice germplasm in Bangladesh.

While the functionality and healthy food value of red rice have increased its popularity, such that market demand for it is expected to rise, most strains suffer from low grain yield. To perform diversity and population structure analyses of red rice germplasm, therefore, becomes essential for improving yields for commercial production. In this study, fifty red rice germplasm from the Bangladesh Rice Research Institute (BRRI) genebank were characterized both morphologically and genetically using fifty simple sequence repeat (SSR) markers. Overall, 162 alleles were detected by the markers with the detected allele number varying from two to seven. Additionally, 22 unique alleles were identified for use as a germplasm diagnostic tool. The highest and lowest polymorphic information content (PIC) indices were 0.75 and 0.04 found in markers RM282 and RM304, respectively, and genetic diversity was moderate, varying from 0.05 to 0.78 (average: 0.35). While phylogenetic cluster analysis of the fifteen distance-based agro-morphological traits divided the germplasm into five clusters (I, II, III, IV and V), a similar SSR analysis yielded only three major groups (I, II, and III), and a model-based population structure analysis yielded four (A, B, C and D). Both principal component and neighbors joining tree analysis from the population structure method showed the tested germplasm as highly diverse in structure. Moreover, an analysis of molecular variance (AMOVA), as well as a pairwise FST analysis, both indicated significant differentiation (ranging from 0.108 to 0.207) among all pairs of populations, suggesting that all four population structure groups differed significantly. Populations A and D were the most differentiated from each other by FST. Findings from this study suggest that the diverse germplasm and polymorphic trait-linked SSR markers of red rice are suitable for the detection of economically desirable trait loci/genes for use in future molecular breeding programs.

Comparison of the proliferation and excretion of Bartonella quintana between body and head lice following oral challenge.

Human body and head lice are highly related haematophagous ectoparasites but only the body louse has been shown to transmit Bartonella quintana, the causative agent of trench fever. The mechanisms by which body lice became a vector for B. quintana, however, are poorly understood. Following oral challenge, green fluorescent protein-expressing B. quintana proliferated over 9 days postchallenge with the number of bacteria being significantly higher in whole body vs. head lice. The numbers of B. quintana detected in faeces from infected lice, however, were approximately the same in both lice. Nevertheless, the viability of B. quintana was significantly higher in body louse faeces. Comparison of immune responses in alimentary tract tissues revealed that basal transcription levels of peptidoglycan recognition protein and defensins were lower in body lice and the transcription of defensin 1 was up-regulated by oral challenge with wild-type B. quintana in head but not in body lice. In addition, the level of cytotoxic reactive oxygen species generated by epithelial cells was significantly lower in body lice. Although speculative at this time, the reduced immune response is consistent with the higher vector competence seen in body vs. head lice in terms of B. quintana infection.

Differential gene expression in laboratory strains of human head and body lice when challenged with Bartonella quintana, a pathogenic bacterium.

Human head and body lice are obligatory hematophagous ectoparasites that belong to a single species, Pediculus humanus. Only body lice, however, are vectors of the infectious Gram-negative bacterium Bartonella quintana. Because of their near identical genomes, yet differential vector competence, head and body lice provide a unique model system to study the gain or loss of vector competence. Using our in vitro louse-rearing system, we infected head and body lice with blood containing B. quintana in order to detect both differences in the proliferation of B. quintana and transcriptional differences of immune-related genes in the lice. B. quintana proliferated rapidly in body lice at 6 days post-infection, but plateaued in head lice at 4 days post-infection. RNAseq and quantitative real-time PCR validation analyses determined gene expression differences. Eight immunoresponse genes were observed to be significantly different with many associated with the Toll pathway: Fibrinogen-like protein, Spaetzle, Defensin 1, Serpin, Scavenger receptor A and Apolipoporhrin 2. Our findings support the hypothesis that body lice, unlike head lice, fight infection from B. quintana only at the later stages of its proliferation.

Modeling evolution of resistance by Maruca vitrata (Lepidoptera: Crambidae) to transgenic insecticidal cowpea in Africa.

We created a detailed model of the Maruca vitrata (F.) and cowpea [Vigna unguiculata (L.) Walp] system to study the possible evolution of resistance by the insect to transgenic insecticidal cowpea, which is under development. We focused on population dynamics and genetics in a region of west Africa. We simulated single-toxin and pyramided (two-toxin) cowpea and emphasized conservative, worst-case scenarios in our analysis. The results indicate that as long as a pyramided, transgenic cowpea can be developed, seed saving by farmers and reliance on natural refuge are not major problems for resistance management. Furthermore, it is possible that one or both toxins in the pyramid may not need to be high dose for evolution to be delayed significantly (>20 yr or 80 generations for resistance to become a concern if transgenic cowpea is deployed in areas where M. vitrata is endemic). If efforts are made to deploy transgenic cowpea only into the regions where M. vitrata is not endemic, then there is little to no concern with resistance emerging in the M. vitrata population.

The spatial genetic differentiation of the legume pod borer, Maruca vitrata F. (Lepidoptera: Crambidae) populations in West Africa.

The legume pod borer, Maruca vitrata, is an endemic insect pest that causes significant yield loss to the cowpea crop in West Africa. The application of population genetic tools is important in the management of insect pests but such data on M. vitrata is lacking. We applied a set of six microsatellite markers to assess the population structure of M. vitrata collected at five sites from Burkina Faso, Niger and Nigeria. Observed polymorphisms ranged from one (marker 3393) to eight (marker 32008) alleles per locus. Observed and expected heterozygosities ranged from 0.0 to 0.8 and 0.0 to 0.6, respectively. Three of the loci in samples from Nigeria and Burkina Faso deviated significantly from Hardy-Weinberg Equilibrium (HWE), whereas no loci deviated significantly in samples from Niger. Analysis of molecular variance (AMOVA) indicated that 67.3% level of the genetic variation was within individuals compared to 17.3% among populations. A global estimate of F ST=0.1 (ENA corrected F ST=0.1) was significant (P⩽0.05) and corroborated by pairwise F ST values that were significant among all possible comparisons. A significant correlation was predicted between genetic divergence and geographic distance between subpopulations (R2=0.6, P=0.04), and cluster analysis by the program STRUCTURE predicted that co-ancestry of genotypes were indicative of three distinct populations. The spatial genetic variance among M. vitrata in West Africa may be due to limited gene flow, south-north seasonal movement pattern or other reproductive barriers. This information is important for the cultural, chemical and biological control strategies for managing M. vitrata.

Comparison of the transcriptional profiles of head and body lice.

Head and body lice are both blood-feeding parasites of humans although only the body louse is a potent disease vector. In spite of numerous morphological and life history differences, head and body lice have recently been hypothesized to be ecotypes of the same species. We took a comparative genomics approach to measure nucleotide diversity by comparing expressed sequence tag data sets from head and body lice. A total of 10 771 body louse and 10 770 head louse transcripts were predicted from a combined assembly of Roche 454 and Illumina sequenced cDNAs from whole body tissues collected at all life stages and during pesticide exposure and bacterial infection treatments. Illumina reads mapped to the 10 775 draft body louse gene models from the whole genome assembly predicted nine presence/absence differences, but PCR confirmation resulted in a single gene difference. Read per million base pair estimates indicated that 14 genes showed significant differential expression between head and body lice under our treatment conditions. One novel microRNA was predicted in both lice species and 99% of the 544 transcripts from Candidatus riesia indicate that they share the same endosymbiont. Overall, few differences exist, which supports the hypothesis that these two organisms are ecotypes of the same species.

Brief exposures of human body lice to sublethal amounts of ivermectin over-transcribes detoxification genes involved in tolerance.

Transcriptional profiling results, using our non-invasive induction assay {short exposure intervals (2-5 h) to sublethal amounts of insecticides [< lethal concentration 3% (LC(3)) at 24 h] administered by stress-reducing means (contact vs. immersion screen) and with induction assessed in a time frame when tolerance is still present [~lethal concentration 90% (LC(90)) in 2-4 h]}, showed that ivermectin-induced detoxification genes from body lice are identified by quantitative real-time PCR analyses. Of the cytochrome P450 monooxygenase and ATP binding cassette transporter genes induced by ivermectin, CYP6CJ1, CYP9AG1, CYP9AG2 and PhABCC4 were respectively most significantly over-expressed, had high basal expression levels and were most closely related to genes from other organisms that metabolized insecticides, including ivermectin. Injection of double-stranded RNAs (dsRNAs) against either CYP9AG2 or PhABCC4 into non-induced female lice reduced their respective transcript level and resulted in increased sensitivity to ivermectin, indicating that these two genes are involved in the xenobiotic metabolism of ivermectin and in the production of tolerance.

Decreased detoxification genes and genome size make the human body louse an efficient model to study xenobiotic metabolism.

The human body louse, Pediculus humanus humanus, has one of the smallest insect genomes, containing ∼10 775 annotated genes. Annotation of detoxification [cytochrome P450 monooxygenase (P450), glutathione-S-transferase (GST), esterase (Est) and ATP-binding cassette transporter (ABC transporter)] genes revealed that they are dramatically reduced in P. h. humanus compared to other insects except for Apis mellifera. There are 37 P450, 13 GST and 17 Est genes present in P. h. humanus, approximately half the number found in Drosophila melanogaster and Anopheles gambiae. The number of putatively functional ABC transporter genes in P. h. humanus and Ap. mellifera are the same (36) but both have fewer than An. gambiae (44) or Dr. melanogaster (65). The reduction of detoxification genes in P. h. humanus may be a result of this louse's simple life history, in which it does not encounter a wide variety of xenobiotics. Neuronal component genes are highly conserved across different insect species as expected because of their critical function. Although reduced in number, P. h. humanus still retains at least a minimum repertoire of genes known to confer metabolic or toxicokinetic resistance to xenobiotics (eg Cyp3 clade P450s, Delta GSTs, B clade Ests and B/C subfamily ABC transporters), suggestive of its high potential for resistance development.

Bowman-Birk inhibitor affects pathways associated with energy metabolism in Drosophila melanogaster.

Bowman-Birk inhibitor (BBI) is toxic when fed to certain insects, including the fruit fly, Drosophila melanogaster. Dietary BBI has been demonstrated to slow growth and increase insect mortality by inhibiting the digestive enzymes trypsin and chymotrypsin, resulting in a reduced supply of amino acids. In mammals, BBI influences cellular energy metabolism. Therefore, we tested the hypothesis that dietary BBI affects energy-associated pathways in the D. melanogaster midgut. Through microarray and metabolomic analyses, we show that dietary BBI affects energy utilization pathways in the midgut cells of D. melanogaster. In addition, ultrastructure studies indicate that microvilli are significantly shortened in BBI-fed larvae. These data provide further insights into the complex cellular response of insects to dietary protease inhibitors.

Transcriptional signatures in response to wheat germ agglutinin and starvation in Drosophila melanogaster larval midgut.

One function of plant lectins such as wheat germ agglutinin is to serve as defences against herbivorous insects. The midgut is one critical site affected by dietary lectins. We observed marked cellular, structural and gene expression changes in the midguts of Drosophila melanogaster third instar larvae that were fed wheat germ agglutinin. Some of these changes were similar to those observed in the midguts of starved D. melanogaster. Dietary wheat germ agglutinin caused shortening, branching, swelling, distortion and in some cases disintegration of the midgut microvilli. Starvation was accompanied primarily by shortening of the microvilli. Microarray analyses revealed that dietary wheat germ agglutinin evoked differential expression of 61 transcripts; seven of these were also differentially expressed in starved D. melanogaster. The differentially transcribed gene clusters in wheat germ agglutinin-fed larvae were associated with (1) cytoskeleton organization; (2) digestive enzymes; (3) detoxification reactions; and (4) energy metabolism. Four possible transcription factor binding motifs were associated with the differentially expressed genes. One of these exhibited substantial similarity to MyoD, a transcription factor binding motif associated with cellular structures in mammals. These results are consistent with the hypothesis that wheat germ agglutinin caused a starvation-like effect and structural changes of midgut cells of D. melanogaster third-instar larvae.

A mixture model approach for the analysis of small exploratory microarray experiments.

The microarray is an important and powerful tool for prescreening of genes for further research. However, alternative solutions are needed to increase power in small microarray experiments. Use of traditional parametric and even non-parametric tests for such small experiments lack power and have distributional problems. A mixture model is described that is performed directly on expression differences assuming that genes in alternative treatments are expressed or not in all combinations (i) not expressed in either condition, (ii) expressed only under the first condition, (iii) expressed only under the second condition, and (iv) expressed under both conditions, giving rise to 4 possible clusters with two treatments. The approach is termed a Mean-Difference-Mixture-Model (MD-MM) method. Accuracy and power of the MD-MM was compared to other commonly used methods, using both simulations, microarray data, and quantitative real time PCR (qRT-PCR). The MD-MM was found to be generally superior to other methods in most situations. The advantage was greatest in situations where there were few replicates, poor signal to noise ratios, or non-homogenous variances.

Transcriptomic profiles of Drosophila melanogaster third instar larval midgut and responses to oxidative stress.

Oligoarray analysis was used to determine the number and nature of genes expressed in third instar Drosophila melanogaster larval midguts. The majority of transcripts were associated with protein synthesis and metabolism. Serine proteases were the main proteolytic enzymes detected. Some 40% of the cytochrome P450 genes and 74% of the glutathione S transferases (GSTs) in the genome of D. melanogaster were observed to be expressed in the midgut by oligoarray analysis. We also identified potential transcription factor binding motifs (TFBMs) of P450s, GSTs and carboxylesterases. Many of the midgut-expressed GST genes contained candidate TFBMs homologous to TFBMs in mammals that have been associated with responses to oxidative stress. We also investigated the response of GSTs in the midgut to dietary H2O2, which showed a dosage-based differential response.

Changes in Drosophila melanogaster midgut proteins in response to dietary Bowman-Birk inhibitor.

The midgut proteome of Drosophila melanogaster was compared in larvae fed dietary Bowman-Birk inhibitor (BBI) vs. larvae fed a control diet. By using two-dimensional gel electrophoresis, nine differentially expressed proteins were observed, which were associated with enzymes or transport functions such as sterol carrier protein X (SCPX), ubiquitin-conjugating enzyme, endopeptidase, receptor signalling protein kinase, ATP-dependent RNA helicase and alpha-tocopherol transport. Quantitative real-time PCR verified differential expression of transcripts coding for six of the proteins observed from the proteomic analysis. BBI evidently affects expression of proteins associated with protein degradation, transport and fatty acid catabolism. We then tested the hypothesis that SCPX was critical for the Drosophila third instars' response to BBI treatment. Inhibition of SCPX caused the third instars to become more susceptible to dietary BBI.

Sequencing of a new target genome: the Pediculus humanus humanus (Phthiraptera: Pediculidae) genome project.

The human body louse, Pediculus humanus humanus (L.), and the human head louse, Pediculus humanus capitis, belong to the hemimetabolous order Phthiraptera. The body louse is the primary vector that transmits the bacterial agents of louse-borne relapsing fever, trench fever, and epidemic typhus. The genomes of the bacterial causative agents of several of these aforementioned diseases have been sequenced. Thus, determining the body louse genome will enhance studies of host-vector-pathogen interactions. Although not important as a major disease vector, head lice are of major social concern. Resistance to traditional pesticides used to control head and body lice have developed. It is imperative that new molecular targets be discovered for the development of novel compounds to control these insects. No complete genome sequence exists for a hemimetabolous insect species primarily because hemimetabolous insects often have large (2000 Mb) to very large (up to 16,300 Mb) genomes. Fortuitously, we determined that the human body louse has one of the smallest genome sizes known in insects, suggesting it may be a suitable choice as a minimal hemimetabolous genome in which many genes have been eliminated during its adaptation to human parasitism. Because many louse species infest birds and mammals, the body louse genome-sequencing project will facilitate studies of their comparative genomics. A 6-8X coverage of the body louse genome, plus sequenced expressed sequence tags, should provide the entomological, evolutionary biology, medical, and public health communities with useful genetic information.

Genome-wide analysis of phenobarbital-inducible genes in Drosophila melanogaster.

An oligoarray analysis was conducted to determine the differential expression of genes due to phenobarbital exposure in Drosophila melanogaster (w(1118) strain) third instar larvae. Seventeen genes were observed to be induced with increased expression by a statistical analysis of microarrays approach with a q < or = 0.05. At q < or = 0.12, four more genes (Cyp12d1, DmGstd4, and two genes with unknown function) were found to be up-regulated, and 11 genes with unknown function were found to be down-regulated. Fifteen of these genes, Cyp4d14, Cyp6a2, Cyp6a8, Cyp12d1, Cyp6d5, Cyp6w1, CG2065, DmGstd6, DmGstd7, Amy-p/Amy-d, Ugt86Dd, GC5724, Jheh1, Jheh2 and CG11893, were verified using quantitative real time polymerase chain reaction. Some of these genes have been shown to be over-transcribed in metabolically DDT-resistant Drosophila strains.

Gene expression profiles among immature and adult reproductive castes of the termite Reticulitermes flavipes.

Array-based genomic studies were conducted with the goal of identifying immature (i.e. nymph) and adult reproductive caste-biased gene expression in the termite Reticulitermes flavipes. Using cDNA macro-arrays, we identified thirty-four nymph-biased genes falling into eight ontogenic categories. Based on gene expression profiles among diverse castes and developmental stages (determined by quantitative PCR), several important trends emerged. These findings highlight the importance of several developmental and survival-based factors among immature and adult termite reproductives, including: vitellogenesis, nutrient storage, juvenile hormone sequestration, ribosomal translational and filtering mechanisms, fatty acid biosynthesis, apoptosis inhibition, and both endogenous and symbiont cellulase-assisted nutrition. These findings are highly significant as they are the first to elucidate the molecular biology underlying termite reproductive caste differentiation and reproductive caste-specific biology. Other gene expression results are in agreement with previous findings that suggest roles for vitellogenin-like haemolymph proteins in soldier caste differentiation.

Expression of Cyp6g1 and Cyp12d1 in DDT resistant and susceptible strains of Drosophila melanogaster.

The Rst(2)DDT locus (loci) in Drosophila is associated with the over-expression of two cytochrome P450 genes, Cyp6g1 and Cyp12d1. Using northern and western blot analysis we observed the expression pattern of these two genes in two DDT susceptible (Canton-S and 91-C) and three DDT resistant strains (Wisconsin, 91-R and Hikone-R). In Canton-S and 91-R, the CYP6G1 protein was constitutively expressed throughout development. In the Wisconsin strain, CYP6G1 was not expressed in third instar larvae unless the larvae are exposed to DDT. CYP12D1 protein was only expressed in adults. Cyp12d1 mRNA is induced in DDT resistant strains post-exposure to DDT and the expression patterns of Cyp12d1 mRNA varied across DDT resistant strains. Our data support the hypothesis that there is evolutionary plasticity in the expression patterns of P450s associated with metabolic pesticide resistance.

Transcriptome analysis of the cowpea weevil bruchid: identification of putative proteinases and alpha-amylases associated with food breakdown.

We describe here the first systematic work to discover insect genes involved in food breakdown using a cDNA library enriched for gut-expressed transcripts from Callosobruchus maculatus. A total of 1056 clones were screened for cDNA insert-containing plasmids, and 503 nonredundant open reading frames were discovered. Twenty-three inferred genes potentially involved in digestive processes in cowpea weevil were identified, including proteinases and amylases. The predicted catalytic sites were identified in the inferred cysteine and aspartic acid proteinases, and in alpha-amylases. Transcriptome analysis of the cowpea bruchid will potentially permit gene discovery in other beetles, an insect order of major economic and ecological importance that is poorly represented in genomic databases.