HPV16 E7 Genetic Conservation Is Critical to Carcinogenesis.

Although most cervical human papillomavirus type 16 (HPV16) infections become undetectable within 1-2 years, persistent HPV16 causes half of all cervical cancers. We used a novel HPV whole-genome sequencing technique to evaluate an exceptionally large collection of 5,570 HPV16-infected case-control samples to determine whether viral genetic variation influences risk of cervical precancer and cancer. We observed thousands of unique HPV16 genomes; very few women shared the identical HPV16 sequence, which should stimulate a careful re-evaluation of the clinical implications of HPV mutation rates, transmission, clearance, and persistence. In case-control analyses, HPV16 in the controls had significantly more amino acid changing variants throughout the genome. Strikingly, E7 was devoid of variants in precancers/cancers compared to higher levels in the controls; we confirmed this in cancers from around the world. Strict conservation of the 98 amino acids of E7, which disrupts Rb function, is critical for HPV16 carcinogenesis, presenting a highly specific target for etiologic and therapeutic research.

The combined finding of HPV 16, 18, or 45 and cytologic Atypical Glandular Cells (AGC) indicates a greatly elevated risk of in situ and invasive cervical adenocarcinoma.

BACKGROUND: Cervical screening has not effectively controlled cervical adenocarcinoma (AC). Human papillomavirus (HPV) testing is recommended for cervical screening but the optimal management of HPV-positive individuals to prevent AC remains a question. Cytology and HPV typing are two triage options to predict the risk of AC. We combined two potential biomarkers (atypical glandular cell, AGC, cytology and HPV-types 16, 18, or 45) to assess their joint effect on detecting AC. METHODS: Kaiser Permanente Northern California (KPNC) used triennial co-testing with cytology and HPV testing (positive/negative) for routine cervical screening between 2003 and 2020. HPV typing of a sample of residual HPV test specimens was performed on a separate cohort selected from KPNC (Persistence and Progression, PaP, cohort). We compared risk of prevalent and incident histologic AC/AIS (adenocarcinoma in situ) associated with preceding combinations of cytologic results and HPV typing. Risk of squamous cell cancer (SCC)/cervical intraepithelial neoplasia grade 3 (CIN3) (SCC/CIN3) was also included for comparison. RESULTS: Among HPV-positive individuals in PaP cohort, 99% of prevalent AC and 96% of AIS were linked to HPV-types 16, 18, or 45 (denoted HPV 16/18/45). Although rare (0.09% of screening population), the concurrent detection of HPV 16/18/45 with AGC cytology predicted a highly elevated relative risk of underlying histologic AC/AIS; the absolute risk of diagnosing AC/AIS was 12% and odds ratio (OR) was 1341 (95%CI:495-3630) compared to patients with other high-risk HPV types and normal cytology. Cumulatively (allowing non-concurrent results), approximately one-third of the AC/AIS cases ever had HPV 16/18/45 and AGC cytology (OR = 1785; 95%CI:872-3656). AGC was not as strongly associated with SCC/CIN3. CONCLUSION: Detection of HPV 16/18/45 positivity elevates risk of adenocarcinoma, particularly if AGC cytology is also found.

Association of HPV35 with cervical carcinogenesis among women of African ancestry: Evidence of viral-host interaction with implications for disease intervention.

HPV35 has been found in only ∼2% of invasive cervical cancers (ICC) worldwide but up to 10% in Sub-Saharan Africa, warranting further investigation and consideration of impact on preventive strategies. We studied HPV35 and ethnicity, in relation to the known steps in cervical carcinogenesis, using multiple large epidemiologic studies in the U.S. and internationally. Combining five U.S. studies, we measured HPV35 positivity and, in Northern California, observed HPV35 type-specific population prevalence and estimated 5-year risk of developing precancer when HPV35-positive. HPV35 genetic variation was examined for differences in carcinogenicity in 1053 HPV35+ cervical specimens from a U.S. cohort and an international collection. African-American women had more HPV35 (12.1% vs 5.1%, P < .001) and more HPV35-associated precancers (7.4% vs 2.1%, P < .001) compared to other ethnicities. Precancer risks after HPV35 infection did not vary by ethnicity (global P = .52). The HPV35 A2 sublineage showed an increased association with precancer/cancer in African-Americans (OR = 5.6 vs A1, 95% CI = 1.3-24.8) and A2 was more prevalent among ICC in Africa than other world regions (41.9% vs 10.4%, P < .01). Our analyses support a strong link between HPV35 and cervical carcinogenesis in women of African ancestry. Current HPV vaccines cover the majority of cervical precancer/cancer across all ethnic groups; additional analyses are required to determine whether the addition of HPV35 to the already highly effective nine-valent HPV vaccine would provide better protection for women in Africa or of African ancestry.

