Coding Sequence Variants Research Articles

Within species expressed genetic variability and gene expression response to different temperatures in the rotifer Brachionus calyciflorus sensu stricto

Genetic divergence is impacted by many factors, including phylogenetic history, gene flow, genetic drift, and divergent selection. Rotifers are an important component of aquatic ecosystems, and genetic variation is essential to their ongoing adaptive diversification and local adaptation. In addition to coding sequence divergence, variation in gene expression may relate to variable heat tolerance, and can impose ecological barriers within species. Temperature plays a significant role in aquatic ecosystems by affecting species abundance, spatio-temporal distribution, and habitat colonization. Recently described (formerly cryptic) species of the Brachionus calyciflorus complex exhibit different temperature tolerance both in natural and in laboratory studies, and show that B. calyciflorus sensu stricto (s.s.) is a thermotolerant species. Even within B. calyciflorus s.s., there is a tendency for further temperature specializations. Comparison of expressed genes allows us to assess the impact of stressors on both expression and sequence divergence among disparate populations within a single species. Here, we have used RNA-seq to explore expressed genetic diversity in B. calyciflorus s.s. in two mitochondrial DNA lineages with different phylogenetic histories and differences in thermotolerance. We identify a suite of candidate genes that may underlie local adaptation, with a particular focus on the response to sustained high or low temperatures. We do not find adaptive divergence in established candidate genes for thermal adaptation. Rather, we detect divergent selection among our two lineages in genes related to metabolism (lipid metabolism, metabolism of xenobiotics).

Variation in the BrHMA3 coding region controls natural variation in cadmium accumulation in Brassica rapa vegetables.

Brassica rapa includes several important leafy vegetable crops with the potential for high cadmium (Cd) accumulation, posing a risk to human health. This study aims to understand the genetic basis underlying the variation in Cd accumulation among B. rapa vegetables. Cd uptake and translocation in 64 B. rapa accessions were compared. The role of the heavy metal ATPase gene BrHMA3 in the variation of Cd accumulation was investigated. BrHMA3 encodes a tonoplast-localized Cd transporter. Five full-length and four truncated haplotypes of the BrHMA3 coding sequence were identified, explaining >80% of the variation in the Cd root to shoot translocation among the 64 accessions and in F2 progeny. Truncated BrHMA3 haplotypes had a 2.3 and 9.3 times higher shoot Cd concentration and Cd translocation ratio, respectively, than full-length haplotypes. When expressed in yeast and Arabidopsis thaliana, full-length BrHMA3 showed activity consistent with a Cd transport function, whereas truncated BrHMA3 did not. Variation in the BrHMA3 promoter sequence had little effect on Cd translocation. Variation in the BrHMA3 coding sequence is a key determinant of Cd translocation to and accumulation in the leaves of B. rapa. Strong alleles of BrHMA3 can be used to breed for B. rapa vegetables that are low in Cd in their edible portions.

Role of SCN5A coding and non-coding sequences in Brugada syndrome onset: What's behind the scenes?

BackgroundBrugada syndrome (BrS) is a rare inherited cardiac arrhythmia associated with a high risk of sudden cardiac death (SCD) due to ventricular fibrillation (VF). BrS is characterized by coved-type ST-segment elevation in the right precordial leads (V1-V3). Mutations in SCN5A gene coding for the α-subunit of the NaV1.5 cardiac sodium channel are identified in 15–30% of BrS cases. Genetic testing of BrS patients generally involves sequencing of the protein-coding portions and flanking intronic regions of SCN5A. This excludes the 5′UTR and 3′UTR from the routine genetic testing. MethodsWe here screened the coding sequence, the flanking intronic regions as well as the 5′ and 3′UTR regions of SCN5A gene and further five candidate genes (GPD1L, SCN1B, KCNE3, SCN4B, and MOG1) in a Tunisian family diagnosed with BrS. ResultsA new SCN5A-Q1000K mutation was identified along with two common polymorphisms (H558R and D1819). Multiple genetic variants were identified on the SCN5A 3′UTR, one of which is predicted to create additional microRNA binding site for miR-1270. Additionally, we identified the hsa-miR-219a-rs107822. No relevant coding sequence variant was identified in the remaining studied candidate genes. ConclusionsThe absence of genotype-phenotype concordance within all the identified genetic variants in this family gives extra evidences about the complexity of the disease and suggests that the occurrence and prognosis of BrS is most likely controlled by a combination of multiple genetic factors, rather than a single variant. Most SCN5A variants were localized in non-coding regions hypothesizing an impact on the miRNA-target complementarities.

