Chromatin three-dimensional interactions mediate genetic effects on gene expression.

Studying the genetic basis of gene expression and chromatin organization is key to characterizing the effect of genetic variability on the function and structure of the human genome. Here we unravel how genetic variation perturbs gene regulation using a dataset combining activity of regulatory elements, gene expression, and genetic variants across 317 individuals and two cell types. We show that variability in regulatory activity is structured at the intra- and interchromosomal levels within 12,583 cis-regulatory domains and 30 trans-regulatory hubs that highly reflect the local (that is, topologically associating domains) and global (that is, open and closed chromatin compartments) nuclear chromatin organization. These structures delimit cell type-specific regulatory networks that control gene expression and coexpression and mediate the genetic effects of cis- and trans-acting regulatory variants on genes.

MRPL53, a New Candidate Gene for Orofacial Clefting, Identified Using an eQTL Approach.

A valuable approach to understand how individual and population genetic differences can predispose to disease is to assess the impact of genetic variants on cellular functions (e.g., gene expression) of cell and tissue types related to pathological states. To understand the genetic basis of nonsyndromic cleft lip with or without cleft palate (NSCL/P) susceptibility, a complex and highly prevalent congenital malformation, we searched for genetic variants with a regulatory role in a disease-related tissue, the lip muscle (orbicularis oris muscle [OOM]), of affected individuals. From 46 OOM samples, which are frequently discarded during routine corrective surgeries on patients with orofacial clefts, we derived mesenchymal stem cells and correlated the individual genetic variants with gene expression from these cultured cells. Through this strategy, we detected significant cis-eQTLs (i.e., DNA variants affecting gene expression) and selected a few candidates to conduct an association study in a large Brazilian cohort (624 patients and 668 controls). This resulted in the discovery of a novel susceptibility locus for NSCL/P, rs1063588, the best eQTL for the MRPL53 gene, where evidence for association was mostly driven by the Native American ancestry component of our Brazilian sample. MRPL53 (2p13.1) encodes a 39S protein subunit of mitochondrial ribosomes and interacts with MYC, a transcription factor required for normal facial morphogenesis. Our study illustrates not only the importance of sampling admixed populations but also the relevance of measuring the functional effects of genetic variants over gene expression to dissect the complexity of disease phenotypes.

Chromosomal contacts connect loci associated with autism, BMI and head circumference phenotypes.

Copy number variants (CNVs) are major contributors to genomic imbalance disorders. Phenotyping of 137 unrelated deletion and reciprocal duplication carriers of the distal 16p11.2 220 kb BP2-BP3 interval showed that these rearrangements are associated with autism spectrum disorders and mirror phenotypes of obesity/underweight and macrocephaly/microcephaly. Such phenotypes were previously associated with rearrangements of the non-overlapping proximal 16p11.2 600 kb BP4-BP5 interval. These two CNV-prone regions at 16p11.2 are reciprocally engaged in complex chromatin looping, as successfully confirmed by 4C-seq, fluorescence in situ hybridization and Hi-C, as well as coordinated expression and regulation of encompassed genes. We observed that genes differentially expressed in 16p11.2 BP4-BP5 CNV carriers are concomitantly modified in their chromatin interactions, suggesting that disruption of chromatin interplays could participate in the observed phenotypes. We also identified cis- and trans-acting chromatin contacts to other genomic regions previously associated with analogous phenotypes. For example, we uncovered that individuals with reciprocal rearrangements of the trans-contacted 2p15 locus similarly display mirror phenotypes on head circumference and weight. Our results indicate that chromosomal contacts' maps could uncover functionally and clinically related genes.

A functional circadian clock is required for proper insulin secretion by human pancreatic islet cells.

AIM: To determine the impact of a functional human islet clock on insulin secretion and gene transcription. METHODS: Efficient circadian clock disruption was achieved in human pancreatic islet cells by small interfering RNA-mediated knockdown of CLOCK. Human islet secretory function was assessed in the presence or absence of a functional circadian clock by stimulated insulin secretion assays, and by continuous around-the-clock monitoring of basal insulin secretion. Large-scale transcription analysis was accomplished by RNA sequencing, followed by quantitative RT-PCR analysis of selected targets. RESULTS: Circadian clock disruption resulted in a significant decrease in both acute and chronic glucose-stimulated insulin secretion. Moreover, basal insulin secretion by human islet cells synchronized in vitro exhibited a circadian pattern, which was perturbed upon clock disruption. RNA sequencing analysis suggested alterations in 352 transcript levels upon circadian clock disruption. Among them, key regulators of the insulin secretion pathway (GNAQ, ATP1A1, ATP5G2, KCNJ11) and transcripts required for granule maturation and release (VAMP3, STX6, SLC30A8) were affected. CONCLUSIONS: Using our newly developed experimental approach for efficient clock disruption in human pancreatic islet cells, we show for the first time that a functional ß-cell clock is required for proper basal and stimulated insulin secretion. Moreover, clock disruption has a profound impact on the human islet transcriptome, in particular, on the genes involved in insulin secretion.

