Copy Number Polymorphisms Research Articles

Differentiation in detoxification gene complements, including neofunctionalization of duplicated cytochrome P450 genes, between lineages of cotton bollworm, Helicoverpa armigera.

Here we investigate the evolutionary dynamics of five enzyme superfamilies (CYPs, GSTs, UGTs, CCEs and ABCs) involved in detoxification in Helicoverpa armigera. The reference assembly for an African isolate of the major lineages, H. a. armigera, has 373 genes in the five superfamilies. Most of its CYPs, GSTs, UGTs and CCEs and a few of its ABCs occur in blocks and most of the clustered genes are in subfamilies specifically implicated in detoxification. Most of the genes have orthologues in the reference genome for the Oceania lineage, H. a. conferta. However, clustered orthologues and subfamilies specifically implicated in detoxification show greater sequence divergence and less constraint on non-synonymous differences between the two assemblies than do other members of the five superfamilies. Two duplicated CYPs, which were found in the H. a. armigera but not H. a. conferta reference genome, were also missing in 16 Chinese populations spanning two different lineages of H. a. armigera. The enzyme produced by one of these duplicates has higher activity against esfenvalerate than a previously described chimeric CYP mutant conferring pyrethroid resistance. Various transposable elements were found in the introns of most detoxification genes, generating diverse gene structures. Extensive resequencing data for the Chinese H. a. armigera and H. a. conferta lineages also revealed complex copy number polymorphisms in 17 CCE001s in a cluster also implicated in pyrethroid metabolism, with substantial haplotype differences between all three lineages. Our results suggest that cotton bollworm has a versatile complement of detoxification genes which are evolving in diverse ways across its range.

Ethnic and functional differentiation of copy number polymorphisms in Tunisian and HapMap population unveils insights on genome organizational plasticity

Admixture mapping has been useful in identifying genetic variations linked to phenotypes, adaptation and diseases. Copy number variations (CNVs) represents genomic structural variants spanning large regions of chromosomes reaching several megabases. In this investigation, the “Canary” algorithm was applied to 102 Tunisian samples and 991 individuals from eleven HapMap III populations to genotype 1279 copy number polymorphisms (CNPs). In this present work, we investigate the Tunisian population structure using the CNP makers previously identified among Tunisian. The study revealed that Sub-Saharan African populations exhibited the highest diversity with the highest proportions of allelic CNPs. Among all the African populations, Tunisia showed the least diversity. Individual ancestry proportions computed using STRUCTURE analysis revealed a major European component among Tunisians with lesser contribution from Sub-Saharan Africa and Asia. Population structure analysis indicated the genetic proximity with Europeans and noticeable distance from the Sub-Saharan African and East Asian clusters. Seven genes harbouring Tunisian high-frequent CNPs were identified known to be associated with 9 Mendelian diseases and/or phenotypes. Functional annotation of genes under selection highlighted a noteworthy enrichment of biological processes to receptor pathway and activity as well as glutathione metabolism. Additionally, pathways of potential concern for health such as drug metabolism, infectious diseases and cancers exhibited significant enrichment. The distinctive genetic makeup of the Tunisians might have been influenced by various factors including natural selection and genetic drift, resulting in the development of distinct genetic variations playing roles in specific biological processes. Our research provides a justification for focusing on the exclusive genome organization of this population and uncovers previously overlooked elements of the genome.

Genome-wide identification of candidate copy number polymorphism genes associated with complex traits of Tibetan-sheep

Copy number variation (CNV) is a genetic structural polymorphism important for phenotypic diversity and important economic traits of livestock breeds, and it plays an important role in the desired genetic variation. This study used whole genome sequencing to detect the CNV variation in the genome of 6 local Tibetan sheep groups. We detected 69,166 CNV events and 7230 copy number variable regions (CNVRs) after merging the overlapping CNVs, accounting for 2.72% of the reference genome. The CNVR length detected ranged from 1.1 to 1693.5 Kb, with a total length of 118.69 Mb and an average length of 16.42 Kb per CNVR. Functional GO cluster analysis showed that the CNVR genes were mainly involved in sensory perception systems, response to stimulus, and signal transduction. Through CNVR-based Vst analysis, we found that the CACNA2D3 and CTBP1 genes related to hypoxia adaptation, the HTR1A gene related to coat color, and the TRNAS-GGA and PIK3C3 genes related to body weight were all strongly selected. The findings of our study will contribute novel insights into the genetic structural variation underlying hypoxia adaptation and economically important traits in Tibetan sheep.

