Exome sequencing in multiplex families with left-sided cardiac defects has high yield for disease gene discovery.

Congenital heart disease (CHD) is a common group of birth defects with a strong genetic contribution to their etiology, but historically the diagnostic yield from exome studies of isolated CHD has been low. Pleiotropy, variable expressivity, and the difficulty of accurately phenotyping newborns contribute to this problem. We hypothesized that performing exome sequencing on selected individuals in families with multiple members affected by left-sided CHD, then filtering variants by population frequency, in silico predictive algorithms, and phenotypic annotations from publicly available databases would increase this yield and generate a list of candidate disease-causing variants that would show a high validation rate. In eight of the nineteen families in our study (42%), we established a well-known gene/phenotype link for a candidate variant or performed confirmation of a candidate variant's effect on protein function, including variants in genes not previously described or firmly established as disease genes in the body of CHD literature: BMP10, CASZ1, ROCK1 and SMYD1. Two plausible variants in different genes were found to segregate in the same family in two instances suggesting oligogenic inheritance. These results highlight the need for functional validation and demonstrate that in the era of next-generation sequencing, multiplex families with isolated CHD can still bring high yield to the discovery of novel disease genes.

Unveiling the immune dynamics of Neisseria persistent oral colonization.

Commensal bacteria are crucial in maintaining host physiological homeostasis, immune system development, and protection against pathogens. Despite their significance, the factors influencing persistent bacterial colonization and their impact on the host still need to be fully understood. Animal models have served as valuable tools to investigate these interactions, but most have limitations. The bacterial genus Neisseria, which includes both commensal and pathogenic species, has been studied from a pathogenicity to humans perspective but lacks models that study immune responses in the context of long-term persistence. Neisseria musculi, a recently described natural commensal of mice, offers a unique opportunity to study long-term host-commensal interactions. In this study, for the first time, we have used this model to study the transcriptional, phenotypic, and functional dynamics of immune cell signatures in the mucosal and systemic tissue of mice in response to N. musculi colonization. We found key genes and pathways vital for immune homeostasis in palate tissue, validated by flow cytometry of immune cells from the lung, blood, and spleen. This study offers a novel avenue for advancing our understanding of host-bacteria dynamics and may provide a platform for developing efficacious interventions against mucosal persistence by pathogenic Neisseria.

Full-length isoform concatenation sequencing to resolve cancer transcriptome complexity.

BACKGROUND: Cancers exhibit complex transcriptomes with aberrant splicing that induces isoform-level differential expression compared to non-diseased tissues. Transcriptomic profiling using short-read sequencing has utility in providing a cost-effective approach for evaluating isoform expression, although short-read assembly displays limitations in the accurate inference of full-length transcripts. Long-read RNA sequencing (Iso-Seq), using the Pacific Biosciences (PacBio) platform, can overcome such limitations by providing full-length isoform sequence resolution which requires no read assembly and represents native expressed transcripts. A constraint of the Iso-Seq protocol is due to fewer reads output per instrument run, which, as an example, can consequently affect the detection of lowly expressed transcripts. To address these deficiencies, we developed a concatenation workflow, PacBio Full-Length Isoform Concatemer Sequencing (PB_FLIC-Seq), designed to increase the number of unique, sequenced PacBio long-reads thereby improving overall detection of unique isoforms. In addition, we anticipate that the increase in read depth will help improve the detection of moderate to low-level expressed isoforms. RESULTS: In sequencing a commercial reference (Spike-In RNA Variants; SIRV) with known isoform complexity we demonstrated a 3.4-fold increase in read output per run and improved SIRV recall when using the PB_FLIC-Seq method compared to the same samples processed with the Iso-Seq protocol. We applied this protocol to a translational cancer case, also demonstrating the utility of the PB_FLIC-Seq method for identifying differential full-length isoform expression in a pediatric diffuse midline glioma compared to its adjacent non-malignant tissue. Our data analysis revealed increased expression of extracellular matrix (ECM) genes within the tumor sample, including an isoform of the Secreted Protein Acidic and Cysteine Rich (SPARC) gene that was expressed 11,676-fold higher than in the adjacent non-malignant tissue. Finally, by using the PB_FLIC-Seq method, we detected several cancer-specific novel isoforms. CONCLUSION: This work describes a concatenation-based methodology for increasing the number of sequenced full-length isoform reads on the PacBio platform, yielding improved discovery of expressed isoforms. We applied this workflow to profile the transcriptome of a pediatric diffuse midline glioma and adjacent non-malignant tissue. Our findings of cancer-specific novel isoform expression further highlight the importance of long-read sequencing for characterization of complex tumor transcriptomes.

