A structural variation reference for medical and population genetics.

Structural variants (SVs) rearrange large segments of DNA1 and can have profound consequences in evolution and human disease2,3. As national biobanks, disease-association studies, and clinical genetic testing have grown increasingly reliant on genome sequencing, population references such as the Genome Aggregation Database (gnomAD)4 have become integral in the interpretation of single-nucleotide variants (SNVs)5. However, there are no reference maps of SVs from high-coverage genome sequencing comparable to those for SNVs. Here we present a reference of sequence-resolved SVs constructed from 14,891 genomes across diverse global populations (54% non-European) in gnomAD. We discovered a rich and complex landscape of 433,371 SVs, from which we estimate that SVs are responsible for 25-29% of all rare protein-truncating events per genome. We found strong correlations between natural selection against damaging SNVs and rare SVs that disrupt or duplicate protein-coding sequence, which suggests that genes that are highly intolerant to loss-of-function are also sensitive to increased dosage6. We also uncovered modest selection against noncoding SVs in cis-regulatory elements, although selection against protein-truncating SVs was stronger than all noncoding effects. Finally, we identified very large (over one megabase), rare SVs in 3.9% of samples, and estimate that 0.13% of individuals may carry an SV that meets the existing criteria for clinically important incidental findings7. This SV resource is freely distributed via the gnomAD browser8 and will have broad utility in population genetics, disease-association studies, and diagnostic screening.

Paired-Duplication Signatures Mark Cryptic Inversions and Other Complex Structural Variation.

Copy-number variants (CNVs) have been the predominant focus of genetic studies of structural variation, and chromosomal microarray (CMA) for genome-wide CNV detection is the recommended first-tier genetic diagnostic screen in neurodevelopmental disorders. We compared CNVs observed by CMA to the structural variation detected by whole-genome large-insert sequencing in 259 individuals diagnosed with autism spectrum disorder (ASD) from the Simons Simplex Collection. These analyses revealed a diverse landscape of complex duplications in the human genome. One remarkably common class of complex rearrangement, which we term dupINVdup, involves two closely located duplications ("paired duplications") that flank the breakpoints of an inversion. This complex variant class is cryptic to CMA, but we observed it in 8.1% of all subjects. We also detected other paired-duplication signatures and duplication-mediated complex rearrangements in 15.8% of all ASD subjects. Breakpoint analysis showed that the predominant mechanism of formation of these complex duplication-associated variants was microhomology-mediated repair. On the basis of the striking prevalence of dupINVdups in this cohort, we explored the landscape of all inversion variation among the 235 highest-quality libraries and found abundant complexity among these variants: only 39.3% of inversions were canonical, or simple, inversions without additional rearrangement. Collectively, these findings indicate that dupINVdups, as well as other complex duplication-associated rearrangements, represent relatively common sources of genomic variation that is cryptic to population-based microarray and low-depth whole-genome sequencing. They also suggest that paired-duplication signatures detected by CMA warrant further scrutiny in genetic diagnostic testing given that they might mark complex rearrangements of potential clinical relevance.

Association Between Titin Loss-of-Function Variants and Early-Onset Atrial Fibrillation.

Importance: Atrial fibrillation (AF) is the most common arrhythmia affecting 1% of the population. Young individuals with AF have a strong genetic association with the disease, but the mechanisms remain incompletely understood. Objective: To perform large-scale whole-genome sequencing to identify genetic variants related to AF. Design, Setting, and Participants: The National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine Program includes longitudinal and cohort studies that underwent high-depth whole-genome sequencing between 2014 and 2017 in 18 526 individuals from the United States, Mexico, Puerto Rico, Costa Rica, Barbados, and Samoa. This case-control study included 2781 patients with early-onset AF from 9 studies and identified 4959 controls of European ancestry from the remaining participants. Results were replicated in the UK Biobank (346 546 participants) and the MyCode Study (42 782 participants). Exposures: Loss-of-function (LOF) variants in genes at AF loci and common genetic variation across the whole genome. Main Outcomes and Measures: Early-onset AF (defined as AF onset in persons <66 years of age). Due to multiple testing, the significance threshold for the rare variant analysis was P = 4.55 × 10-3. Results: Among 2781 participants with early-onset AF (the case group), 72.1% were men, and the mean (SD) age of AF onset was 48.7 (10.2) years. Participants underwent whole-genome sequencing at a mean depth of 37.8 fold and mean genome coverage of 99.1%. At least 1 LOF variant in TTN, the gene encoding the sarcomeric protein titin, was present in 2.1% of case participants compared with 1.1% in control participants (odds ratio [OR], 1.76 [95% CI, 1.04-2.97]). The proportion of individuals with early-onset AF who carried a LOF variant in TTN increased with an earlier age of AF onset (P value for trend, 4.92 × 10-4), and 6.5% of individuals with AF onset prior to age 30 carried a TTN LOF variant (OR, 5.94 [95% CI, 2.64-13.35]; P = 1.65 × 10-5). The association between TTN LOF variants and AF was replicated in an independent study of 1582 patients with early-onset AF (cases) and 41 200 control participants (OR, 2.16 [95% CI, 1.19-3.92]; P = .01). Conclusions and Relevance: In a case-control study, there was a statistically significant association between an LOF variant in the TTN gene and early-onset AF, with the variant present in a small percentage of participants with early-onset AF (the case group). Further research is necessary to understand whether this is a causal relationship.