Genomic insights into linezolid-resistant Enterococci revealed its evolutionary diversity and poxtA copy number heterogeneity

ObjectivesThis study aimed to determine the molecular mechanisms of linezolid-resistant enterococci (LRE) in swine slaughterhouses in China and apply the “One Health” perspective to analyse the evolutionary dynamics of poxtA-positive E. faecium in clinical and non-clinical settings worldwide. MethodsThe phenotypic and genomic characteristics of multiple LRE isolates were systematically investigated using antimicrobial susceptibility testing, transfer assays, evolutionary experiments, quantitative RT-PCR assays, whole-genome sequencing, and bioinformatics analyses. ResultsSwine faeces served as a significant reservoir for LRE isolates, and optrA and poxtA were the primary contributors to linezolid resistance. Co-occurrence network analysis revealed a significant interconnection between optrA and several other ARGs. The poxtA copy number heterogeneity and polymorphism were initially observed in E. faecium parental and evolved isolates. The poxtA-carrying tandem repeat region exhibits high mobility and has undergone extensive duplication owing to linezolid pressure. The poxtA copy number varies from four copies on the plasmid of E. faecium IC25 to 11 copies on the plasmid and six copies on the chromosome in the evolved isolate IC25–50_poxtA. Furthermore, phylogenetic analysis of 185 poxtA-positive E. faecium strains worldwide found that one isolate from a French patient in 2018 shared only two SNPs with CC17 E. faecium isolates IC25 and IC7–2 from this study, highlighting the potential global transmission of CC17 poxtA-positive E. faecium between humans and animals. ConclusionThis study identified amplification of poxtA as a response of E. faecium to linezolid pressure. Phylogenetic analysis shed light on the potential global transmission of hospital-associated CC17 poxtA-positive E. faecium in clinical and non-clinical settings.

Y-Linked Copy Number Polymorphism of Target of Rapamycin Is Associated with Sexual Size Dimorphism in Seed Beetles.

The Y chromosome is theorized to facilitate evolution of sexual dimorphism by accumulating sexually antagonistic loci, but empirical support is scarce. Due to the lack of recombination, Y chromosomes are prone to degenerative processes, which poses a constraint on their adaptive potential. Yet, in the seed beetle, Callosobruchus maculatus segregating Y linked variation affects male body size and thereby sexual size dimorphism (SSD). Here, we assemble C. maculatus sex chromosome sequences and identify molecular differences associated with Y-linked SSD variation. The assembled Y chromosome is largely euchromatic and contains over 400 genes, many of which are ampliconic with a mixed autosomal and X chromosome ancestry. Functional annotation suggests that the Y chromosome plays important roles in males beyond primary reproductive functions. Crucially, we find that, besides an autosomal copy of the gene target of rapamycin (TOR), males carry an additional TOR copy on the Y chromosome. TOR is a conserved regulator of growth across taxa, and our results suggest that a Y-linked TOR provides a male specific opportunity to alter body size. A comparison of Y haplotypes associated with male size difference uncovers a copy number variation for TOR, where the haplotype associated with decreased male size, and thereby increased sexual dimorphism, has two additional TOR copies. This suggests that sexual conflict over growth has been mitigated by autosome to Y translocation of TOR followed by gene duplications. Our results reveal that despite of suppressed recombination, the Y chromosome can harbor adaptive potential as a male-limited supergene.

Complement C4 gene copy number variations and polymorphisms, autoantibodies, and clinical manifestations of juvenile dermatomyositis– a multi-center study

