Pan-genome analysis of 33 genetically diverse rice accessions reveals hidden genomic variations.

Structural variations (SVs) and gene copy number variations (gCNVs) have contributed to crop evolution, domestication, and improvement. Here, we assembled 31 high-quality genomes of genetically diverse rice accessions. Coupling with two existing assemblies, we developed pan-genome-scale genomic resources including a graph-based genome, providing access to rice genomic variations. Specifically, we discovered 171,072 SVs and 25,549 gCNVs and used an Oryza glaberrima assembly to infer the derived states of SVs in the Oryza sativa population. Our analyses of SV formation mechanisms, impacts on gene expression, and distributions among subpopulations illustrate the utility of these resources for understanding how SVs and gCNVs shaped rice environmental adaptation and domestication. Our graph-based genome enabled genome-wide association study (GWAS)-based identification of phenotype-associated genetic variations undetectable when using only SNPs and a single reference assembly. Our work provides rich population-scale resources paired with easy-to-access tools to facilitate rice breeding as well as plant functional genomics and evolutionary biology research.

Features of the Copy Number Variation of Certain Genes in Tumor Cells in Patients with Serous Ovarian Adenocarcinoma.

We analyzed the peculiarities of the copy number variation of genes that regulate apoptosis, DNA repair, cell proliferation, metabolism, and estrogen reception in tumor and normal cells of high-grade and low-grade serous adenocarcinoma of the ovaries. Using real-time qPCR method, the relative copy number of 34 genes (BAX, BCL2, TP53, MDM2, CASP9, CASP3, CASP7, CASP8, PRKCI, SOX2, OCT4, PIK3, PTEN, C-MYC, SOX18, AKT1, NOTCH1, BRCA1, BRCA2, EXO1, SCNN1A, KRAS, EGFR, BRAF, CYP1A1, CYP1A2, CYP1B1, CYP19A, ESR1, ESR2, GPER, STS, SULT1A, and SULT1E1) was determined in normal and tumor cells of the ovaries extracted by contactless capture laser microsection from FFPE-blocks of 200 patients. The most typical molecular markers of ovarian serous adenocarcinoma cells were identified: copy number of PIK3CA, BCL2, BAX, CASP3, and CASP8 genes. Based on the differences in the gene copy number variation, two molecular subtypes of serous adenocarcinoma were identified, corresponding to two histological subtypes: high-grade (MDM2, SOX2, ESR1, CYP1B1, SULT1E1, TP53, BRCA2) and low-grade (PIK3CA, PTEN, BCL2, BAX, and CASP3). Each of these subtypes was also characterized by molecular heterogeneity and can be subdivided into several subgroups: 3 subgroups for high-grade and 4 subgroups for low-grade serous adenocarcinoma. These findings extend our understanding of the molecular mechanisms of carcinogenesis in the ovarian tissue, confirm molecular difference between the two histological subtypes of serous adenocarcinoma probably underlying their different clinical course.

[Whole genomic copy number variation score: an indicator of potential diagnostic and prognostic value for lung adenocarcinoma].

