PTEN deletion drives aberrations of DNA methylome and transcriptome in different stages of prostate cancer.

Phosphatase and tensin homolog located on chromosome 10 (PTEN) is a tumor suppressor gene and one of the most frequently mutated/deleted genes in human prostate cancer (PCa). However, how PTEN deletion would impact the epigenome and transcriptome alterations remain unknown. This hypothesis was tested in a prostate-specific PTEN-/- (KO) mouse prostatic adenocarcinoma model through DNA methyl-Seq and RNA-Seq analyses. Examination of cancer genomic datasets revealed that PTEN is expressed at lower levels in PTEN-deleted tumor samples than in normal solid tissue samples. Methylome and transcriptome profiling identified several inflammatory responses and immune response signaling pathways, including NF-kB signaling, IL-6 signaling, LPS/IL-1-mediated inhibition of RXR Function, PI3K in B lymphocytes, iCOS-iCOSL in T helper cells, and the role of NFAT in regulating the immune response, were affected by PTEN deletion. Importantly, a small subset of genes that showed DNA hypermethylation or hypomethylation was correlated with decreased or increased gene expression including CXCL1. quantitative polymerase chain reaction analyses of representative genes validated the RNA-Seq results. Histopathological examinations showed that the severity of prostatic intraepithelial neoplasia and inflammation development gradually increased as PTEN null mice aged. Collectively, these findings suggest that loss of PTEN drives global changes in DNA CpG methylation and transcriptomic gene expression and highly associated with several inflammatory and immune molecular pathways during PCa development. These biomarkers could be valuable molecular targets for cancer drug discovery and development against PCa.

Candidate gene identification of existing or induced mutations with pipelines applicable to large genomes.

Bulked segregant analysis (BSA) is used to identify existing or induced variants that are linked to phenotypes. Although it is widely used in Arabidopsis and rice, it remains challenging for crops with large genomes, such as maize. Moreover, analysis of huge data sets can present a bottleneck linking phenotypes to their molecular basis, especially for geneticists without programming experience. Here, we identified two genes of maize defective kernel mutants with newly developed analysis pipelines that require no programing skills and should be applicable to any large genome. In the 1970s, Neuffer and Sheridan generated a chemically induced defective kernel (dek) mutant collection with the potential to uncover critical genes for seed development. To locate such mutations, the dek phenotypes were introgressed into two inbred lines to take advantage of maize haplotype variations and their sequenced genomes. We generated two pipelines that take fastq files derived from next-generation (nextGen) paired-end DNA and cDNA sequencing as input, call on several well established and freely available genomic analysis tools to call SNPs and INDELs, and generate lists of the most likely causal mutations together with variant index plots to locate the mutation to a specific sequence position on a chromosome. The pipelines were validated with a known strawberry mutation before cloning the dek mutants, thereby enabling phenotypic analysis of large genomes by next-generation sequencing.

Cell-to-cell movement of mitochondria in plants.

We report cell-to-cell movement of mitochondria through a graft junction. Mitochondrial movement was discovered in an experiment designed to select for chloroplast transfer from Nicotiana sylvestris into Nicotiana tabacum cells. The alloplasmic N. tabacum line we used carries Nicotiana undulata cytoplasmic genomes, and its flowers are male sterile due to the foreign mitochondrial genome. Thus, rare mitochondrial DNA transfer from N. sylvestris to N. tabacum could be recognized by restoration of fertile flower anatomy. Analyses of the mitochondrial genomes revealed extensive recombination, tentatively linking male sterility to orf293, a mitochondrial gene causing homeotic conversion of anthers into petals. Demonstrating cell-to-cell movement of mitochondria reconstructs the evolutionary process of horizontal mitochondrial DNA transfer and enables modification of the mitochondrial genome by DNA transmitted from a sexually incompatible species. Conversion of anthers into petals is a visual marker that can be useful for mitochondrial transformation.

Analysis of tandem gene copies in maize chromosomal regions reconstructed from long sequence reads.