Abstract Juvenile dermatomyositis (JDM) is an autoimmune myopathy characterized by rash and muscle weakness. Complement C4 gene copy number (GCN) variation is a known risk factor for JDM. JDM patients often develop autoantibodies. We investigated the relationship between C4 GCN, clinical features, and antibody development. METHODS: Subjects were recruited (n=255) from the US, Sweden, and the Czech Republic. C4 GCN was determined by real time PCR. Comparative analyses were performed via student’s t test or Mann Whitney U test. Correlation was assessed via linear regression or Fisher’s exact test. All studies were IRB approved, and informed consent was obtained. RESULTS: C4 GCN correlated with muscle pathology and extra-muscular disease at diagnosis. Lower C4A and C4L GCN were associated with a higher number of abnormal muscle enzymes (p=0.0062 and p=0.0029), while the opposite was true for C4S (p=0.024). Higher C4S GCN correlated with lower muscle strength scores (p=0.039). Lower C4L correlated with abnormal MRI (p= 0.042). C4B GCN correlated with dysphagia (p=0.035). Higher GCN of C4S was a risk factor for arthritis (p=0.032), systemic symptoms (p=0.020), and dysphagia (p=0.0078). C4 GCN was associated with myositis-specific (MSA) and myositis-associated (MAA) antibodies. Patients with homozygous C4A deficiency were more likely to test positive for anti-NXP2 (OR 20.0, p=0.011), and C4A GCN was lower in anti-NXP2 positive subjects (p=0.030). A higher C4B GCN correlated with positive MAA (p=0.043). DISCUSSION: Here, we show that low C4A/C4L GCN and high C4B/C4S GCN correlate with worse muscle disease, extra-muscular pathology, and autoantibodies. Our data suggest that complement C4 may play a key role in the ongoing pathogenesis of JDM. Supported by grants from NIH (R21 AR070509) and the CureJM Foundation

Radiation-Induced Chromosome Instability in WTK1 and TK6 Human Lymphoblastoid Cells

The effects of ionizing radiation on single nucleotide polymorphism (SNP) copy number variations between TK6 and WTK1 cell lines are described herein. Specifically, the integrity of the chromosomes for two WIL2-derived lymphoblastoid cell lines (TK6 and WTK1) was analyzed in the presence and absence of ionizing radiation. WTK1 cells contain a p53 mutation, whereas the TK6 cell line has the native p53 tumor-suppressor gene. Each cell line was isolated post-irradiation for SNP analysis, which showed significant, genome-wide impacts on both cell lines; for the mutant WTK1 sample, there were a total of 48 gene deletions and no gene amplifications, whereas for the wild-type TK6 sample, there were 217 gene deletions and 9 gene amplifications. It appears that both cell lines are affected in the areas of cell-cycle control, but that other affected areas differ significantly between the two.

Transposable element-mediated rearrangements are prevalent in human genomes

Transposable elements constitute about half of human genomes, and their role in generating human variation through retrotransposition is broadly studied and appreciated. Structural variants mediated by transposons, which we call transposable element-mediated rearrangements (TEMRs), are less well studied, and the mechanisms leading to their formation as well as their broader impact on human diversity are poorly understood. Here, we identify 493 unique TEMRs across the genomes of three individuals. While homology directed repair is the dominant driver of TEMRs, our sequence-resolved TEMR resource allows us to identify complex inversion breakpoints, triplications or other high copy number polymorphisms, and additional complexities. TEMRs are enriched in genic loci and can create potentially important risk alleles such as a deletion in TRIM65, a known cancer biomarker and therapeutic target. These findings expand our understanding of this important class of structural variation, the mechanisms responsible for their formation, and establish them as an important driver of human diversity.

Complement component C4 structural variation and quantitative traits contribute to sex-biased vulnerability in systemic sclerosis

Copy number (CN) polymorphisms of complement C4 play distinct roles in many conditions, including immune-mediated diseases. We investigated the association of C4 CN with systemic sclerosis (SSc) risk. Imputed total C4, C4A, C4B, and HERV-K CN were analyzed in 26,633 individuals and validated in an independent cohort. Our results showed that higher C4 CN confers protection to SSc, and deviations from CN parity of C4A and C4B augmented risk. The protection contributed per copy of C4A and C4B differed by sex. Stronger protection was afforded by C4A in men and by C4B in women. C4 CN correlated well with its gene expression and serum protein levels, and less C4 was detected for both in SSc patients. Conditioned analysis suggests that C4 genetics strongly contributes to the SSc association within the major histocompatibility complex locus and highlights classical alleles and amino acid variants of HLA-DRB1 and HLA-DPB1 as C4-independent signals.

Killer-Cell Immunoglobulin-like Receptor Diversity in an Admixed South American Population.