Case-control analysis identifies shared properties of rare germline variation in cancer predisposing genes.

Along with traditional effects of aging and carcinogen exposure-inherited DNA variation has substantial contribution to cancer risk. Extraordinary progress made in analysis of common variation with GWAS methodology does not provide sufficient resolution to understand rare variation. To fulfill missing classification for rare germline variation we assembled dataset of whole exome sequences from>2000 patients (selected cases tested negative for candidate genes and unselected cases) with different types of cancers (breast cancer, colon cancer, and cutaneous and ocular melanomas) matched to more than 7000 non-cancer controls and analyzed germline variation in known cancer predisposing genes to identify common properties of disease-associated DNA variation and aid the future searches for new cancer susceptibility genes. Cancer predisposing genes were divided into non-overlapping classes according to the mode of inheritance of the related cancer syndrome or known tumor suppressor activity. Out of all classes only genes linked to dominant syndromes presented significant rare germline variants enrichment in cases. Separate analysis of protein-truncating and missense variation in this list of genes confirmed significant prevalence of protein-truncating variants in cases only in loss-of-function tolerant genes (pLI < 0.1), while ultra-rare missense variants were significantly overrepresented in cases only in constrained genes (pLI > 0.9). In addition to findings in genetically enriched cases, we observed significant burden of rare variation in unselected cases, suggesting substantial role of inherited variation even in relatively late cancer manifestation. Taken together, our findings provide reference for distribution and types of DNA variation underlying inherited predisposition to some common cancer types.

A Multinational Arab Genome-Wide Association Study Identifies New Genetic Associations for Rheumatoid Arthritis.

OBJECTIVE: Genetic factors underlying susceptibility to rheumatoid arthritis (RA) in Arab populations are largely unknown. This genome-wide association study (GWAS) was undertaken to explore the generalizability of previously reported RA loci to Arab subjects and to discover new Arab-specific genetic loci. METHODS: The Genetics of Rheumatoid Arthritis in Some Arab States Study was designed to examine the genetics and clinical features of RA patients from Jordan, the Kingdom of Saudi Arabia, Lebanon, Qatar, and the United Arab Emirates. In total, >7 million single-nucleotide polymorphisms (SNPs) were tested for association with RA overall and with seropositive or seronegative RA in 511 RA cases and 352 healthy controls. In addition, replication of 15 signals was attempted in 283 RA cases and 221 healthy controls. A genetic risk score of 68 known RA SNPs was also examined in this study population. RESULTS: Three loci (HLA region, intergenic 5q13, and 17p13 at SMTNL2/GGT6) reached genome-wide significance in the analyses of association with RA and with seropositive RA, and for all 3 loci, evidence of independent replication was demonstrated. Consistent with the findings in European and East Asian populations, the association of RA with HLA-DRB1 amino acid position 11 conferred the strongest effect (P = 4.8 × 10-16 ), and a weighted genetic risk score of previously associated RA loci was found to be associated with RA (P = 3.41 × 10-5 ) and with seropositive RA (P = 1.48 × 10-6 ) in this population. In addition, 2 novel associations specific to Arab populations were found at the 5q13 and 17p13 loci. CONCLUSION: This first RA GWAS in Arab populations confirms that established HLA-region and known RA risk alleles contribute strongly to the risk and severity of disease in some Arab groups, suggesting that the genetic architecture of RA is similar across ethnic groups. Moreover, this study identified 2 novel RA risk loci in Arabs, offering further population-specific insights into the pathophysiology of RA.

Defining the diverse spectrum of inversions, complex structural variation, and chromothripsis in the morbid human genome.

