Correction: Mendelian and Non-Mendelian Regulation of Gene Expression in Maize.

[This corrects the article DOI: 10.1371/journal.pgen.1003202.].

Evaluation of variant calling tools for large plant genome re-sequencing.

BACKGROUND: Discovering single nucleotide polymorphisms (SNPs) from agriculture crop genome sequences has been a widely used strategy for developing genetic markers for several applications including marker-assisted breeding, population diversity studies for eco-geographical adaption, genotyping crop germplasm collections, and others. Accurately detecting SNPs from large polyploid crop genomes such as wheat is crucial and challenging. A few variant calling methods have been previously developed but they show a low concordance between their variant calls. A gold standard of variant sets generated from one human individual sample was established for variant calling tool evaluations, however hitherto no gold standard of crop variant set is available for wheat use. The intent of this study was to evaluate seven SNP variant calling tools (FreeBayes, GATK, Platypus, Samtools/mpileup, SNVer, VarScan, VarDict) with the two most popular mapping tools (BWA-mem and Bowtie2) on wheat whole exome capture (WEC) re-sequencing data from allohexaploid wheat. RESULTS: We found the BWA-mem mapping tool had both a higher mapping rate and a higher accuracy rate than Bowtie2. With the same mapping quality (MQ) cutoff, BWA-mem detected more variant bases in mapping reads than Bowtie2. The reads preprocessed with quality trimming or duplicate removal did not significantly affect the final mapping performance in terms of mapped reads. Based on the concordance and receiver operating characteristic (ROC), the Samtools/mpileup variant calling tool with BWA-mem mapping of raw sequence reads outperformed other tests followed by FreeBayes and GATK in terms of specificity and sensitivity. VarDict and VarScan were the poorest performing variant calling tools with the wheat WEC sequence data. CONCLUSION: The BWA-mem and Samtools/mpileup pipeline, with no need to preprocess the raw read data before mapping onto the reference genome, was ascertained the optimum for SNP calling for the complex wheat genome re-sequencing. These results also provide useful guidelines for reliable variant identification from deep sequencing of other large polyploid crop genomes.

Atypical chromatin structure of immune-related genes expressed in chicken erythrocytes.

The major biological role of red blood cells is to carry oxygen to the tissues in the body. However, another role of the erythroid cell is to participate in the immune response. Mature erythrocytes from chickens express Toll-like receptors and several cytokines in response to stimulation of the immune system. We previously reported the application of a biochemical fractionation protocol to isolate highly enriched transcribed DNA from polychromatic erythrocytes from chickens. In conjunction with next-generation DNA, RNA sequencing, chromatin immunoprecipitation-DNA sequencing, and formaldehyde-assisted isolation of regulatory elements (FAIRE) sequencing, we identified the active chromosomal compartments and determined their structural signatures in relation to expression levels. Here, we present the detailed chromatin characteristics of erythroid genes participating in the innate immune response. Our studies revealed an atypical chromatin structure for several genes coding for Toll-like receptors, interleukins, and interferon regulatory factors. The body of these genes had nucleosome-free regions intermingled with nucleosomes modified with H3K4me3 and H3K27ac, suggesting a dynamic unstable chromatin structure. We further show that human genes involved in cell identity have gene bodies with the same chromatin-instability features as the chicken polychromatic erythrocyte genes participating in the innate immune response.

Visible-Light Photocatalytic Activity of Fe and/or Ni Doped Ilmenite Derived-Titanium Dioxide Nanoparticles.