Natural Killer (NK) cells are innate immune cells that mediate antiviral and antitumor responses. NK cell activation and induction of effector functions are tightly regulated by the integration of activating and inhibitory receptors such as killer immunoglobulin-like receptors (KIR). KIR genes are characterized by a high degree of diversity due to presence or absence, gene copy number and allelic polymorphism. The aim of this study was to establish the distribution of KIR genes and genotypes, to infer the most common haplotypes in an admixed Colombian population and to compare these KIR gene frequencies with some Central and South American populations and worldwide. A total of 161 individuals from Medellin, Colombia were included in the study. Genomic DNA was used for KIR and HLA genotyping. We analyzed only KIR gene-content (presence or absence) based on PCR-SSO. The KIR genotype, most common haplotypes and combinations of KIR and HLA ligands frequencies were estimated according to the presence or absence of KIR and HLA genes. Dendrograms, principal component (PC) analysis and Heatmap analysis based on genetic distance were constructed to compare KIR gene frequencies among Central and South American, worldwide and Amerindian populations. The 16 KIR genes analyzed were distributed in 37 different genotypes and the 7 most frequent KIR inferred haplotypes. Importantly, we found three new genotypes not previously reported in any other ethnic group. Our genetic distance, PC and Heatmap analysis revealed marked differences in the distribution of KIR gene frequencies in the Medellin population compared to worldwide populations. These differences occurred mainly in the activating KIR isoforms, which are more frequent in our population, particularly KIR3DS1. Finally, we observed unique structural patterns of genotypes, which evidences the potential diversity and variability of this gene family in our population, and the need for exhaustive genetic studies to expand our understanding of the KIR gene complex in Colombian populations.

Genome-Wide Detection of Copy Number Variations and Evaluation of Candidate Copy Number Polymorphism Genes Associated With Complex Traits of Pigs

Copy number variation (CNV) has been considered to be an important source of genetic variation for important phenotypic traits of livestock. In this study, we performed whole-genome CNV detection on Suhuai (SH) (n = 23), Chinese Min Zhu (MZ) (n = 11), and Large White (LW) (n = 12) pigs based on next-generation sequencing data. The copy number variation regions (CNVRs) were annotated and analyzed, and 10,885, 10,836, and 10,917 CNVRs were detected in LW, MZ, and SH pigs, respectively. Some CNVRs have been randomly selected for verification of the variation type by real-time PCR. We found that SH and LW pigs are closely related, while MZ pigs are distantly related to the SH and LW pigs by CNVR-based genetic structure, PCA, VST, and QTL analyses. A total of 14 known genes annotated in CNVRs were unique for LW pigs. Among them, the cyclin T2 (CCNT2) is involved in cell proliferation and the cell cycle. The FA Complementation Group M (FANCM) is involved in defective DNA repair and reproductive cell development. Ten known genes annotated in 47 CNVRs were unique for MZ pigs. The genes included glycerol-3-phosphate acyltransferase 3 (GPAT3) is involved in fat synthesis and is essential to forming the glycerol triphosphate. Glutathione S-transferase mu 4 (GSTM4) gene plays an important role in detoxification. Eleven known genes annotated in 23 CNVRs were unique for SH pigs. Neuroligin 4 X-linked (NLGN4X) and Neuroligin 4 Y-linked (NLGN4Y) are involved with nerve disorders and nerve signal transmission. IgLON family member 5 (IGLON5) is related to autoimmunity and neural activities. The unique characteristics of LW, MZ, and SH pigs are related to these genes with CNV polymorphisms. These findings provide important information for the identification of candidate genes in the molecular breeding of pigs.

Epitope characterization of a monoclonal antibody that selectively recognizes KIR2DL1 allotypes.

Killer immunoglobulin‐like receptor (KIR) genes code for a family of inhibitory and activating receptors, finely tuning NK cell function. Numerous studies reported the relevance of KIR allelic polymorphism on KIR expression, ligand affinity, and strength in signal transduction. Although KIR variability, including gene copy number and allelic polymorphism, in combination with HLA class I polymorphism, impacts both KIR expression and NK cell education, only a precise phenotypic analysis can define the size of the different KIRpos NK cell subsets. In this context, reagents recognizing a limited number of KIRs is essential. In this study, we have characterized the specificity of an anti‐KIR mAb termed HP‐DM1. Testing its binding to HEK‐293T cells transfected with plasmids coding for different KIRs, we demonstrated that HP‐DM1 mAb exclusively reacts with KIR2DL1. Using site‐directed mutagenesis, we identified the four amino acids relevant for HP‐DM1 recognition: M44, S67, R68, and T70. HP‐DM1 mAb binds to a conformational epitope including M44, the residue crucial for HLA‐C K80 recognition by KIR2DL1. Based on the HP‐DM1 epitope characterization, we could extend its reactivity to all KIR2DL1 allotypes identified except for KIR2DL1*022 and, most likely, KIR2DL1*020, predicting that it does not recognize any other KIR with the only exception of KIR2DS1*013. Moreover, by identifying the residues relevant for HP‐DM1 binding, continuously updating of its reactivity will be facilitated.