Haplotype variation not only involves SNPs but also insertions and deletions, in particular gene copy number variations. However, comparisons of individual genomes have been difficult because traditional sequencing methods give too short reads to unambiguously reconstruct chromosomal regions containing repetitive DNA sequences. An example of such a case is the protein gene family in maize that acts as a sink for reduced nitrogen in the seed. Previously, 41-48 gene copies of the alpha zein gene family that spread over six loci spanning between 30- and 500-kb chromosomal regions have been described in two Iowa Stiff Stalk (SS) inbreds. Analyses of those regions were possible because of overlapping BAC clones, generated by an expensive and labor-intensive approach. Here we used single-molecule real-time (Pacific Biosciences) shotgun sequencing to assemble the six chromosomal regions from the Non-Stiff Stalk maize inbred W22 from a single DNA sequence dataset. To validate the reconstructed regions, we developed an optical map (BioNano genome map; BioNano Genomics) of W22 and found agreement between the two datasets. Using the sequences of full-length cDNAs from W22, we found that the error rate of PacBio sequencing seemed to be less than 0.1% after autocorrection and assembly. Expressed genes, some with premature stop codons, are interspersed with nonexpressed genes, giving rise to genotype-specific expression differences. Alignment of these regions with those from the previous analyzed regions of SS lines exhibits in part dramatic differences between these two heterotic groups.

Using whole-exome sequencing to investigate the genetic bases of lysosomal storage diseases of unknown etiology.

Lysosomes are membrane-bound, acidic eukaryotic cellular organelles that play important roles in the degradation of macromolecules. Mutations that cause the loss of lysosomal protein function can lead to a group of disorders categorized as the lysosomal storage diseases (LSDs). Suspicion of LSD is frequently based on clinical and pathologic findings, but in some cases, the underlying genetic and biochemical defects remain unknown. Here, we performed whole-exome sequencing (WES) on 14 suspected LSD cases to evaluate the feasibility of using WES for identifying causal mutations. By examining 2,157 candidate genes potentially associated with lysosomal function, we identified eight variants in five genes as candidate disease-causing variants in four individuals. These included both known and novel mutations. Variants were corroborated by targeted sequencing and, when possible, functional assays. In addition, we identified nonsense mutations in two individuals in genes that are not known to have lysosomal function. However, mutations in these genes could have resulted in phenotypes that were diagnosed as LSDs. This study demonstrates that WES can be used to identify causal mutations in suspected LSD cases. We also demonstrate cases where a confounding clinical phenotype may potentially reflect more than one lysosomal protein defect.

The complete genome sequence of 'Candidatus Liberibacter americanus', associated with Citrus huanglongbing.

Liberibacter spp. form a Rhizobiaceae clade of phloem-limited pathogens of limited host range. Two obligately parasitic species have been sequenced: 'Candidatus Liberibacter asiaticus', which causes citrus huanglongbing (HLB) worldwide, and 'Ca. L. solanacearum', which causes potato "zebra chip" disease. A third (proposed) species, Liberibacter crescens, was isolated from mountain papaya, grown in axenic culture, and sequenced. In an effort to identify common host determinants, the complete genomic DNA sequence of a second HLB species, 'Ca. L. americanus' strain 'São Paulo' was determined. The circular genome of 1,195,201 bp had an average 31.12% GC content and 983 predicted protein encoding genes, 800 (81.4%) of which had a predicted function. There were 658 genes common to all sequenced Liberibacter spp. and only 8 genes common to 'Ca. L. americanus' and 'Ca. L. asiaticus' but not found in 'Ca. L. solanacearum'. Surprisingly, most of the lipopolysaccharide biosynthetic genes were missing from the 'Ca. L. americanus' genome, as well as OmpA and a key regulator of flagellin, all indicating a 'Ca. L. americanus' strategy of avoiding production of major pathogen-associated molecular patterns present in 'Ca. L. asiaticus' and 'Ca. L. solanacearum'. As with 'Ca. L. asiaticus', one of two 'Ca. L. americanus' prophages replicated as an excision plasmid and carried potential lysogenic conversion genes that appeared fragmentary or degenerated in 'Ca. L. solanacearum'.