Common Variants in the Glycerol Kinase Gene Reduce Tuberculosis Drug Efficacy.

Despite the administration of multiple drugs that are highly effective in vitro, tuberculosis (TB) treatment requires prolonged drug administration and is confounded by the emergence of drug-resistant strains. To understand the mechanisms that limit antibiotic efficacy, we performed a comprehensive genetic study to identify Mycobacterium tuberculosis genes that alter the rate of bacterial clearance in drug-treated mice. Several functionally distinct bacterial genes were found to alter bacterial clearance, and prominent among these was the glpK gene that encodes the glycerol-3-kinase enzyme that is necessary for glycerol catabolism. Growth on glycerol generally increased the sensitivity of M. tuberculosis to antibiotics in vitro, and glpK-deficient bacteria persisted during antibiotic treatment in vivo, particularly during exposure to pyrazinamide-containing regimens. Frameshift mutations in a hypervariable homopolymeric region of the glpK gene were found to be a specific marker of multidrug resistance in clinical M. tuberculosis isolates, and these loss-of-function alleles were also enriched in extensively drug-resistant clones. These data indicate that frequently observed variation in the glpK coding sequence produces a drug-tolerant phenotype that can reduce antibiotic efficacy and may contribute to the evolution of resistance.IMPORTANCE TB control is limited in part by the length of antibiotic treatment needed to prevent recurrent disease. To probe mechanisms underlying survival under antibiotic pressure, we performed a genetic screen for M. tuberculosis mutants with altered susceptibility to treatment using the mouse model of TB. We identified multiple genes involved in a range of functions which alter sensitivity to antibiotics. In particular, we found glycerol catabolism mutants were less susceptible to treatment and that common variation in a homopolymeric region in the glpK gene was associated with drug resistance in clinical isolates. These studies indicate that reversible high-frequency variation in carbon metabolic pathways can produce phenotypically drug-tolerant clones and have a role in the development of resistance.

QTL analysis of natural Saccharomyces cerevisiae isolates reveals unique alleles involved in lignocellulosic inhibitor tolerance.

Decoding the genetic basis of lignocellulosic inhibitor tolerance in Saccharomyces cerevisiae is crucial for rational engineering of bioethanol strains with enhanced robustness. The genetic diversity of natural strains present an invaluable resource for the exploration of complex traits of industrial importance from a pan-genomic perspective to complement the limited range of specialised, tolerant industrial strains. Natural S. cerevisiae isolates have lately garnered interest as a promising toolbox for engineering novel, genetically encoded tolerance phenotypes into commercial strains. To this end, we investigated the genetic basis for lignocellulosic inhibitor tolerance of natural S. cerevisiae isolates. A total of 12 quantitative trait loci underpinning tolerance were identified by next-generation sequencing linked bulk-segregant analysis of superior interbred pools. Our findings corroborate the current perspective of lignocellulosic inhibitor tolerance as a multigenic, complex trait. Apart from a core set of genetic variants required for inhibitor tolerance, an additional genetic background-specific response was observed. Functional analyses of the identified genetic loci revealed the uncharacterised ORF, YGL176C and the bud-site selection XRN1/BUD13 as potentially beneficial alleles contributing to tolerance to a complex lignocellulosic inhibitor mixture. We present evidence for the consideration of both regulatory and coding sequence variants for strain improvement.

