Author Correction: Characterization of Nigerian breast cancer reveals prevalent homologous recombination deficiency and aggressive molecular features

The original version of this Article contained an error in the author affiliations. The affiliation of Kevin P. White with Tempus Labs, Inc. Chicago, IL, USA was inadvertently omitted. This has now been corrected in both the PDF and HTML versions of the Article.

Characterization of Nigerian breast cancer reveals prevalent homologous recombination deficiency and aggressive molecular features

Racial/ethnic disparities in breast cancer mortality continue to widen but genomic studies rarely interrogate breast cancer in diverse populations. Through genome, exome, and RNA sequencing, we examined the molecular features of breast cancers using 194 patients from Nigeria and 1037 patients from The Cancer Genome Atlas (TCGA). Relative to Black and White cohorts in TCGA, Nigerian HR + /HER2 - tumors are characterized by increased homologous recombination deficiency signature, pervasive TP53 mutations, and greater structural variation-indicating aggressive biology. GATA3 mutations are also more frequent in Nigerians regardless of subtype. Higher proportions of APOBEC-mediated substitutions strongly associate with PIK3CA and CDH1 mutations, which are underrepresented in Nigerians and Blacks. PLK2, KDM6A, and B2M are also identified as previously unreported significantly mutated genes in breast cancer. This dataset provides novel insights into potential molecular mechanisms underlying outcome disparities and lay a foundation for deployment of precision therapeutics in underserved populations.

Evolution of transcript modification by N6-methyladenosine in primates

Phenotypic differences within populations and between closely related species are often driven by variation and evolution of gene expression. However, most analyses have focused on the effects of genomic variation at cis-regulatory elements such as promoters and enhancers that control transcriptional activity, and little is understood about the influence of post-transcriptional processes on transcript evolution. Post-transcriptional modification of RNA by N⁶-methyladenosine (m⁶A) has been shown to be widespread throughout the transcriptome, and this reversible mark can affect transcript stability and translation dynamics. Here we analyze m⁶A mRNA modifications in lymphoblastoid cell lines (LCLs) from human, chimpanzee and rhesus, and we identify patterns of m⁶A evolution among species. We find that m⁶A evolution occurs in parallel with evolution of consensus RNA sequence motifs known to be associated with the enzymatic complexes that regulate m⁶A dynamics, and expression evolution of m⁶A-modified genes occurs in parallel with m⁶A evolution.

Clinical Validation of a Next-Generation Sequencing Genomic Oncology Panel via Cross-Platform Benchmarking against Established Amplicon Sequencing Assays

Next-generation sequencing (NGS) genomic oncology profiling assays have emerged as key drivers of personalized cancer care and translational research. However, validation of these assays to meet strict clinical standards has been historically problematic because of both significant assay complexity and a scarcity of optimal validation samples. Herein, we present the clinical validation of 76 genes from a novel 1212-gene large-scale hybrid capture cancer sequencing assay (University of Chicago Medicine OncoPlus) using full-data comparisons against multiple clinical NGS amplicon-based assays to yield dramatic increases in per-sample data comparison efficiency compared with previously published validations. Using a sample set of 104 normal, solid tumor, and hematopoietic malignancy specimens, head-to-head NGS data analyses allowed for 6.8 million individual clinical base call comparisons, including 2729 previously confirmed variants, with 100% sensitivity and specificity. University of Chicago Medicine OncoPlus showed excellent performance for detection of single-nucleotide variants, insertions/deletions up to 52 bp, and FLT3 internal tandem duplications of up to 102 bp or larger. Highly concordant copy number variant and ALK/RET/ROS1 gene fusion detection were also observed. In addition to underlining the efficiency of NGS validation via full-data benchmarking against existing clinical NGS assays, this study also highlights the degree of performance similarity between hybrid capture and amplicon assays that is attainable with the application of strict quality control parameters and optimized computational analytics.

Divergent patterns of endogenous small RNA populations from seed and vegetative tissues of Glycine max

BACKGROUND

Small non-coding RNAs (smRNAs) are known to have major roles in gene regulation in eukaryotes. In plants, knowledge of the biogenesis and mechanisms of action of smRNA classes including microRNAs (miRNAs), short interfering RNAs (siRNAs), and trans-acting siRNAs (tasiRNAs) has been gained mostly through studies with Arabidopsis. In recent years, high throughput sequencing of smRNA populations has enabled extension of knowledge from model systems to plants with larger, more complex genomes. Soybean (Glycine max) now has many genomics resources available including a complete genome sequence and predicted gene models. Relatively little is known, however, about the full complement of its endogenous smRNAs populations and the silenced genes.

RESULTS

Using Illumina sequencing and computational analysis, we characterized eight smRNA populations from multiple tissues and organs of soybean including developing seed and vegetative tissues. A total of 41 million raw sequence reads collapsed into 135,055 unique reads were mapped to the soybean genome and its predicted cDNA gene models. Bioinformatic analyses were used to distinguish miRNAs and siRNAs and to determine their genomic origins and potential target genes. In addition, we identified two soybean TAS3 gene homologs, the miRNAs that putatively guide cleavage of their transcripts, and the derived tasiRNAs that could target soybean genes annotated as auxin response factors. Tissue-differential expression based on the flux of normalized miRNA and siRNA abundances in the eight smRNA libraries was evident, some of which was confirmed by smRNA blotting. Our global view of these smRNA populations also revealed that the size classes of smRNAs varied amongst different tissues, with the developing seed and seed coat having greater numbers of unique smRNAs of the 24-nt class compared to the vegetative tissues of germinating seedlings. The 24-nt class is known to be derived from repetitive elements including transposons. Detailed analysis of the size classes associated with ribosomal RNAs and transposable element families showed greater diversity of smRNAs in the 22- and 24-nt size classes.