Differential selection after duplication in mammalian developmental genes.

Gene duplication provides the opportunity for subsequent refinement of distinct functions of the duplicated copies. Either through changes in coding sequence or changes in regulatory regions, duplicate copies appear to obtain new or tissue-specific functions. If this divergence were driven by natural selection, we would expect duplicated copies to have differentiated patterns of substitutions. We tested this hypothesis using genes that duplicated before the human/mouse split and whose orthologous relations were clear. The null hypothesis is that the number of amino acid changes between humans and mice was distributed similarly across different paralogs. We used a method modified from Tang and Lewontin to detect heterogeneity in the amino acid substitution pattern between those different paralogs. Our results show that many of the paralogous gene pairs appear to be under differential selection in the human/mouse comparison. The properties that led to diversification appear to have arisen before the split of the human and mouse lineages. Further study of the diverged genes revealed insights regarding the patterns of amino acid substitution that resulted in differences in function and/or expression of these genes. This approach has utility in the study of newly identified members of gene families in genomewide data mining and for contrasting the merits of alternative hypotheses for the evolutionary divergence of function of duplicated genes.

Nontransitivity of sperm precedence in Drosophila.

Sperm competition is an important component of fitness in Drosophila, but we still do not have a clear understanding of the unit of selection that is relevant to sperm competition. Here we demonstrate that sperm competitive ability is not a property of the sperm haplotype, but rather of the diploid male's genotype. Then we test whether the relative sperm competitive ability of males can be ranked on a linear array or whether competitive ability instead depends on particular pairwise contests among males. Sperm precedence of six chromosome-extracted lines was tested against three different visible marker lines (cn bw, bwD, and Cy), and the rank order of the six lines differed markedly among the mutant lines. Population genetic theory has shown that departures from transitivity of sperm precedence may be important to the maintenance of polymorphism for genes that influence sperm competitive ability. The nontransitivity seen in sperm precedence should theoretically increase the opportunity for polymorphism in genes that influence this phenotype.

Non-Mendelian segregation of sex chromosomes in heterospecific Drosophila males.

Interspecific hybrids and backcrossed organisms generally suffer from reduced viability and/or fertility. To identify and genetically map these defects, we introgressed regions of the Drosophila sechellia genome into the D. simulans genome. A female-biased sex ratio was observed in 24 of the 221 recombinant inbred lines, and subsequent tests attributed the skew to failure of Y-bearing sperm to fertilize the eggs. Apparently these introgressed lines fail to suppress a normally silent meiotic drive system. Using molecular markers we mapped two regions of the Drosophila genome that appear to exhibit differences between D. simulans and D. sechellia in their regulation of sex chromosome segregation distortion. The data indicate that the sex ratio phenotype results from an epistatic interaction between at least two factors. We discuss whether this observation is relevant to the meiotic drive theory of hybrid male sterility.

Negative covariance suggests mutation bias in a two-locus microsatellite system in the fish Sparus aurata.

Constraints on microsatellite length appear to vary in a species-specific manner. We know very little about the nature of these constraints and why they should vary among species. While surveying microsatellite variation in the Mediterranean gilthead sea bream, Sparus aurata, we discovered an unusual pattern of covariation between two closely linked microsatellite loci. One- and two-locus haplotypes were scored from PCR amplification products of each locus separately and both loci together. In a sample of 211 fish, there was a strong negative covariance in repeat number between the two loci, which suggests a mechanism that maintains the combined length below a constrained size. In addition, there were two clusters of the same combined haplotype length, one consisting of a long repeat array at one locus and a short array at the other and vice versa. We demonstrate that several models of biased mutation or natural selection, in theory, could generate this pattern of covariance. The common feature of all the models is the idea that tightly linked microsatellites do not evolve in complete independence, and that whatever size dependence there is to the process, it appears to "read" the combined size of the two loci.