Pure TiO2 nanoparticles and ones doped with Fe and/or Ni were successfully prepared by a co-precipitation method from ilmenite. The samples were structurally characterized by XRD, XPS, FT-IR, UV-vis, SEM, EDX, AAS and BET measurement. The XRD results showed that all samples were anatase TiO2, and no characteristic peaks of dopants were observed. The crystallite sizes of all doped TiO2 nanoparticles were less than 20 nm and doping TiO2 with metal ions can suppress the crystal growth of the particles. The XRD and XPS results indicated that TiO2 was uniformly doped and its crystalline phase was not changed by doping. The specific surface area of Fe-Ni/TiO2 is bigger than that of the un-doped TiO2. The pore size and pore volume of Fe-Ni/TiO2 is smaller than that of the un-doped. UV-vis spectra of the samples showed that the absorption edge red shifted with increasing doped metal content. The photocatalytic activity was evaluated in oxidative degradation of methylene blue (MB) with H2O2 under visible light irradiation. When doped with a single type of transition metal, the photocatalytic performance of Ni-doped samples was lower than that of Fe-doped ones. For the co-doped catalysts, the catalytic efficiency of 0.5%Fe4%Ni/TiO2 was the highest, reaching 93.34% after 250 min. Metal doping enhanced the photocatalytic decomposition of methylene blue compared with that of pure TiO2 by up to 1.5 times. The synergistic effects of the two metal ions improved the photocatalytic performance. The particles exhibited pronounced activity in degradation of MB as well as efficient recyclability. The photocatalytic degradation mechanism of methylene blue was analyzed.

CuO modified vanadium-based SCR catalysts for Hg0 oxidation and NO reduction.

The catalytic performance of Hg0 oxidation over vanadium-based SCR catalysts modified by different addition amounts of CuO was investigated. All catalysts were prepared by impregnation method and characterized. The 7% Cu/VWTi exhibited high Hg0 oxidation as well as a desired NO removal efficiency at 280-360 °C. The characterization revealed the enhancement of redox properties and well-dispersed active species results in the high catalytic performance after modification. The incorporation model showed that CuO in 7% Cu/VWTi was present in the monolayer dispersion, leading to the highest performance. Moreover, the effects of O2, NO, SO2, NH3 and HCl were explored. It showed all flue gas except NH3 could promote Hg0 oxidation. Fortunately, the inhibiting effect of NH3 could be scavenged if the catalyst is installed at the downstream of the SCR reactor. In addition, the mechanism of Hg0 oxidation over Cu/VWTi was discussed.

Genes and functions from breast cancer signatures.

BACKGROUND: Breast cancer is a heterogeneous disease and personalized medicine is the hope for the improvement of the clinical outcome. Multi-gene signatures for breast cancer stratification have been extensively studied in the past decades and more than 30 different signatures have been reported. A major concern is the minimal overlap of genes among the reported signatures. We investigated the breast cancer signature genes to address our hypothesis that the genes of different signature may share common functions, as well as to use these previously reported signature genes to build better prognostic models. METHODS: A total of 33 signatures and the corresponding gene lists were investigated. We first examined the gene frequency and the gene overlap in these signatures. Then the gene functions of each signature gene list were analysed and compared by the KEGG pathways and gene ontology (GO) terms. A classifier built using the common genes was tested using the METABRIC (Molecular Taxonomy of Breast Cancer International Consortium) data. The common genes were also tested for building the Yin Yang gene mean expression ratio (YMR) signature using public datasets (GSE1456 and GSE2034). RESULTS: Among a total of 2239 genes collected from the 33 breast cancer signatures, only 238 genes overlapped in at least two signatures; while from a total of 1979 function terms enriched in the 33 signature gene lists, 429 terms were common in at least two signatures. Most of the common function terms were involved in cell cycle processes. While there is almost no common overlapping genes between signatures developed for ER-positive (e.g. 21-gene signature) and those developed for ER-negative (e.g. basal signatures) tumours, they have common function terms such as cell death, regulation of cell proliferation. We used the 62 genes that were common in at least three signatures as a classifier and subtyped 1141 METABRIC cases including 144 normal samples into nine subgroups. These subgroups showed different clinical outcome. Among the 238 common genes, we selected those genes that are more highly expressed in normal breast tissue than in tumours as Yang genes and those more highly expressed in tumours than in normal as Yin genes and built a YMR model signature. This YMR showed significance in risk stratification in two datasets (GSE1456 and GSE2034). CONCLUSIONS: The lack of significant numbers of overlapping genes among most breast cancer signatures can be partially explained by our discovery that these signature genes represent groups with similar functions. The genes collected from these previously reported signatures are valuable resources for new model development. The subtype classifier and YMR signature built from the common genes showed promising results.

The exploration of contrasting pathways in Triple Negative Breast Cancer (TNBC).