Picarbutrazox effectiveness added to a seed treatment mixture for management of Oomycetes that impact soybean in Ohio.

None of the current oomycota fungicides are effective towards all species of Phytophthora, Phytopythium, Globisporangium, and Pythium that affect soybean seed and seedlings in Ohio. Picarbutrazox is a new oomyceticide with a novel mode of action towards Oomycete pathogens. Our objectives were to evaluate picarbutrazox to determine i) baseline sensitivity (EC50) to 189 isolates of 29 species, ii) the efficacy with a base seed treatment with three cultivars with different levels of resistance in 14 field environments; and iii) if the rhizosphere microbiome was affected by the addition of the seed treatment on a moderately susceptible cultivar. The mycelial growth of all isolates was inhibited beginning at 0.001µg and the EC50 ranged from 0.0013 to 0.0483 µg a.i. ml-1. The effect of seed treatment was significantly different for plant population and yield in 8 of 14 and 6 of 12 environments, respectively. The addition of picarbutrazox at 1 and 2.5 g a.i. 100 kg seed-1 to the base seed treatment compared to the base alone was associated with higher plant populations and yield in 3 and 1 environment, respectively. There was limited impact of the seed treatment mefenoxam 7.5 g a.i. plus picarbutrazox 1 g a.i. per 100 kg seed-1 on the oomycetes detected in the rhizosphere of soybean seedlings collected at the V1 growth stage. Picarbutrazox has efficacy towards a wider range of oomycetes that cause disease on soybean and this will be another oomyceticide tool to combat early season damping-off in areas where environmental conditions highly favor disease development.

PTEN somatic mutations contribute to spectrum of cerebral overgrowth.

Phosphatase and tensin homologue (PTEN) regulates cell growth and survival through inhibition of the mammalian target of rapamycin (MTOR) signalling pathway. Germline genetic variation of PTEN is associated with autism, macrocephaly and PTEN hamartoma tumour syndromes. The effect of developmental PTEN somatic mutations on nervous system phenotypes is not well understood, although brain somatic mosaicism of MTOR pathway genes is an emerging cause of cortical dysplasia and epilepsy in the paediatric population. Here we report two somatic variants of PTEN affecting a single patient presenting with intractable epilepsy and hemimegalencephaly that varied in clinical severity throughout the left cerebral hemisphere. High-throughput sequencing analysis of affected brain tissue identified two somatic variants in PTEN. The first variant was present in multiple cell lineages throughout the entire hemisphere and associated with mild cerebral overgrowth. The second variant was restricted to posterior brain regions and affected the opposite PTEN allele, resulting in a segmental region of more severe malformation, and the only neurons in which it was found by single-nuclei RNA-sequencing had a unique disease-related expression profile. This study reveals brain mosaicism of PTEN as a disease mechanism of hemimegalencephaly and furthermore demonstrates the varying effects of single- or bi-allelic disruption of PTEN on cortical phenotypes.

The macronuclear genome of anaerobic ciliate Entodinium caudatum reveals its biological features adapted to the distinct rumen environment.