Tet-dependent 5-hydroxymethyl-Cytosine modification of mRNA regulates axon guidance genes in Drosophila.

Modifications of mRNA, especially methylation of adenosine, have recently drawn much attention. The much rarer modification, 5-hydroxymethylation of cytosine (5hmC), is not well understood and is the subject of this study. Vertebrate Tet proteins are 5-methylcytosine (5mC) hydroxylases and catalyze the transition of 5mC to 5hmC in DNA. These enzymes have recently been shown to have the same function in messenger RNAs in both vertebrates and in Drosophila. The Tet gene is essential in Drosophila as Tet knock-out animals do not reach adulthood. We describe the identification of Tet-target genes in the embryo and larval brain by mapping one, Tet DNA-binding sites throughout the genome and two, the Tet-dependent 5hmrC modifications transcriptome-wide. 5hmrC modifications are distributed along the entire transcript, while Tet DNA-binding sites are preferentially located at the promoter where they overlap with histone H3K4me3 peaks. The identified mRNAs are preferentially involved in neuron and axon development and Tet knock-out led to a reduction of 5hmrC marks on specific mRNAs. Among the Tet-target genes were the robo2 receptor and its slit ligand that function in axon guidance in Drosophila and in vertebrates. Tet knock-out embryos show overlapping phenotypes with robo2 and both Robo2 and Slit protein levels were markedly reduced in Tet KO larval brains. Our results establish a role for Tet-dependent 5hmrC in facilitating the translation of modified mRNAs primarily in cells of the nervous system.

Comparative genomic and transcriptome analyses of pathotypes of Xanthomonas citri subsp. citri provide insights into mechanisms of bacterial virulence and host range.

BACKGROUND: Citrus bacterial canker is a disease that has severe economic impact on citrus industries worldwide and is caused by a few species and pathotypes of Xanthomonas. X. citri subsp. citri strain 306 (XccA306) is a type A (Asiatic) strain with a wide host range, whereas its variant X. citri subsp. citri strain A(w)12879 (Xcaw12879, Wellington strain) is restricted to Mexican lime. RESULTS: To characterize the mechanism for the differences in host range of XccA and Xcaw, the genome of Xcaw12879 that was completed recently was compared with XccA306 genome. Effectors xopAF and avrGf1 are present in Xcaw12879, but were absent in XccA306. AvrGf1 was shown previously for Xcaw to cause hypersensitive response in Duncan grapefruit. Mutation analysis of xopAF indicates that the gene contributes to Xcaw growth in Mexican lime but does not contribute to the limited host range of Xcaw. RNA-Seq analysis was conducted to compare the expression profiles of Xcaw12879 and XccA306 in Nutrient Broth (NB) medium and XVM2 medium, which induces hrp gene expression. Two hundred ninety two and 281 genes showed differential expression in XVM2 compared to in NB for XccA306 and Xcaw12879, respectively. Twenty-five type 3 secretion system genes were up-regulated in XVM2 for both XccA and Xcaw. Among the 4,370 common genes of Xcaw12879 compared to XccA306, 603 genes in NB and 450 genes in XVM2 conditions were differentially regulated. Xcaw12879 showed higher protease activity than XccA306 whereas Xcaw12879 showed lower pectate lyase activity in comparison to XccA306. CONCLUSIONS: Comparative genomic analysis of XccA306 and Xcaw12879 identified strain specific genes. Our study indicated that AvrGf1 contributes to the host range limitation of Xcaw12879 whereas XopAF contributes to virulence. Transcriptome analyses of XccA306 and Xcaw12879 presented insights into the expression of the two closely related strains of X. citri subsp. citri. Virulence genes including genes encoding T3SS components and effectors are induced in XVM2 medium. Numerous genes with differential expression in Xcaw12879 and XccA306 were identified. This study provided the foundation to further characterize the mechanisms for virulence and host range of pathotypes of X. citri subsp. citri.