Sequence polymorphisms of PR39 cathelicidins and extensive copy variations in commercial pig breeds

Copy number polymorphisms (CNPs) of antimicrobial peptides (AMPs) in livestock can influence the innate immune response of individuals. We conducted a high-resolution analysis of the genomic variations of porcine cathelicidin PR39 using cloned PR39 amplicons corresponding to the 5′ untranslated region (UTR) to 3′ UTR from four individuals of three different pig breeds. We identified 15 different sequences corresponding to 9 different coding domain sequences (CDSs), encoding 7 different protein sequences consisting of 3 functional and 4 non-functional forms. Subsequently, we developed a PR39 CNP typing method using real-time polymerase chain reaction (PCR) and analyzed the PR39 copy numbers from 44 pigs of six breeds. Significant variations in PR39 copies ranging from 2 to 10 copies, with a mean copy number of 5, were observed among all commercial breeds, except the wild boar. Among the different breeds, the PR39 copy number was highest (10) in Korean native pigs. Gene expression analysis showed that PR39 expression was correlated with the copy number. Moreover, the comparative analysis of the cathelicidin cluster-containing region among eight mammalian species showed the complete evolutionary conservation of the region, except for differences in the degree of cathelicidin expansion in each species. Therefore, characterization of CNPs in AMP genes could aid in improving the genetic potential of innate immune responses in livestock animals.

WI12Rhg1 interacts with DELLAs and mediates soybean cyst nematode resistance through hormone pathways.

SummaryThe soybean cyst nematode (SCN) is one of the most important causes of soybean yield loss. The major source of genetic resistance to SCN is the Rhg1 repeat, a tandem copy number polymorphism of three genes. The roles of these genes are only partially understood. Moreover, nematode populations virulent on Rhg1‐carrying soybeans are becoming more common, increasing the need to understand the most successful genetic resistance mechanism. Here, we show that a Rhg1‐locus gene (Glyma.18G02270) encoding a wound‐inducible protein (WI12 Rhg1 ) is needed for SCN resistance. Furthermore, knockout of WI12Rhg1 reduces the expression of DELLA18, and the expression of WI12Rhg1 is itself induced by either JA, SA or GA. The content of the defence hormone SA is significantly lower whilst GA12 and GA53 are increased in WI12Rhg1 knockout roots compared with unedited hairy roots. We find that WI12 Rhg1 directly interacts with DELLA18 (Glyma.18G040000) in yeast and plants and that double knockout of DELLA18 and its homeolog DELLA11 (Glyma.11G216500) significantly reduces SCN resistance and alters the root morphology. As DELLA proteins are implicated in hormone signalling, we explored the content of defence hormones (JA and SA) in DELLA knockout and unedited roots, finding reduced levels of JA and SA after the knockout of DELLA. Additionally, the treatment of DELLA‐knockout roots with JA or SA rescues SCN resistance lost by the knockout. Meanwhile, the SCN resistance of unedited roots decreases after the treatment with GA, but increases with JA or SA. Our findings highlight the critical roles of WI12 Rhg1 and DELLA proteins in SCN resistance through interconnection with hormone signalling.

Etiology of Acute Leukemia: A Review.