Entodinium caudatum is an anaerobic binucleated ciliate representing the most dominant protozoal species in the rumen. However, its biological features are largely unknown due to the inability to establish an axenic culture. In this study, we primally sequenced its macronucleus (MAC) genome to aid the understanding of its metabolism, physiology, ecology. We isolated the MAC of E. caudatum strain MZG-1 and sequenced the MAC genome using Illumina MiSeq, MinION, and PacBio RSII systems. De novo assembly of the MiSeq sequence reads followed with subsequent scaffolding with MinION and PacBio reads resulted in a draft MAC genome about 117 Mbp. A large number of carbohydrate-active enzymes were likely acquired through horizontal gene transfer. About 8.74% of the E. caudatum predicted proteome was predicted as proteases. The MAC genome of E. caudatum will help better understand its important roles in rumen carbohydrate metabolism, and interaction with other members of the rumen microbiome.

Discovery of clinically relevant fusions in pediatric cancer.

BACKGROUND: Pediatric cancers typically have a distinct genomic landscape when compared to adult cancers and frequently carry somatic gene fusion events that alter gene expression and drive tumorigenesis. Sensitive and specific detection of gene fusions through the analysis of next-generation-based RNA sequencing (RNA-Seq) data is computationally challenging and may be confounded by low tumor cellularity or underlying genomic complexity. Furthermore, numerous computational tools are available to identify fusions from supporting RNA-Seq reads, yet each algorithm demonstrates unique variability in sensitivity and precision, and no clearly superior approach currently exists. To overcome these challenges, we have developed an ensemble fusion calling approach to increase the accuracy of identifying fusions. RESULTS: Our Ensemble Fusion (EnFusion) approach utilizes seven fusion calling algorithms: Arriba, CICERO, FusionMap, FusionCatcher, JAFFA, MapSplice, and STAR-Fusion, which are packaged as a fully automated pipeline using Docker and Amazon Web Services (AWS) serverless technology. This method uses paired end RNA-Seq sequence reads as input, and the output from each algorithm is examined to identify fusions detected by a consensus of at least three algorithms. These consensus fusion results are filtered by comparison to an internal database to remove likely artifactual fusions occurring at high frequencies in our internal cohort, while a "known fusion list" prevents failure to report known pathogenic events. We have employed the EnFusion pipeline on RNA-Seq data from 229 patients with pediatric cancer or blood disorders studied under an IRB-approved protocol. The samples consist of 138 central nervous system tumors, 73 solid tumors, and 18 hematologic malignancies or disorders. The combination of an ensemble fusion-calling pipeline and a knowledge-based filtering strategy identified 67 clinically relevant fusions among our cohort (diagnostic yield of 29.3%), including RBPMS-MET, BCAN-NTRK1, and TRIM22-BRAF fusions. Following clinical confirmation and reporting in the patient's medical record, both known and novel fusions provided medically meaningful information. CONCLUSIONS: The EnFusion pipeline offers a streamlined approach to discover fusions in cancer, at higher levels of sensitivity and accuracy than single algorithm methods. Furthermore, this method accurately identifies driver fusions in pediatric cancer, providing clinical impact by contributing evidence to diagnosis and, when appropriate, indicating targeted therapies.

Hypomorphic alleles pose challenges in rare disease genomic variant interpretation.

Exon skipping associated with an ATP7B intronic variant in a patient with Wilson's disease. (A) Sashimi plot visualization of aligned RNA sequencing data from proband liver tissue at ATP7B exons 14-13-12. The red track shows traditional RNA-seq data; the blue track shows RNA-seq enriched with exon capture (cDNA-cap) which achieves higher depth of protein-coding transcripts. The histogram indicates overall sequencing depth while arcs tabulate the number of junction-spanning reads supporting exon pairs. (B) The domain structure (top) and exon structure (bottom) of ATP7B. Loss of exon 13 (dashed box) would remove a transmembrane domain and disrupt the first phosphorylation domain.

Farm Stage, Bird Age, and Body Site Dominantly Affect the Quantity, Taxonomic Composition, and Dynamics of Respiratory and Gut Microbiota of Commercial Layer Chickens.