Impaired auxin biosynthesis in the defective endosperm18 mutant is due to mutational loss of expression in the ZmYuc1 gene encoding endosperm-specific YUCCA1 protein in maize.

The phytohormone auxin (indole-3-acetic acid [IAA]) plays a fundamental role in vegetative and reproductive plant development. Here, we characterized a seed-specific viable maize (Zea mays) mutant, defective endosperm18 (de18) that is impaired in IAA biosynthesis. de18 endosperm showed large reductions of free IAA levels and is known to have approximately 40% less dry mass, compared with De18. Cellular analyses showed lower total cell number, smaller cell volume, and reduced level of endoreduplication in the mutant endosperm. Gene expression analyses of seed-specific tryptophan-dependent IAA pathway genes, maize Yucca1 (ZmYuc1), and two tryptophan-aminotransferase co-orthologs were performed to understand the molecular basis of the IAA deficiency in the mutant. Temporally, all three genes showed high expression coincident with high IAA levels; however, only ZmYuc1 correlated with the reduced IAA levels in the mutant throughout endosperm development. Furthermore, sequence analyses of ZmYuc1 complementary DNA and genomic clones revealed many changes specific to the mutant, including a 2-bp insertion that generated a premature stop codon and a truncated YUC1 protein of 212 amino acids, compared with the 400 amino acids in the De18. The putative, approximately 1.5-kb, Yuc1 promoter region also showed many rearrangements, including a 151-bp deletion in the mutant. Our concurrent high-density mapping and annotation studies of chromosome 10, contig 395, showed that the De18 locus was tightly linked to the gene ZmYuc1. Collectively, the data suggest that the molecular changes in the ZmYuc1 gene encoding the YUC1 protein are the causal basis of impairment in a critical step in IAA biosynthesis, essential for normal endosperm development in maize.

Wolf-Hirschhorn Syndrome with Hyperparathyroidism: A Case Report and a Narrative Review of the Literature.

Wolf-Hirschhorn syndrome (WHS) is a contiguous gene deletion condition. The WHS core phenotype includes developmental delays, intellectual disabilities, seizures, and distinctive facial features. Various other comorbidities have also been reported, such as hearing loss, heart defects, as well as eye problems and kidney problems. In this report, we present a case of WHS accompanied by hyperparathyroidism and hypercalcemia, which has not been previously reported. A girl was born at 37 weeks of gestation by vaginal delivery. She was small for the gestational age (2,045 g) and admitted to neonatal intensive care unit. She had typical WHS facial features and was found to have bilateral small kidneys associated with transient metabolic acidosis and renal insufficiency. She had right-sided sensorineural hearing loss, a small atrial septal defect, and colpocephaly and hypoplasia of corpus callosum. She had a single seizure which was well controlled with an oral antiepileptic medication. Cytogenetic studies demonstrated a large terminal chromosome 4p deletion (21.4 Mb) and 4p duplication (2.1 Mb) adjacent to the deletion. A unique finding in this patient is her consistently elevated levels of parathyroid hormone and serum calcium, suggesting hyperparathyroidism. We present this rare case along with a review of the literature and hope to draw an attention to a potential relationship between WHS and hyperparathyroidism.

Tet-dependent 5-hydroxymethyl-Cytosine modification of mRNA regulates axon guidance genes in Drosophila.

