A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Common variants in 94 loci have been associated with breast cancer including 15 loci with genome-wide significant associations (P<5 × 10(-8)) with oestrogen receptor (ER)-negative breast cancer and BRCA1-associated breast cancer risk. In this study, to identify new ER-negative susceptibility loci, we performed a meta-analysis of 11 genome-wide association studies (GWAS) consisting of 4,939 ER-negative cases and 14,352 controls, combined with 7,333 ER-negative cases and 42,468 controls and 15,252 BRCA1 mutation carriers genotyped on the iCOGS array. We identify four previously unidentified loci including two loci at 13q22 near KLF5, a 2p23.2 locus near WDR43 and a 2q33 locus near PPIL3 that display genome-wide significant associations with ER-negative breast cancer. In addition, 19 known breast cancer risk loci have genome-wide significant associations and 40 had moderate associations (P<0.05) with ER-negative disease. Using functional and eQTL studies we implicate TRMT61B and WDR43 at 2p23.2 and PPIL3 at 2q33 in ER-negative breast cancer aetiology. All ER-negative loci combined account for ∼11% of familial relative risk for ER-negative disease and may contribute to improved ER-negative and BRCA1 breast cancer risk prediction.

Genome-wide association studies have reported 11 regions conferring risk of high-grade serous epithelial ovarian cancer (HGSOC). Expression quantitative trait locus (eQTL) analyses can identify candidate susceptibility genes at risk loci. Here we evaluate cis-eQTL associations at 47 regions associated with HGSOC risk (P≤10(-5)). For three cis-eQTL associations (P<1.4 × 10(-3), FDR<0.05) at 1p36 (CDC42), 1p34 (CDCA8) and 2q31 (HOXD9), we evaluate the functional role of each candidate by perturbing expression of each gene in HGSOC precursor cells. Overexpression of HOXD9 increases anchorage-independent growth, shortens population-doubling time and reduces contact inhibition. Chromosome conformation capture identifies an interaction between rs2857532 and the HOXD9 promoter, suggesting this SNP is a leading causal variant. Transcriptomic profiling after HOXD9 overexpression reveals enrichment of HGSOC risk variants within HOXD9 target genes (P=6 × 10(-10) for risk variants (P<10(-4)) within 10 kb of a HOXD9 target gene in ovarian cells), suggesting a broader role for this network in genetic susceptibility to HGSOC.

Retinoic acid inducible gene-I (RIG-I) ensures immune surveillance of viral RNAs bearing a 5'-triphosphate (5'ppp) moiety. Mutations in RIG-I (C268F and E373A) lead to impaired ATPase activity, thereby driving hyperactive signaling associated with autoimmune diseases. Here we report, using hydrogen/deuterium exchange, mechanistic models for dysregulated RIG-I proofreading that ultimately result in the improper recognition of cellular RNAs bearing 7-methylguanosine and N1-2'-O-methylation (Cap1) on the 5' end. Cap1-RNA compromises its ability to stabilize RIG-I helicase and blunts caspase activation and recruitment domains (CARD) partial opening by threefold. RIG-I H830A mutation restores Cap1-helicase engagement as well as CARDs partial opening event to a level comparable to that of 5'ppp. However, E373A RIG-I locks the receptor in an ATP-bound state, resulting in enhanced Cap1-helicase engagement and a sequential CARDs stimulation. C268F mutation renders a more tethered ring architecture and results in constitutive CARDs signaling in an ATP-independent manner.

The analysis of whole-genome sequencing studies is challenging due to the large number of rare variants in noncoding regions and the lack of natural units for testing. We propose a statistical method to detect and localize rare and common risk variants in whole-genome sequencing studies based on a recently developed knockoff framework. It can (1) prioritize causal variants over associations due to linkage disequilibrium thereby improving interpretability; (2) help distinguish the signal due to rare variants from shadow effects of significant common variants nearby; (3) integrate multiple knockoffs for improved power, stability, and reproducibility; and (4) flexibly incorporate state-of-the-art and future association tests to achieve the benefits proposed here. In applications to whole-genome sequencing data from the Alzheimer's Disease Sequencing Project (ADSP) and COPDGene samples from NHLBI Trans-Omics for Precision Medicine (TOPMed) Program we show that our method compared with conventional association tests can lead to substantially more discoveries.