Damaging Variants in Proangiogenic Genes Impair Growth in Fetuses with Cardiac Defects

To determine the impact of damaging genetic variation in proangiogenic pathways on placental function, complications of pregnancy, fetal growth, and clinical outcomes in pregnancies with fetal congenital heart defect. Families delivering a baby with a congenital heart defect requiring surgical repair in infancy were recruited. The placenta and neonate were weighed and measured. Hemodynamic variables were recorded from a third trimester (36.4±1.7weeks) fetal echocardiogram. Exome sequencing was performed on the probands (N=133) and consented parents (114 parent-child trios, and 15 parent-child duos) and the GeneVetter analysis tool used to identify damaging coding sequence variants in 163 genes associated with the positive regulation of angiogenesis (PRA) (GO:0045766). In total, 117 damaging variants were identified in PRA genes in 133 congenital heart defect probands with 73 subjects having at least 1 variant. Presence of a damaging PRA variant was associated with increased umbilical artery pulsatility index (mean 1.11 with variant vs 1.00 without; P=.01). The presence of a damaging PRA variant was also associated with lower neonatal length and head circumference for age z score at birth (mean -0.44 and -0.47 with variant vs 0.23 and -0.05 without; P=.01 and .04, respectively). During median 3.1years (IQR 2.0-4.1years) of follow-up, deaths occurred in 2 of 60 (3.3%) subjects with no PRA variant and in 9 of 73 (12.3%) subjects with 1 or more PRA variants (P=.06). Damaging variants in proangiogenic genes may impact placental function and are associated with impaired fetal growth in pregnancies involving a fetus with congenital heart defect.

Calcyphosine-like (CAPSL) is regulated in Multiple Symmetric Lipomatosis and is involved in Adipogenesis

Little is known on the causes and pathogenesis of the adipose tissue disorder (familial) Multiple Symmetric Lipomatosis (MSL). In a four-generation MSL-family, we performed whole exome sequencing (WES) in 3 affected individuals and 1 obligate carrier and identified Calcyphosine-like (CAPSL) as the most promising candidate gene for this family. Screening of 21 independent patients excluded CAPSL coding sequence variants as a common monogenic cause, but using immunohistochemistry we found that CAPSL was down-regulated in adipose tissue not only from the index patient but also in 10 independent sporadic MSL-patients. This suggests that CAPSL is regulated in sporadic MSL irrespective of the underlying genetic/multifactorial cause. Furthermore, we cultivated pre-adipocytes from MSL-patients and generated 3T3-L1-based Capsl knockout and overexpressing cell models showing altered autophagy, adipogenesis, lipogenesis and Sirtuin-1 (SIRT1) expression. CAPSL seems to be involved in adipocyte biology and perturbation of autophagy is a potential mechanism in the pathogenesis of MSL. Downregulation of CAPSL and upregulation of UCP1 were common features in MSL fat while the known MSL genes MFN2 and LIPE did not show consistent alterations. CAPSL immunostainings could serve as first diagnostic tools in MSL clinical care with a potential to improve time to diagnosis and healthcare options.

Molecular cloning of WIF1 and HMGA2 reveals ear-preferential expression while uncovering a missense mutation associated with porcine ear size in WIF1.

