A novel framework for functional annotation of variants of uncertain significance in ID/ASD risk gene CC2D1A.

Intellectual disability (ID) and autism spectrum disorder (ASD) are genetically heterogeneous with hundreds of identified risk genes, most affecting only a few patients. Novel missense variants in these genes are being discovered as clinical exome sequencing is now routinely integrated into diagnosis, yet most of them are annotated as variants of uncertain significance (VUS). VUSs are a major roadblock in using patient genetics to inform clinical action. We developed a framework to characterize VUSs in Coiled-coil and C2 domain containing 1A (CC2D1A), a gene causing autosomal recessive ID with comorbid ASD in 40% of cases. We analyzed seven VUSs (p.Pro319Leu, p.Ser327Leu, p.Gly441Val, p.Val449Met, p.Thr580Ile, p.Arg886His and p.Glu910Lys) from four cases of individuals with ID and ASD. Variants were cloned and overexpressed in HEK293 individually and in their respective heterozygous combination. CC2D1A is a signaling scaffold that positively regulates PKA-CREB signaling by repressing phosphodiesterase 4D (PDE4D) to prevent cAMP degradation. After testing multiple parameters including direct interaction between PDE4D and CC2D1A, cAMP levels and CREB activation, we found that the most sensitive readout was CREB transcriptional activity using a luciferase assay. Compared to WT CC2D1A, five VUSs (p.Pro319Leu, p.Gly441Val, p.Val449Met, p.Thr580Ile, and p.Arg886His) led to significantly blunted response to forskolin induced CREB activation. This luciferase assay approach can be scaled up to annotate ~150 CC2D1A VUSs that are currently listed in ClinVar. Since CREB activation is a common denominator for multiple ASD/ID genes, our paradigm can also be adapted for their VUSs.

Fine structures of intrinsically disordered proteins.

We report simulation studies of 33 single intrinsically disordered proteins (IDPs) using coarse-grained bead-spring models where interactions among different amino acids are introduced through a hydropathy matrix and additional screened Coulomb interaction for the charged amino acid beads. Our simulation studies of two different hydropathy scales (HPS1, HPS2) [Dignon et al., PLoS Comput. Biol. 14, e1005941 (2018); Tesei et al. Proc. Natl. Acad. Sci. U. S. A. 118, e2111696118 (2021)] and the comparison with the existing experimental data indicate an optimal interaction parameter ϵ = 0.1 and 0.2 kcal/mol for the HPS1 and HPS2 hydropathy scales. We use these best-fit parameters to investigate both the universal aspects as well as the fine structures of the individual IDPs by introducing additional characteristics. (i) First, we investigate the polymer-specific scaling relations of the IDPs in comparison to the universal scaling relations [Bair et al., J. Chem. Phys. 158, 204902 (2023)] for the homopolymers. By studying the scaled end-to-end distances ⟨RN2⟩/(2Lℓp) and the scaled transverse fluctuations l̃⊥2=⟨l⊥2⟩/L, we demonstrate that IDPs are broadly characterized with a Flory exponent of ν ≃ 0.56 with the conclusion that conformations of the IDPs interpolate between Gaussian and self-avoiding random walk chains. Then, we introduce (ii) Wilson charge index (W) that captures the essential features of charge interactions and distribution in the sequence space and (iii) a skewness index (S) that captures the finer shape variation of the gyration radii distributions as a function of the net charge per residue and charge asymmetry parameter. Finally, our study of the (iv) variation of ⟨Rg⟩ as a function of salt concentration provides another important metric to bring out finer characteristics of the IDPs, which may carry relevant information for the origin of life.

Universality in conformations and transverse fluctuations of a semi-flexible polymer in a crowded environment.