Rare germline variants in known melanoma susceptibility genes in familial melanoma.

Known high-risk cutaneous malignant melanoma (CMM) genes account for melanoma risk in <40% of melanoma-prone families, suggesting the existence of additional high-risk genes or perhaps a polygenic mechanism involving multiple genetic modifiers. The goal of this study was to systematically characterize rare germline variants in 42 established melanoma genes among 144 CMM patients in 76 American CMM families without known mutations using data from whole-exome sequencing. We identified 68 rare (<0.1% in public and in-house control datasets) nonsynonymous variants in 25 genes. We technically validated all loss-of-function, inframe insertion/deletion, and missense variants predicted as deleterious, and followed them up in 1, 559 population-based CMM cases and 1, 633 controls. Several of these variants showed disease co-segregation within families. Of particular interest, a stopgain variant in TYR was present in five of six CMM cases/obligate gene carriers in one family and a single population-based CMM case. A start gain variant in the 5'UTR region of PLA2G6 and a missense variant in ATM were each seen in all three affected people in a single family, respectively. Results from rare variant burden tests showed that familial and population-based CMM patients tended to have higher frequencies of rare germline variants in albinism genes such as TYR, TYRP1, and OCA2 (P < 0.05). Our results suggest that rare nonsynonymous variants in low- or intermediate-risk CMM genes may influence familial CMM predisposition, warranting further investigation of both common and rare variants in genes affecting functionally important pathways (such as melanogenesis) in melanoma risk assessment.

Evidence for SNP-SNP interaction identified through targeted sequencing of cleft case-parent trios.

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is the most common craniofacial birth defect in humans, affecting 1 in 700 live births. This malformation has a complex etiology where multiple genes and several environmental factors influence risk. At least a dozen different genes have been confirmed to be associated with risk of NSCL/P in previous studies. However, all the known genetic risk factors cannot fully explain the observed heritability of NSCL/P, and several authors have suggested gene-gene (G × G) interaction may be important in the etiology of this complex and heterogeneous malformation. We tested for G × G interactions using common single nucleotide polymorphic (SNP) markers from targeted sequencing in 13 regions identified by previous studies spanning 6.3 Mb of the genome in a study of 1,498 NSCL/P case-parent trios. We used the R-package trio to assess interactions between polymorphic markers in different genes, using a 1 degree of freedom (1df) test for screening, and a 4 degree of freedom (4df) test to assess statistical significance of epistatic interactions. To adjust for multiple comparisons, we performed permutation tests. The most significant interaction was observed between rs6029315 in MAFB and rs6681355 in IRF6 (4df P = 3.8 × 10-8 ) in case-parent trios of European ancestry, which remained significant after correcting for multiple comparisons. However, no significant interaction was detected in trios of Asian ancestry.

Prevalence and spectrum of germline rare variants in BRCA1/2 and PALB2 among breast cancer cases in Sarawak, Malaysia.

PURPOSE: To characterize the spectrum of germline mutations in BRCA1, BRCA2, and PALB2 in population-based unselected breast cancer cases in an Asian population. METHODS: Germline DNA from 467 breast cancer patients in Sarawak General Hospital, Malaysia, where 93% of the breast cancer patients in Sarawak are treated, was sequenced for the entire coding region of BRCA1; BRCA2; PALB2; Exons 6, 7, and 8 of TP53; and Exons 7 and 8 of PTEN. Pathogenic variants included known pathogenic variants in ClinVar, loss of function variants, and variants that disrupt splice site. RESULTS: We found 27 pathogenic variants (11 BRCA1, 10 BRCA2, 4 PALB2, and 2 TP53) in 34 patients, which gave a prevalence of germline mutations of 2.8, 3.23, and 0.86% for BRCA1, BRCA2, and PALB2, respectively. Compared to mutation non-carriers, BRCA1 mutation carriers were more likely to have an earlier age at onset, triple-negative subtype, and lower body mass index, whereas BRCA2 mutation carriers were more likely to have a positive family history. Mutation carrier cases had worse survival compared to non-carriers; however, the association was mostly driven by stage and tumor subtype. We also identified 19 variants of unknown significance, and some of them were predicted to alter splicing or transcription factor binding sites. CONCLUSION: Our data provide insight into the genetics of breast cancer in this understudied group and suggest the need for modifying genetic testing guidelines for this population with a much younger age at diagnosis and more limited resources compared with Caucasian populations.