HNF1B is overexpressed in clear cell epithelial ovarian cancer, and we observed epigenetic silencing in serous epithelial ovarian cancer, leading us to hypothesize that variation in this gene differentially associates with epithelial ovarian cancer risk according to histological subtype. Here we comprehensively map variation in HNF1B with respect to epithelial ovarian cancer risk and analyse DNA methylation and expression profiles across histological subtypes. Different single-nucleotide polymorphisms associate with invasive serous (rs7405776 odds ratio (OR)=1.13, P=3.1 × 10(-10)) and clear cell (rs11651755 OR=0.77, P=1.6 × 10(-8)) epithelial ovarian cancer. Risk alleles for the serous subtype associate with higher HNF1B-promoter methylation in these tumours. Unmethylated, expressed HNF1B, primarily present in clear cell tumours, coincides with a CpG island methylator phenotype affecting numerous other promoters throughout the genome. Different variants in HNF1B associate with risk of serous and clear cell epithelial ovarian cancer; DNA methylation and expression patterns are also notably distinct between these subtypes. These findings underscore distinct mechanisms driving different epithelial ovarian cancer histological subtypes.

Understanding the molecular mechanisms regulating the maintenance and destruction of intervertebral disc may lead to the development of new therapies for intervertebral disc degeneration (IDD). Here we present evidence from miRNA microarray analyses of clinical data sets along with in vitro and in vivo experiments that miR-141 is a key regulator of IDD. Gain- and loss-of-function studies show that miR-141 drives IDD by inducing nucleus pulposus (NP) apoptosis. Furthermore, miR-141 KO in mice attenuated spontaneous and surgically induced IDD. Mechanistically, miR-141 promotes IDD development by targeting and depleting SIRT1, a negative regulator of NF-κB pathway. Therapeutically, upregulation or downregulation of miR-141 by nanoparticle delivery in IDD model aggravated or alleviated experimental IDD, respectively. Our findings reveal a novel mechanism by which miR-141, in part, promotes IDD progression by interacting with SIRT1/NF-κB pathway. Blockade of miR-141 in vivo may serve as a potential therapeutic approach in the treatment of IDD.

Mucinous ovarian carcinoma (MOC) is a unique subtype of ovarian cancer with an uncertain etiology, including whether it genuinely arises at the ovary or is metastatic disease from other organs. In addition, the molecular drivers of invasive progression, high-grade and metastatic disease are poorly defined. We perform genetic analysis of MOC across all histological grades, including benign and borderline mucinous ovarian tumors, and compare these to tumors from other potential extra-ovarian sites of origin. Here we show that MOC is distinct from tumors from other sites and supports a progressive model of evolution from borderline precursors to high-grade invasive MOC. Key drivers of progression identified are TP53 mutation and copy number aberrations, including a notable amplicon on 9p13. High copy number aberration burden is associated with worse prognosis in MOC. Our data conclusively demonstrate that MOC arise from benign and borderline precursors at the ovary and are not extra-ovarian metastases.

Plant viruses adopt diverse virulence strategies to inhibit host antiviral defense. However, general antiviral defense directly targeted by different types of plant viruses have rarely been studied. Here, we show that the single rice DELLA protein, SLENDER RICE 1 (SLR1), a master negative regulator in Gibberellin (GA) signaling pathway, is targeted by several different viral effectors for facilitating viral infection. Viral proteins encoded by different types of rice viruses all directly trigger the rapid degradation of SLR1 by promoting association with the GA receptor OsGID1. SLR1-mediated broad-spectrum resistance was subverted by these independently evolved viral proteins, which all interrupted the functional crosstalk between SLR1 and jasmonic acid (JA) signaling. This decline of JA antiviral further created the advantage of viral infection. Our study reveals a common viral counter-defense strategy in which different types of viruses convergently target SLR1-mediated broad-spectrum resistance to benefit viral infection in the monocotyledonous crop rice.

Systemic lupus erythematosus is a heritable autoimmune disease that predominantly affects young women. To improve our understanding of genetic etiology, we conduct multi-ancestry and multi-trait meta-analysis of genome-wide association studies, encompassing 12 systemic lupus erythematosus cohorts from 3 different ancestries and 10 genetically correlated autoimmune diseases, and identify 16 novel loci. We also perform transcriptome-wide association studies, computational drug repurposing analysis, and cell type enrichment analysis. We discover putative drug classes, including a histone deacetylase inhibitor that could be repurposed to treat lupus. We also identify multiple cell types enriched with putative target genes, such as non-classical monocytes and B cells, which may be targeted for future therapeutics. Using this newly assembled result, we further construct polygenic risk score models and demonstrate that integrating polygenic risk score with clinical lab biomarkers improves the diagnostic accuracy of systemic lupus erythematosus using the Vanderbilt BioVU and Michigan Genomics Initiative biobanks.