We study the universal aspects of polymer conformations and transverse fluctuations for a single swollen chain characterized by a contour length L and a persistence length ℓp in two dimensions (2D) and three dimensions (3D) in the bulk, as well as in the presence of excluded volume (EV) particles of different sizes occupying different area/volume fractions. In the absence of the EV particles, we extend the previously established universal scaling relations in 2D [Huang et al., J. Chem. 140, 214902 (2014)] to include 3D and demonstrate that the scaled end-to-end distance ⟨RN2⟩/(2Lℓp) and the scaled transverse fluctuation ⟨l⊥2⟩/L as a function of L/ℓp collapse onto the same master curve, where ⟨RN2⟩ and ⟨l⊥2⟩ are the mean-square end-to-end distance and transverse fluctuations. However, unlike in 2D, where the Gaussian regime is absent due to the extreme dominance of the EV interaction, we find that the Gaussian regime is present, albeit very narrow in 3D. The scaled transverse fluctuation in the limit L/ℓp ≪ 1 is independent of the physical dimension and scales as ⟨l⊥2⟩/L∼(L/ℓp)ζ-1, where ζ = 1.5 is the roughening exponent. For L/ℓp ≫ 1, the scaled fluctuation scales as ⟨l⊥2⟩/L∼(L/ℓp)ν-1, where ν is the Flory exponent for the corresponding spatial dimension (ν2D = 0.75 and ν3D = 0.58). When EV particles of different sizes for different area or volume fractions are added into 2D and 3D systems, our results indicate that the crowding density either does not or does only weakly affect the universal scaling relations. We discuss the implications of these results in living matter by showing the experimental result for a dsDNA on the master plot.

How capture affects polymer translocation in a solitary nanopore.

DNA capture with high fidelity is an essential part of nanopore translocation. We report several important aspects of the capture process and subsequent translocation of a model DNA polymer through a solid-state nanopore in the presence of an extended electric field using the Brownian dynamics simulation that enables us to record statistics of the conformations at every stage of the translocation process. By releasing the equilibrated DNAs from different equipotentials, we observe that the capture time distribution depends on the initial starting point and follows a Poisson process. The field gradient elongates the DNA on its way toward the nanopore and favors a successful translocation even after multiple failed threading attempts. Even in the limit of an extremely narrow pore, a fully flexible chain has a finite probability of hairpin-loop capture, while this probability decreases for a stiffer chain and promotes single file translocation. Our in silico studies identify and differentiate characteristic distributions of the mean first passage time due to single file translocation from those due to translocation of different types of folds and provide direct evidence of the interpretation of the experimentally observed folds [M. Gershow and J. A. Golovchenko, Nat. Nanotechnol. 2, 775 (2007) and Mihovilovic et al., Phys. Rev. Lett. 110, 028102 (2013)] in a solitary nanopore. Finally, we show a new finding-that a charged tag attached at the 5' end of the DNA enhances both the multi-scan rate and the uni-directional translocation (5' → 3') probability that would benefit the genomic barcoding and sequencing experiments.

Polymer escape through a three dimensional double-nanopore system.

We study the escape dynamics of a double-stranded DNA (dsDNA) through an idealized double nanopore geometry subject to two equal and opposite forces (tug-of-war) using Brownian dynamics (BD) simulation. In addition to the geometrical restrictions imposed on the cocaptured dsDNA segment in between the pores, the presence of tug-of-war forces at each pore results in a variation of the local chain stiffness for the segment of the chain in between the pores, which increases the overall stiffness of the chain. We use the BD simulation results to understand how the intrinsic chain stiffness and the tug-of-war forces affect the escape dynamics by monitoring the local chain persistence length ℓp, the residence time of the individual monomers W(m) in the nanopores, and the chain length dependence of the escape time ⟨τ⟩ and its distribution. Finally, we generalize the scaling theory for the unbiased single nanopore translocation for a fully flexible chain for the escape of a semi-flexible chain through a double nanopore in the presence of tug-of-war forces. We establish that the stiffness dependent part of the escape time is approximately independent of the translocation mechanism so that ⟨τ⟩∼ℓp 2/D+2, and therefore, the generalized escape time for a semi-flexible chain can be written as ⟨τ⟩=ANαℓp 2/D+2. We use the BD simulation results to compare the predictions of the scaling theory. Our numerical studies supplemented by scaling analysis provide fundamental insights to design new experiments where a dsDNA moves slowly through a series of graphene nanopores.

Population diversity and adaptive evolution in keratinization genes: impact of environment in shaping skin phenotypes.