BACKGROUND: Structural variation (SV) influences genome organization and contributes to human disease. However, the complete mutational spectrum of SV has not been routinely captured in disease association studies. RESULTS: We sequenced 689 participants with autism spectrum disorder (ASD) and other developmental abnormalities to construct a genome-wide map of large SV. Using long-insert jumping libraries at 105X mean physical coverage and linked-read whole-genome sequencing from 10X Genomics, we document seven major SV classes at ~5 kb SV resolution. Our results encompass 11,735 distinct large SV sites, 38.1% of which are novel and 16.8% of which are balanced or complex. We characterize 16 recurrent subclasses of complex SV (cxSV), revealing that: (1) cxSV are larger and rarer than canonical SV; (2) each genome harbors 14 large cxSV on average; (3) 84.4% of large cxSVs involve inversion; and (4) most large cxSV (93.8%) have not been delineated in previous studies. Rare SVs are more likely to disrupt coding and regulatory non-coding loci, particularly when truncating constrained and disease-associated genes. We also identify multiple cases of catastrophic chromosomal rearrangements known as chromoanagenesis, including somatic chromoanasynthesis, and extreme balanced germline chromothripsis events involving up to 65 breakpoints and 60.6 Mb across four chromosomes, further defining rare categories of extreme cxSV. CONCLUSIONS: These data provide a foundational map of large SV in the morbid human genome and demonstrate a previously underappreciated abundance and diversity of cxSV that should be considered in genomic studies of human disease.

HIV-1 pol diversity among female bar and hotel workers in Northern Tanzania.

A national ART program was launched in Tanzania in October 2004. Due to the existence of multiple HIV-1 subtypes and recombinant viruses co-circulating in Tanzania, it is important to monitor rates of drug resistance. The present study determined the prevalence of HIV-1 drug resistance mutations among ART-naive female bar and hotel workers, a high-risk population for HIV-1 infection in Moshi, Tanzania. A partial HIV-1 pol gene was analyzed by single-genome amplification and sequencing in 45 subjects (622 pol sequences total; median number of sequences per subject, 13; IQR 5-20) in samples collected in 2005. The prevalence of HIV-1 subtypes A1, C, and D, and inter-subtype recombinant viruses, was 36%, 29%, 9% and 27%, respectively. Thirteen different recombination patterns included D/A1/D, C/A1, A1/C/A1, A1/U/A1, C/U/A1, C/A1, U/D/U, D/A1/D, A1/C, A1/C, A2/C/A2, CRF10_CD/C/CRF10_CD and CRF35_AD/A1/CRF35_AD. CRF35_AD was identified in Tanzania for the first time. All recombinant viruses in this study were unique, suggesting ongoing recombination processes among circulating HIV-1 variants. The prevalence of multiple infections in this population was 16% (n = 7). Primary HIV-1 drug resistance mutations to RT inhibitors were identified in three (7%) subjects (K65R plus Y181C; N60D; and V106M). In some subjects, polymorphisms were observed at the RT positions 41, 69, 75, 98, 101, 179, 190, and 215. Secondary mutations associated with NNRTIs were observed at the RT positions 90 (7%) and 138 (6%). In the protease gene, three subjects (7%) had M46I/L mutations. All subjects in this study had HIV-1 subtype-specific natural polymorphisms at positions 36, 69, 89 and 93 that are associated with drug resistance in HIV-1 subtype B. These results suggested that HIV-1 drug resistance mutations and natural polymorphisms existed in this population before the initiation of the national ART program. With increasing use of ARV, these results highlight the importance of drug resistance monitoring in Tanzania.

Frequent intra-subtype recombination among HIV-1 circulating in Tanzania.

The study estimated the prevalence of HIV-1 intra-subtype recombinant variants among female bar and hotel workers in Tanzania. While intra-subtype recombination occurs in HIV-1, it is generally underestimated. HIV-1 env gp120 V1-C5 quasispecies from 45 subjects were generated by single-genome amplification and sequencing (median (IQR) of 38 (28-50) sequences per subject). Recombination analysis was performed using seven methods implemented within the recombination detection program version 3, RDP3. HIV-1 sequences were considered recombinant if recombination signals were detected by at least three methods with p-values of ≤0.05 after Bonferroni correction for multiple comparisons. HIV-1 in 38 (84%) subjects showed evidence for intra-subtype recombination including 22 with HIV-1 subtype A1, 13 with HIV-1 subtype C, and 3 with HIV-1 subtype D. The distribution of intra-patient recombination breakpoints suggested ongoing recombination and showed selective enrichment of recombinant variants in 23 (60%) subjects. The number of subjects with evidence of intra-subtype recombination increased from 29 (69%) to 36 (82%) over one year of follow-up, although the increase did not reach statistical significance. Adjustment for intra-subtype recombination is important for the analysis of multiplicity of HIV infection. This is the first report of high prevalence of intra-subtype recombination in the HIV/AIDS epidemic in Tanzania, a region where multiple HIV-1 subtypes co-circulate. HIV-1 intra-subtype recombination increases viral diversity and presents additional challenges for HIV-1 vaccine design.