Multiple rare variants in high-risk pancreatic cancer-related genes may increase risk for pancreatic cancer in a subset of patients with and without germline CDKN2A mutations.

The risk of pancreatic cancer (PC) is increased in melanoma-prone families but the causal relationship between germline CDKN2A mutations and PC risk is uncertain, suggesting the existence of non-CDKN2A factors. One genetic possibility involves patients having mutations in multiple high-risk PC-related genes; however, no systematic examination has yet been conducted. We used next-generation sequencing data to examine 24 putative PC-related genes in 43 PC patients with and 23 PC patients without germline CDKN2A mutations and 1001 controls. For each gene and the four pathways in which they occurred, we tested whether PC patients (overall or CDKN2A+ and CDKN2A- cases separately) had an increased number of rare nonsynonymous variants. Overall, we identified 35 missense variants in PC patients, 14 in CDKN2A+ and 21 in CDKN2A- PC cases. We found nominally significant associations for mismatch repair genes (MLH1, MSH2, MSH6, PMS2) in all PC patients and for ATM, CPA1, and PMS2 in CDKN2A- PC patients. Further, nine CDKN2A+ and four CDKN2A- PC patients had rare potentially deleterious variants in multiple PC-related genes. Loss-of-function variants were only observed in CDKN2A- PC patients, with ATM having the most pathogenic variants. Also, ATM variants (n = 5) were only observed in CDKN2A- PC patients with a family history that included digestive system tumors. Our results suggest that a subset of PC patients may have increased risk because of germline mutations in multiple PC-related genes.

Physicochemical characteristics and microbial community succession during oat silage prepared without or with Lactiplantibacillus plantarum or Lentilactobacillus buchneri.

IMPORTANCE: Ensiled whole-plant oats are an important feedstuff for ruminants in large parts of the world. Oat silage is rich in dietary fibers, minerals, vitamins, and phytochemicals beneficial to animal health. The fermentation of oat silage is a complex biochemical process that includes interactions between various microorganisms. The activity of many microbes in silage may cause an extensive breakdown of nutrition and lead to undesirable fermentation. Moreover, it is difficult to make high-quality oat silage because the number of epiphytic lactic acid bacterium microflora was lower than the requirement. Understanding the complex microbial community during the fermentation process and its relationship with community functions is therefore important in the context of developing improved fermentation biotechnology systems. These results suggested that the addition of Lactobacillus plantarum or Lactobacillus buchneri regulated the ensiling performance and microbial community in oat silage by shaping the metabolic pathways.

Variety-driven rhizosphere microbiome bestows differential salt tolerance to alfalfa for coping with salinity stress.

Soil salinization is a global environmental issue and a significant abiotic stress that threatens crop production. Root-associated rhizosphere microbiota play a pivotal role in enhancing plant tolerance to abiotic stresses. However, limited information is available concerning the specific variations in rhizosphere microbiota driven by different plant genotypes (varieties) in response to varying levels of salinity stress. In this study, we compared the growth performance of three alfalfa varieties with varying salt tolerance levels in soils with different degrees of salinization. High-throughput 16S rRNA and ITS sequencing were employed to analyze the rhizosphere microbial communities. Undoubtedly, the increasing salinity significantly inhibited alfalfa growth and reduced rhizosphere microbial diversity. However, intriguingly, salt-tolerant varieties exhibited relatively lower susceptibility to salinity, maintaining more stable rhizosphere bacterial community structure, whereas the reverse was observed for salt-sensitive varieties. Bacillus emerged as the dominant species in alfalfa's adaptation to salinity stress, constituting 21.20% of the shared bacterial genera among the three varieties. The higher abundance of Bacillus, Ensifer, and Pseudomonas in the rhizosphere of salt-tolerant alfalfa varieties is crucial in determining their elevated salt tolerance. As salinity levels increased, salt-sensitive varieties gradually accumulated a substantial population of pathogenic fungi, such as Fusarium and Rhizoctonia. Furthermore, rhizosphere bacteria of salt-tolerant varieties exhibited increased activity in various metabolic pathways, including biosynthesis of secondary metabolites, carbon metabolism, and biosynthesis of amino acids. It is suggested that salt-tolerant alfalfa varieties can provide more carbon sources to the rhizosphere, enriching more effective plant growth-promoting bacteria (PGPB) such as Pseudomonas to mitigate salinity stress. In conclusion, our results highlight the variety-mediated enrichment of rhizosphere microbiota in response to salinity stress, confirming that the high-abundance enrichment of specific dominant rhizosphere microbes and their vital roles play a significant role in conferring high salt adaptability to these varieties.