The digestive and respiratory tracts of chickens are colonized by bacteria that are believed to play important roles in the overall health and performance of the birds. Most of the current research on the commensal bacteria (microbiota) of chickens has focused on broilers and gut microbiota, and less attention has been given to layers and respiratory microbiota. This research bias has left significant gaps in our knowledge of the layer microbiome. This study was conducted to define the core microbiota colonizing the upper respiratory tract (URT) and lower intestinal tract (LIT) in commercial layers under field conditions. One hundred eighty-one chickens were sampled from a flock of >80,000 birds at nine times to collect samples for 16S rRNA gene-based bacterial metabarcoding. Generally, the body site and age/farm stage had very dominant effects on the quantity, taxonomic composition, and dynamics of core bacteria. Remarkably, ileal and URT microbiota were compositionally more related to each other than to that from the cecum. Unique taxa dominated in each body site yet some taxa overlapped between URT and LIT sites, demonstrating a common core. The overlapping bacteria also contained various levels of several genera with well-recognized avian pathogens. Our findings suggest that significant interaction exists between gut and respiratory microbiota, including potential pathogens, in all stages of the farm sequence. The baseline data generated in this study can be useful for the development of effective microbiome-based interventions to enhance production performance and to prevent and control disease in commercial chicken layers.IMPORTANCE The poultry industry is faced with numerous challenges associated with infectious diseases and suboptimal performance of flocks. As microbiome research continues to grow, it is becoming clear that poultry health and production performance are partly influenced by nonpathogenic symbionts that occupy different habitats within the bird. This study has defined the baseline composition and overlaps between respiratory and gut bacteria in healthy, optimally performing chicken layers across all stages of the commercial farm sequence. Consequently, the study has set the groundwork for the development of interventions that seek to enhance production performance and to prevent and control infectious diseases through the modulation of gut and respiratory bacteria.

Differentiated transcriptional signatures in the maize landraces of Chiapas, Mexico.

BACKGROUND: Landrace farmers are the keepers of crops locally adapted to the environments where they are cultivated. Patterns of diversity across the genome can provide signals of past evolution in the face of abiotic and biotic change. Understanding this rich genetic resource is imperative especially since diversity can provide agricultural security as climate continues to shift. RESULTS: Here we employ RNA sequencing (RNA-seq) to understand the role that conditions that vary across a landscape may have played in shaping genetic diversity in the maize landraces of Chiapas, Mexico. We collected landraces from three distinct elevational zones and planted them in a midland common garden. Early season leaf tissue was collected for RNA-seq and we performed weighted gene co-expression network analysis (WGCNA). We then used association analysis between landrace co-expression module expression values and environmental parameters of landrace origin to elucidate genes and gene networks potentially shaped by environmental factors along our study gradient. Elevation of landrace origin affected the transcriptome profiles. Two co-expression modules were highly correlated with temperature parameters of landrace origin and queries into their 'hub' genes suggested that temperature may have led to differentiation among landraces in hormone biosynthesis/signaling and abiotic and biotic stress responses. We identified several 'hub' transcription factors and kinases as candidates for the regulation of these responses. CONCLUSIONS: These findings indicate that natural selection may influence the transcriptomes of crop landraces along an elevational gradient in a major diversity center, and provide a foundation for exploring the genetic basis of local adaptation. While we cannot rule out the role of neutral evolutionary forces in the patterns we have identified, combining whole transcriptome sequencing technologies, established bioinformatics techniques, and common garden experimentation can powerfully elucidate structure of adaptive diversity across a varied landscape. Ultimately, gaining such understanding can facilitate the conservation and strategic utilization of crop genetic diversity in a time of climate change.

Johnsongrass mosaic virus Contributes to Maize Lethal Necrosis in East Africa.