Objective: To verify the value of whole genomic copy number variation (WGCNV) detection and scoring system in the diagnosis and prognosis of lung adenocarcinoma. Methods: Seventy-six lung adenocarcinoma specimens including ninety-one tumor samples and twenty adjacent non-tumor lung tissue samples were collected using Laser capture microdissection (LCM). Whole genomic amplification (WGA) was used to enrich DNA and construct a sequencing library for next generation sequencing (NGS). Changes of larger than 5Mb CNV in this study were analyzed and scored. The nuclear grading and score of tumor cells in the surgery and pleural effusion cytology of lung adenocarcinoma specimens were evaluated separately. For each case, we evaluated (1) nuclear size, (2) mitotic counts, (3) nuclear atypia, (4) atypical mitoses. The data of disease-free survive (DFS) and overall survive (OS) were collected for assessing the prognostic value of WGCNV score. Meanwhile, receiver operating characteristic (ROC) and area under curve (AUC) were used to define a cut-off value and evaluate the diagnostic significance in lung adenocarcinoma. Results: The WGCNV scores of twenty adjacent non-tumor lung tissue samples were treated as normal control and all of WGCNV scores of tumor samples range from 0 to 9.95, the median score was 2.7. The WGCNV scores were divided into three groups: low score group <1.74, medium score grade 1.74~4.23, high score grade >4.23. The WGCNV score was positively associated with the nuclear grade scoring (r=0.780 90, P<0.001). The result for evaluation of prognostic value of the WGCNV scores showed that comparing with low WGCNV score group, Hazard Ratio (HR) of medium score group was 4.11 (95%CI=0.72~23.57) and high score group was 2.07 (95%CI=0.30~14.12). These results suggested that the risks of the medium and high WGCNV score group elevated. According to the analysis results of ROC curve, when the cut off value was 0.01, the sensitivity and specificity for lung adenocarcinoma diagnosis were 97.8% and 95.0% respectively, the positive predictive value (PPV) and negative predictive value (NPV) were 99.0% and 90.1%, respectively, the AUC was 0.981. In the differentiation of adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma (MIA) group and invasive adenocarcinoma group, when the cut off value was 1.8, the sensitivity and specificity between the two groups were 78.1% and 94.4%, and the PPV and NPV were 98.0% and 52.0%, respectively, the AUC was 0.896. Conclusion: This study verifies that WGCNV scoring system has a potential diagnostic and prognostic value in lung adenocarcinoma.

Glyphosate Resistance and EPSPS Gene Duplication: Convergent Evolution in Multiple Plant Species.

One of the increasingly widespread mechanisms of resistance to the herbicide glyphosate is copy number variation (CNV) of the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene. EPSPS gene duplication has been reported in 8 weed species, ranging from 3 to 5 extra copies to more than 150 extra copies. In the case of Palmer amaranth (Amaranthus palmeri), a section of >300 kb containing EPSPS and many other genes has been replicated and inserted at new loci throughout the genome, resulting in significant increase in total genome size. The replicated sequence contains several classes of mobile genetic elements including helitrons, raising the intriguing possibility of extra-chromosomal replication of the EPSPS-containing sequence. In kochia (Kochia scoparia), from 3 to more than 10 extra EPSPS copies are arranged as a tandem gene duplication at one locus. In the remaining 6 weed species that exhibit EPSPS gene duplication, little is known about the underlying mechanisms of gene duplication or their entire sequence. There is mounting evidence that adaptive gene amplification is an important mode of evolution in the face of intense human-mediated selection pressure. The convergent evolution of CNVs for glyphosate resistance in weeds, through at least 2 different mechanisms, may be indicative of a more general importance for this mechanism of adaptation in plants. CNVs warrant further investigation across plant functional genomics for adaptation to biotic and abiotic stresses, particularly for adaptive evolution on rapid time scales.

Comparative analysis of the EGFR, HER2, c-MYC, and MET variations in colorectal cancer determined by three different measures: gene copy number gain, amplification status and the 2013 ASCO/CAP guideline criterion for HER2 testing of breast cancer.

BACKGROUND: The purpose of this study was to explore gene copy number (GCN) variation of EGFR, HER2, c-MYC, and MET in patients with primary colorectal cancer (CRC). METHODS: Dual-colour silver-enhanced in situ hybridization was performed in tissue samples of 334 primary CRC patients. The amplification status (GCN ratio ≥2) and GCN gain (average GCN ≥4) data for the EGFR, HER2, c-MYC and MET genes were obtained. GCN variation was also assessed by the criterion of the 2013 ASCO/CAP guidelines for HER2 testing. RESULTS: Amplification of EGFR, HER2, c-MYC and MET was detected in 8 (2.4%), 20 (6.0%), 29 (8.7%), and 14 (4.2%) patients, respectively. Of 66 patients with at least one amplified gene, five exhibited co-amplification of genes studied (HER2-MET co-amplification: two patients; HER2-c-MYC co-amplification: two patients; EGFR-c-MYC co-amplification: one patient). There were 109 patients with GCN gains of one or more genes (EGFR: 11/334, HER2: 29/334, c-MYC; 60/334, MET: 48/334) and 32.1% (35/109) had multiple GCN gains. When each GCN was assessed by the criterion of the ASCO/CAP 2013 guideline for HER2 testing, 116 people showed positive or equivocal results for one or more genes. The cumulative amplification status had no association with patients' outcome. However, the cumulative results of the GCN gain and GCN status determined according to the ASCO/CAP guideline had a significant prognostic correlation in the univariate analysis (P values of 0.006 and 0.022, respectively). In the multivariate analysis, GCN gain and GCN status were independent prognostic factors (P values of 0.010 and 0.017, respectively). CONCLUSIONS: In this study, we evaluated GCN variation of four genes in a large sample of Korean CRC patients. The amplification status was not related to patient outcome. However, the GCN gain and GCN status according to the ASCO/CAP 2013 guideline were independent prognostic factors.