BACKGROUND: Triple Negative Breast Cancers (TNBCs) lack the appropriate targets for currently used breast cancer therapies, conferring an aggressive phenotype, more frequent relapse and poorer survival rates. The biological heterogeneity of TNBC complicates the clinical treatment further. We have explored and compared the biological pathways in TNBC and other subtypes of breast cancers, using an in silico approach and the hypothesis that two opposing effects (Yin and Yang) pathways in cancer cells determine the fate of cancer cells. Identifying breast subgroup specific components of these opposing pathways may aid in selecting potential therapeutic targets as well as further classifying the heterogeneous TNBC subtype. METHODS: Gene expression and patient clinical data from The Cancer Genome Atlas (TCGA) and the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) were used for this study. Gene Set Enrichment Analysis (GSEA) was used to identify the more active pathways in cancer (Yin) than in normal and the more active pathways in normal (Yang) than in cancer. The clustering analysis was performed to compare pathways of TNBC with other types of breast cancers. The association of pathway classified TNBC sub-groups to clinical outcomes was tested using Cox regression model. RESULTS: Among 4729 curated canonical pathways in GSEA database, 133 Yin pathways (FDR < 0.05) and 71 Yang pathways (p-value <0.05) were discovered in TNBC. The FOXM1 is the top Yin pathway while PPARα is the top Yang pathway in TNBC. The TNBC and other types of breast cancers showed different pathways enrichment significance profiles. Using top Yin and Yang pathways as classifier, the TNBC can be further subtyped into six sub-groups each having different clinical outcomes. CONCLUSION: We first reported that the FOMX1 pathway is the most upregulated and the PPARα pathway is the most downregulated pathway in TNBC. These two pathways could be simultaneously targeted in further studies. Also the pathway classifier we performed in this study provided insight into the TNBC heterogeneity.

A 16 Yin Yang gene expression ratio signature for ER+/node- breast cancer.

Breast cancer is one of the leading causes of cancer death in women. It is a complex and heterogeneous disease with different clinical outcomes. Stratifying patients into subgroups with different outcomes could help guide clinical decision making. In this study, we used two opposing groups of genes, Yin and Yang, to develop a prognostic expression ratio signature. Using the METABRIC cohort we identified a16-gene signature capable of stratifying breast cancer patients into four risk levels with intention that low-risk patients would not undergo adjuvant systemic therapy, intermediate-low-risk patients will be treated with hormonal therapy only, and intermediate-high- and high-risk groups will be treated by chemotherapy in addition to the hormonal therapy. The 16-gene signature for four risk level stratifications of breast cancer patients has been validated using 14 independent datasets. Notably, the low-risk group (n = 51) of 205 estrogen receptor-positive and node negative (ER+/node-) patients from three different datasets who had not had any systemic adjuvant therapy had 100% 15-year disease-specific survival rate. The Concordance Index of YMR for ER+/node negative patients is close to the commercially available signatures. However, YMR showed more significance (HR = 3.7, p = 8.7e-12) in stratifying ER+/node- subgroup than OncotypeDx (HR = 2.7, p = 1.3e-7), MammaPrint (HR = 2.5, p = 5.8e-7), rorS (HR = 2.4, p = 1.4e-6), and NPI (HR = 2.6, p = 1.2e-6). YMR signature may be developed as a clinical tool to select a subgroup of low-risk ER+/node- patients who do not require any adjuvant hormonal therapy (AHT).

Histone H3K4 trimethylation: dynamic interplay with pre-mRNA splicing.

Histone H3 lysine 4 trimethylation (H3K4me3) is often stated as a mark of transcriptionally active promoters. However, closer study of the positioning of H3K4me3 shows the mark locating primarily after the first exon at the 5' splice site and overlapping with a CpG island in mammalian cells. There are several enzyme complexes that are involved in the placement of the H3K4me3 mark, including multiple protein complexes containing SETD1A, SETD1B, and MLL1 enzymes (writers). CXXC1, which is associated with SETD1A and SETD1B, target these enzymes to unmethylated CpG islands. Lysine demethylases (KDM5 family members, erasers) demethylate H3K4me3. The H3K4me3 mark is recognized by several proteins (readers), including lysine acetyltransferase complexes, chromatin remodelers, and RNA bound proteins involved in pre-mRNA splicing. Interestingly, attenuation of H3K4me3 impacts pre-mRNA splicing, and inhibition of pre-mRNA splicing attenuates H3K4me3.