Reduced BRCA1 expression causes homologous recombination (HR) repair defects in high-grade serous ovarian cancers (HGSOCs). Here, we demonstrate that BRCA1 is transcriptionally activated by a previously unknown function of ZC3H18. We show that ZC3H18 is a DNA-binding protein that interacts with an E2F site in the BRCA1 promoter where it facilitates recruitment of E2F4 to an adjacent E2F site to promote BRCA1 transcription. Consistent with ZC3H18 role in activating BRCA1 expression, ZC3H18 depletion induces BRCA1 promoter methylation, reduces BRCA1 expression, disrupts HR, and sensitizes cells to DNA crosslinkers and poly(ADP-ribose) polymerase inhibitors. Moreover, in patient-derived xenografts and primary HGSOC tumors, ZC3H18 and E2F4 mRNA levels are positively correlated with BRCA1 mRNA levels, further supporting ZC3H18 role in regulating BRCA1. Given that ZC3H18 lies within 16q24.2, a region with frequent copy number loss in HGSOC, these findings suggest that ZC3H18 copy number losses could contribute to HR defects in HGSOC.

The synovium is an important component of any synovial joint and is the major target tissue of inflammatory arthritis. However, the multi-omics landscape of synovium required for functional inference is absent from large-scale resources. Here we integrate genomics with transcriptomics and chromatin accessibility features of human synovium in up to 245 arthritic patients, to characterize the landscape of genetic regulation on gene expression and the regulatory mechanisms mediating arthritic diseases predisposition. We identify 4765 independent primary and 616 secondary cis-expression quantitative trait loci (cis-eQTLs) in the synovium and find that the eQTLs with multiple independent signals have stronger effects and heritability than single independent eQTLs. Integration of genome-wide association studies (GWASs) and eQTLs identifies 84 arthritis related genes, revealing 38 novel genes which have not been reported by previous studies using eQTL data from the GTEx project or immune cells. We further develop a method called eQTac to identify variants that could affect gene expression by affecting chromatin accessibility and identify 1517 regions with potential regulatory function of chromatin accessibility. Altogether, our study provides a comprehensive synovium multi-omics resource for arthritic diseases and gains new insights into the regulation of gene expression.

CRISPR advances genome engineering by directing endonuclease sequence specificity with a guide RNA molecule (gRNA). For precisely targeting a gene for modification, each genetic construct requires a unique gRNA. By generating a gRNA against the flippase recognition target (FRT) site, a common genetic element shared by multiple genetic collections, CRISPR-FRT circumvents this design constraint to provide a broad platform for fast, scarless, off-the-shelf genome engineering.

GWAS have identified a breast cancer susceptibility locus on 2q35. Here we report the fine mapping of this locus using data from 101,943 subjects from 50 case-control studies. We genotype 276 SNPs using the 'iCOGS' genotyping array and impute genotypes for a further 1,284 using 1000 Genomes Project data. All but two, strongly correlated SNPs (rs4442975 G/T and rs6721996 G/A) are excluded as candidate causal variants at odds against >100:1. The best functional candidate, rs4442975, is associated with oestrogen receptor positive (ER+) disease with an odds ratio (OR) in Europeans of 0.85 (95% confidence interval=0.84-0.87; P=1.7 × 10(-43)) per t-allele. This SNP flanks a transcriptional enhancer that physically interacts with the promoter of IGFBP5 (encoding insulin-like growth factor-binding protein 5) and displays allele-specific gene expression, FOXA1 binding and chromatin looping. Evidence suggests that the g-allele confers increased breast cancer susceptibility through relative downregulation of IGFBP5, a gene with known roles in breast cell biology.

Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (rg = 0.57, p = 4.6 × 10-8), breast and ovarian cancer (rg = 0.24, p = 7 × 10-5), breast and lung cancer (rg = 0.18, p =1.5 × 10-6) and breast and colorectal cancer (rg = 0.15, p = 1.1 × 10-4). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.