Effect of isolated lactic acid bacteria on the quality and bacterial diversity of native grass silage.

Objective: The objective of this study was to isolate lactic acid bacteria (LAB) from native grasses and naturally fermented silages, determine their identity, and assess their effects on silage quality and bacterial communities of the native grasses of three steppe types fermented for 60 days. Methods: Among the 58 isolated LAB strains, Limosilactobacillus fermentum (BL1) and Latilactobacillus graminis (BL5) were identified using 16S rRNA sequences. Both strains showed normal growth at 15- 45°C temperature, 3-6.5% NaCl concentration, and pH 4-9. Two isolated LAB strains (labeled L1 and L5) and two commercial additives (Lactiplantibacillus plantarum and Lentilactobacillus buchneri; designated as LP and LB, respectively) were added individually to native grasses of three steppe types (meadow steppe, MS; typical steppe, TS; desert steppe, DS), and measured after 60 d of fermentation. The fresh material (FM) of different steppe types was treated with LAB (1 × 105 colony forming units/g fresh weight) or distilled water (control treatment [CK]). Results: Compared with CK, the LAB treatment showed favorable effects on all three steppe types, i.e., reduced pH and increased water-soluble carbohydrate content, by modulating the microbiota. The lowest pH was found in the L5 treatment of three steppe types, at the same time, the markedly (p < 0.05) elevated acetic acid (AA) concentration was detected in the L1 and LB treatment. The composition of bacterial community in native grass silage shifted from Pantoea agglomerans and Rosenbergiella nectarea to Lentilactobacillus buchneri at the species level. The abundance of Lentilactobacillus buchneri and Lactiplantibacillus plantarum increased significantly in L1, L5, LP, and LB treatments, respectively, compared with CK (p < 0.05). Conclusion: In summary, the addition of LAB led to the shifted of microbiota and modified the quality of silage, and L. fermentum and L. graminis improved the performance of native grass silage.

Community Synergy of Lactic Acid Bacteria and Cleaner Fermentation of Oat Silage Prepared with a Multispecies Microbial Inoculant.

To investigate community synergy of lactic acid bacteria (LAB) and cleaner fermentation of oat silage, oat silages were prepared with or without (control) commercial LAB inoculants LI1 (containing Lactiplantibacillus plantarum, Lentilactobacillus buchneri, Lacticaseibacillus paracasei, and Pediococcus acidilactici) and LI2 (containing Lactiplantibacillus plantarum and Lentilactobacillus buchneri). The microbial community, LAB synergy, and cleaner fermentation were analyzed at 1, 3, 6, 15, 35, and 90 days of ensiling. The LAB inoculant improved fermentation quality, with significantly (P < 0.05) lower pH, ammonia nitrogen content, and gas production and higher lactic acid and acetic acid contents than those of the control. Enterobacteriaceae was the main bacterial community in early stage of fermentation, which utilizes sugar to produce CO2 gas, causing dry matter (DM) and energy loss. As fermentation progressed, the microbial diversity decreased, and the microbial community shifted from Gram-negative to Gram-positive bacteria. The inoculation of multispecies LAB displayed community synergy; Pediococcus acidilactici formed a dominant community in the early stage of fermentation, which produced an acid and anaerobic environment for the subsequent growth of Lentilactobacillus and Lacticaseibacillus species, thus forming a LAB-dominated microbial community. The predicted functional profile indicated that the silage inoculated with LI1 enhanced the carbohydrate metabolism pathway but inhibited the amino acid metabolism pathway, which played a role in promoting faster lactic acid production, reducing the decomposition of protein to ammonia nitrogen, and improving the fermentation quality of silage. Therefore, oat silage can be processed to high-quality and cleaner fermented feed by using an LAB inoculant, and LI1 showed better efficiency than LI2. IMPORTANCE Oat natural silage is rich in Enterobacteriaceae, increasing gas production and fermentation loss. Lactic acid bacteria interact synergistically to form a dominant community during ensiling. Pediococci grow vigorously in the early stage of fermentation and create an anaerobic environment. Lactobacilli inhibit the harmful microorganisms and result in cleaner fermentation of oat silage.