Several studies have demonstrated the role of climatic factors in shaping skin phenotypes, particularly pigmentation. Keratinization is another well-designed feature of human skin, which is involved in modulating transepidermal water loss (TEWL). Although this physiological process is closely linked to climate, presently it is not clear whether genetic diversity is observed in keratinization and whether this process also responds to the environmental pressure. To address this, we adopted a multipronged approach, which involved analysis of 1) copy number variations in diverse Indian and HapMap populations from varied geographical regions; 2) genetic association with geoclimatic parameters in 61 populations of dbCLINE database in a set of 549 genes from four processes namely keratinization, pigmentation, epidermal differentiation, and housekeeping functions; 3) sequence divergence in 4,316 orthologous promoters and corresponding exonic regions of human and chimpanzee with macaque as outgroup, and 4) protein sequence divergence (Ka/Ks) across nine vertebrate classes, which differ in their extent of TEWL. Our analyses demonstrate that keratinization and epidermal differentiation genes are under accelerated evolution in the human lineage, relative to pigmentation and housekeeping genes. We show that this entire pathway may have been driven by environmental selection pressure through concordant functional polymorphisms across several genes involved in skin keratinization. Remarkably, this underappreciated function of skin may be a crucial determinant of adaptation to diverse environmental pressures across world populations.

Semiflexible macromolecules in quasi-one-dimensional confinement: Discrete versus continuous bond angles.

The conformations of semiflexible polymers in two dimensions confined in a strip of width D are studied by computer simulations, investigating two different models for the mechanism by which chain stiffness is realized. One model (studied by molecular dynamics) is a bead-spring model in the continuum, where stiffness is controlled by a bond angle potential allowing for arbitrary bond angles. The other model (studied by Monte Carlo) is a self-avoiding walk chain on the square lattice, where only discrete bond angles (0° and ±90°) are possible, and the bond angle potential then controls the density of kinks along the chain contour. The first model is a crude description of DNA-like biopolymers, while the second model (roughly) describes synthetic polymers like alkane chains. It is first demonstrated that in the bulk the crossover from rods to self-avoiding walks for both models is very similar, when one studies average chain linear dimensions, transverse fluctuations, etc., despite their differences in local conformations. However, in quasi-one-dimensional confinement two significant differences between both models occur: (i) The persistence length (extracted from the average cosine of the bond angle) gets renormalized for the lattice model when D gets less than the bulk persistence length, while in the continuum model it stays unchanged. (ii) The monomer density near the repulsive walls for semiflexible polymers is compatible with a power law predicted for the Kratky-Porod model in the case of the bead-spring model, while for the lattice case it tends to a nonzero constant across the strip. However, for the density of chain ends, such a constant behavior seems to occur for both models, unlike the power law observed for flexible polymers. In the regime where the bulk persistence length ℓp is comparable to D, hairpin conformations are detected, and the chain linear dimensions are discussed in terms of a crossover from the Daoud/De Gennes "string of blobs"-picture to the flexible rod picture when D decreases and/or the chain stiffness increases. Introducing a suitable further coarse-graining of the chain contours of the continuum model, direct estimates for the deflection length and its distribution could be obtained.

Conformations, transverse fluctuations, and crossover dynamics of a semi-flexible chain in two dimensions.

We present a unified scaling description for the dynamics of monomers of a semiflexible chain under good solvent condition in the free draining limit. We consider both the cases where the contour length L is comparable to the persistence length ℓ(p) and the case L ≫ ℓ(p). Our theory captures the early time monomer dynamics of a stiff chain characterized by t(3/4) dependence for the mean square displacement of the monomers, but predicts a first crossover to the Rouse regime of t(2ν/1 + 2ν) for τ¹ ~ ℓ(p)³, and a second crossover to the purely diffusive dynamics for the entire chain at τ2 ∼ L(5/2). We confirm the predictions of this scaling description by studying monomer dynamics of dilute solution of semi-flexible chains under good solvent conditions obtained from our Brownian dynamics (BD) simulation studies for a large choice of chain lengths with number of monomers per chain N = 16-2048 and persistence length ℓ(p) = 1-500 Lennard-Jones units. These BD simulation results further confirm the absence of Gaussian regime for a two-dimensional (2D) swollen chain from the slope of the plot of ⟨R(N)²⟩/2Lℓ(p) ~ L/ℓ(p) which around L/ℓ(p) ∼ 1 changes suddenly from (L/ℓ(p)) → (L/ℓ(p))(0.5), also manifested in the power law decay for the bond autocorrelation function disproving the validity of the worm-like-chain in 2D. We further observe that the normalized transverse fluctuations of the semiflexible chains for different stiffness √(⟨l(⊥)²⟩)/L as a function of renormalized contour length L/ℓ(p) collapse on the same master plot and exhibits power law scaling √(⟨l(⊥)²⟩)/L ~ (L/ℓ(p))(n) at extreme limits, where η = 0.5 for extremely stiff chains (L/ℓ(p) ≫ 1), and η = -0.25 for fully flexible chains. Finally, we compare the radial distribution functions obtained from our simulation studies with those obtained analytically.