Alpha-tryptase gene variation is associated with levels of circulating IgE and lung function in asthma.

BACKGROUND: Tryptase, a major secretory product of human mast cells has been implicated as a key mediator of allergic inflammation. Genetic variation in the tryptases is extensive, and α-tryptase, an allelic variant of the more extensively studied ß-tryptase, is absent in substantial numbers of the general population. The degree to which α-tryptase expression may be associated with asthma has not been studied. We have investigated the α-tryptase gene copy number variation and its potential associations with phenotypes of asthma. OBJECTIVES: Caucasian families (n = 341) with at least two asthmatic siblings (n = 1350) were genotyped for the α-tryptase alleles, using high-resolution melting assays. Standards for the possible α-/ß-tryptase ratios were constructed by cloning α-and ß-tryptase PCR products to generate artificial templates. Association analysis of asthma affection status and related phenotypes [total and allergen-specific serum IgE, bronchial hyperresponsiveness to methacholine, forced expiratory volume in 1s (FEV1 ) and atopy and asthma severity scores] was undertaken using family-based association tests (FBAT). RESULTS: Four consistent melting patterns for the α-tryptase genotype were identified with alleles carrying null, one or two copies of the α-tryptase allele. Possessing one copy of α-tryptase was significantly associated with lower serum levels of total and dust mite-specific IgE levels and higher FEV1 measurements, while two copies were related to higher serum concentrations of total and dust mite-specific IgE and greater atopy severity scores. CONCLUSIONS AND CLINICAL RELEVANCE: Associations of α-tryptase copy number with serum IgE levels, atopy scores and bronchial function may reflect roles for tryptases in regulating IgE production and other processes in asthma.

A review of the molecular mechanisms of acaricide resistance in mites and ticks.

The Arachnida subclass of Acari comprises many harmful pests that threaten agriculture as well as animal health, including herbivorous spider mites, the bee parasite Varroa, the poultry mite Dermanyssus and several species of ticks. Especially in agriculture, acaricides are often used intensively to minimize the damage they inflict, promoting the development of resistance. Beneficial predatory mites used in biological control are also subjected to acaricide selection in the field. The development and use of new genetic and genomic tools such as genome and transcriptome sequencing, bulked segregant analysis (QTL mapping), and reverse genetics via RNAi or CRISPR/Cas9, have greatly increased our understanding of the molecular genetic mechanisms of resistance in Acari, especially in the spider mite Tetranychus urticae which emerged as a model species. These new techniques allowed to uncover and validate new resistance mutations in a larger range of species. In addition, they provided an impetus to start elucidating more challenging questions on mechanisms of gene regulation of detoxification associated with resistance.

Association of PKLR gene copy number, expression levels and enzyme activity with 2,3,7,8-TCDD exposure in individuals exposed to Agent Orange/Dioxin in Vietnam.