Integrated Analysis of Coexpression and Exome Sequencing to Prioritize Susceptibility Genes for Familial Cutaneous Melanoma.

The application of whole-exome sequencing has led to the identification of high- and moderate-risk variants that contribute to cutaneous melanoma susceptibility. However, confirming disease-causing variants remains challenging. We applied a gene coexpression network analysis to prioritize the candidate genes identified from whole-exome sequencing of 34 melanoma-prone families, with at least three affected members sequenced per family (N = 119 cases). A coexpression network was constructed from genotype-tissue expression project, skin melanoma from the cancer genome atlas, and primary melanocyte cultures. We performed module-specific enrichment and focused on modules associated with pigmentation processes because they are the best-studied and most well-known risk factors for melanoma susceptibility. We found that pigmentation-associated modules across the four expression datasets examined were enriched for well-known melanoma susceptibility genes plus genes associated with pigmentation. We also used network properties to prioritize genes within pigmentation modules as candidate susceptibility genes. Integrating information from coexpression network analysis and variant prioritization, we identified 36 genes (such as DCT, TPCN2, TRPM1, ATP10A, and EPHA5) as potential melanoma risk genes in the families. Our approach also allowed us to link families with private gene mutations on the basis of gene coexpression patterns and thereby may provide an innovative perspective in gene identification in high-risk families.

Phylogenomic Analysis of Human Papillomavirus Type 31 and Cervical Carcinogenesis: A Study of 2093 Viral Genomes.

Human papillomavirus (HPV) type 31 (HPV31) is closely related to the most carcinogenic type, HPV16, but only accounts for 4% of cervical cancer cases worldwide. Viral genetic and epigenetic variations have been associated with carcinogenesis for other high-risk HPV types, but little is known about HPV31. We sequenced 2093 HPV31 viral whole genomes from two large studies, one from the U.S. and one international. In addition, we investigated CpG methylation in a subset of 175 samples. We evaluated the association of HPV31 lineages/sublineages, single nucleotide polymorphisms (SNPs) and viral methylation with cervical carcinogenesis. HPV31 A/B clade was >1.8-fold more associated with cervical intraepithelial neoplasia grade 3 and cancer (CIN3+) compared to the most common C lineage. Lineage/sublineage distribution varied by race/ethnicity and geographic region. A viral genome-wide association analysis identified SNPs within the A/B clade associated with CIN3+, including H23Y (C626T) (odds ratio = 1.60, confidence intervals = 1.17-2.19) located in the pRb CR2 binding-site within the E7 oncogene. Viral CpG methylation was higher in lineage B, compared to the other lineages, and was most elevated in CIN3+. In conclusion, these data support the increased oncogenicity of the A/B lineages and suggest variation of E7 as a contributing risk factor.

Rare Germline Variants in Chordoma-Related Genes and Chordoma Susceptibility.

BACKGROUND: Chordoma is a rare bone cancer with an unknown etiology. TBXT is the only chordoma susceptibility gene identified to date; germline single nucleotide variants and copy number variants in TBXT have been associated with chordoma susceptibility in familial and sporadic chordoma. However, the genetic susceptibility of chordoma remains largely unknown. In this study, we investigated rare germline genetic variants in genes involved in TBXT/chordoma-related signaling pathways and other biological processes in chordoma patients from North America and China. METHODS: We identified variants that were very rare in general population and internal control datasets and showed evidence for pathogenicity in 265 genes in a whole exome sequencing (WES) dataset of 138 chordoma patients of European ancestry and in a whole genome sequencing (WGS) dataset of 80 Chinese patients with skull base chordoma. RESULTS: Rare and likely pathogenic variants were identified in 32 of 138 European ancestry patients (23%), including genes that are part of notochord development, PI3K/AKT/mTOR, Sonic Hedgehog, SWI/SNF complex and mesoderm development pathways. Rare pathogenic variants in COL2A1, EXT1, PDK1, LRP2, TBXT and TSC2, among others, were also observed in Chinese patients. CONCLUSION: We identified several rare loss-of-function and predicted deleterious missense variants in germline DNA from patients with chordoma, which may influence chordoma predisposition and reflect a complex susceptibility, warranting further investigation in large studies.