Gene expression identifies heterogeneity of metastatic behavior among gastrointestinal stromal tumors.

BACKGROUND: Adjuvant imatinib is useful in patients with gastrointestinal stromal tumors (GIST) at high risk of recurrence. At present, the risk of recurrence is determined based on tumor size, mitotic rate, tumor site, and tumor rupture. Previous studies using various biochemical pathways identified gene expression patterns that distinguish two subsets of aggressive fibromatosis (AF), serous ovarian carcinoma (OVCA), and clear cell renal cell carcinoma (RCC). These gene sets separated soft tissue sarcomas into two groups with different probabilities of developing metastatic disease. The present study used these gene sets to identify GIST subgroups with different probabilities of developing metastatic disease. METHODS: We utilized these three gene sets, hierarchical clustering, and Kaplan-Meier analysis, to examine 60 primary resected GIST samples using Agilent chip expression profiling. RESULTS: Hierarchical clustering using both the combined and individual AF-, OVCA-, and RCC- gene sets identified differences in probabilities of developing metastatic disease between the clusters defined by the first branch point of the clustering dendrograms (p = 0.029 for the combined gene set, p = 0.003 for the AF-gene set, p < 0.001 for the OVCA-gene set, and p = 0.003 for the RCC-gene set). CONCLUSIONS: Hierarchical clustering using these gene sets identified at least two subsets of GIST with distinct clinical behavior and risk of metastatic disease. The use of gene expression analysis along with other known prognostic factors may better predict the long-term outcome following surgery, and thus restrict the use of adjuvant therapy to high-risk GIST, and reduce heterogeneity among groups in clinical trials of new drugs.

The transcriptome of nitrofen-induced pulmonary hypoplasia in the rat model of congenital diaphragmatic hernia.

BACKGROUND: We currently do not know how the herbicide nitrofen induces lung hypoplasia and congenital diaphragmatic hernia in rats. Our aim was to compare the differentially expressed transcriptome of nitrofen-induced hypoplastic lungs to control lungs in embryonic day 13 rat embryos before the development of embryonic diaphragmatic defects. METHODS: Using next-generation sequencing technology, we identified the expression profile of microRNA (miRNA) and mRNA genes. Once the dataset was validated by both RT-qPCR and digital-PCR, we conducted gene ontology, miRNA target analysis, and orthologous miRNA sequence matching for the deregulated miRNAs in silico. RESULTS: Our study identified 186 known mRNA and 100 miRNAs which were differentially expressed in nitrofen-induced hypoplastic lungs. Sixty-four rat miRNAs homologous to known human miRNAs were identified. A subset of these genes may promote lung hypoplasia in rat and/or human, and we discuss their associations. Potential miRNA pathways relevant to nitrofen-induced lung hypoplasia include PI3K, TGF-ß, and cell cycle kinases. CONCLUSION: Nitrofen-induced hypoplastic lungs have an abnormal transcriptome that may lead to impaired development.

RNA-dependent dynamic histone acetylation regulates MCL1 alternative splicing.

Histone deacetylases (HDACs) and lysine acetyltransferases (KATs) catalyze dynamic histone acetylation at regulatory and coding regions of transcribed genes. Highly phosphorylated HDAC2 is recruited within corepressor complexes to regulatory regions, while the nonphosphorylated form is associated with the gene body. In this study, we characterized the nonphosphorylated HDAC2 complexes recruited to the transcribed gene body and explored the function of HDAC-complex-mediated dynamic histone acetylation. HDAC1 and 2 were coimmunoprecipitated with several splicing factors, including serine/arginine-rich splicing factor 1 (SRSF1) which has roles in alternative splicing. The co-chromatin immunoprecipitation of HDAC1/2 and SRSF1 to the gene body was RNA-dependent. Inhibition of HDAC activity and knockdown of HDAC1, HDAC2 or SRSF1 showed that these proteins were involved in alternative splicing of MCL1. HDAC1/2 and KAT2B were associated with nascent pre-mRNA in general and with MCL1 pre-mRNA specifically. Inhibition of HDAC activity increased the occupancy of KAT2B and acetylation of H3 and H4 of the H3K4 methylated alternative MCL1 exon 2 nucleosome. Thus, nonphosphorylated HDAC1/2 is recruited to pre-mRNA by splicing factors to act at the RNA level with KAT2B and other KATs to catalyze dynamic histone acetylation of the MCL1 alternative exon and alter the splicing of MCL1 pre-mRNA.