2,3,7,8-tetrachlorodibenzo-p-dioxin (2,3,7,8-TCDD) is the most toxic congener of dioxin and has serious long-term effects on the environment and human health. Pyruvate Kinase L/R (PKLR) gene expression levels and gene variants are associated with pyruvate kinase enzyme deficiency, which has been identified as the cause of several diseases linked to dioxin exposure. In this study, we estimated PKLR gene copy number and gene expression levels using real-time quantitative PCR (RT-qPCR) assays, genotyped PKLR SNP rs3020781 by Sanger sequencing, and quantified plasma pyruvate kinase enzyme activity in 100 individuals exposed to Agent Orange/Dioxin near Bien Hoa and Da Nang airfields in Vietnam and 100 healthy controls. The means of PKLR copy numbers and PKLR gene expression levels were significantly higher, while pyruvate kinase enzyme activity was significantly decreased in Agent Orange/Dioxin-exposed individuals compared to healthy controls (P < 0.0001). Positive correlations of PKLR gene copy number and gene expression with 2,3,7,8-TCDD concentrations were observed (r = 0.2, P = 0.045 and r = 0.54, P < 0.0001, respectively). In contrast, pyruvate kinase enzyme activity was inversely correlated with 2,3,7,8-TCDD concentrations (r = -0.52, P < 0.0001). PKLR gene copy number and gene expression levels were also inversely correlated with pyruvate kinase enzyme activity. Additionally, PKLR SNP rs3020781 was found to be associated with 2,3,7,8-TCDD concentrations and PKLR gene expression. In conclusion, PKLR copy number, gene expression levels, and pyruvate kinase enzyme activity are associated with 2,3,7,8-TCDD exposure in individuals living in Agent Orange/Dioxin-contaminated areas.

Trypanosoma cruzi Genomic Variability: Array Comparative Genomic Hybridization Analysis of Clone and Parental Strain.

Trypanosoma cruzi, the etiological agent of Chagas disease, exhibits extensive inter- and intrastrain genetic diversity. As we have previously described, there are some genetic differences between the parental G strain and its clone D11, which was isolated by the limiting dilution method and infection of cultured mammalian cells. Electrophoretic karyotyping and Southern blot hybridization of chromosomal bands with specific markers revealed chromosome length polymorphisms of small size with additional chromosomal bands in clone D11 and the maintenance of large syntenic groups. Both G strain and clone D11 belong to the T. cruzi lineage TcI. Here, we designed intraspecific array-based comparative genomic hybridization (aCGH) to identify chromosomal regions harboring copy-number variations between clone D11 and the G strain. DNA losses were more extensive than DNA gains in clone D11. Most alterations were flanked by repeated sequences from multigene families that could be involved in the duplication and deletion events. Several rearrangements were detected by chromoblot hybridization and confirmed by aCGH. We have integrated the information of genomic sequence data obtained by aCGH to the electrophoretic karyotype, allowing the reconstruction of possible recombination events that could have generated the karyotype of clone D11. These rearrangements may be explained by unequal crossing over between sister or homologous chromatids mediated by flanking repeated sequences and unequal homologous recombination via break-induced replication. The genomic changes detected by aCGH suggest the presence of a dynamic genome that responds to environmental stress by varying the number of gene copies and generating segmental aneuploidy.

Clinical azole cross-resistance in Candida parapsilosis is related to a novel MRR1 gain-of-function mutation.

OBJECTIVES: Hereby, we describe the molecular mechanisms underlying the acquisition of azole resistance by a Candida parapsilosis isolate following fluconazole treatment due to candiduria. METHODS: A set of three consecutive C. parapsilosis isolates were recovered from the urine samples of a patient with candiduria. Whole-genome sequencing and antifungal susceptibility assays were performed. The expression of MRR1, MDR1, ERG11 and CDR1B (CPAR2_304370) was quantified by RT-qPCR. RESULTS: The initial isolate CPS-A was susceptible to all three azoles tested (fluconazole, voriconazole and posaconazole); isolate CPS-B, collected after the second cycle of treatment, exhibited a susceptible-dose-dependent phenotype to fluconazole and isolate CPS-C, recovered after the third cycle, exhibited a cross-resistance profile to fluconazole and voriconazole. Whole-genome sequencing revealed a putative resistance mechanism in isolate CPS-C, associated with a G1810A nucleotide substitution, leading to a G604R change in the Mrr1p transcription factor. Introducing this mutation into the susceptible CPS-A isolate (MRR1RI) resulted in resistance to fluconazole and voriconazole, as well as up-regulation of MRR1 and MDR1. Interestingly, the susceptible-dose-dependent phenotype exhibited by isolate CPS-B was associated with an increased copy number of the CDR1B gene. The expression of CDR1B was increased in both isolates CPS-B and CPS-C and in the MRR1RI strain, harbouring the gain-of-function mutation. CONCLUSIONS: Our results describe clinical azole cross-resistance acquisition in C. parapsilosis due to a G1810A (G604R) gain-of-function mutation, resulting in MRR1 hyperactivation and consequently, MDR1 efflux pump overexpression. We also associated amplification of the CDR1B gene with decreased fluconazole susceptibility and showed that it is a putative target of the MRR1 gain-of-function mutation.