Using whole-exome sequencing and protein interaction networks to prioritize candidate genes for germline cutaneous melanoma susceptibility.

Although next-generation sequencing has demonstrated great potential for novel gene discovery, confirming disease-causing genes after initial discovery remains challenging. Here, we applied a network analysis approach to prioritize candidate genes identified from whole-exome sequencing analysis of 98 cutaneous melanoma patients from 27 families. Using a network propagation method, we ranked candidate genes by their similarity to known disease genes in protein-protein interaction networks and identified gene clusters with functional connectivity. Using this approach, we identified several new candidate susceptibility genes that warrant future investigations such as NGLY1, IL1RN, FABP2, PRKDC, and PROSER2. The propagated network analysis also allowed us to link families that did not have common underlying genes but that carried variants in genes that interact on protein-protein interaction networks. In conclusion, our study provided an analysis perspective for gene prioritization in the context of genetic heterogeneity across families and prioritized top potential candidate susceptibility genes in our dataset.

Mutations in the HPV16 genome induced by APOBEC3 are associated with viral clearance.

HPV16 causes half of cervical cancers worldwide; for unknown reasons, most infections resolve within two years. Here, we analyze the viral genomes of 5,328 HPV16-positive case-control samples to investigate mutational signatures and the role of human APOBEC3-induced mutations in viral clearance and cervical carcinogenesis. We identify four de novo mutational signatures, one of which matches the COSMIC APOBEC-associated signature 2. The viral genomes of the precancer/cancer cases are less likely to contain within-host somatic HPV16 APOBEC3-induced mutations (Fisher's exact test, P = 6.2 x 10-14), and have a 30% lower nonsynonymous APOBEC3 mutation burden compared to controls. We replicate the low prevalence of HPV16 APOBEC3-induced mutations in 1,749 additional cases. APOBEC3 mutations also historically contribute to the evolution of HPV16 lineages. We demonstrate that cervical infections with a greater burden of somatic HPV16 APOBEC3-induced mutations are more likely to be benign or subsequently clear, suggesting they may reduce persistence, and thus progression, within the host.

HPV16 Sublineage Associations With Histology-Specific Cancer Risk Using HPV Whole-Genome Sequences in 3200 Women.

BACKGROUND: HPV16 is a common sexually transmitted infection although few infections lead to cervical precancer/cancer; we cannot distinguish nor mechanistically explain why only certain infections progress. HPV16 can be classified into four main evolutionary-derived variant lineages (A, B, C, D) that have been previously suggested to have varying disease risks. METHODS: We used a high-throughput HPV16 whole-genome sequencing assay to investigate variant lineage risk among 3215 HPV16-infected women. Using sublineages A1/A2 as the reference, we assessed all variant lineage associations with infection outcome over three or more years of follow-up: 1107 control subjects (

Rare Germline Copy Number Variations and Disease Susceptibility in Familial Melanoma.

Mounting evidence suggests that copy number variations (CNVs) can contribute to cancer susceptibility. The main goal of this study was to evaluate the role of germline CNVs in melanoma predisposition in high-risk melanoma families. We used genome-wide tiling comparative genomic hybridization and single nucleotide polymorphism arrays to characterize CNVs in 335 individuals (240 melanoma cases) from American melanoma-prone families (22 with germline CDKN2A or CDK4 mutations). We found that the global burden of overall CNVs (or deletions or duplications separately) was not significantly associated with case-control or CDKN2A/CDK4 mutation status after accounting for the familial dependence. However, we identified several rare CNVs that either involved known melanoma genes (e.g., PARP1, CDKN2A) or cosegregated with melanoma (duplication on 10q23.23, 3p12.2 and deletions on 8q424.3, 2q22.1) in families without mutations in known melanoma high-risk genes. Some of these CNVs were correlated with expression changes in disrupted genes based on RNASeq data from a subset of melanoma cases included in the CNV study. These results suggest that rare cosegregating CNVs may influence melanoma susceptibility in some melanoma-prone families and genes found in our study warrant further evaluation in future genetic analyses of melanoma.