DNA Barcodes Using a Dual Nanopore Device.

We report a novel method based on the current blockade (CB) characteristics obtained from a dual nanopore device that can determine DNA barcodes with near-perfect accuracy using a Brownian dynamics simulation strategy. The method supersedes our previously reported velocity correction algorithm (S. Seth and A. Bhattacharya, RSC Advances, 11:20781-20787, 2021), taking advantage of the better measurement of the time-of-flight (TOF) protocol offered by the dual nanopore setup. We demonstrate the efficacy of the method by comparing our simulation data from a coarse-grained model of a polymer chain consisting of 2048 excluded volume beads of diameter 𝜎 = 24 bp using with those obtained from experimental CB data from a 48,500 bp λ-phage DNA, providing a 48500 2400 ≅ 24 base pair resolution in simulation. The simulation time scale is compared to the experimental time scale by matching the simulated time-of-flight (TOF) velocity distributions with those obtained experimentally (Rand et al., ACS Nano, 16:5258-5273, 2022). We then use the evolving coordinates of the dsDNA and the molecular features to reconstruct the current blockade characteristics on the fly using a volumetric model based on the effective van der Waal radii of the species inside and in the immediate vicinity of the pore. Our BD simulation mimics the control-zoom-in-logic to understand the origin of the TOF distributions due to the relaxation of the out-of-equilibrium conformations followed by a reversal of the electric fields. The simulation algorithm is quite general and can be applied to differentiate DNA barcodes from different species.

Driven translocation of a semi-flexible chain through a nanopore: a Brownian dynamics simulation study in two dimensions.

We study translocation dynamics of a semi-flexible polymer chain through a nanoscopic pore in two dimensions using Langevin dynamics simulation in presence of an external bias F inside the pore. For chain length N and stiffness parameter κb considered in this paper, we observe that the mean first passage time <τ> increases as <τ(κb)>~<τ(κb=0)>lp(aN) , where κb and lp are the stiffness parameter and persistence length, respectively, and aN is a constant that has a weak N dependence. We monitor the time dependence of the last monomer xN(t) at the cis compartment and calculate the tension propagation time (TP) ttp directly from simulation data for ~ t as alluded in recent nonequlibrium TP theory [T. Sakaue, Phys. Rev. E 76, 021803 (2007)] and its modifications to Brownian dynamics tension propagation theory [T. Ikonen, A. Bhattacharya, T. Ala-Nissila, and W. Sung, Phys. Rev. E 85, 051803 (2012); and J. Chem. Phys. 137, 085101 (2012)] originally developed to study translocation of a fully flexible chain. We also measure ttp from peak position of the waiting time distribution W(s) of the translocation coordinate s (i.e., the monomer inside the pore), and explicitly demonstrate the underlying TP picture along the chain backbone of a translocating chain to be valid for semi-flexible chains as well. From the simulation data, we determine the dependence of ttp on chain persistence length lp and show that the ratio ttp∕<τ> is independent of the bias F.

Knot formation of dsDNA pushed inside a nanochannel.

Recent experiments demonstrated that knots in single molecule dsDNA can be formed by compression in a nanochannel. In this manuscript, we further elucidate the underlying molecular mechanisms by carrying out a compression experiment in silico, where an equilibrated coarse-grained double-stranded DNA confined in a square channel is pushed by a piston. The probability of forming knots is a non-monotonic function of the persistence length and can be enhanced significantly by increasing the piston speed. Under compression knots are abundant and delocalized due to a backfolding mechanism from which chain-spanning loops emerge, while knots are less frequent and only weakly localized in equilibrium. Our in silico study thus provides insights into the formation, origin and control of DNA knots in nanopores.