Considerable diversity exists in porcine ear size, which is an important morphological feature of pig breeds. Previously, we localized four crucial candidate genes-high mobility group AT-hook 2 (HMGA2), LEM domain-containing 3 (LEMD3), methionine sulfoxide reductase B3 (MSRB3) and Wnt inhibitory factor 1 (WIF1)-on Sus Scrofa chromosome 5 affecting porcine ear size, then cloned LEMD3 and MSBR3. In this study, we performed rapid amplification of cDNA ends to obtain full-length cDNA sequences of 2338-bp WIF1 and 2998-bp HMGA2. Using quantitative real-time PCR, we revealed that WIF1 expression was highest in ear cartilage of 60-day-old pigs and that this is therefore a better candidate gene for ear size than HMGA2. We further screened coding sequence variants in both genes and identified only one missense mutation (WIF1:c.1167C>G) in a conserved epidermal growth factor-like domain from the mammalian WIF1 protein. The protein-altering mutation was significantly associated with ear size across the Large White×Minzhu hybrid and Beijing Black pig populations. When WIF1:c.1167C>G was included as fixed effect in the model to re-run a genome-wide association study in the Large White×Minzhu intercross population the P-value of the peak SNP on SSC5 from re-running the genome-wide association study dropped from 2.45E-12 to 7.33E-05. Taken together, the WIF1:c.1167C>G could be an important mutation associated with ear size. Our findings provide helpful information for further studies of the molecular mechanisms controlling porcine ear size.

The Oculome Panel Test: Next-Generation Sequencing to Diagnose a Diverse Range of Genetic Developmental Eye Disorders

To develop a comprehensive next-generation sequencing panel assay that screens genes known to cause developmental eye disorders and inherited eye disease and to evaluate its diagnostic yield in a pediatric cohort with malformations of the globe, anterior segment anomalies, childhood glaucoma, or a combination thereof. Evaluation of diagnostic test. Two hundred seventy-seven children, 0 to 16 years of age, diagnosed with nonsyndromic or syndromic developmental eye defects without a genetic diagnosis. We developed a new oculome panel using a custom-designed Agilent SureSelect QXT target capture method (Agilent Technologies, Santa Clara, CA) to capture and perform parallel high-throughput sequencing analysis of 429 genes associated with eye disorders. Bidirectional Sanger sequencing confirmed suspected pathogenic variants. Collated clinical details and oculome molecular genetic results. The oculome design covers 429 known eye disease genes; these are subdivided into 5 overlapping virtual subpanels for anterior segment developmental anomalies including glaucoma (ASDA; 59 genes), microphthalmia-anophthalmia-coloboma (MAC; 86 genes), congenital cataracts and lens-associated conditions (70 genes), retinal dystrophies (RET; 235 genes), and albinism (15 genes), as well as additional genes implicated in optic atrophy and complex strabismus (10 genes). Panel development and testing included analyzing 277 clinical samples and 3 positive control samples using Illumina sequencing platforms; more than30× read depth was achieved for 99.5% of the targeted 1.77-Mb region. Bioinformatics analysis performed using a pipeline based on Freebayes and ExomeDepth to identify coding sequence and copy number variants, respectively, resulted in a definitive diagnosis in 68 of 277 samples, with variability in diagnostic yield between phenotypic subgroups: MAC, 8.2% (8 of 98 cases solved); ASDA, 24.8% (28 of 113 cases solved); other or syndromic, 37.5% (3 of 8 cases solved); RET, 42.8% (21 of 49 cases solved); and congenital cataracts and lens-associated conditions, 88.9% (8of 9 cases solved). The oculome test diagnoses a comprehensive range of genetic conditions affecting the development of the eye, potentially replacing protracted and costly multidisciplinary assessments and allowing for faster targeted management. The oculome enabled molecular diagnosis of a significant number of cases in our sample cohort of varied ocular birth defects.

Mutations in EPAS1 in congenital heart disease in Tibetans.