Origin and Evolution of the Cannabinoid Oxidocyclase Gene Family.

Cannabis is an ancient crop representing a rapidly increasing legal market, especially for medicinal purposes. Medicinal and psychoactive effects of Cannabis rely on specific terpenophenolic ligands named cannabinoids. Recent whole-genome sequencing efforts have uncovered variation in multiple genes encoding the final steps in cannabinoid biosynthesis. However, the origin, evolution, and phylogenetic relationships of these cannabinoid oxidocyclase genes remain unclear. To elucidate these aspects, we performed comparative genomic analyses of Cannabis, related genera within the Cannabaceae family, and selected outgroup species. Results show that cannabinoid oxidocyclase genes originated in the Cannabis lineage from within a larger gene expansion in the Cannabaceae family. Localization and divergence of oxidocyclase genes in the Cannabis genome revealed two main syntenic blocks, each comprising tandemly repeated cannabinoid oxidocyclase genes. By comparing these blocks with those in genomes from closely related species, we propose an evolutionary model for the origin, neofunctionalization, duplication, and diversification of cannabinoid oxidocycloase genes. Based on phylogenetic analyses, we propose a comprehensive classification of three main clades and seven subclades that are intended to aid unequivocal referencing and identification of cannabinoid oxidocyclase genes. Our data suggest that cannabinoid phenotype is primarily determined by the presence/absence of single-copy genes. Although wild populations of Cannabis are still unknown, increased sampling of landraces and wild/feral populations across its native geographic range is likely to uncover additional cannabinoid oxidocyclase sequence variants.

Prognostic value and characterization of NTRK1 variation by fluorescence in situ hybridization in esophageal squamous cell carcinoma.

PURPOSE: Rearrangement of the neurotrophic tyrosine kinase receptor (NTRK) 1 gene is a target of tropomyosin receptor kinase A (TRKA) inhibitors, and its targeted drug (larotrectinib) has been approved by the US Food and Drug Administration. We investigated the existence and prognostic importance of NTRK1 variation in esophageal squamous cell carcinoma (ESCC). METHODS: Fluorescence in situ hybridization of a NTRK1 rearrangement was conducted on 523 ESCC samples through tissue microarrays. Kaplan-Meier curves with log-rank tests were used to evaluate survival. RESULTS: We identified 8 (1.5%), 35(6.7%) and 109 (20.8%) cases with a NTRK1 rearrangement using 15%, 10% and 5% as cut-off values, respectively. We observed copy number (CN) variation of NTRK1 in some cases: 79 (15.1%) cases had a gain in NTRK1 CN ≥ 3, and 24 (4.6%) cases had NTRK1 CN ≥ 4. A NTRK1 rearrangement at the above-mentioned thresholds was not related to disease-free survival (DFS, P = 0.45, 0.47, 0.87) and overall survival (OS, P = 0.80, 0.74, 0.57), respectively. Gain in NTRK1 CN was associated with a poor prognosis irrespective of whether NTRK1 CN ≥ 4 (DFS, P = 0.015; OS, P = 0.035) or NTRK1 CN ≥ 3 (DFS, P = 0.039; OS, P = 0.025). CONCLUSION: A NTRK1 rearrangement occurred rarely in ESCC. The increased CN of NTRK1 might be a prognostic indicator for DFS and OS in patients with ESCC.

Comparative genomics of Leishmania isolates from Brazil confirms the presence of Leishmania major in the Americas.

Leishmania (Leishmania) major is an important agent of cutaneous leishmaniasis, having as a vector sandflies belonging to the genus Phlebotomus. Although this species has been described as restricted to the Old World, parasites similar to L. major have been isolated from South American patients who have never travelled abroad. These parasites were named "L. major-like", and several studies have been carried out to characterise them biochemically, molecularly, and biologically. However, the phylogenetic origin of these isolates is still unknown. In the present study we characterised three L. major-like isolates, named BH49, BH121 and BH129, using comparative genomics approaches. We evaluated the presence of gene and segmental duplications/deletions and the presence of aneuploidies that could explain the differences in infectivity observed in the BH49 and BH121 isolates. All isolates presented a pattern of mosaic aneuploidy and gene copy number variation, which are common in the genus Leishmania. Virulence factors such as phosphatases and peptidases were found to have increased gene copy numbers in the infective isolate, which could explain the difference in infectivity previously observed between BH121 and BH49. Phylogenetic analyses revealed that BH49, BH121 and BH129 L. major-like grouped with L. major isolates, and suggest they were imported from the Old World in at least two independent events. We suggest that new epidemiological inquiries should also evaluate L. major infections in South America, to assess the epidemiological importance of this species in the New World.