Simple SummaryAcute leukemias are some of the most common cancers affecting all age groups. Despite a significant improvement made in the treatment of acute leukemias, their cause remains unknown. A number of genetic and environmental factors for the development of acute leukemias have been proposed, but none have been proven. Undoubtedly, genetics have a major role in the development of these diseases. The effects of a variety of environmental factors, occupations and hobbies have been explored. A recent “two-hit” theory” for the development of acute lymphoblastic leukemia has been proposed. This combines genetic factors and exposure to infections for the development of this disease. Several genetic factors are suggested. Most recently, for the infection portion, exposure to a virus containing Aspergillus Flavus has been proposed. This review summarizes what is currently known about the factors that are proposed for the development of acute leukemias.Acute leukemias constitute some of the most common malignant disorders. Despite significant progress made in the treatment of these disorders, their etiology remains unknown. A large and diverse group of genetic and environmental variables have been proposed. The role of a variety of factors, including pre-existing and acquired genetic mutations, exposure to radiation and various chemicals during preconception, pregnancy and throughout life, have been explored. The effects of inherited genetic variations and disorders, pre-existing diseases, infectious agents, hobbies, occupations, prior treatments, and a host of other factors have been proposed, but none is universally applicable to all cases. Variation in the incidence and prognosis based on the age, sex, race, type of the disease, geographic area of residence and other factors are intriguing but remain unexplained. Advances in genomic profiling, including genome-wide gene expression, DNA copy number and single nucleotide polymorphism (SNP) genotype, may shed some light on the role of genetics in these disparities. Separate two-hit hypotheses for the development of acute myeloblastic and lymphoblastic leukemia have been proposed. The latter combines genetics and infection factors resulting in leukemogenesis. A number of pre- and post-natal environmental conditions and exposure to infections, including a mycovirus infected Aspergillus flavus, have been suggested. The exact nature, timing, sequence of the events and mechanisms resulting in the occurrence of leukemia requires further investigations. This review summarizes some of the above factors in acute lymphoblastic and myeloblastic leukemias and the direction for future research on the etiology of these disorders.

The structure, function and evolution of a complete human chromosome 8

The complete assembly of each human chromosome is essential for understanding human biology and evolution1,2. Here we use complementary long-read sequencing technologies to complete the linear assembly of human chromosome 8. Our assembly resolves the sequence of five previously long-standing gaps, including a 2.08-Mb centromeric α-satellite array, a 644-kb copy number polymorphism in the β-defensin gene cluster that is important for disease risk, and an 863-kb variable number tandem repeat at chromosome 8q21.2 that can function as a neocentromere. We show that the centromeric α-satellite array is generally methylated except for a 73-kb hypomethylated region of diverse higher-order α-satellites enriched with CENP-A nucleosomes, consistent with the location of the kinetochore. In addition, we confirm the overall organization and methylation pattern of the centromere in a diploid human genome. Using a dual long-read sequencing approach, we complete high-quality draft assemblies of the orthologous centromere from chromosome 8 in chimpanzee, orangutan and macaque to reconstruct its evolutionary history. Comparative and phylogenetic analyses show that the higher-order α-satellite structure evolved in the great ape ancestor with a layered symmetry, in which more ancient higher-order repeats locate peripherally to monomeric α-satellites. We estimate that the mutation rate of centromeric satellite DNA is accelerated by more than 2.2-fold compared to the unique portions of the genome, and this acceleration extends into the flanking sequence.

DNA replication, transcription, and H3K56 acetylation regulate copy number and stability at tandem repeats

Tandem repeats are inherently unstable and exhibit extensive copy number polymorphisms. Despite mounting evidence for their adaptive potential, the mechanisms associated with regulation of the stability and copy number of tandem repeats remain largely unclear. To study copy number variation at tandem repeats, we used two well-studied repetitive arrays in the budding yeast genome, the ribosomal DNA (rDNA) locus, and the copper-inducible CUP1 gene array. We developed powerful, highly sensitive, and quantitative assays to measure repeat instability and copy number and used them in multiple high-throughput genetic screens to define pathways involved in regulating copy number variation. These screens revealed that rDNA stability and copy number are regulated by DNA replication, transcription, and histone acetylation. Through parallel studies of both arrays, we demonstrate that instability can be induced by DNA replication stress and transcription. Importantly, while changes in stability in response to stress are observed within a few cell divisions, a change in steady state repeat copy number requires selection over time. Further, H3K56 acetylation is required for regulating transcription and transcription-induced instability at the CUP1 array, and restricts transcription-induced amplification. Our work suggests that the modulation of replication and transcription is a direct, reversible strategy to alter stability at tandem repeats in response to environmental stimuli, which provides cells rapid adaptability through copy number variation. Additionally, histone acetylation may function to promote the normal adaptive program in response to transcriptional stress. Given the omnipresence of DNA replication, transcription, and chromatin marks like histone acetylation, the fundamental mechanisms we have uncovered significantly advance our understanding of the plasticity of tandem repeats more generally.