EPAS1 encodes HIF2 and is closely related to high altitude chronic hypoxia. Mutations in the EPAS1 coding sequence are associated with several kinds of human diseases, including syndromic congenital heart disease (CHD). However, whether there are rare EPAS1 coding variants related to Tibetan non-syndromic CHD have not been fully investigated. A group of 286 Tibetan patients with non-syndromic CHD and 250 unrelated Tibetan healthy controls were recruited from Qinghai, China. Sanger sequencing was performed to identify variations in the EPAS1 coding sequence. The novelty of identified variants was confirmed by the examination of 1000G and ExAC databases. Control samples were screened to establish that the rare candidate variants were specific to the Tibetan patients with non-syndromic CHD. Bioinformatics software was used to assess the conservation of the mutations and to predict their effects. The effect of EPAS1 mutations on the transcription of its target gene, VEGF, was assessed by dual-luciferase reporter assay. The mammalian two-hybrid assay was used to study the protein interactions between HIF2 and PHD2 or pVHL. We identified two novel EPAS1 mutations (NM_001430: c.607A>C, p.N203H; c.2170G>T, p.G724W) in two patients. The N203H mutation significantly affected the transcription activity of the VEGF promoter, especially in conditions of hypoxia. The N203H mutation also showed enhanced protein–protein interactions between HIF2 and PHD2, and HIF2 and pVHL, especially in conditions of hypoxia. However, the G724W mutation did not demonstrate the same effects. Our results indicate that EPAS1 mutations might have a potential causative effect on the development of Tibetan non-syndromic CHD.

Genetic and transcriptional variations in NRAMP-2 and OPAQUE1 genes are associated with salt stress response in wheat

Key messageSNP alleles on chromosomes 4BL and 6AL are associated with sensitivity to salt tolerance in wheat and upon validation can be exploited in the development of salt-tolerant wheat varieties.The dissection of the genetic and molecular components of salt stress response offers strong opportunities toward understanding and improving salt tolerance in crops. In this study, GWAS was employed to identify a total of 106 SNP loci (R2 = 0.12–63.44%) linked to salt stress response in wheat using leaf chlorophyll fluorescence, grain quality and shoot ionic (Na+ and K+ ions) attributes. Among them, 14 SNP loci individually conferred pleiotropic effects on multiple independent salinity tolerance traits including loci at 99.04 cM (R2 ≥ 14.7%) and 68.45 cM (R2 ≥ 4.10%) on chromosomes 6AL and 4BL, respectively, that influenced shoot Na+-uptake, shoot K+/Na+ ratio, and specific energy fluxes for absorption (ABS/RC) and dissipation (DIo/RC). Analysis of the open reading frame (ORF) containing the SNP markers revealed that they are orthologous to genes involved in photosynthesis and plant stress (salt) response. Further transcript abundance and qRT-PCR analyses indicated that the genes are mostly up-regulated in salt-tolerant and down-regulated in salt-sensitive wheat genotypes including NRAMP-2 and OPAQUE1 genes on 4BL and 6AL, respectively. Both genes showed highest differential expression between contrasting genotypes when expressions of all the genes within their genetic intervals were analyzed. Possible cis-acting regulatory elements and coding sequence variation that may be involved in salt stress response were also identified in both genes. This study identified genetic and molecular components of salt stress response that are associated with Na+-uptake, shoot Na+/K+ ratio, ABS/RC, DIo/RC, and grain quality traits and upon functional validation would facilitate the development of gene-specific markers that could be deployed to improve salinity tolerance in wheat.

Low-Frequency and Rare-Coding Variation Contributes to Multiple Sclerosis Risk

SummaryMultiple sclerosis is a complex neurological disease, with ∼20% of risk heritability attributable to common genetic variants, including >230 identified by genome-wide association studies. Multiple strands of evidence suggest that much of the remaining heritability is also due to additive effects of common variants rather than epistasis between these variants or mutations exclusive to individual families. Here, we show in 68,379 cases and controls that up to 5% of this heritability is explained by low-frequency variation in gene coding sequence. We identify four novel genes driving MS risk independently of common-variant signals, highlighting key pathogenic roles for regulatory T cell homeostasis and regulation, IFNγ biology, and NFκB signaling. As low-frequency variants do not show substantial linkage disequilibrium with other variants, and as coding variants are more interpretable and experimentally tractable than non-coding variation, our discoveries constitute a rich resource for dissecting the pathobiology of MS.