Human Complement C4B Allotypes and Deficiencies in Selected Cases With Autoimmune Diseases.

Human complement C4 is one of the most diverse but heritable effectors for humoral immunity. To help understand the roles of C4 in the defense and pathogenesis of autoimmune and inflammatory diseases, we determined the bases of polymorphisms including the frequent genetic deficiency of C4A and/or C4B isotypes. We demonstrated the diversities of C4A and C4B proteins and their gene copy number variations (CNVs) in healthy subjects and patients with autoimmune disease, such as type 1 diabetes, systemic lupus erythematosus (SLE) and encephalitis. We identified subjects with (a) the fastest migrating C4B allotype, B7, or (b) a deficiency of C4B protein caused by genetic mutation in addition to gene copy-number variation. Those variants and mutants were characterized, sequenced and specific techniques for detection developed. Novel findings were made in four case series. First, the amino acid sequence determinant for C4B7 was likely the R729Q variation at the anaphylatoxin-like region. Second, in healthy White subject MS630, a C-nucleotide deletion at codon-755 led to frameshift mutations in his single C4B gene, which was a private mutation. Third, in European family E94 with multiplex lupus-related mortality and low serum C4 levels, the culprit was a recurrent haplotype with HLA-A30, B18 and DR7 that segregated with two defective C4B genes and identical mutations at the donor splice site of intron-28. Fourth, in East-Asian subject E133P with anti-NMDA receptor encephalitis, the C4B gene had a mutation that changed tryptophan-660 to a stop-codon (W660x), which was present in a haplotype with HLA-DRB1*04:06 and B*15:27. The W660x mutation is recurrent among East-Asians with a frequency of 1.5% but not detectable among patients with SLE. A meticulous annotation of C4 sequences revealed clusters of variations proximal to sites for protein processing, activation and inactivation, and binding of interacting molecules.

Both Methylation and Copy Number Variation Participated in the Varied Expression of PRAME in Multiple Myeloma.

PURPOSE: The cancer-testis antigen, which is a preferentially expressed antigen of melanoma (PRAME), is an ideal target for immunotherapy and cancer vaccines. Since the expression of this antigen is relevant to therapy responses, the heterogeneity in its expression and the underlying mechanism need to be investigated. PATIENTS AND METHODS: Plasma cell sorting was performed in 48 newly diagnosed multiple myeloma (MM) patients. Real-time quantitative PCR was performed to examine the PRAME transcript levels and gene copy numbers. Bisulfate clone sequencing of the PRAME promoter and exon 1b regions was performed in 4 patients. Quantitative methylation-specific PCR of the +287 CpG site was performed for all patients. The human MM cell lines RPMI8226, LP-1 and MOLP-2 were treated with 5-azacytidine. RESULTS: The median PRAME transcript level was 3.1% (range: 0-298.3%) in the plasma cells sorted from the 48 MM patients. Eleven (22.9%) and 37 (77.1%) patients were individually categorized into the PRAME low- and high-expression groups according to the cut-off value of 0.05%. The methylation ratios of the promoter and the 3' region of exon 1b region were both negatively related to the transcript levels. The degrees of methylation at the +287 CpG site were significantly negatively related to the transcript levels in all 48 patients (r=-0.44, P=0.0018), and those in the high-expression group (r=-0.69, P<0.0001) but not those in the low-expression group (r=-0.27, P=0.43). All 5 patients with homozygous deletions were categorized into the low-expression group. There were no significant differences in the PRAME transcript levels between the hemizygous deletion (n=8) and no deletion (n=35) groups (P=0.40). Furthermore, the PRAME transcript levels significantly increased in the MM cell lines after treatment with 5-azacytidine. CONCLUSION: Both methylation and copy number variation may participate in the regulation of PRAME expression in MM; in patients with no homozygous deletion, PRAME expression is mainly controlled by methylation, and a proportion of fairly low expression is caused by homozygous deletion.