A genomic amplification affecting a carboxylesterase gene cluster confers organophosphate resistance in the mosquito Aedes aegypti: From genomic characterization to high-throughput field detection.

By altering gene expression and creating paralogs, genomic amplifications represent a key component of short‐term adaptive processes. In insects, the use of insecticides can select gene amplifications causing an increased expression of detoxification enzymes, supporting the usefulness of these DNA markers for monitoring the dynamics of resistance alleles in the field. In this context, the present study aims to characterize a genomic amplification event associated with resistance to organophosphate insecticides in the mosquito Aedes aegypti and to develop a molecular assay to monitor the associated resistance alleles in the field. An experimental evolution experiment using a composite population from Laos supported the association between the over‐transcription of multiple contiguous carboxylesterase genes on chromosome 2 and resistance to multiple organophosphate insecticides. Combining whole genome sequencing and qPCR on specific genes confirmed the presence of a ~100‐Kb amplification spanning at least five carboxylesterase genes at this locus with the co‐existence of multiple structural duplication haplotypes. Field data confirmed their circulation in South‐East Asia and revealed high copy number polymorphism among and within populations suggesting a trade‐off between this resistance mechanism and associated fitness costs. A dual‐color multiplex TaqMan assay allowing the rapid detection and copy number quantification of this amplification event in Ae. aegypti was developed and validated on field populations. The routine use of this novel assay will improve the tracking of resistance alleles in this major arbovirus vector.

The combinatorial diversity of KIR and HLA class I allotypes in Peninsular Malaysia.

Killer cell immunoglobulin-like receptors (KIRs) interact with polymorphic human leucocyte antigen (HLA) class I molecules, modulating natural killer (NK) cell functions and affecting both the susceptibility and outcome of immune-mediated diseases. The KIR locus is highly diverse in gene content, copy number and allelic polymorphism within individuals and across geographical populations. To analyse currently under-represented Asian and Pacific populations, we investigated the combinatorial diversity of KIR and HLA class I in 92 unrelated Malay and 75 Malaysian Chinese individuals from the Malay Peninsula. We identified substantial allelic and structural diversity of the KIR locus in both populations and characterized novel variations at each analysis level. The Malay population is more diverse than Malay Chinese, likely representing a unique history including admixture with immigrating populations spanning several thousand years. Characterizing the Malay population are KIR haplotypes with large structural variants present in 10% individuals, and KIR and HLA alleles previously identified in Austronesian populations. Despite the differences in ancestries, the proportion of HLA allotypes that serve as KIR ligands is similar in each population. The exception is a significantly reduced frequency of interactions of KIR2DL1 with C2+ HLA-C in the Malaysian Chinese group, caused by the low frequency of C2+ HLA. One likely implication is a greater protection from preeclampsia, a pregnancy disorder associated with KIR2DL1, which shows higher incidence in the Malay than in the Malaysian Chinese. This first complete, high-resolution, characterization of combinatorial diversity of KIR and HLA in Malaysians will form a valuable reference for future clinical and population studies.

The Multilocus Multispecies Coalescent: A Flexible New Model of Gene Family Evolution.

Incomplete lineage sorting (ILS), the interaction between coalescence and speciation, can generate incongruence between gene trees and species trees, as can gene duplication (D), transfer (T), and loss (L). These processes are usually modeled independently, but in reality, ILS can affect gene copy number polymorphism, that is, interfere with DTL. This has been previously recognized, but not treated in a satisfactory way, mainly because DTL events are naturally modeled forward-in-time, while ILS is naturally modeled backward-in-time with the coalescent. Here, we consider the joint action of ILS and DTL on the gene tree/species tree problem in all its complexity. In particular, we show that the interaction between ILS and duplications/transfers (without losses) can result in patterns usually interpreted as resulting from gene loss, and that the realized rate of D, T, and L becomes nonhomogeneous in time when ILS is taken into account. We introduce algorithmic solutions to these problems. Our new model, the multilocus multispecies coalescent, which also accounts for any level of linkage between loci, generalizes the multispecies coalescent (MSC) model and offers a versatile, powerful framework for proper simulation, and inference of gene family evolution. [Gene duplication; gene loss; horizontal gene transfer; incomplete lineage sorting; multispecies coalescent; hemiplasy; recombination.].