5\u2019-flanking variants of equine casein genes (CSN1S1, CSN1S2, CSN2, CSN3) and their relationship with gene expression and milk composition

Genes encoding casein proteins are important candidates for milk composition traits in mammals. In the case of the domestic horse, our knowledge of casein genes is limited mainly to coding sequence variants. This study involved screening for polymorphism in 5’-flanking regions of four genes encoding equine caseins (CSN1S1, CSN1S2, CSN2, and CSN3) and making a preliminary assessment of their effect on the gene expression (on the mRNA and protein levels) and milk composition traits in selected horse breeds. Altogether, 23 polymorphisms (21 described previously SNPs and two novel InDels) were found in the studied sequences, the majority of which are common in various horse breeds. Statistical analysis revealed that some are putatively associated with gene expression or milk composition — for example, the c.-2047_-2048insAT polymorphism (CSN1S1) turns out to be related to the total milk protein content in Polish Primitive Horse (p < 0.05), whereas c.-2105C>G SNP (CSN2) is related to beta-casein relative mRNA level and milk lactose concentration in the Polish Coldblood Horse breed (p < 0.05). We have also found significant effects of horse breed and lactation time-point on gene expression and mare’s milk composition. Our study indicates that the 5’-regulatory regions of genes encoding casein proteins are interesting targets for functional studies of their expression and the composition traits of mare’s milk.

Positive selection in coding regions and motif duplication in regulatory regions of bottlenose dolphin MHC class II genes.

The vertebrate immune response is mediated through highly adaptive, quickly evolving cell surface receptors, the major histocompatibility complex (MHC). MHC molecules bind and present a diverse array of pathogenic molecules and trigger a cascade of defenses. Use of MHC variation as a marker for population health has also evolved quickly following advances in sequencing methods. We applied a combination of traditional and next generation sequencing methodology to characterize coding (peptide binding region) and regulatory (proximal promoter) sequence variation in MHC Class II DQA and DQB genes between estuarine and coastal populations of the bottlenose dolphin, Tursiops truncatus, an apex predator whose health status is indicative of anthropogenic impacts on the ecosystem. The coding regions had 10 alleles each at DQA and DQB; the promoters had 6 and 7 alleles at DQA and DQB, respectively with variation within key regulatory motifs. Positive selection was observed for the coding regions of both genes while both coding and promoter regions exhibited geographic differences in allele composition that likely indicates diversifying selection across habitats. Most notable was the discovery of a complete duplication of a 14-bp T-box motif in the DQA promoter. Four class II promoter regions (DQA, DQB, DRA, DRB) were characterized in species from four cetacean families (Delphinidae, Monodontidae, Lipotidae, and Physeteridae) and revealed substantial promoter structural diversity across this order. Peptide binding regions may not be the only source of adaptive potential within cetacean MHC for responding to pathogenic threats. These findings are the first analysis of cetacean MHC regulatory motifs, which may divulge unique immunogenetic strategies among cetaceans and reveal how MHC transcriptional control continues to evolve. The combined MHC regulatory and coding data provide new genetic context for distinct vulnerability profiles between coastal and estuarine populations, which are key concerns for health and risk management.

APOL1 risk variants cause podocytes injury through enhancing endoplasmic reticulum stress.