Mendelian and non-Mendelian regulation of gene expression in maize.

Transcriptome variation plays an important role in affecting the phenotype of an organism. However, an understanding of the underlying mechanisms regulating transcriptome variation in segregating populations is still largely unknown. We sought to assess and map variation in transcript abundance in maize shoot apices in the intermated B73 × Mo17 recombinant inbred line population. RNA-based sequencing (RNA-seq) allowed for the detection and quantification of the transcript abundance derived from 28,603 genes. For a majority of these genes, the population mean, coefficient of variation, and segregation patterns could be predicted by the parental expression levels. Expression quantitative trait loci (eQTL) mapping identified 30,774 eQTL including 96 trans-eQTL "hotspots," each of which regulates the expression of a large number of genes. Interestingly, genes regulated by a trans-eQTL hotspot tend to be enriched for a specific function or act in the same genetic pathway. Also, genomic structural variation appeared to contribute to cis-regulation of gene expression. Besides genes showing Mendelian inheritance in the RIL population, we also found genes whose expression level and variation in the progeny could not be predicted based on parental difference, indicating that non-Mendelian factors also contribute to expression variation. Specifically, we found 145 genes that show patterns of expression reminiscent of paramutation such that all the progeny had expression levels similar to one of the two parents. Furthermore, we identified another 210 genes that exhibited unexpected patterns of transcript presence/absence. Many of these genes are likely to be gene fragments resulting from transposition, and the presence/absence of their transcripts could influence expression levels of their ancestral syntenic genes. Overall, our results contribute to the identification of novel expression patterns and broaden the understanding of transcriptional variation in plants.

Triclosan can select for an AdeIJK-overexpressing mutant of Acinetobacter baumannii ATCC 17978 that displays reduced susceptibility to multiple antibiotics.

In order to determine if triclosan can select for mutants of Acinetobacter baumannii ATCC 17978 that display reduced susceptibilities to antibiotics, we isolated a triclosan-resistant mutant, A. baumannii AB042, by serial passaging of A. baumannii ATCC 17978 in growth medium supplemented with triclosan. The antimicrobial susceptibility of AB042 was analyzed by the 2-fold serial dilution method. Expression of five different resistance-nodulation-division (RND) pump-encoding genes (adeB, adeG, adeJ, A1S_2818, and A1S_3217), two outer membrane porin-encoding genes (carO and oprD), and the MATE family pump-encoding gene abeM was analyzed using quantitative reverse transcriptase (qRT) PCR. A. baumannii AB042 exhibited elevated resistance to multiple antibiotics, including piperacillin-tazobactam, doxycycline, moxifloxacin, ceftriaxone, cefepime, meropenem, doripenem, ertapenem, ciprofloxacin, aztreonam, tigecycline, and trimethoprim-sulfamethoxazole, in addition to triclosan. Genome sequencing of A. baumannii AB042 revealed a (116)G→V mutation in fabI, the gene encoding the target enzyme for triclosan. Expression analysis of efflux pumps showed overexpression of the AdeIJK pump, and sequencing of adeN, the gene that encodes the repressor of the adeIJK operon, revealed a 73-bp deletion which would cause a premature termination of translation, resulting in an inactive truncated AdeN protein. This work shows that triclosan can select for mutants of A. baumannii that display reduced susceptibilities to multiple antibiotics from chemically distinct classes in addition to triclosan resistance. This multidrug resistance can be explained by the overexpression of the AdeIJK efflux pump.

Gene expression identifies heterogeneity of metastatic behavior among high-grade non-translocation associated soft tissue sarcomas.