Maize lethal necrosis (MLN), a severe virus disease of maize, has emerged in East Africa in recent years with devastating effects on production and food security where maize is a staple subsistence crop. In extensive surveys of MLN-symptomatic plants in East Africa, sequences of Johnsongrass mosaic virus (JGMV) were identified in Uganda, Kenya, Rwanda, and Tanzania. The East African JGMV is distinct from previously reported isolates and infects maize, sorghum, and Johnsongrass but not wheat or oat. This isolate causes MLN in coinfection with Maize chlorotic mottle virus (MCMV), as reported for other potyviruses, and was present in MLN-symptomatic plants in which the major East African potyvirus, Sugarcane mosaic virus (SCMV), was not detected. Virus titers were compared in single and coinfections by quantitative reverse transcription-polymerase chain reaction. MCMV titer increased in coinfected plants whereas SCMV, Maize dwarf mosaic virus, and JGMV titers were unchanged compared with single infections at 11 days postinoculation. Together, these results demonstrate the presence of an East African JGMV that contributes to MLN in the region.

Association Analysis for Bacterial Spot Resistance in a Directionally Selected Complex Breeding Population of Tomato.

Bacterial spot of tomato is caused by at least four species of Xanthomonas with multiple physiological races. We developed a complex breeding population for simultaneous discovery of marker-trait linkage, validation of existing quantitative trait loci (QTL), and pyramiding of resistance. Six advanced accessions with resistance from distinct sources were crossed in all combinations and their F1 hybrids were intercrossed. Over 1,100 segregating progeny were evaluated in the field following inoculation with X. euvesicatoria race T1 strains. We selected 5% of the most resistant and 5% of the most susceptible progeny for evaluation as plots in two subsequent replicated field trials inoculated with T1 and T3 (X. perforans) strains. The estimated heritability of T1 resistance was 0.32. In order to detect previously reported resistance genes, as well as novel QTL, we explored methods to correct for population structure and analysis based on single markers or haplotypes. Both single-point and haplotype analyses identified strong associations in the genomic regions known to carry Rx-3 (chromosome 5) and Rx-4/Xv3 (chromosome 11). Accounting for kinship and structure generally improved the fit of statistical models. Detection of known loci was improved by adding kinship or a combination of kinship and structure using a Q matrix from model-based clustering. Additional QTL were detected on chromosomes 1, 4, 6, and 7 for T1 resistance and chromosomes 2, 4, and 6 for T3 resistance (P < 0.01). Haplotype analysis improved our ability to trace the origin of positive alleles. These results demonstrate that both known and novel associations can be identified using complex breeding populations that have experienced directional selection.

Virus-independent and common transcriptome responses of leafhopper vectors feeding on maize infected with semi-persistently and persistent propagatively transmitted viruses.

BACKGROUND: Insects are the most important epidemiological factors for plant virus disease spread, with >75% of viruses being dependent on insects for transmission to new hosts. The black-faced leafhopper (Graminella nigrifrons Forbes) transmits two viruses that use different strategies for transmission: Maize chlorotic dwarf virus (MCDV) which is semi-persistently transmitted and Maize fine streak virus (MFSV) which is persistently and propagatively transmitted. To date, little is known regarding the molecular and cellular mechanisms in insects that regulate the process and efficiency of transmission, or how these mechanisms differ based on virus transmission strategy. RESULTS: RNA-Seq was used to examine transcript changes in leafhoppers after feeding on MCDV-infected, MFSV-infected and healthy maize for 4 h and 7 d. After sequencing cDNA libraries constructed from whole individuals using Illumina next generation sequencing, the Rnnotator pipeline in Galaxy was used to reassemble the G. nigrifrons transcriptome. Using differential expression analyses, we identified significant changes in transcript abundance in G. nigrifrons. In particular, transcripts implicated in the innate immune response and energy production were more highly expressed in insects fed on virus-infected maize. Leafhoppers fed on MFSV-infected maize also showed an induction of transcripts involved in hemocoel and cell-membrane linked immune responses within four hours of feeding. Patterns of transcript expression were validated for a subset of transcripts by quantitative real-time reverse transcription polymerase chain reaction using RNA samples collected from insects fed on healthy or virus-infected maize for between a 4 h and seven week period. CONCLUSIONS: We expected, and found, changes in transcript expression in G. nigrifrons feeding of maize infected with a virus (MFSV) that also infects the leafhopper, including induction of immune responses in the hemocoel and at the cell membrane. The significant induction of the innate immune system in G. nigrifrons fed on a foregut-borne virus (MCDV) that does not infect leafhoppers was less expected. The changes in transcript accumulation that occur independent of the mode of pathogen transmission could be key for identifying insect factors that disrupt vector-mediated plant virus transmission.