Discriminating protein tags on a dsDNA construct using a Dual Nanopore Device.

We report Brownian dynamics simulation results with the specific goal to identify key parameters controlling the experimentally measurable characteristics of protein tags on a dsDNA construct translocating through a double nanopore setup. First, we validate the simulation scheme in silico by reproducing and explaining the physical origin of the asymmetric experimental dwell time distributions of the oligonucleotide flap markers on a 48 kbp long dsDNA at the left and the right pore. We study the effect of the electric field inside and beyond the pores, critical to discriminate the protein tags based on their effective charges and masses revealed through a generic power-law dependence of the average dwell time at each pore. The simulation protocols monitor piecewise dynamics at a sub-nanometer length scale and explain the disparate velocity using the concepts of nonequilibrium tension propagation theory. We further justify the model and the chosen simulation parameters by calculating the Péclet number which is in close agreement with the experiment. We demonstrate that our carefully chosen simulation strategies can serve as a powerful tool to discriminate different types of neutral and charged tags of different origins on a dsDNA construct in terms of their physical characteristics and can provide insights to increase both the efficiency and accuracy of an experimental dual-nanopore setup.

Increased TGFß1 and SMAD3 Contribute to Age-Related Aortic Valve Calcification.

Aims: Calcific aortic valve disease (CAVD) is a progressive heart disease that is particularly prevalent in elderly patients. The current treatment of CAVD is surgical valve replacement, but this is not a permanent solution, and it is very challenging for elderly patients. Thus, a pharmacological intervention for CAVD may be beneficial. In this study, we intended to rescue aortic valve (AV) calcification through inhibition of TGFß1 and SMAD3 signaling pathways. Methods and Results: The klotho gene, which was discovered as an aging-suppressor gene, has been observed to play a crucial role in AV calcification. The klotho knockout (Kl -/-) mice have shorter life span (8-12 weeks) and develop severe AV calcification. Here, we showed that increased TGFß1 and TGFß-dependent SMAD3 signaling were associated with AV calcification in Kl -/- mice. Next, we generated Tgfb1- and Smad3-haploinsufficient Kl -/- mice to determine the contribution of TGFß1 and SMAD3 to the AV calcification in Kl -/- mice. The histological and morphometric evaluation suggested a significant reduction of AV calcification in Kl -/-; Tgfb1 ± mice compared to Kl -/- mice. Smad3 heterozygous deletion was observed to be more potent in reducing AV calcification in Kl -/- mice compared to the Kl -/-; Tgfb1 ± mice. We observed significant inhibition of Tgfb1, Pai1, Bmp2, Alk2, Spp1, and Runx2 mRNA expression in Kl -/-; Tgfb1 ± and Kl -/-; Smad3 ± mice compared to Kl -/- mice. Western blot analysis confirmed that the inhibition of TGFß canonical and non-canonical signaling pathways were associated with the rescue of AV calcification of both Kl -/-; Tgfb1 ± and Kl -/-; Smad3 ± mice. Conclusion: Overall, inhibition of the TGFß1-dependent SMAD3 signaling pathway significantly blocks the development of AV calcification in Kl -/- mice. This information is useful in understanding the signaling mechanisms involved in CAVD.

Unraveling the mysteries of MYT1L: From reprogramming factor to multifaceted regulator of neuronal differentiation.

In this issue of Neuron, Chen et al. (2021) generated a mouse model for haploinsufficiency of MYT1L. MYT1L is widely used in neuronal reprogramming, and de novo mutations have been linked to a neurodevelopmental syndrome. Extensive characterization in this study better delineates MYT1L's role in transcriptional regulation and neuronal differentiation.

DNA barcodes using a double nanopore system.

The potential of a double nanopore system to determine DNA barcodes has been demonstrated experimentally. By carrying out Brownian dynamics simulation on a coarse-grained model DNA with protein tag (barcodes) at known locations along the chain backbone, we demonstrate that due to large variation of velocities of the chain segments between the tags, it is inevitable to under/overestimate the genetic lengths from the experimental current blockade and time of flight data. We demonstrate that it is the tension propagation along the chain's backbone that governs the motion of the entire chain and is the key element to explain the non uniformity and disparate velocities of the tags and DNA monomers under translocation that introduce errors in measurement of the length segments between protein tags. Using simulation data we further demonstrate that it is important to consider the dynamics of the entire chain and suggest methods to accurately decipher barcodes. We introduce and validate an interpolation scheme using simulation data for a broad distribution of tag separations and suggest how to implement the scheme experimentally.