CONCLUSIONS

The flux of endogenous smRNAs within multiple stages and tissues of seed development was contrasted with vegetative tissues of soybean, one of the dominant sources of protein and oil in world markets. The smRNAs varied in size class, complexity of origins, and possible targets. Sequencing revealed tissue-preferential expression for certain smRNAs and expression differences among closely related miRNA family members.

Endogenous, tissue-specific short interfering RNAs silence the chalcone synthase gene family in glycine max seed coats

Two dominant alleles of the I locus in Glycine max silence nine chalcone synthase (CHS) genes to inhibit function of the flavonoid pathway in the seed coat. We describe here the intricacies of this naturally occurring silencing mechanism based on results from small RNA gel blots and high-throughput sequencing of small RNA populations. The two dominant alleles of the I locus encompass a 27-kb region containing two perfectly repeated and inverted clusters of three chalcone synthase genes (CHS1, CHS3, and CHS4). This structure silences the expression of all CHS genes, including CHS7 and CHS8, located on other chromosomes. The CHS short interfering RNAs (siRNAs) sequenced support a mechanism by which RNAs transcribed from the CHS inverted repeat form aberrant double-stranded RNAs that become substrates for dicer-like ribonuclease. The resulting primary siRNAs become guides that target the mRNAs of the nonlinked, highly expressed CHS7 and CHS8 genes, followed by subsequent amplification of CHS7 and CHS8 secondary siRNAs by RNA-dependent RNA polymerase. Most remarkably, this silencing mechanism occurs only in one tissue, the seed coat, as shown by the lack of CHS siRNAs in cotyledons and vegetative tissues. Thus, production of the trigger double-stranded RNA that initiates the process occurs in a specific tissue and represents an example of naturally occurring inhibition of a metabolic pathway by siRNAs in one tissue while allowing expression of the pathway and synthesis of valuable secondary metabolites in all other organs/tissues of the plant.

Features of a 103-kb gene-rich region in soybean include an inverted perfect repeat cluster of CHS genes comprising the I locus

The I locus in soybean (Glycine max) corresponds to a region of chalcone synthase (CHS) gene duplications affecting seed pigmentation. We sequenced and annotated BAC clone 104J7, which harbors a dominant i(i) allele from Glycine max 'Williams 82', to gain insight into the genetic structure of this multigenic region in addition to examining its flanking regions. The 103-kb BAC encompasses a gene-rich region with 11 putatively expressed genes. In addition to six copies of CHS, these genes include: a geranylgeranyltransferase type II beta subunit (E.C.2.5.1.60), a beta-galactosidase, a putative spermine and (or) spermidine synthase (E.C.2.5.1.16), and an unknown expressed gene. Strikingly, sequencing data revealed that the 10.91-kb CHS1, CHS3, CHS4 cluster is present as a perfect inverted repeat separated by 5.87 kb. Contiguous arrangement of CHS paralogs could lead to folding into multiple secondary structures, hypothesized to induce deletions that have previously been shown to effect CHS expression. BAC104J7 also contains several gene fragments representing a cation/hydrogen exchanger, a 40S ribosomal protein, a CBL-interacting protein kinase, and the amino terminus of a subtilisin. Chimeric ESTs were identified that may represent read-through transcription from a flanking truncated gene into a CHS cluster, generating aberrant CHS RNA molecules that could play a role in CHS gene silencing.

Tissue-specific gene silencing mediated by a naturally occurring chalcone synthase gene cluster in Glycine max

Chalcone synthase, a key regulatory enzyme in the flavonoid pathway, constitutes an eight-member gene family in Glycine max (soybean). Three of the chalcone synthase (CHS) gene family members are arranged as inverted repeats in a 10-kb region, corresponding to the I locus (inhibitor). Spontaneous mutations of a dominant allele (I or i(i)) to a recessive allele (i) have been shown to delete promoter sequences, paradoxically increasing total CHS transcript levels and resulting in black seed coats. However, it is not known which of the gene family members contribute toward pigmentation and how this locus affects CHS expression in other tissues. We investigated the unusual nature of the I locus using four pairs of isogenic lines differing with respect to alleles of the I locus. RNA gel blots using a generic open reading frame CHS probe detected similar CHS transcript levels in stems, roots, leaves, young pods, and cotyledons of the yellow and black isolines but not in the seed coats, which is consistent with the dominant I and i(i) alleles mediating CHS gene silencing in a tissue-specific manner. Using real-time RT-PCR, a variable pattern of expression of CHS genes in different tissues was demonstrated. However, increase in pigmentation in the black seed coats was associated with release of the silencing effect specifically on CHS7/CHS8, which occurred at all stages of seed coat development. These expression changes were linked to structural changes taking place at the I locus, shown to encompass a much wider region of at least 27 kb, comprising two identical 10.91-kb stretches of CHS gene duplications. The suppressive effect of this 27-kb I locus in a specific tissue of the G. max plant represents a unique endogenous gene silencing mechanism.