Two coding sequence variants (G1 and G2) of Apolipoprotein L1 (APOL1) gene have been implicated as a higher risk factor for chronic kidney diseases (CKD) in African Americans when compared with European Americans. Previous studies have suggested that the APOL1 G1 and G2 variant proteins are more toxic to kidney cells than the wild-type APOL1 G0, but the underlying mechanisms are poorly understood. To determine whether endoplasmic reticulum (ER) stress contributes to podocyte toxicity, we generated human podocytes (HPs) that stably overexpressed APOL1 G0, G1, or G2 (Vec/HPs, G0/HPs, G1/HPs, and G2/HPs). Propidium iodide staining showed that HP overexpressing the APOL1 G1 or G2 variant exhibited a higher rate of necrosis when compared with those overexpressing the wild-type G0 counterpart. Consistently, the expression levels of nephrin and podocin proteins were significantly decreased in the G1- or G2-overexpressing cells despite the maintenance of their mRNA expressions levels. In contrast, the expression of the 78-kDa glucose-regulated protein ((GRP78), also known as the binding Ig protein, BiP) and the phosphorylation of the eukaryotic translation initiation factor 1 (eIF1) were significantly elevated in the G1/HPs and G2/HPs, suggesting a possible occurrence of ER stress in these cells. Furthermore, ER stress inhibitors not only restored nephrin protein expression, but also provided protection against necrosis in G1/HPs and G2/HPs, suggesting that APOL1 risk variants cause podocyte injury partly through enhancing ER stress.

The transcriptional landscape of polyploid wheat.

The coordinated expression of highly related homoeologous genes in polyploid species underlies the phenotypes of many of the world's major crops. Here we combine extensive gene expression datasets to produce a comprehensive, genome-wide analysis of homoeolog expression patterns in hexaploid bread wheat. Bias in homoeolog expression varies between tissues, with ~30% of wheat homoeologs showing nonbalanced expression. We found expression asymmetries along wheat chromosomes, with homoeologs showing the largest inter-tissue, inter-cultivar, and coding sequence variation, most often located in high-recombination distal ends of chromosomes. These transcriptionally dynamic genes potentially represent the first steps toward neo- or subfunctionalization of wheat homoeologs. Coexpression networks reveal extensive coordination of homoeologs throughout development and, alongside a detailed expression atlas, provide a framework to target candidate genes underpinning agronomic traits in wheat.

Noncoding RET variants explain the strong association with Hirschsprung disease in patients without rare coding sequence variant

The pathogenesis of Hirschsprung disease is complex. Although the RET proto-oncogene is the most frequently affected gene in Hirschsprung disease, rare coding sequence variants explain only a small part of Hirschsprung disease cases. We aimed to assess the genetic background of Hirschsprung disease using a genome-wide association analysis combined with sequencing all RET exons in samples from 105 Hirschsprung disease cases (30 familial and 75 sporadic) and 386 controls.As expected, variants in or near RET showed the strongest overall association with Hirschsprung disease and the most statistically significant association was observed when using a recessive genetic model (rs2435357, NC_000010.10:g.43582056T > C; genotype TT, OR = 17.31, P = 1.462 × 10−21). Previously published associations in variants in SEMA (rs11766001, NC_000007.13:g.84145202A > C; allele C, OR = 2.268, P = 0.009533) and NRG1 (rs4541858, NC_000008.10:g.32410309A > G; allele G, OR = 1.567, P = 0.015; rs7835688, NC_000008.10:g.32411499G > C; allele C, OR = 1.567, P = 0.015) were also replicated in the genome-wide association analysis. Sequencing revealed a total of 12 exonic RET rare variants. Of these, eight amino acid changing rare variants and two frameshift variants caused or possibly caused Hirschsprung disease.Only a minority of the Hirschsprung disease cases (9/30 familial; 7/75 sporadic) carried one of the rare variants. Excluding the rare variant carriers from the genome-wide association analysis did not appreciably change the association of rs2435357 with Hirschsprung disease. We estimate that approximately two thirds of the sporadic cases may be statistically attributed to the recessive action of the common non-coding RET variants. Thus, even though most cases do not carry rare RET variants, combinations of rare variants and the common non-coding RET variant cause the majority of the cases in our population.

Inter and Intraspecific Genomic Divergence in Drosophila montana Shows Evidence for Cold Adaptation.