Elevated serum complement levels and higher gene copy number of complement C4B are associated with hypertension and effective response to statin therapy in childhood-onset systemic lupus erythematosus (SLE).

OBJECTIVE: Systemic lupus erythematosus (SLE) features high frequency of cardiovascular disease (CVD) and fluctuating complement levels. The clinical trial Atherosclerosis Prevention in Pediatric Lupus Erythematosus (APPLE) aimed to evaluate whether atorvastatin treatment reduced the progression of atherosclerosis in 221 patients with childhood-onset SLE (cSLE), using carotid intima media thickness (CIMT) as surrogates. We leveraged APPLE biorepository and trial data to investigate the relationship between complement and CVD in cSLE. METHODS: Gene copy numbers (GCNs) for total C4, C4A and C4B were measured by TaqMan-based real-time PCR and Southern blotting, and analysed with laboratory and clinical parameters through Student's t-test and χ2 analyses. Effects of total C4, C4A and C4B GCNs on the response to placebo or atorvastatin treatment and progression of CIMT were examined by regression analyses. RESULTS: At baseline, C4 protein levels strongly correlated with GCNs of total C4 (p=1.8×10-6). Each copy of C4 gene increased mean serum C4 by 3.28 mg/dL. Compared with those without hypertension (N=142), individuals with hypertension demonstrated significantly elevated serum levels for C4 and C3 at baseline and serially (C4: P=5.0×10-25; C3: P=5.84×10-20). Individuals with ≥2 C4B genes had 2.5 times the odds of having hypertension (p=0.016) and higher diastolic blood pressure (p=0.015) compared with those with C4B deficiency. At the study end, subjects with ≥2 C4B and atorvastatin treatment had significantly slower increase in CIMT compared with those treated with placebo (p=0.018). CONCLUSIONS: cSLE with hypertension had elevated serum levels of C4 and C3 and higher GCN of C4B; cSLE with ≥2 C4B genes would benefit from statins therapy to prevent atherosclerosis.

Impact of C4, C4A and C4B gene copy number variation in the susceptibility, phenotype and progression of systemic lupus erythematosus.

BACKGROUND: Complement component 4 (C4) gene copy number (GCN) affects the susceptibility to systemic lupus erythematosus (SLE) in different populations, however the possible phenotype significance remains to be determined. This study aimed to associate C4A, C4B and total C4 GCN and SLE, focusing on the clinical phenotype and disease progression. METHODS: C4, C4A and C4B GCN were determined by real-time PCR in 427 SLE patients and 301 healthy controls, which underwent a detailed clinical evaluation according to a pre-established protocol. RESULTS: The risk of developing SLE was 2.62 times higher in subjects with low total C4 GCN (< 4 copies, OR = 2.62, CI = 1.77 to 3.87, p < 0.001) and 3.59 times higher in subjects with low C4A GCN (< 2 copies; OR = 3.59, CI = 2.15 to 5.99, p < 0.001) compared to those subjects with normal or high GCN of total C4 (≥4) and C4A (≥2), respectively. An increased risk was also observed regarding low C4B GCN, albeit to a lesser degree (OR = 1.46, CI = 1.03 to 2.08, p = 0.03). Furthermore, subjects with low C4A GCN had higher permanent disease damage as assessed by the Systemic Lupus International Collaborating Clinics - Damage Index (SLICC-DI; median = 1.5, 95% CI = 1.2-1.9) than patients with normal or high copy number of C4A (median = 1.0, 95% CI = 0.8-1.1; p = 0.004). There was a negative association between low C4A GCN and serositis (p = 0.02) as well as between low C4B GCN and arthritis (p = 0.02). CONCLUSIONS: This study confirms the association between low C4 GCN and SLE susceptibility, and originally demonstrates an association between low C4A GCN and disease severity.

Leishmania Genome Dynamics during Environmental Adaptation Reveal Strain-Specific Differences in Gene Copy Number Variation, Karyotype Instability, and Telomeric Amplification.