RNA-sequencing reveals early, dynamic transcriptome changes in the corollas of pollinated petunias.

BACKGROUND: Pollination reduces flower longevity in many angiosperms by accelerating corolla senescence. This response requires hormone signaling between the floral organs and results in the degradation of macromolecules and organelles within the petals to allow for nutrient remobilization to developing seeds. To investigate early pollination-induced changes in petal gene expression, we utilized high-throughput sequencing to identify transcripts that were differentially expressed between corollas of pollinated Petunia × hybrida flowers and their unpollinated controls at 12, 18, and 24 hours after opening. RESULTS: In total, close to 0.5 billion Illumina 101 bp reads were generated, de novo assembled, and annotated, resulting in an EST library of approximately 33 K genes. Over 4,700 unique, differentially expressed genes were identified using comparisons between the pollinated and unpollinated libraries followed by pairwise comparisons of pollinated libraries to unpollinated libraries from the same time point (i.e. 12-P/U, 18-P/U, and 24-P/U) in the Bioconductor R package DESeq2. Over 500 gene ontology terms were enriched. The response to auxin stimulus and response to 1-aminocyclopropane-1-carboxylic acid terms were enriched by 12 hours after pollination (hap). Using weighted gene correlation network analysis (WGCNA), three pollination-specific modules were identified. Module I had increased expression across pollinated corollas at 12, 18, and 24 h, and modules II and III had a peak of expression in pollinated corollas at 18 h. A total of 15 enriched KEGG pathways were identified. Many of the genes from these pathways were involved in metabolic processes or signaling. More than 300 differentially expressed transcription factors were identified. CONCLUSIONS: Gene expression changes in corollas were detected within 12 hap, well before fertilization and corolla wilting or ethylene evolution. Significant changes in gene expression occurred at 18 hap, including the up-regulation of autophagy and down-regulation of ribosomal genes and genes involved in carbon fixation. This transcriptomic database will greatly expand the genetic resources available in petunia. Additionally, it will guide future research aimed at identifying the best targets for increasing flower longevity by delaying corolla senescence.

Increase in Full-Length Dystrophin by Exon Skipping in Duchenne Muscular Dystrophy Patients with Single Exon Duplications: An Open-label Study.

Single exon duplications account for disease in a minority of Duchenne muscular dystrophy patients. Exon skipping in these patients has the potential to be highly therapeutic through restoration of full-length dystrophin expression. We conducted a 48-week open label study of casimersen and golodirsen in 3 subjects with an exon 45 or 53 duplication. Two subjects (aged 18 and 23 years) were non-ambulatory at baseline. Upper limb, pulmonary, and cardiac function appeared stable in the 2 subjects in whom they could be evaluated. Dystrophin expression increased from 0.94 % ±0.59% (mean±SD) of normal to 5.1% ±2.9% by western blot. Percent dystrophin positive fibers also rose from 14% ±17% at baseline to 50% ±42% . Our results provide initial evidence that the use of exon-skipping drugs may increase dystrophin levels in patients with single-exon duplications.

Glucocorticoid- and pioglitazone-induced proteinuria reduction in experimental NS both correlate with glomerular ECM modulation.