BACKGROUND: The biologic heterogeneity of soft tissue sarcomas (STS), even within histological subtypes, complicates treatment. In earlier studies, gene expression patterns that distinguish two subsets of clear cell renal carcinoma (RCC), serous ovarian carcinoma (OVCA), and aggressive fibromatosis (AF) were used to separate 73 STS into two or four groups with different probabilities of developing metastatic disease (PrMet). This study was designed to confirm our earlier observations in a larger independent data set. METHODS: We utilized these gene sets, hierarchical clustering (HC), and Kaplan-Meier analysis, to examine 309 STS, using Affymetrix chip expression profiling. RESULTS: HC using the combined AF-, RCC-, and OVCA-gene sets identified subsets of the STS samples. Analysis revealed differences in PrMet between the clusters defined by the first branch point of the clustering dendrogram (p = 0.048), and also among the four different clusters defined by the second branch points (p < 0.0001). Analysis also revealed differences in PrMet between the leiomyosarcomas (LMS), dedifferentiated liposarcomas (LipoD), and undifferentiated pleomorphic sarcomas (UPS) (p = 0.0004). HC of both the LipoD and UPS sample sets divided the samples into two groups with different PrMet (p = 0.0128, and 0.0002, respectively). HC of the UPS samples also showed four groups with different PrMet (p = 0.0007). HC found no subgroups of the LMS samples. CONCLUSIONS: These data confirm our earlier studies, and suggest that this approach may allow the identification of more than two subsets of STS, each with distinct clinical behavior, and may be useful to stratify STS in clinical trials and in patient management.

Parent-of-origin effects on gene expression and DNA methylation in the maize endosperm.

Imprinting describes the differential expression of alleles based on their parent of origin. Deep sequencing of RNAs from maize (Zea mays) endosperm and embryo tissue 14 d after pollination was used to identify imprinted genes among a set of ~12,000 genes that were expressed and contained sequence polymorphisms between the B73 and Mo17 genotypes. The analysis of parent-of-origin patterns of expression resulted in the identification of 100 putative imprinted genes in maize endosperm, including 54 maternally expressed genes (MEGs) and 46 paternally expressed genes (PEGs). Three of these genes have been previously identified as imprinted, while the remaining 97 genes represent novel imprinted maize genes. A genome-wide analysis of DNA methylation identified regions with reduced endosperm DNA methylation in, or near, 19 of the 100 imprinted genes. The reduced levels of DNA methylation in endosperm are caused by hypomethylation of the maternal allele for both MEGs and PEGs in all cases tested. Many of the imprinted genes with reduced DNA methylation levels also show endosperm-specific expression patterns. The imprinted maize genes were compared with imprinted genes identified in genome-wide screens of rice (Oryza sativa) and Arabidopsis thaliana, and at least 10 examples of conserved imprinting between maize and each of the other species were identified.

Host-Mycobacterium avium subsp. paratuberculosis interactome reveals a novel iron assimilation mechanism linked to nitric oxide stress during early infection.

BACKGROUND: The initial interaction between host cell and pathogen sets the stage for the ensuing infection and ultimately determine the course of disease. However, there is limited knowledge of the transcripts utilized by host and pathogen and how they may impact one another during this critical step. The purpose of this study was to create a host-Mycobacterium avium subsp. paratuberculosis (MAP) interactome for early infection in an epithelium-macrophage co-culture system using RNA-seq. RESULTS: Establishment of the host-MAP interactome revealed a novel iron assimilation system for carboxymycobactin. Iron assimilation is linked to nitric oxide synthase-2 production by the host and subsequent nitric oxide buildup. Iron limitation as well as nitric oxide is a prompt for MAP to enter into an iron sequestration program. This new iron sequestration program provides an explanation for mycobactin independence in some MAP strains grown in vitro as well as during infection within the host cell. Utilization of such a pathway is likely to aid MAP establishment and long-term survival within the host. CONCLUSIONS: The host-MAP interactome identified a number of metabolic, DNA repair and virulence genes worthy for consideration as novel drug targets as well as future pathogenesis studies. Reported interactome data may also be utilized to conduct focused, hypothesis-driven research. Co-culture of uninfected bovine epithelial cells (MAC-T) and primary bovine macrophages creates a tolerant genotype as demonstrated by downregulation of inflammatory pathways. This co-culture system may serve as a model to investigate other bovine enteric pathogens.