DNA barcode by flossing through a cylindrical nanopore.

We report an accurate method to determine DNA barcodes from the dwell time measurement of protein tags (barcodes) along the DNA backbone using Brownian dynamics simulation of a model DNA and use a recursive theoretical scheme which improves the measurements to almost 100% accuracy. The heavier protein tags along the DNA backbone introduce a large speed variation in the chain that can be understood using the idea of non-equilibrium tension propagation theory. However, from an initial rough characterization of velocities into "fast" (nucleotides) and "slow" (protein tags) domains, we introduce a physically motivated interpolation scheme that enables us to determine the barcode velocities rather accurately. Our theoretical analysis of the motion of the DNA through a cylindrical nanopore opens up the possibility of its experimental realization and carries over to multi-nanopore devices used for barcoding.

Myocardial TGFß2 Is Required for Atrioventricular Cushion Remodeling and Myocardial Development.

Among the three transforming growth factor beta (TGFß) ligands, TGFß2 is essential for heart development and is produced by multiple cell types, including myocardium. Heterozygous mutations in TGFB2 in patients of connective tissue disorders result in congenital heart defects and adult valve malformations, including mitral valve prolapse (MVP) with or without regurgitation. Tgfb2 germline knockout fetuses exhibit multiple cardiac defects but the role of myocardial-TGFß2 in heart development is yet to be elucidated. Here, myocardial Tgfb2 conditional knockout (CKO) embryos were generated by crossing Tgfb2flox mice with Tgfb2+/-; cTntCre mice. Tgfb2flox/- embryos were normal, viable. Cell fate mapping was done using dual-fluorescent mT/mG+/- mice. Cre-mediated Tgfb2 deletion was assessed by genomic PCR. RNAscope in situ hybridization was used to detect the loss of myocardial Tgfb2 expression. Histological, morphometric, immunohistochemical, and in situ hybridization analyses of CKOs and littermate controls at different stages of heart development (E12.5-E18.5) were used to determine the role of myocardium-derived TGFß2 in atrioventricular (AV) cushion remodeling and myocardial development. CKOs exhibit a thin ventricular myocardium, AV cushion remodeling defects and developed incomplete AV septation defects. The loss of myocardial Tgfb2 resulted in impaired cushion maturation and dysregulated cell death. Phosphorylated SMAD2, a surrogate for TGFß signaling, was "paradoxically" increased in both AV cushion mesenchyme and ventricular myocardium in the CKOs. Our results indicate that TGFß2 produced by cardiomyocytes acting as cells autonomously on myocardium and via paracrine signaling on AV cushions are required for heart development.

Multiple Alu Exonization in 3'UTR of a Primate-Specific Isoform of CYP20A1 Creates a Potential miRNA Sponge.