Idiopathic nephrotic syndrome (NS) is a common glomerular disease. Although glucocorticoids (GC) are the primary treatment, the PPARγ agonist pioglitazone (Pio) also reduces proteinuria in patients with NS and directly protects podocytes from injury. Because both drugs reduce proteinuria, we hypothesized these effects result from overlapping transcriptional patterns. Systems biology approaches compared glomerular transcriptomes from rats with PAN-induced NS treated with GC vs. Pio and identified 29 commonly regulated genes-of-interest, primarily involved in extracellular matrix (ECM) remodeling. Correlation with clinical idiopathic NS patient datasets confirmed glomerular ECM dysregulation as a potential mechanism of injury. Cellular deconvolution in silico revealed GC- and Pio-induced amelioration of altered genes primarily within podocytes and mesangial cells. While validation studies are indicated, these analyses identified molecular pathways involved in the early stages of NS (prior to scarring), suggesting that targeting glomerular ECM dysregulation may enable a future non-immunosuppressive approach for proteinuria reduction in idiopathic NS.

Molecular mechanisms underpinning quantitative resistance to Phytophthora sojae in Glycine max using a systems genomics approach.

Expression of quantitative disease resistance in many host-pathogen systems is controlled by genes at multiple loci, each contributing a small effect to the overall response. We used a systems genomics approach to study the molecular underpinnings of quantitative disease resistance in the soybean-Phytophthora sojae pathosystem, incorporating expression quantitative trait loci (eQTL) mapping and gene co-expression network analysis to identify the genes putatively regulating transcriptional changes in response to inoculation. These findings were compared to previously mapped phenotypic (phQTL) to identify the molecular mechanisms contributing to the expression of this resistance. A subset of 93 recombinant inbred lines (RILs) from a Conrad × Sloan population were inoculated with P. sojae isolate 1.S.1.1 using the tray-test method; RNA was extracted, sequenced, and the normalized read counts were genetically mapped from tissue collected at the inoculation site 24 h after inoculation from both mock and inoculated samples. In total, more than 100,000 eQTLs were mapped. There was a switch from predominantly cis-eQTLs in the mock treatment to an almost entirely nonoverlapping set of predominantly trans-eQTLs in the inoculated treatment, where greater than 100-fold more eQTLs were mapped relative to mock, indicating vast transcriptional reprogramming due to P. sojae infection occurred. The eQTLs were organized into 36 hotspots, with the four largest hotspots from the inoculated treatment corresponding to more than 70% of the eQTLs, each enriched for genes within plant-pathogen interaction pathways. Genetic regulation of trans-eQTLs in response to the pathogen was predicted to occur through transcription factors and signaling molecules involved in plant-pathogen interactions, plant hormone signal transduction, and MAPK pathways. Network analysis identified three co-expression modules that were correlated with susceptibility to P. sojae and associated with three eQTL hotspots. Among the eQTLs co-localized with phQTLs, two cis-eQTLs with putative functions in the regulation of root architecture or jasmonic acid, as well as the putative master regulators of an eQTL hotspot nearby a phQTL, represent candidates potentially underpinning the molecular control of these phQTLs for resistance.

Synergistic Role of NK Cells and Monocytes in Promoting Atherogenesis in Severe COVID-19 Patients.

Clinical data demonstrate an increased predisposition to cardiovascular disease (CVD) following severe COVID-19 infection. This may be driven by a dysregulated immune response associated with severe disease. Monocytes and vascular tissue resident macrophages play a critical role in atherosclerosis, the main pathology leading to ischemic CVD. Natural killer (NK) cells are a heterogenous group of cells that are critical during viral pathogenesis and are known to be dysregulated during severe COVID-19 infection. Their role in atherosclerotic cardiovascular disease has recently been described. However, the contribution of their altered phenotypes to atherogenesis following severe COVID-19 infection is unknown. We demonstrate for the first time that during and after severe COVID-19, circulating proinflammatory monocytes and activated NK cells act synergistically to increase uptake of oxidized low-density lipoprotein (Ox-LDL) into vascular tissue with subsequent foam cell generation leading to atherogenesis despite recovery from acute infection. Our data provide new insights, revealing the roles of monocytes/macrophages, and NK cells in COVID-19-related atherogenesis.