HIV-1 subtypes and recombinants in Northern Tanzania: distribution of viral quasispecies.

This study analyzed the distribution and prevalence of HIV-1 subtypes, multiplicity of HIV-1 infection, and frequency of inter-subtype recombination among HIV-1-infected female bar and hotel workers in Moshi, Kilimanjaro Region, Tanzania, from 2004 to 2007. The HIV-1 viral sequences spanning the V1-C5 region of HIV-1 env gp120 were analyzed from 50 subjects by single genome amplification and sequencing (SGA/S) technique. A total of 1740 sequences were amplified and sequenced from the HIV-1 proviral DNA template. The median env sequences analyzed per subject per two time points was 38 (IQR 28-50) over one year of HIV infection. In a subset of 14 subjects, a total of 239 sequences were obtained from HIV-1 RNA template at the baseline visit. The most prevalent HIV-1 subtypes were A1 (56%) and C (30%), while HIV-1 subtype D and inter-subtype recombinant viruses were found in 6% and 8% of subjects respectively. Transmission of multiple HIV-1 variants was evident in 27% of the subjects infected with pure HIV-1 subtypes A1, C, or D. The HIV-1 inter-subtype recombinants were found in 8% including HIV-1 C/A, D/A, and complex mosaic recombinants. Multiple viral variants were found in two subjects infected with inter-subtype recombinants. One subject harbored quasispecies of both pure HIV-1 A1 and C/A recombinant. The other subject was infected with two complex mosaic inter-subtype recombinant variants belonging to subtype D. HIV-1 multiple infections and ongoing recombination contribute significantly to the genetic diversity of circulating HIV-1 in Tanzania and have important implications for vaccine design and the development of therapeutic strategies.

Transmission of single and multiple viral variants in primary HIV-1 subtype C infection.

To address whether sequences of viral gag and env quasispecies collected during the early post-acute period can be utilized to determine multiplicity of transmitted HIV's, recently developed approaches for analysis of viral evolution in acute HIV-1 infection [1,2] were applied. Specifically, phylogenetic reconstruction, inter- and intra-patient distribution of maximum and mean genetic distances, analysis of Poisson fitness, shape of highlighter plots, recombination analysis, and estimation of time to the most recent common ancestor (tMRCA) were utilized for resolving multiplicity of HIV-1 transmission in a set of viral quasispecies collected within 50 days post-seroconversion (p/s) in 25 HIV-infected individuals with estimated time of seroconversion. The decision on multiplicity of HIV infection was made based on the model's fit with, or failure to explain, the observed extent of viral sequence heterogeneity. The initial analysis was based on phylogeny, inter-patient distribution of maximum and mean distances, and Poisson fitness, and was able to resolve multiplicity of HIV transmission in 20 of 25 (80%) cases. Additional analysis involved distribution of individual viral distances, highlighter plots, recombination analysis, and estimation of tMRCA, and resolved 4 of the 5 remaining cases. Overall, transmission of a single viral variant was identified in 16 of 25 (64%) cases, and transmission of multiple variants was evident in 8 of 25 (32%) cases. In one case multiplicity of HIV-1 transmission could not be determined. In primary HIV-1 subtype C infection, samples collected within 50 days p/s and analyzed by a single-genome amplification/sequencing technique can provide reliable identification of transmission multiplicity in 24 of 25 (96%) cases. Observed transmission frequency of a single viral variant and multiple viral variants were within the ranges of 64% to 68%, and 32% to 36%, respectively.

Evolutionary gamut of in vivo Gag substitutions during early HIV-1 subtype C infection.

Two analyses of HIV-1 subtype C Gag quasispecies were performed in a prospective cohort of 42 acutely and recently infected individuals by SGA on viral RNA/proviral DNA templates. First, in vivo Gag substitutions were assessed in relation to the HIV-1C consensus sequence, which revealed that 29.3% of detected amino acid substitutions can be classified as reversions to subtype consensus, 61.3% as forward substitutions from subtype consensus, and 9.3% as polymorphisms not associated with the subtype consensus sequence. Second, the proportion, dynamics, and relationships within individual pools of viral quasispecies were analyzed. Among reverse substitutions, 16.1% were minor, 11.0% transient, 13.6% dominant, and 59.2% fixed. In contrast, 31.6% of forward substitutions were minor, 59.3% transient, 3.8% dominant, and 5.3% fixed. The distinct patterns in the spectrum and dynamics of reverse and forward Gag substitutions suggest that these differences should be considered in HIV-1 evolutionary studies and analyses of viral mutational pathways.