The barley UNICULM2 gene resides in a centromeric region and may be associated with signaling and stress responses.

Vegetative axillary meristem (AXM) activity results in the production of branches. In barley (Hordeum vulgare L.), vegetative AXM develop in the crown and give rise to modified branches, referred to as tillers. Mutations in the barley low-tillering mutant uniculm2 block vegetative AXM development and prevent tiller development. The objectives of this work were to examine gene expression in wild-type and cul2 mutant plants, fine map the CUL2 gene, and to examine synteny in the CUL2 region in barley with rice. RNA profiling experiments using two near-isogenic line pairs carrying either the cul2 mutant allele or wild-type CUL2 allele in different genetic backgrounds detected 28 unique gene transcripts exhibiting similar patterns of differential accumulation in both genetic backgrounds, indicating that we have identified key genes impacted by the CUL2 gene. Twenty-four genes had higher abundance in uniculm2 mutant tissues, and nearly half of the annotated genes likely function in stress-response or signal transduction pathways. Genetic mapping identified five co-segregating markers in 1,088 F2 individuals. These markers spanned the centromere region on chromosome 6H, and coincided with a 50-cM region on rice chromosome 2, indicating that it may be difficult to positionally clone CUL2. Taken together, the results revealed stress response and signal transduction pathways that are associated with the CUL2 gene, isolating CUL2 via positional cloning approaches that may be difficult, and the remnants of barley-rice synteny in the CUL2 region.

Suppression of the interferon and NF-κB responses by severe fever with thrombocytopenia syndrome virus.

Severe fever with thrombocytopenia syndrome (SFTS) is an emerging infectious disease characterized by high fever, thrombocytopenia, multiorgan dysfunction, and a high fatality rate between 12 and 30%. It is caused by SFTS virus (SFTSV), a novel Phlebovirus in family Bunyaviridae. Although the viral pathogenesis remains largely unknown, hemopoietic cells appear to be targeted by the virus. In this study we report that human monocytes were susceptible to SFTSV, which replicated efficiently, as shown by an immunofluorescence assay and real-time reverse transcription-PCR. We examined host responses in the infected cells and found that antiviral interferon (IFN) and IFN-inducible proteins were induced upon infection. However, our data also indicated that downregulation of key molecules such as mitochondrial antiviral signaling protein (MAVS) or weakened activation of interferon regulatory factor (IRF) and NF-κB responses may contribute to a restricted innate immunity against the infection. NSs, the nonstructural protein encoded by the S segment, suppressed the beta interferon (IFN-ß) and NF-κB promoter activities, although NF-κB activation appears to facilitate SFTSV replication in human monocytes. NSs was found to be associated with TBK1 and may inhibit the activation of downstream IRF and NF-κB signaling through this interaction. Interestingly, we demonstrated that the nucleoprotein (N), also encoded by the S segment, exhibited a suppressive effect on the activation of IFN-ß and NF-κB signaling as well. Infected monocytes, mainly intact and free of apoptosis, may likely be implicated in persistent viral infection, spreading the virus to the circulation and causing primary viremia. Our findings provide the first evidence in dissecting the host responses in monocytes and understanding viral pathogenesis in humans infected with a novel deadly Bunyavirus.

Structural variants in the soybean genome localize to clusters of biotic stress-response genes.

Genome-wide structural and gene content variations are hypothesized to drive important phenotypic variation within a species. Structural and gene content variations were assessed among four soybean (Glycine max) genotypes using array hybridization and targeted resequencing. Many chromosomes exhibited relatively low rates of structural variation (SV) among genotypes. However, several regions exhibited both copy number and presence-absence variation, the most prominent found on chromosomes 3, 6, 7, 16, and 18. Interestingly, the regions most enriched for SV were specifically localized to gene-rich regions that harbor clustered multigene families. The most abundant classes of gene families associated with these regions were the nucleotide-binding and receptor-like protein classes, both of which are important for plant biotic defense. The colocalization of SV with plant defense response signal transduction pathways provides insight into the mechanisms of soybean resistance gene evolution and may inform the development of new approaches to resistance gene cloning.