Protozoan parasites of the genus Leishmania adapt to environmental change through chromosome and gene copy number variations. Only little is known about external or intrinsic factors that govern Leishmania genomic adaptation. Here, by conducting longitudinal genome analyses of 10 new Leishmania clinical isolates, we uncovered important differences in gene copy number among genetically highly related strains and revealed gain and loss of gene copies as potential drivers of long-term environmental adaptation in the field. In contrast, chromosome rather than gene amplification was associated with short-term environmental adaptation to in vitro culture. Karyotypic solutions were highly reproducible but unique for a given strain, suggesting that chromosome amplification is under positive selection and dependent on species- and strain-specific intrinsic factors. We revealed a progressive increase in read depth towards the chromosome ends for various Leishmania isolates, which may represent a nonclassical mechanism of telomere maintenance that can preserve integrity of chromosome ends during selection for fast in vitro growth. Together our data draw a complex picture of Leishmania genomic adaptation in the field and in culture, which is driven by a combination of intrinsic genetic factors that generate strain-specific phenotypic variations, which are under environmental selection and allow for fitness gain.IMPORTANCE Protozoan parasites of the genus Leishmania cause severe human and veterinary diseases worldwide, termed leishmaniases. A hallmark of Leishmania biology is its capacity to adapt to a variety of unpredictable fluctuations inside its human host, notably pharmacological interventions, thus, causing drug resistance. Here we investigated mechanisms of environmental adaptation using a comparative genomics approach by sequencing 10 new clinical isolates of the L. donovani, L. major, and L. tropica complexes that were sampled across eight distinct geographical regions. Our data provide new evidence that parasites adapt to environmental change in the field and in culture through a combination of chromosome and gene amplification that likely causes phenotypic variation and drives parasite fitness gains in response to environmental constraints. This novel form of gene expression regulation through genomic change compensates for the absence of classical transcriptional control in these early-branching eukaryotes and opens new venues for biomarker discovery.

PCR-based quantification of taxa-specific abundances in microbial communities: Quantifying and avoiding common pitfalls.

The quantification of relative and absolute taxa-specific abundances in complex microbial communities is crucial for understanding and modeling natural and engineered ecosystems. Many errors inherent to this quantification are, though well-known, still insufficiently addressed and can potentially lead to a completely different interpretation of experimental results. This review provides a critical assessment of next generation sequencing (NGS) of amplicons and quantitative real-time PCR for the quantification of relative and absolute taxa-specific genome abundances. Starting from DNA extraction, the following error sources were considered: DNA extraction efficiency, PCR-associated bias, variance of strain-specific 16S rRNA operon copy number per genome, and analysis of quantitative real-time PCR and NGS data. Tools and methods for estimating and minimizing these errors are presented and demonstrated on published data. In conclusion, amplicon sequencing and qPCR of 16S rRNA genes are valuable tools to determine relative and absolute taxa-specific genome abundances, but results can deviate by several orders of magnitudes from the true values if the reviewed error sources are ignored. Many of these errors can be minimized in a cost-efficient manner and large errors can be easily identified by plausibility checks as shown in this review. Finally, the accurate conversion of genome abundances to cell numbers and microbial biomasses was pointed out as an important future research topic for the integration of PCR-based abundances into mathematical models.

A Driving Bioinformatics Approach to Explore Co-regulation of AOX Gene Family Members During Growth and Development.

The alternative oxidase (AOX) gene family is a hot candidate for functional marker development that could help plant breeding on yield stability through more robust plants based on multi-stress tolerance. However, there is missing knowledge on the interplay between gene family members that might interfere with the efficiency of marker development. It is common view that AOX1 and AOX2 have different physiological roles. Nevertheless, both family member groups act in terms of molecular-biochemical function as "typical" alternative oxidases and co-regulation of AOX1 and AOX2 had been reported. Although conserved sequence differences had been identified, the basis for differential effects on physiology regulation is not sufficiently explored.This protocol gives instructions for a bioinformatics approach that supports discovering potential interaction of AOX family members in regulating growth and development. It further provides a strategy to elucidate the relevance of gene sequence diversity and copy number variation for final functionality in target tissues and finally the whole plant. Thus, overall this protocol provides the means for efficiently identifying plant AOX variants as functional marker candidates related to growth and development.