Modifications of mRNA, especially methylation of adenosine, have recently drawn much attention. The much rarer modification, 5-hydroxymethylation of cytosine (5hmC), is not well understood and is the subject of this study. Vertebrate Tet proteins are 5-methylcytosine (5mC) hydroxylases and catalyze the transition of 5mC to 5hmC in DNA. These enzymes have recently been shown to have the same function in messenger RNAs in both vertebrates and in Drosophila. The Tet gene is essential in Drosophila as Tet knock-out animals do not reach adulthood. We describe the identification of Tet-target genes in the embryo and larval brain by mapping one, Tet DNA-binding sites throughout the genome and two, the Tet-dependent 5hmrC modifications transcriptome-wide. 5hmrC modifications are distributed along the entire transcript, while Tet DNA-binding sites are preferentially located at the promoter where they overlap with histone H3K4me3 peaks. The identified mRNAs are preferentially involved in neuron and axon development and Tet knock-out led to a reduction of 5hmrC marks on specific mRNAs. Among the Tet-target genes were the robo2 receptor and its slit ligand that function in axon guidance in Drosophila and in vertebrates. Tet knock-out embryos show overlapping phenotypes with robo2 and both Robo2 and Slit protein levels were markedly reduced in Tet KO larval brains. Our results establish a role for Tet-dependent 5hmrC in facilitating the translation of modified mRNAs primarily in cells of the nervous system.

Tet-dependent 5-hydroxymethyl-Cytosine modification of mRNA regulates the axon guidance genes robo2 and slit in Drosophila.

Modifications of mRNA, especially methylation of adenosine, have recently drawn much attention. The much rarer modification, 5-hydroxymethylation of cytosine (5hmC), is not well understood and is the subject of this study. Vertebrate Tet proteins are 5-methylcytosine (5mC) hydroxylases enzymes catalyzing the transition of 5mC to 5hmC in DNA and have recently been shown to have the same function in messenger RNAs in both vertebrates and in Drosophila. The Tet gene is essential in Drosophila because Tet knock-out animals do not reach adulthood. We describe the identification of Tet-target genes in the embryo and larval brain by determining Tet DNA-binding sites throughout the genome and by mapping the Tet-dependent 5hmrC modifications transcriptome-wide. 5hmrC-modified sites can be found along the entire transcript and are preferentially located at the promoter where they overlap with histone H3K4me3 peaks. The identified mRNAs are frequently involved in neuron and axon development and Tet knock-out led to a reduction of 5hmrC marks on specific mRNAs. Among the Tet-target genes were the robo2 receptor and its slit ligand that function in axon guidance in Drosophila and in vertebrates. Tet knock-out embryos show overlapping phenotypes with robo2 and are sensitized to reduced levels of slit. Both Robo2 and Slit protein levels were markedly reduced in Tet KO larval brains. Our results establish a role for Tet-dependent 5hmrC in facilitating the translation of modified mRNAs, primarily in developing nerve cells.

The Xylella fastidiosa biocontrol strain EB92-1 genome is very similar and syntenic to Pierce's disease strains.

Xylella fastidiosa infects a wide range of plant hosts and causes economically serious diseases, including Pierce's disease (PD) of grapevines. X. fastidiosa biocontrol strain EB92-1 is infectious to grapevines but does not cause symptoms. The draft genome of EB92-1 reveals that it may be missing 10 potential pathogenicity effectors.

Comparative genomic analysis of Xanthomonas axonopodis pv. citrumelo F1, which causes citrus bacterial spot disease, and related strains provides insights into virulence and host specificity.

Xanthomonas axonopodis pv. citrumelo is a citrus pathogen causing citrus bacterial spot disease that is geographically restricted within the state of Florida. Illumina, 454 sequencing, and optical mapping were used to obtain a complete genome sequence of X. axonopodis pv. citrumelo strain F1, 4.9 Mb in size. The strain lacks plasmids, in contrast to other citrus Xanthomonas pathogens. Phylogenetic analysis revealed that this pathogen is very close to the tomato bacterial spot pathogen X. campestris pv. vesicatoria 85-10, with a completely different host range. We also compared X. axonopodis pv. citrumelo to the genome of citrus canker pathogen X. axonopodis pv. citri 306. Comparative genomic analysis showed differences in several gene clusters, like those for type III effectors, the type IV secretion system, lipopolysaccharide synthesis, and others. In addition to pthA, effectors such as xopE3, xopAI, and hrpW were absent from X. axonopodis pv. citrumelo while present in X. axonopodis pv. citri. These effectors might be responsible for survival and the low virulence of this pathogen on citrus compared to that of X. axonopodis pv. citri. We also identified unique effectors in X. axonopodis pv. citrumelo that may be related to the different host range as compared to that of X. axonopodis pv. citri. X. axonopodis pv. citrumelo also lacks various genes, such as syrE1, syrE2, and RTX toxin family genes, which were present in X. axonopodis pv. citri. These may be associated with the distinct virulences of X. axonopodis pv. citrumelo and X. axonopodis pv. citri. Comparison of the complete genome sequence of X. axonopodis pv. citrumelo to those of X. axonopodis pv. citri and X. campestris pv. vesicatoria provides valuable insights into the mechanism of bacterial virulence and host specificity.

