ABSTRACT
Insect crop pests threaten global food security. This threat is amplified through the spread of nonnative species and through adaptation of native pests to control measures. Adaptations such as pesticide resistance can result from selection on variation within a population, or through gene flow from another population. We investigate these processes in an economically important noctuid crop pest, Helicoverpa zea, which has evolved resistance to a wide range of pesticides. Its sister species Helicoverpa armigera, first detected as an invasive species in Brazil in 2013, introduced the pyrethroid-resistance gene CYP337B3 to South American H. zea via adaptive introgression. To understand whether this could contribute to pesticide resistance in North America, we sequenced 237 H. zea genomes across 10 sample sites. We report H. armigera introgression into the North American H. zea population. Two individuals sampled in Texas in 2019 carry H. armigera haplotypes in a 4 Mbp region containing CYP337B3. Next, we identify signatures of selection in the panmictic population of nonadmixed H. zea, identifying a selective sweep at a second cytochrome P450 gene: CYP333B3. We estimate that its derived allele conferred a â¼5% fitness advantage and show that this estimate explains independently observed rare nonsynonymous CYP333B3 mutations approaching fixation over a â¼20-year period. We also detect putative signatures of selection at a kinesin gene associated with Bt resistance. Overall, we document two mechanisms of rapid adaptation: the introduction of fitness-enhancing alleles through interspecific introgression, and selection on intraspecific variation.
Subject(s)
Genetic Introgression , Insecticide Resistance , Moths , Animals , Moths/genetics , Insecticide Resistance/genetics , Cytochrome P-450 Enzyme System/genetics , North America , Adaptation, Biological/genetics , Adaptation, Physiological/genetics , Selection, Genetic , Introduced SpeciesABSTRACT
Males of some species, from horses to humans, require medical help for subfertility problems. There is an urgent need for novel molecular assays that reflect spermatozoal function. In the last 25 years, studies examined RNAs in spermatozoa as a window into gene expression during their development and, more recently, for their functions in early embryo development. In clinics, more dense spermatozoa are isolated by density gradient centrifugation before use in artificial insemination to increase pregnancy rates. The objectives of the current study were to discover and quantify the microRNAs in stallion spermatozoa and identify those with differential expression levels in more dense versus less dense spermatozoa. First, spermatozoa from seven stallions were separated into more dense and less dense populations by density gradient centrifugation. Next, small RNAs were sequenced from each of the 14 RNA samples. We identified 287 different mature microRNAs within the 11,824,720 total mature miRNA reads from stallion spermatozoa. The most prevalent was miR-10a/b-5p. The less dense spermatozoa had fewer mature microRNAs and more microRNA precursor sequences than more dense spermatozoa, perhaps indicating that less dense spermatozoa are less mature. Two of the most prevalent microRNAs in more dense stallion spermatozoa were predicted to target mRNAs that encode proteins that accelerate mRNA decay. Nine microRNAs were more highly expressed in more dense spermatozoa. Three of those microRNAs were predicted to target mRNAs that encode proteins involved in protein decay. Both mRNA and protein decay are very active in late spermiogenesis but not in mature spermatozoa. The identified microRNAs may be part of the mechanism to shut down those processes. The microRNAs with greater expression in more dense spermatozoa may be useful biomarkers for spermatozoa with greater functional capabilities.
Subject(s)
Biomarkers , MicroRNAs , Spermatozoa , Horses , Male , Animals , Spermatozoa/physiology , Spermatozoa/metabolism , MicroRNAs/genetics , MicroRNAs/metabolism , Biomarkers/metabolismABSTRACT
BACKGROUND: Rice is one of the most important cereals consumed worldwide. Two major abiotic factors affecting rice plants in different growth stages are flooding stress and cold stress. These abiotic stresses can take place independently or simultaneously and significantly affect rice plants during germination and seedling growth. Fortunately, a wide array of phenotypic responses conferring flooding stress and chilling stress tolerance exist within the rice germplasm, indicating the presence of different molecular mechanisms underlying tolerance to these stresses. Understanding these differences may assist in developing improved rice cultivars having higher tolerance to both stresses. In this study, we conducted a comparative global gene expression analysis of two rice genotypes with contrasting phenotypes under cold stress, anaerobic stress, and combined cold and anaerobic stress during germination. RESULTS: The differential gene expression analysis revealed that 5571 differentially expressed genes (DEGs), 7206 DEGs, and 13279 DEGs were identified under anaerobic stress, cold stress, and combined stress, respectively. Genes involved in the carbohydrate metabolic process, glucosyltransferase activity, regulation of nitrogen compound metabolic process, protein metabolic process, lipid metabolic process, cellular nitrogen compound biosynthetic process, lipid biosynthetic process, and a microtubule-based process were enriched across all stresses. Notably, the common Gene Ontology (GO) analysis identified three hub genes, namely Os08g0176800 (similar to mRNA-associated protein mrnp 41), Os11g0454200 (dehydrin), and OS10g0505900 (expressed protein). CONCLUSION: A large number of differentially expressed genes were identified under anaerobic, cold conditions during germination and the combination of the two stress conditions in rice. These results will assist in the identification of promising candidate genes for possible manipulation toward rice crops that are more tolerant under flooding and cold during germination, both independently and concurrently.
Subject(s)
Oryza , Transcriptome , Seedlings , Cold-Shock Response/genetics , Anaerobiosis , Gene Expression Profiling , Nitrogen/metabolism , Lipids , Gene Expression Regulation, Plant , Cold TemperatureABSTRACT
KEY MESSAGE: A novel locus was discovered on chromosome 7 associated with a lesion mimic in maize; this lesion mimic had a quantitative and heritable phenotype and was predicted better via subset genomic markers than whole genome markers across diverse environments. Lesion mimics are a phenotype of leaf micro-spotting in maize (Zea mays L.), which can be early signs of biotic or abiotic stresses. Dissecting its inheritance is helpful to understand how these loci behave across different genetic backgrounds. Here, 538 maize recombinant inbred lines (RILs) segregating for a novel lesion mimic were quantitatively phenotyped in Georgia, Texas, and Wisconsin. These RILs were derived from three bi-parental crosses using a tropical pollinator (Tx773) as the common parent crossed with three inbreds (LH195, LH82, and PB80). While this lesion mimic was heritable across three environments based on phenotypic ([Formula: see text] = 0.68) and genomic ([Formula: see text] = 0.91) data, transgressive segregation was observed. A genome-wide association study identified a single novel locus on chromosome 7 (at 70.6 Mb) also covered by a quantitative trait locus interval (69.3-71.0 Mb), explaining 11-15% of the variation, depending on the environment. One candidate gene identified in this region, Zm00001eb308070, is related to the abscisic acid pathway involving in cell death. Genomic predictions were applied to genome-wide markers (39,611 markers) contrasted with a marker subset (51 markers). Population structure explained more variation than environment in genomic prediction, but other substantial genetic background effects were additionally detected. Subset markers explained substantially less genetic variation (24.9%) for the lesion mimic than whole genome markers (55.4%) in the model, yet predicted the lesion mimic better (0.56-0.66 vs. 0.26-0.29). These results indicate this lesion mimic phenotype was less affected by environment than by epistasis and genetic background effects, which explain its transgressive segregation.
Subject(s)
Genome-Wide Association Study , Zea mays , Zea mays/genetics , Epistasis, Genetic , Chromosome Mapping , Phenotype , Genetic Background , Polymorphism, Single NucleotideABSTRACT
Using imbalanced historical yield data to predict performance and select new lines is an arduous breeding task. Genome-wide association studies (GWAS) and high throughput genotyping based on sequencing techniques can increase prediction accuracy. An association mapping panel of 227 Texas elite (TXE) wheat breeding lines was used for GWAS and a training population to develop prediction models for grain yield selection. An imbalanced set of yield data collected from 102 environments (year-by-location) over 10 years, through testing yield in 40-66 lines each year at 6-14 locations with 38-41 lines repeated in the test in any two consecutive years, was used. Based on correlations among data from different environments within two adjacent years and heritability estimated in each environment, yield data from 87 environments were selected and assigned to two correlation-based groups. The yield best linear unbiased estimation (BLUE) from each group, along with reaction to greenbug and Hessian fly in each line, was used for GWAS to reveal genomic regions associated with yield and insect resistance. A total of 74 genomic regions were associated with grain yield and two of them were commonly detected in both correlation-based groups. Greenbug resistance in TXE lines was mainly controlled by Gb3 on chromosome 7DL in addition to two novel regions on 3DL and 6DS, and Hessian fly resistance was conferred by the region on 1AS. Genomic prediction models developed in two correlation-based groups were validated using a set of 105 new advanced breeding lines and the model from correlation-based group G2 was more reliable for prediction. This research not only identified genomic regions associated with yield and insect resistance but also established the method of using historical imbalanced breeding data to develop a genomic prediction model for crop improvement. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01287-8.
ABSTRACT
Wheat flour's end-use quality is tightly linked to the quantity and composition of storage proteins in the endosperm. TAM 111 and TAM 112 are two popular cultivars grown in the Southern US Great Plains with significantly different protein content. To investigate regulatory differences, transcriptome data were analyzed from developing grains at early- and mid-filling stages. At the mid-filling stage, TAM 111 preferentially upregulated starch metabolism-related pathways compared to TAM 112, whereas amino acid metabolism and transporter-related pathways were over-represented in TAM 112. Elemental analyses also indicated a higher N percentage in TAM 112 at the mid-filling stage. To explore the regulatory variation, weighted correlation gene network was constructed from publicly available RNAseq datasets to identify the modules differentially regulated in TAM 111 and TAM 112. Further, the potential transcription factors (TFs) regulating those modules were identified using graphical least absolute shrinkage and selection operator (GLASSO). Homologs of the OsNF-Y family members with known starch metabolism-related functions showed higher connectivities in TAM 111. Multiple TFs with high connectivity in TAM 112 had predicted functions associated with ABA response in grain. These results will provide novel targets for breeders to explore and further our understanding in mechanisms regulating grain development.
Subject(s)
Plant Proteins , Triticum , Triticum/metabolism , Plant Proteins/metabolism , Flour , Gene Expression Profiling , Edible Grain/metabolism , Transcriptome , Transcription Factors/metabolism , Starch/metabolism , Amino Acids/metabolism , Gene Expression Regulation, PlantABSTRACT
Translational control during cell division determines when cells start a new cell cycle, how fast they complete it, the number of successive divisions, and how cells coordinate proliferation with available nutrients. The translational efficiencies of mRNAs in cells progressing synchronously through the mitotic cell cycle, while preserving the coupling of cell division with cell growth, remain uninvestigated. We now report comprehensive ribosome profiling of a yeast cell size series from the time of cell birth, to identify mRNAs under periodic translational control. The data reveal coordinate translational activation of mRNAs encoding lipogenic enzymes late in the cell cycle including Acc1p, the rate-limiting enzyme acetyl-CoA carboxylase. An upstream open reading frame (uORF) confers the translational control of ACC1 and adjusts Acc1p protein levels in different nutrients. The ACC1 uORF is relevant for cell division because its ablation delays cell cycle progression, reduces cell size, and suppresses the replicative longevity of cells lacking the Sch9p protein kinase regulator of ribosome biogenesis. These findings establish an unexpected relationship between lipogenesis and protein synthesis in mitotic cell divisions.
Subject(s)
Acetyl-CoA Carboxylase/biosynthesis , Gene Expression Regulation, Fungal , Mitosis , Protein Biosynthesis , Yeasts/growth & development , Yeasts/genetics , Acetyl-CoA Carboxylase/genetics , Lipid Metabolism , Open Reading Frames , Ribosomes/metabolism , Yeasts/metabolismABSTRACT
Toll-like receptor-4 (TLR4) is a transmembrane receptor that initiates an immune response following a bacterial infection or host derived molecules associated with cellular distress. Beyond triggering inflammation, TLR4 has been implicated in modulating behavioral and cognitive processes in a physiologically normal state, as young adult TLR4 deficient mice show learning enhancements in select tasks. Currently unknown is whether these benefits are present in both sexes and persist with aging. The present study evaluated spatial memory, anxiety-like behavior, and central levels of pro- and anti-inflammatory molecules in young (4-5â¯months) and aged (18-19â¯months) TLR4 deficient (TLR4-/-) and wild-type (WT) male and female mice. Results confirmed that TLR4-/- mice show enhanced spatial memory compared to WT mice. These effects were age- and sex-specific, as memory retention was superior in the TLR4-/- young males and aged females. While TLR4-/- mice showed age-related changes in behavior, these changes were attenuated relative to aged WT mice. Further, aged TLR4-/- mice showed differential expression of molecules involved in interleukin (IL)-1 signaling in the hippocampus. For instance, aged TLR4-/- females showed heightened expression of IL-1 receptor antagonist (IL-1ra) and the IL-1 accessory proteins AcP and AcPb. Collectively, these data provide the initial evidence that TLR4 deficiency enhances cognitive function and modulates the inflammatory profile of the hippocampus in a sex- and age-dependent manner.
Subject(s)
Interleukin-1/genetics , Spatial Memory/physiology , Toll-Like Receptor 4/metabolism , Age Factors , Animals , Anxiety/genetics , Anxiety/metabolism , Cognition/physiology , Female , Hippocampus/metabolism , Inflammation/metabolism , Interleukin 1 Receptor Antagonist Protein/metabolism , Interleukin-1/metabolism , Lipopolysaccharides/metabolism , Male , Maze Learning/physiology , Mice , Mice, Inbred C57BL , Mice, Knockout , Receptors, Interleukin-1/metabolism , Sex Factors , Signal Transduction/physiology , Toll-Like Receptor 4/geneticsABSTRACT
BACKGROUND: The information content of genomes plays a crucial role in the existence and proper development of living organisms. Thus, tremendous effort has been dedicated to developing DNA sequencing technologies that provide a better understanding of the underlying mechanisms of cellular processes. Advances in the development of sequencing technology have made it possible to sequence genomes in a relatively fast and inexpensive way. However, as with any measurement technology, there is noise involved and this needs to be addressed to reach conclusions based on the resulting data. In addition, there are multiple intermediate steps and degrees of freedom when constructing genome assemblies that lead to ambiguous and inconsistent results among assemblers. METHODS: Here we introduce HiMMe, an HMM-based tool that relies on genetic patterns to score genome assemblies. Through a Markov chain, the model is able to detect characteristic genetic patterns, while, by introducing emission probabilities, the noise involved in the process is taken into account. Prior knowledge can be used by training the model to fit a given organism or sequencing technology. RESULTS: Our results show that the method presented is able to recognize patterns even with relatively small k-mer size choices and limited computational resources. CONCLUSIONS: Our methodology provides an individual quality metric per contig in addition to an overall genome assembly score, with a time complexity well below that of an aligner. Ultimately, HiMMe provides meaningful statistical insights that can be leveraged by researchers to better select contigs and genome assemblies for downstream analysis.
Subject(s)
Genomics/methods , Markov Chains , Algorithms , Bayes Theorem , Reproducibility of ResultsABSTRACT
BACKGROUND: The soybean cyst nematode (SCN), Heterodera glycines, is one of the most devastating diseases limiting soybean production worldwide. It is known that small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), play important roles in regulating plant growth and development, defense against pathogens, and responses to environmental changes. RESULTS: In order to understand the role of soybean miRNAs during SCN infection, we analyzed 24 small RNA libraries including three biological replicates from two soybean cultivars (SCN susceptible KS4607, and SCN HG Type 7 resistant KS4313N) that were grown under SCN-infested and -noninfested soil at two different time points (SCN feeding establishment and egg production). In total, 537 known and 70 putative novel miRNAs in soybean were identified from a total of 0.3 billion reads (average about 13.5 million reads for each sample) with the programs of Bowtie and miRDeep2 mapper. Differential expression analyses were carried out using edgeR to identify miRNAs involved in the soybean-SCN interaction. Comparative analysis of miRNA profiling indicated a total of 60 miRNAs belonging to 25 families that might be specifically related to cultivar responses to SCN. Quantitative RT-PCR validated similar miRNA interaction patterns as sequencing results. CONCLUSION: These findings suggest that miRNAs are likely to play key roles in soybean response to SCN. The present work could provide a framework for miRNA functional identification and the development of novel approaches for improving soybean SCN resistance in future studies.
Subject(s)
Genomics , Glycine max/genetics , High-Throughput Nucleotide Sequencing , MicroRNAs/genetics , Plant Diseases , Sequence Analysis, RNA , Tylenchoidea/physiology , Animals , Glycine max/physiologyABSTRACT
BACKGROUND: Colletotrichum graminicola is a hemibiotrophic fungal pathogen that causes maize anthracnose disease. It progresses through three recognizable phases of pathogenic development in planta: melanized appressoria on the host surface prior to penetration; biotrophy, characterized by intracellular colonization of living host cells; and necrotrophy, characterized by host cell death and symptom development. A "Mixed Effects" Generalized Linear Model (GLM) was developed and applied to an existing Illumina transcriptome dataset, substantially increasing the statistical power of the analysis of C. graminicola gene expression during infection and colonization. Additionally, the in planta transcriptome of the wild-type was compared with that of a mutant strain impaired in the establishment of biotrophy, allowing detailed dissection of events occurring specifically during penetration, and during early versus late biotrophy. RESULTS: More than 2000 fungal genes were differentially transcribed during appressorial maturation, penetration, and colonization. Secreted proteins, secondary metabolism genes, and membrane receptors were over-represented among the differentially expressed genes, suggesting that the fungus engages in an intimate and dynamic conversation with the host, beginning prior to penetration. This communication process probably involves reception of plant signals triggering subsequent developmental progress in the fungus, as well as production of signals that induce responses in the host. Later phases of biotrophy were more similar to necrotrophy, with increased production of secreted proteases, inducers of plant cell death, hydrolases, and membrane bound transporters for the uptake and egress of potential toxins, signals, and nutrients. CONCLUSIONS: This approach revealed, in unprecedented detail, fungal genes specifically expressed during critical phases of host penetration and biotrophic establishment. Many encoded secreted proteins, secondary metabolism enzymes, and receptors that may play roles in host-pathogen communication necessary to promote susceptibility, and thus may provide targets for chemical or biological controls to manage this important disease. The differentially expressed genes could be used as 'landmarks' to more accurately identify developmental progress in compatible versus incompatible interactions involving genetic variants of both host and pathogen.
Subject(s)
Colletotrichum/genetics , Plant Diseases/microbiology , Transcriptome , Colletotrichum/pathogenicity , Gene Expression Regulation, Fungal , Genes, Fungal , Host-Pathogen Interactions , Linear Models , RNA, Fungal/genetics , Secondary Metabolism , Sequence Analysis, RNA , Zea mays/microbiologyABSTRACT
Folate is a vitamin required for cell growth and is present in fortified foods in the form of folic acid to prevent congenital abnormalities. The impact of low folate status on life-long health is poorly understood. We found that limiting folate levels with the folate antagonist methotrexate increased the lifespan of yeast and worms. We then restricted folate intake in aged mice and measured various health metrics, metabolites, and gene expression signatures. Limiting folate intake decreased anabolic biosynthetic processes in mice and enhanced metabolic plasticity. Despite reduced serum folate levels in mice with limited folic acid intake, these animals maintained their weight and adiposity late in life, and we did not observe adverse health outcomes. These results argue that the effectiveness of folate dietary interventions may vary depending on an individual's age and sex. A higher folate intake is advantageous during the early stages of life to support cell divisions needed for proper development. However, a lower folate intake later in life may result in healthier aging.
ABSTRACT
Folate is a vitamin required for cell growth and is present in fortified foods in the form of folic acid to prevent congenital abnormalities. The impact of low-folate status on life-long health is poorly understood. We found that limiting folate levels with the folate antagonist methotrexate increased the lifespan of yeast and worms. We then restricted folate intake in aged mice and measured various health metrics, metabolites, and gene expression signatures. Limiting folate intake decreased anabolic biosynthetic processes in mice and enhanced metabolic plasticity. Despite reduced serum folate levels in mice with limited folic acid intake, these animals maintained their weight and adiposity late in life, and we did not observe adverse health outcomes. These results argue that the effectiveness of folate dietary interventions may vary depending on an individual's age and sex. A higher folate intake is advantageous during the early stages of life to support cell divisions needed for proper development. However, a lower folate intake later in life may result in healthier aging.
Subject(s)
Folic Acid , Longevity , Animals , Folic Acid/administration & dosage , Folic Acid/metabolism , Mice , Male , Female , Aging/metabolism , Diet/methods , Mice, Inbred C57BL , Methotrexate/pharmacology , Folic Acid Deficiency/metabolism , Caenorhabditis elegans , Saccharomyces cerevisiae/metabolismABSTRACT
BACKGROUND: A key goal of systems biology and translational genomics is to utilize high-throughput measurements of cellular states to develop expression-based classifiers for discriminating among different phenotypes. Recent developments of Next Generation Sequencing (NGS) technologies can facilitate classifier design by providing expression measurements for tens of thousands of genes simultaneously via the abundance of their mRNA transcripts. Because NGS technologies result in a nonlinear transformation of the actual expression distributions, their application can result in data that are less discriminative than would be the actual expression levels themselves, were they directly observable. RESULTS: Using state-of-the-art distributional modeling for the NGS processing pipeline, this paper studies how that pipeline, via the resulting nonlinear transformation, affects classification and feature selection. The effects of different factors are considered and NGS-based classification is compared to SAGE-based classification and classification directly on the raw expression data, which is represented by a very high-dimensional model previously developed for gene expression. As expected, the nonlinear transformation resulting from NGS processing diminishes classification accuracy; however, owing to a larger number of reads, NGS-based classification outperforms SAGE-based classification. CONCLUSIONS: Having high numbers of reads can mitigate the degradation in classification performance resulting from the effects of NGS technologies. Hence, when performing a RNA-Seq analysis, using the highest possible coverage of the genome is recommended for the purposes of classification.
Subject(s)
Genomics/methods , High-Throughput Nucleotide Sequencing/methods , Systems Biology/methods , Gene Expression Profiling , Genome/genetics , Models, Genetic , RNA, Messenger/analysis , RNA, Messenger/genetics , RNA, Messenger/metabolismABSTRACT
Haematobia irritans irritans (Linnaeus, 1758: Diptera: Muscidae), the horn fly, is an external parasite of penned and pastured livestock that causes a major economic impact on cattle production worldwide. Pesticides such as synthetic pyrethroids and organophosphates are routinely used to control horn flies; however, resistance to these chemicals has become a concern in several countries. To further elucidate the molecular mechanisms of resistance in horn fly populations, we sequenced the transcriptomes of ten populations of horn flies from the southern US possessing varying degrees of pesticide resistance levels to pyrethroids, organophosphates, and endosulfans. We employed an Illumina paired end HiSeq approach, followed by de novo assembly of the transcriptomes using CLC Genomics Workbench 8.0.1 De Novo Assembler using multiple kmers, and annotation using Blast2GO PRO version 5.2.5. The Gene Ontology biological process term Response to Insecticide was found in all the populations, but at an increased frequency in the populations with higher levels of insecticide resistance. The raw sequence reads are archived in the Sequence Read Archive (SRA) and assembled population transcriptomes in the Transcriptome Shotgun Assembly (TSA) at the National Center for Biotechnology Information (NCBI).
ABSTRACT
Knowledge of host associations of blood-feeding vectors may afford insights into managing disease systems and protecting public health. However, the ability of methods to distinguish bloodmeal sources varies widely. We used two methods-Sanger sequencing and amplicon deep sequencing-to target a 228 bp region of the vertebrate Cytochrome b gene and determine hosts fed upon by triatomines (n = 115) collected primarily in Texas, USA. Direct Sanger sequencing of PCR amplicons was successful for 36 samples (31%). Sanger sequencing revealed 15 distinct host species, which included humans, domestic animals (Canis lupus familiaris, Ovis aries, Gallus gallus, Bos taurus, Felis catus, and Capra hircus), wildlife (Rattus rattus, Incilius nebulifer, Sciurus carolinensis, Sciurus niger, and Odocoileus virginianus), and captive animals (Panthera tigris, Colobus spp., and Chelonoidis carbonaria). Samples sequenced by the Sanger method were also subjected to Illumina MiSeq amplicon deep sequencing. The amplicon deep sequencing results (average of 302,080 usable reads per sample) replicated the host community revealed using Sanger sequencing, and detected additional hosts in five triatomines (13.9%), including two additional blood sources (Procyon lotor and Bassariscus astutus). Up to four bloodmeal sources were detected in a single triatomine (I. nebulifer, Homo sapiens, C. lupus familiaris, and S. carolinensis). Enhanced understanding of vector-host-parasite networks may allow for integrated vector management programs focusing on highly-utilized and highly-infected host species.
Subject(s)
Chagas Disease , Deer , Trypanosoma cruzi , Animals , Animals, Domestic/genetics , Cats , Cattle , Chagas Disease/parasitology , Deer/genetics , Dogs , High-Throughput Nucleotide Sequencing , Trypanosoma cruzi/geneticsABSTRACT
The phagocytosis and destruction of pathogens in lysosomes constitute central elements of innate immune defense. Here, we show that Brucella, the causative agent of brucellosis, the most prevalent bacterial zoonosis globally, subverts this immune defense pathway by activating regulated IRE1α-dependent decay (RIDD) of Bloc1s1 mRNA encoding BLOS1, a protein that promotes endosome-lysosome fusion. RIDD-deficient cells and mice harboring a RIDD-incompetent variant of IRE1α were resistant to infection. Inactivation of the Bloc1s1 gene impaired the ability to assemble BLOC-1-related complex (BORC), resulting in differential recruitment of BORC-related lysosome trafficking components, perinuclear trafficking of Brucella-containing vacuoles (BCVs), and enhanced susceptibility to infection. The RIDD-resistant Bloc1s1 variant maintains the integrity of BORC and a higher-level association of BORC-related components that promote centrifugal lysosome trafficking, resulting in enhanced BCV peripheral trafficking and lysosomal destruction, and resistance to infection. These findings demonstrate that host RIDD activity on BLOS1 regulates Brucella intracellular parasitism by disrupting BORC-directed lysosomal trafficking. Notably, coronavirus murine hepatitis virus also subverted the RIDD-BLOS1 axis to promote intracellular replication. Our work establishes BLOS1 as a novel immune defense factor whose activity is hijacked by diverse pathogens.
Subject(s)
Brucella , Brucellosis , Animals , Brucellosis/metabolism , Brucellosis/microbiology , Endoribonucleases/metabolism , Endosomes/metabolism , Mice , Protein Serine-Threonine KinasesABSTRACT
BACKGROUND: RNA-Seq is the recently developed high-throughput sequencing technology for profiling the entire transcriptome in any organism. It has several major advantages over current hybridization-based approach such as microarrays. However, the cost per sample by RNA-Seq is still prohibitive for most laboratories. With continued improvement in sequence output, it would be cost-effective if multiple samples are multiplexed and sequenced in a single lane with sufficient transcriptome coverage. The objective of this analysis is to evaluate what sequencing depth might be sufficient to interrogate gene expression profiling in the chicken by RNA-Seq. RESULTS: Two cDNA libraries from chicken lungs were sequenced initially, and 4.9 million (M) and 1.6 M (60 bp) reads were generated, respectively. With significant improvements in sequencing technology, two technical replicate cDNA libraries were re-sequenced. Totals of 29.6 M and 28.7 M (75 bp) reads were obtained with the two samples. More than 90% of annotated genes were detected in the data sets with 28.7-29.6 M reads, while only 68% of genes were detected in the data set with 1.6 M reads. The correlation coefficients of gene expression between technical replicates within the same sample were 0.9458 and 0.8442. To evaluate the appropriate depth needed for mRNA profiling, a random sampling method was used to generate different number of reads from each sample. There was a significant increase in correlation coefficients from a sequencing depth of 1.6 M to 10 M for all genes except highly abundant genes. No significant improvement was observed from the depth of 10 M to 20 M (75 bp) reads. CONCLUSION: The analysis from the current study demonstrated that 30 M (75 bp) reads is sufficient to detect all annotated genes in chicken lungs. Ten million (75 bp) reads could detect about 80% of annotated chicken genes, and RNA-Seq at this depth can serve as a replacement of microarray technology. Furthermore, the depth of sequencing had a significant impact on measuring gene expression of low abundant genes. Finally, the combination of experimental and simulation approaches is a powerful approach to address the relationship between the depth of sequencing and transcriptome coverage.
Subject(s)
Chickens/genetics , Gene Expression Profiling , High-Throughput Nucleotide Sequencing , Lung/metabolism , Sequence Analysis, RNA , Animals , Gene Library , Molecular Sequence Annotation , RNA, Messenger/geneticsABSTRACT
Global climate change has increased the number of severe flooding events that affect agriculture, including rice production in the U.S. and internationally. Heavy rainfall can cause rice plants to be completely submerged, which can significantly affect grain yield or completely destroy the plants. Recently, a major effect submergence tolerance QTL during the vegetative stage, qSub8.1, which originated from Ciherang-Sub1, was identified in a mapping population derived from a cross between Ciherang-Sub1 and IR10F365. Ciherang-Sub1 was, in turn, derived from a cross between Ciherang and IR64-Sub1. Here, we characterize the qSub8.1 region by analyzing the sequence information of Ciherang-Sub1 and its two parents (Ciherang and IR64-Sub1) and compare the whole genome profile of these varieties with the Nipponbare and Minghui 63 (MH63) reference genomes. The three rice varieties were sequenced with 150 bp pair-end whole-genome shotgun sequencing (Illumina HiSeq4000), followed by performing the Trimmomatic-SOAPdenovo2-MUMmer3 pipeline for genome assembly, resulting in approximate genome sizes of 354.4, 343.7, and 344.7 Mb, with N50 values of 25.1, 25.4, and 26.1 kb, respectively. The results showed that the Ciherang-Sub1 genome is composed of 59-63% Ciherang, 22-24% of IR64-Sub1, and 15-17% of unknown sources. The genome profile revealed a more detailed genomic composition than previous marker-assisted breeding and showed that the qSub8.1 region is mostly from Ciherang, with some introgressed segments from IR64-Sub1 and currently unknown source(s).
ABSTRACT
Ticks from the genus Rhipicephalus have enormous global economic impact as ectoparasites of cattle. Rhipicephalus microplus and Rhipicephalus annulatus are known to harbor infectious pathogens such as Babesia bovis, Babesia bigemina, and Anaplasma marginale. Having reference quality genomes of these ticks would advance research to identify druggable targets for chemical entities with acaricidal activity and refine anti-tick vaccine approaches. We sequenced and assembled the genomes of R. microplus and R. annulatus, using Pacific Biosciences and HiSeq 4000 technologies on very high molecular weight genomic DNA. We used 22 and 29 SMRT cells on the Pacific Biosciences Sequel for R. microplus and R. annulatus, respectively, and 3 lanes of the Illumina HiSeq 4000 platform for each tick. The PacBio sequence yields for R. microplus and R. annulatus were 21.0 and 27.9 million subreads, respectively, which were assembled with Canu v. 1.7. The final Canu assemblies consisted of 92,167 and 57,796 contigs with an average contig length of 39,249 and 69,055 bp for R. microplus and R. annulatus, respectively. Annotated genome quality was assessed by BUSCO analysis to provide quantitative measures for each assembled genome. Over 82% and 92% of the 1066 member BUSCO gene set was found in the assembled genomes of R. microplus and R. annulatus, respectively. For R. microplus, only 189 of the 1066 BUSCO genes were missing and only 140 were present in a fragmented condition. For R. annulatus, only 75 of the BUSCO genes were missing and only 109 were present in a fragmented condition. The raw sequencing reads and the assembled contigs/scaffolds are archived at the National Center for Biotechnology Information.