Alu repeats contribute to phylogenetic novelties in conserved regulatory networks in primates. Our study highlights how exonized Alus could nucleate large-scale mRNA-miRNA interactions. Using a functional genomics approach, we characterize a transcript isoform of an orphan gene, CYP20A1 (CYP20A1_Alu-LT) that has exonization of 23 Alus in its 3'UTR. CYP20A1_Alu-LT, confirmed by 3'RACE, is an outlier in length (9 kb 3'UTR) and widely expressed. Using publically available data sets, we demonstrate its expression in higher primates and presence in single nucleus RNA-seq of 15,928 human cortical neurons. miRanda predicts ∼4,700 miRNA recognition elements (MREs) for ∼1,000 miRNAs, primarily originated within these 3'UTR-Alus. CYP20A1_Alu-LT could be a potential multi-miRNA sponge as it harbors ≥10 MREs for 140 miRNAs and has cytosolic localization. We further tested whether expression of CYP20A1_Alu-LT correlates with mRNAs harboring similar MRE targets. RNA-seq with conjoint miRNA-seq analysis was done in primary human neurons where we observed CYP20A1_Alu-LT to be downregulated during heat shock response and upregulated in HIV1-Tat treatment. In total, 380 genes were positively correlated with its expression (significantly downregulated in heat shock and upregulated in Tat) and they harbored MREs for nine expressed miRNAs which were also enriched in CYP20A1_Alu-LT. MREs were significantly enriched in these 380 genes compared with random sets of differentially expressed genes (P = 8.134e-12). Gene ontology suggested involvement of these genes in neuronal development and hemostasis pathways thus proposing a novel component of Alu-miRNA-mediated transcriptional modulation that could govern specific physiological outcomes in higher primates.

Tug of war in a double-nanopore system.

We simulate a tug-of-war (TOW) scenario for a model double-stranded DNA threading through a double nanopore (DNP) system. The DNA, simultaneously captured at both pores, is subject to two equal and opposite forces -f[over ⃗]_{L}=f[over ⃗]_{R} (TOW), where f[over ⃗]_{L} and f[over ⃗]_{R} are the forces applied to the left and the right pore, respectively. Even though the net force on the DNA polymer Δf[over ⃗]_{LR}=f[over ⃗]_{L}+f[over ⃗]_{R}=0, the mean first passage time (MFPT) 〈τ〉 depends on the magnitude of the TOW forces |f_{L}|=|f_{R}|=f_{LR}. We qualitatively explain this dependence of 〈τ〉 on f_{LR} from the known results for the single-pore translocation of a triblock copolymer A-B-A with ℓ_{pB}>ℓ_{pA}, where ℓ_{pA} and ℓ_{pB} are the persistence length of the A and B segments, respectively. We demonstrate that the time of flight of a monomer with index m [〈τ_{LR}(m)〉] from one pore to the other exhibits quasiperiodic structure commensurate with the distance between the pores d_{LR}. Finally, we study the situation where we offset the TOW biases so that Δf[over ⃗]_{LR}=f[over ⃗]_{L}+f[over ⃗]_{R}≠0, and qualitatively reproduce the experimental result of the dependence of the MFPT on Δf[over ⃗]_{LR}. We demonstrate that, for a moderate bias, the MFPT for the DNP system for a chain length N follows the same scaling ansatz as that for the single nanopore, 〈τ〉=(AN^{1+ν}+η_{pore}N)(Δf_{LR})^{-1}, where η_{pore} is the pore friction, which enables us to estimate 〈τ〉 for a long chain. Our Brownian dynamics simulation studies provide fundamental insights and valuable information about the details of the translocation speed obtained from 〈τ_{LR}(m)〉, and accuracy of the translation of the data obtained in the time domain to units of genomic distances.

Transforming Growth Factor Beta3 is Required for Cardiovascular Development.

Transforming growth factor beta3 (TGFB3) gene mutations in patients of arrhythmogenic right ventricular dysplasia/cardiomyopathy (ARVD1) and Loeys-Dietz syndrome-5 (LDS5)/Rienhoff syndrome are associated with cardiomyopathy, cardiac arrhythmia, cardiac fibrosis, cleft palate, aortic aneurysms, and valvular heart disease. Although the developing heart of embryos express Tgfb3, its overarching role remains unclear in cardiovascular development and disease. We used histological, immunohistochemical, and molecular analyses of Tgfb3-/- fetuses and compared them to wildtype littermate controls. The cardiovascular phenotypes were diverse with approximately two thirds of the Tgfb3-/- fetuses having one or more cardiovascular malformations, including abnormal ventricular myocardium (particularly of the right ventricle), outflow tract septal and alignment defects, abnormal aortic and pulmonary trunk walls, and thickening of semilunar and/or atrioventricular valves. Ventricular septal defects (VSD) including the perimembranous VSDs were observed in Tgfb3-/- fetuses with myocardial defects often accompanied by the muscular type VSD. In vitro studies using TGFß3-deficient fibroblasts in 3-D collagen lattice formation assays indicated that TGFß3 was required for collagen matrix reorganization. Biochemical studies indicated the 'paradoxically' increased activation of canonical (SMAD-dependent) and noncanonical (MAP kinase-dependent) pathways. TGFß3 is required for cardiovascular development to maintain a balance of canonical and noncanonical TGFß signaling pathways.