Dynamics and timing of in vivo mutations at Gag residue 242 during primary HIV-1 subtype C infection.

Viral mutations at Gag residue 242 and relevant viral polymorphisms were analyzed in a cohort of 42 individuals with primary HIV-1 subtype C infection using single-genome amplification/sequencing. In HLA-B*57/5801-negative subjects infected with 242N escape variant, reversion to Asn appeared at median (IQR) 103 days (97-213 days) post-seroconversion (p/s) and became dominant at 193 days (170-215 days) p/s. In subjects expressing HLA-B*57/5801 and infected with the wild-type virus, the T242N escape appeared at 203 days (196-231) p/s, reached dominance at 277 days (265-315 days) p/s, and became complete at 323 days (289-373 days) p/s. HLA-B*57/5801-negative subjects infected with 242N escape variant did not show reduced viral load or increased CD4 count. The study highlights the differential selection of T242N escape by HLA-B*57 and B*5801 and suggests that the presence of HLA-B*57/5801-mediated immune pressure is able to control replication of the wild-type virus encoding Thr at Gag residue 242 but fails to suppress the T242N escape variant.

Evolution of proviral gp120 over the first year of HIV-1 subtype C infection.

The evolution of proviral gp120 during the first year after seroconversion in HIV-1 subtype C infection was addressed in a case series of eight subjects. Multiple viral variants were found in two out of eight cases. Slow rate of viral RNA decline and high early viral RNA set point were associated with a higher level of proviral diversity from 0 to 200 days after seroconversion. Proviral divergence from MRCA over the same period also differed between subjects with slow and fast decline of viral RNA, suggesting that evolution of proviral gp120 early in infection may be linked to the level of viral RNA replication. Changes in the length of variable loops were minimal, and length reduction was more common than length increase. Potential N-linked glycosylation sites ranged +/-one site, showing common fluctuations in the V4 and V5 loops. These results highlight the role of proviral gp120 diversity and diversification in the pathogenesis of acute HIV-1 subtype C infection.

Timing constraints of in vivo gag mutations during primary HIV-1 subtype C infection.

BACKGROUND: Aiming to answer the broad question "When does mutation occur?" this study examined the time of appearance, dominance, and completeness of in vivo Gag mutations in primary HIV-1 subtype C infection. METHODS: A primary HIV-1C infection cohort comprised of 8 acutely and 34 recently infected subjects were followed frequently up to 500 days post-seroconversion (p/s). Gag mutations were analyzed by employing single-genome amplification and direct sequencing. Gag mutations were determined in relation to the estimated time of seroconversion. Time of appearance, dominance, and completeness was compared for different types of in vivo Gag mutations. RESULTS: Reverse mutations to the wild type appeared at a median (IQR) of 62 (44;139) days p/s, while escape mutations from the wild type appeared at 234 (169;326) days p/s (p<0.001). Within the subset of mutations that became dominant, reverse and escape mutations appeared at 54 (30;78) days p/s and 104 (47;198) days p/s, respectively (p<0.001). Among the mutations that reached completeness, reverse and escape mutations appeared at 54 (30;78) days p/s and 90 (44;196) days p/s, respectively (p=0.006). Time of dominance for reverse mutations to and escape mutations from the wild type was 58 (44;105) days p/s and 219 (90;326) days p/s, respectively (p<0.001). Time of completeness for reverse and escape mutations was 152 (100;176) days p/s and 243 (101;370) days p/s, respectively (p=0.001). Fitting a Cox proportional hazards model with frailties confirmed a significantly earlier time of appearance (hazard ratio (HR): 2.6; 95% CI: 2.3-3.0), dominance (4.8 (3.4-6.8)), and completeness (3.6 (2.3-5.5)) of reverse mutations to the wild type Gag than escape mutations from the wild type. Some complex mutational pathways in Gag included sequential series of reversions and escapes. CONCLUSIONS: The study identified the timing of different types of in vivo Gag mutations in primary HIV-1 subtype C infection in relation to the estimated time of seroconversion. Overall, the in vivo reverse mutations to the wild type occurred significantly earlier than escape mutations from the wild type. This shorter time to incidence of reverse mutations remained in the subsets of in vivo Gag mutations that reached dominance or completeness.