Dissection of two soybean QTL conferring partial resistance to Phytophthora sojae through sequence and gene expression analysis.

BACKGROUND: Phytophthora sojae is the primary pathogen of soybeans that are grown on poorly drained soils. Race-specific resistance to P. sojae in soybean is gene-for-gene, although in many areas of the US and worldwide there are populations that have adapted to the most commonly deployed resistance to P. sojae ( Rps) genes. Hence, this system has received increased attention towards identifying mechanisms and molecular markers associated with partial resistance to this pathogen. Several quantitative trait loci (QTL) have been identified in the soybean cultivar 'Conrad' that contributes to the expression of partial resistance to multiple P. sojae isolates. RESULTS: In this study, two of the Conrad QTL on chromosome 19 were dissected through sequence and expression analysis of genes in both resistant (Conrad) and susceptible ('Sloan') genotypes. There were 1025 single nucleotide polymorphisms (SNPs) in 87 of 153 genes sequenced from Conrad and Sloan. There were 304 SNPs in 54 genes sequenced from Conrad compared to those from both Sloan and Williams 82, of which 11 genes had SNPs unique to Conrad. Eleven of 19 genes in these regions analyzed with qRT-PCR had significant differences in fold change of transcript abundance in response to infection with P. sojae in lines with QTL haplotype from the resistant parent compared to those with the susceptible parent haplotype. From these, 8 of the 11 genes had SNPs in the upstream, untranslated region, exon, intron, and/or downstream region. These 11 candidate genes encode proteins potentially involved in signal transduction, hormone-mediated pathways, plant cell structural modification, ubiquitination, and basal resistance. CONCLUSIONS: These findings may indicate a complex defense network with multiple mechanisms underlying these two soybean QTL conferring resistance to P. sojae. SNP markers derived from these candidate genes can contribute to fine mapping of QTL and marker assisted breeding for resistance to P. sojae.

RNA-Seq and molecular docking reveal multi-level pesticide resistance in the bed bug.

BACKGROUND: Bed bugs (Cimex lectularius) are hematophagous nocturnal parasites of humans that have attained high impact status due to their worldwide resurgence. The sudden and rampant resurgence of C. lectularius has been attributed to numerous factors including frequent international travel, narrower pest management practices, and insecticide resistance. RESULTS: We performed a next-generation RNA sequencing (RNA-Seq) experiment to find differentially expressed genes between pesticide-resistant (PR) and pesticide-susceptible (PS) strains of C. lectularius. A reference transcriptome database of 51,492 expressed sequence tags (ESTs) was created by combining the databases derived from de novo assembled mRNA-Seq tags (30,404 ESTs) and our previous 454 pyrosequenced database (21,088 ESTs). The two-way GLMseq analysis revealed ~15,000 highly significant differentially expressed ESTs between the PR and PS strains. Among the top 5,000 differentially expressed ESTs, 109 putative defense genes (cuticular proteins, cytochrome P450s, antioxidant genes, ABC transporters, glutathione S-transferases, carboxylesterases and acetyl cholinesterase) involved in penetration resistance and metabolic resistance were identified. Tissue and development-specific expression of P450 CYP3 clan members showed high mRNA levels in the cuticle, Malpighian tubules, and midgut; and in early instar nymphs, respectively. Lastly, molecular modeling and docking of a candidate cytochrome P450 (CYP397A1V2) revealed the flexibility of the deduced protein to metabolize a broad range of insecticide substrates including DDT, deltamethrin, permethrin, and imidacloprid. CONCLUSIONS: We developed significant molecular resources for C. lectularius putatively involved in metabolic resistance as well as those participating in other modes of insecticide resistance. RNA-Seq profiles of PR strains combined with tissue-specific profiles and molecular docking revealed multi-level insecticide resistance in C. lectularius. Future research that is targeted towards RNA interference (RNAi) on the identified metabolic targets such as cytochrome P450s and cuticular proteins could lay the foundation for a better understanding of the genetic basis of insecticide resistance in C. lectularius.