Regulatory Potential of Long Non-Coding RNAs (lncRNAs) in Boar Spermatozoa with Good and Poor Freezability.

Long non-coding RNAs (lncRNAs) are suggested to play an important role in the sperm biological processes. We performed de novo transcriptome assembly to characterize lncRNAs in spermatozoa, and to investigate the role of the potential target genes of the differentially expressed lncRNAs (DElncRNAs) in sperm freezability. We detected approximately 4007 DElncRNAs, which were differentially expressed in spermatozoa from boars classified as having good and poor semen freezability (GSF and PSF, respectively). Most of the DElncRNAs were upregulated in boars of the PSF group and appeared to significantly affect the sperm's response to the cryopreservation conditions. Furthermore, we predicted that the potential target genes were regulated by DElncRNAs in cis or trans. It was found that DElncRNAs of both freezability groups had potential cis- and trans-regulatory effects on different protein-coding genes, such as COX7A2L, TXNDC8 and SOX-7. Gene Ontology (GO) enrichment revealed that the DElncRNA target genes are associated with numerous biological processes, including signal transduction, response to stress, cell death (apoptosis), motility and embryo development. Significant differences in the de novo assembled transcriptome expression profiles of the DElncRNAs between the freezability groups were confirmed by quantitative real-time PCR analysis. This study reveals the potential effects of protein-coding genes of DElncRNAs on sperm functions, which could contribute to further research on their relevance in semen freezability.

Cytogenetic and molecular landscape and its potential clinical significance in Hispanic CMML patients from Puerto Rico.

Chronic myelomonocytic leukemia (CMML) is a clonal hematopoietic neoplasm that exhibits myelodysplastic and myeloproliferative characteristics with heterogeneous clinical and pathological features. There are limited publications on the ethnic and racial disparity of cytogenetics and genomics in CMML patients. This study aims to define the cytogenetic and molecular landscape in Hispanic CMML patients from Puerto Rico and explore its possible clinical significance. One hundred and eleven (111) Hispanic CMML patients from Puerto Rico were diagnosed in our institute from 2009 to 2018. Karyotypes were available in one hundred and seven (107) patients. Seventeen (17) patients had abnormal karyotypes (17/107, 16%). Compared to previously published data, Hispanic CMML patients in this study had significantly lower rates of overall cytogenetic abnormalities (16% vs 27-28%, p < 0.05) and trisomy 8 (2% vs 7%, p < 0.05). Among one hundred and eleven (111) Hispanic CMML patients, 40-gene myeloid molecular profile tests were performed in fifty-six (56) CMML patients. Gene mutations were identified in fifty-four (54) patients (96%). The most frequent mutated genes were: TET2, SRSF2, ASXL1, ZRSR2, DNMT3A, NRAS, CBL, and RUNX1. Twenty-nine (29) out of fifty-six (56) patients (29/56, 52%) had mutated TET2/wild type ASXL1 (muTET2/wtASXL1). Previous studies indicated that mutated ASXL1, DNMT3A, NRAS, RUNX1, and SETBP1 may associate with an unfavorable prognosis and muTET2/wtASXL1 may associate with a favorable prognosis in CMML patients. Compared to previously published data, Hispanic CMML patients from Puerto Rico in this study had significantly lower mutation rates in ASXL1 and SETBP1, and a higher rate of muTET2/wtASXL1. The findings raise the possibility of a favorable prognosis in Hispanic CMML patients.