The genomes of species that are ecological specialists will likely contain signatures of genomic adaptation to their niche. However, distinguishing genes related to ecological specialism from other sources of selection and more random changes is a challenge. Here, we describe the genome of Drosophila montana, which is the most extremely cold-adapted Drosophila species known. We use branch tests to identify genes showing accelerated divergence in contrasts between cold- and warm-adapted species and identify about 250 genes that show differences, possibly driven by a lower synonymous substitution rate in cold-adapted species. We also look for evidence of accelerated divergence between D. montana and D. virilis, a previously sequenced relative, but do not find strong evidence for divergent selection on coding sequence variation. Divergent genes are involved in a variety of functions, including cuticular and olfactory processes. Finally, we also resequenced three populations of D. montana from across its ecological and geographic range. Outlier loci were more likely to be found on the X chromosome and there was a greater than expected overlap between population outliers and those genes implicated in cold adaptation between Drosophila species, implying some continuity of selective process at these different evolutionary scales.

Multiple massive domestication and recent amplification of Kolobok superfamily transposons in the clawed frog Xenopus

BackgroundDNA transposons are generally destroyed by mutations and have short lifespans in hosts, as they are neutral or harmful to the host and therefore not conserved by natural selection. The clawed frog Xenopus harbors many DNA transposons and certain families, such as T2-MITE, have extremely long lives. These have ancient origins, but have shown recent transposition activity. In addition, certain transposase genes may have been “domesticated” by Xenopus and conserved over long time periods by natural selection. The aim of this study was to elucidate the evolutionary interactions between the host and the long-lived DNA transposon family it contains. Here, we investigated the molecular evolution of the Kolobok DNA transposon superfamily. Kolobok is thought to contribute to T2-MITE transposition.ResultsIn the diploid western clawed frog Xenopus tropicalis and the allotetraploid African clawed frog Xenopus laevis, we searched for transposase genes homologous to those in the Kolobok superfamily. To determine the amplification and domestication of these genes, we used molecular phylogenetics and analyses of copy numbers, conserved motifs, orthologous gene synteny, and coding sequence divergence between the orthologs of X. laevis and X. tropicalis, or between those of two distant X. tropicalis lineages. Among 38 X. tropicalis and 24 X. laevis prospective transposase genes, 10 or more in X. tropicalis and 14 or more in X. laevis were apparently domesticated. These genes may have undergone multiple independent domestications from before the divergence of X. laevis and X. tropicalis. In contrast, certain other transposases may have retained catalytic activity required for transposition and could therefore have been recently amplified.ConclusionMultiple domestication of certain transposases and prolonged conservation of the catalytic activity in others suggest that Kolobok superfamily transposons were involved in complex, mutually beneficial relationships with their Xenopus hosts. Some transposases may serve to activate long-lived T2-MITE subfamilies.

TMTC2 variant associated with sensorineural hearing loss and auditory neuropathy spectrum disorder in a family dyad.

BackgroundSensorineural hearing loss (SNHL) is a common form of hearing loss that can be inherited or triggered by environmental insults; auditory neuropathy spectrum disorder (ANSD) is a SNHL subtype with unique diagnostic criteria. The genetic factors associated with these impairments are vast and diverse, but causal genetic factors are rarely characterized.MethodsA family dyad, both cochlear implant recipients, presented with a hearing history of bilateral, progressive SNHL, and ANSD. Whole‐exome sequencing was performed to identify coding sequence variants shared by both family members, and screened against genes relevant to hearing loss and variants known to be associated with SNHL and ANSD.ResultsBoth family members are successful cochlear implant users, demonstrating effective auditory nerve stimulation with their devices. Genetic analyses revealed a mutation (rs35725509) in the TMTC2 gene, which has been reported previously as a likely genetic cause of SNHL in another family of Northern European descent.ConclusionThis study represents the first confirmation of the rs35725509 variant in an independent family as a likely cause for the complex hearing loss phenotype (SNHL and ANSD) observed in this family dyad.