Due to a lack of Mesozoic fossil records, the origins and early evolution of feather-feeding behaviors by insects are obscure. Here, we report ten nymph specimens of a new lineage of insect, Mesophthirus engeli gen et. sp. nov. within Mesophthiridae fam. nov. from the mid-Cretaceous (ca. 100 Mya) Myanmar (Burmese) amber. This new insect clade shows a series of ectoparasitic morphological characters such as tiny wingless body, head with strong chewing mouthparts, robust and short antennae having long setae, legs with only one single tarsal claw associated with two additional long setae, etc. Most significantly, these insects are preserved with partially damaged dinosaur feathers, the damage of which was probably made by these insects' integument-feeding behaviors. This finding demonstrates that feather-feeding behaviors of insects originated at least in mid-Cretaceous, accompanying the radiation of feathered dinosaurs including early birds.

African Americans have a significantly higher risk of developing chronic kidney disease, especially focal segmental glomerulosclerosis -, than European Americans. Two coding variants (G1 and G2) in the APOL1 gene play a major role in this disparity. While 13% of African Americans carry the high-risk recessive genotypes, only a fraction of these individuals develops FSGS or kidney failure, indicating the involvement of additional disease modifiers. Here, we show that the presence of the APOL1 p.N264K missense variant, when co-inherited with the G2 APOL1 risk allele, substantially reduces the penetrance of the G1G2 and G2G2 high-risk genotypes by rendering these genotypes low-risk. These results align with prior functional evidence showing that the p.N264K variant reduces the toxicity of the APOL1 high-risk alleles. These findings have important implications for our understanding of the mechanisms of APOL1-associated nephropathy, as well as for the clinical management of individuals with high-risk genotypes that include the G2 allele.

Protocell fitness under extreme prebiotic conditions is critical in understanding the origin of life. However, little is known about protocell's survival and fitness under prebiotic radiations. Here we present a radioresistant protocell model based on assembly of two types of coacervate droplets, which are formed through interactions of inorganic polyphosphate (polyP) with divalent metal cation and cationic tripeptide, respectively. Among the coacervate droplets, only the polyP-Mn droplet is radiotolerant and provides strong protection for recruited proteins. The radiosensitive polyP-tripeptide droplet sequestered with both proteins and DNA could be encapsulated inside the polyP-Mn droplet, and form into a compartmentalized protocell. The protocell protects the inner nucleoid-like condensate through efficient reactive oxygen species' scavenging capacity of intracellular nonenzymic antioxidants including Mn-phosphate and Mn-peptide. Our results demonstrate a radioresistant protocell model with redox reaction system in response to ionizing radiation, which might enable the protocell fitness to prebiotic radiation on the primitive Earth preceding the emergence of enzyme-based fitness. This protocell might also provide applications in synthetic biology as bioreactor or drug delivery system.

The adaptive immune response is under circadian control, yet, why adaptive immune reactions continue to exhibit circadian changes over long periods of time is unknown. Using a combination of experimental and mathematical modeling approaches, we show here that dendritic cells migrate from the skin to the draining lymph node in a time-of-day-dependent manner, which provides an enhanced likelihood for functional interactions with T cells. Rhythmic expression of TNF in the draining lymph node enhances BMAL1-controlled ICAM-1 expression in high endothelial venules, resulting in lymphocyte infiltration and lymph node expansion. Lymph node cellularity continues to be different for weeks after the initial time-of-day-dependent challenge, which governs the immune response to vaccinations directed against Hepatitis A virus as well as SARS-CoV-2. In this work, we present a mechanistic understanding of the time-of-day dependent development and maintenance of an adaptive immune response, providing a strategy for using time-of-day to optimize vaccination regimes.

Acute kidney injury (AKI) exhibits high morbidity and mortality. Kidney injury molecule-1 (KIM1) is dramatically upregulated in renal tubules upon injury, and acts as a biomarker for various renal diseases. However, the exact role and underlying mechanism of KIM1 in the progression of AKI remain elusive. Herein, we report that renal tubular specific knockout of Kim1 attenuates cisplatin- or ischemia/reperfusion-induced AKI in male mice. Mechanistically, transcription factor Yin Yang 1 (YY1), which is downregulated upon AKI, binds to the promoter of KIM1 and represses its expression. Injury-induced KIM1 binds to the ECD domain of death receptor 5 (DR5), which activates DR5 and the following caspase cascade by promoting its multimerization, thus induces renal cell apoptosis and exacerbates AKI. Blocking the KIM1-DR5 interaction with rationally designed peptides exhibit reno-protective effects against AKI. Here, we reveal a YY1-KIM1-DR5 axis in the progression of AKI, which warrants future exploration as therapeutic targets.