Habitat-Specific Clock Variation and Its Consequence on Reproductive Fitness.

The circadian clock controls daily activities at the cellular and organismic level, allowing an organism to anticipate incoming stresses and to use resources accordingly. The circadian clock has therefore been considered a fitness trait in multiple organisms. However, the mechanism of how circadian clock variation influences organismal reproductive fitness is still not well understood. Here we describe habitat-specific clock variation (HSCV) of asexual reproduction in Neurospora discreta, a species that is adapted to 2 different habitats, under or above tree bark. African (AF) N. discreta strains, whose habitat is above the tree bark in light-dark (LD) conditions, display a higher rhythmicity index compared with North American (NA) strains, whose habitat is under the tree bark in constant dark (DD). Although AF-type strains demonstrated an overall fitness advantage under LD and DD conditions, NA-type strains exhibit a habitat-specific fitness advantage in DD over the LD condition. In addition, we show that allelic variation of the clock-controlled gene, Ubiquinol cytochrome c oxidoreductase (NEUDI_158280), plays a role in HSCV by modulating cellular reactive oxygen species levels. Our results demonstrate a mechanism by which local adaptation involving circadian clock regulation influences reproductive fitness.

Comparative Analysis of the Liver Transcriptome among Cattle Breeds Using RNA-seq.

Global gene expression in liver transcriptome varies among cattle breeds. The present investigation was aimed to identify the differentially expressed genes (DEGs), metabolic gene networks and metabolic pathways in bovine liver transcriptome of young bulls. In this study, we comparatively analyzed the bovine liver transcriptome of dairy (Polish Holstein Friesian (HF); n = 6), beef (Hereford; n = 6), and dual purpose (Polish-Red; n = 6) cattle breeds. This study identified 895, 338, and 571 significant (p < 0.01) differentially expressed (DE) gene-transcripts represented as 745, 265, and 498 hepatic DE genes through the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-HF versus Polish-Red breeds comparisons, respectively. By combining all breeds comparisons, 75 hepatic DE genes (p < 0.01) were identified as commonly shared among all the three breed comparisons; 70, 160, and 38 hepatic DE genes were commonly shared between the following comparisons: (i) Polish-Red versus Hereford and Polish-HF versus Hereford; (ii) Polish-Red versus Hereford and Polish-HF versus Polish-Red; and (iii) Polish-HF versus Hereford and Polish-HF versus Polish-Red, respectively. A total of 440, 82, and 225 hepatic DE genes were uniquely observed for the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-Red versus Polish-HF comparisons, respectively. Gene ontology (GO) analysis identified top-ranked enriched GO terms (p < 0.01) including 17, 16, and 31 functional groups and 151, 61, and 140 gene functions that were DE in all three breed liver transcriptome comparisons. Gene network analysis identified several potential metabolic pathways involved in glutamine family amino-acid, triglyceride synthesis, gluconeogenesis, p38MAPK cascade regulation, cholesterol biosynthesis (Polish-Red versus Hereford); IGF-receptor signaling, catecholamine transport, lipoprotein lipase, tyrosine kinase binding receptor (Polish-HF versus Hereford), and PGF-receptor binding, (Polish-HF versus Polish-Red). Validation results showed that the relative expression values were consistent to those obtained by RNA-seq, and significantly correlated between the quantitative reverse transcription PCR (RT-qPCR) and RNA-seq (Pearson's r > 0.90). Our results provide new insights on bovine liver gene expressions among dairy versus dual versus beef breeds by identifying the large numbers of DEGs markers submitted to NCBI gene expression omnibus (GEO) accession number GSE114233, which can serve as useful genetic tools to develop the gene assays for trait-associated studies as well as, to effectively implement in genomics selection (GS) cattle breeding programs in Poland.

Competitive Ability of Maize Pollen Grains Requires Paralogous Serine Threonine Protein Kinases STK1 and STK2.

serine threonine kinase1 (stk1) and serine threonine kinase2 (stk2) are closely related maize paralogous genes predicted to encode serine/threonine protein kinases. Pollen mutated in stk1 or stk2 competes poorly with normal pollen, pointing to a defect in pollen tube germination or growth. Both genes are expressed in pollen, but not in most other tissues. In germination media, STK1 and STK2 fluorescent fusion proteins localize to the plasma membrane of the vegetative cell. RNA-seq experiments identified 534 differentially expressed genes in stk1 mutant pollen relative to wild type. Gene ontology (GO) molecular functional analysis uncovered several differentially expressed genes with putative ribosome initiation and elongation functions, suggesting that stk1 might affect ribosome function. Of the two paralogs, stk1 may play a more important role in pollen development than stk2, as stk2 mutations have a smaller pollen transmission effect. However, stk2 does act as an enhancer of stk1 because the double mutant combination is only infrequently pollen-transmitted in double heterozygotes. We conclude that the stk paralogs play an essential role in pollen development.

Single nucleotide polymorphism discovery in bovine liver using RNA-seq technology.

BACKGROUND: RNA-seq is a useful next-generation sequencing (NGS) technology that has been widely used to understand mammalian transcriptome architecture and function. In this study, a breed-specific RNA-seq experiment was utilized to detect putative single nucleotide polymorphisms (SNPs) in liver tissue of young bulls of the Polish Red, Polish Holstein-Friesian (HF) and Hereford breeds, and to understand the genomic variation in the three cattle breeds that may reflect differences in production traits. RESULTS: The RNA-seq experiment on bovine liver produced 107,114,4072 raw paired-end reads, with an average of approximately 60 million paired-end reads per library. Breed-wise, a total of 345.06, 290.04 and 436.03 million paired-end reads were obtained from the Polish Red, Polish HF, and Hereford breeds, respectively. Burrows-Wheeler Aligner (BWA) read alignments showed that 81.35%, 82.81% and 84.21% of the mapped sequencing reads were properly paired to the Polish Red, Polish HF, and Hereford breeds, respectively. This study identified 5,641,401 SNPs and insertion and deletion (indel) positions expressed in the bovine liver with an average of 313,411 SNPs and indel per young bull. Following the removal of the indel mutations, a total of 195,3804, 152,7120 and 205,3184 raw SNPs expressed in bovine liver were identified for the Polish Red, Polish HF, and Hereford breeds, respectively. Breed-wise, three highly reliable breed-specific SNP-databases (SNP-dbs) with 31,562, 24,945 and 28,194 SNP records were constructed for the Polish Red, Polish HF, and Hereford breeds, respectively. Using a combination of stringent parameters of a minimum depth of ≥10 mapping reads that support the polymorphic nucleotide base and 100% SNP ratio, 4,368, 3,780 and 3,800 SNP records were detected in the Polish Red, Polish HF, and Hereford breeds, respectively. The SNP detections using RNA-seq data were successfully validated by kompetitive allele-specific PCR (KASPTM) SNP genotyping assay. The comprehensive QTL/CG analysis of 110 QTL/CG with RNA-seq data identified 20 monomorphic SNP hit loci (CARTPT, GAD1, GDF5, GHRH, GHRL, GRB10, IGFBPL1, IGFL1, LEP, LHX4, MC4R, MSTN, NKAIN1, PLAG1, POU1F1, SDR16C5, SH2B2, TOX, UCP3 and WNT10B) in all three cattle breeds. However, six SNP loci (CCSER1, GHR, KCNIP4, MTSS1, EGFR and NSMCE2) were identified as highly polymorphic among the cattle breeds. CONCLUSIONS: This study identified breed-specific SNPs with greater SNP ratio and excellent mapping coverage, as well as monomorphic and highly polymorphic putative SNP loci within QTL/CGs of bovine liver tissue. A breed-specific SNP-db constructed for bovine liver yielded nearly six million SNPs. In addition, a KASPTM SNP genotyping assay, as a reliable cost-effective method, successfully validated the breed-specific putative SNPs originating from the RNA-seq experiments.